zstxb 发表于 2020-12-3 22:22:23

python求解复杂方程

0.978147600733806*(25.14*(1 - 0.039788672413643*sqrt(632.0196/(0.0526732314729343*sin(phi1)**2 + 1) + 1.75351898592098*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**3))/(sqrt(1 + 0.00277446931380131*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**2)*sqrt(0.0526732314729343*sin(phi1)**2 + 1)) + 0.069750158548965*sin(phi1)**2*cos(phi1)**2/(sqrt(1 + 0.00277446931380131*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**2)*(0.0526732314729343*sin(phi1)**2 + 1)**(5/2)))/sqrt((0.0533469650208306*sin(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(3/2) - 0.0533469650208306*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(3/2) + 0.00842987111076221*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(5/2))**2 + 1.02574523451281/(0.0526732314729343*sin(phi1)**2 + 1)) + 0.207911690817759*(1.32420503922957*(1 - 0.039788672413643*sqrt(632.0196/(0.0526732314729343*sin(phi1)**2 + 1) + 1.75351898592098*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**3))*sin(phi1)*cos(phi1)/(sqrt(1 + 0.00277446931380131*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**2)*(0.0526732314729343*sin(phi1)**2 + 1)**(3/2)) - 1.32420503922957*sin(phi1)*cos(phi1)/(sqrt(1 + 0.00277446931380131*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**2)*(0.0526732314729343*sin(phi1)**2 + 1)**(3/2)))/sqrt((0.0533469650208306*sin(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(3/2) - 0.0533469650208306*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(3/2) + 0.00842987111076221*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(5/2))**2 + 1.02574523451281/(0.0526732314729343*sin(phi1)**2 + 1)) + 0.539308918904122*(-0.039788672413643*sqrt(632.0196/(0.0526732314729343*sin(phi1)**2 + 1) + 1.75351898592098*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**3) + 1 + (-0.0526732314729343*sin(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1) + 0.0526732314729343*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1) - 0.00554893862760263*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**2)/(1 + 0.00277446931380131*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**2))/sqrt((0.0533469650208306*sin(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(3/2) - 0.0533469650208306*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(3/2) + 0.00842987111076221*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(5/2))**2 + 1.02574523451281/(0.0526732314729343*sin(phi1)**2 + 1)) + 0.0383663379417315*(0.039788672413643*sqrt(632.0196/(0.0526732314729343*sin(phi1)**2 + 1) + 1.75351898592098*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**3) - 1 - (-0.0526732314729343*sin(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1) + 0.0526732314729343*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1) - 0.00554893862760263*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**2)/(1 + 0.00277446931380131*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**2))/sqrt((0.0533469650208306*sin(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(3/2) - 0.0533469650208306*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(3/2) + 0.00842987111076221*sin(phi1)**2*cos(phi1)**2/(0.0526732314729343*sin(phi1)**2 + 1)**(5/2))**2 + 1.02574523451281/(0.0526732314729343*sin(phi1)**2 + 1))

coderyzl 发表于 2020-12-3 22:54:05

你疯了吧

wp231957 发表于 2020-12-4 08:36:43

这里的phi1 应该有初始值

import math

>>> phi1=60
>>> 0.978147600733806*(25.14*(1 - 0.039788672413643*math.sqrt(632.0196/(0.0526732314729343*math.sin(phi1)**2 + 1) + 1.75351898592098*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(ph
i1)**2 + 1)**3))/(math.sqrt(1 + 0.00277446931380131*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**2)*math.sqrt(0.0526732314729343*math.sin(phi1)**2 + 1)) + 0.06975015
8548965*math.sin(phi1)**2*math.cos(phi1)**2/(math.sqrt(1 + 0.00277446931380131*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**2)*(0.0526732314729343*math.sin(phi1)**2
+ 1)**(5/2)))/math.sqrt((0.0533469650208306*math.sin(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(3/2) - 0.0533469650208306*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(3/2)
+ 0.00842987111076221*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(5/2))**2 + 1.02574523451281/(0.0526732314729343*math.sin(phi1)**2 + 1)) + 0.207911690817759*(1.3
2420503922957*(1 - 0.039788672413643*math.sqrt(632.0196/(0.0526732314729343*math.sin(phi1)**2 + 1) + 1.75351898592098*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**3)
)*math.sin(phi1)*math.cos(phi1)/(math.sqrt(1 + 0.00277446931380131*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**2)*(0.0526732314729343*math.sin(phi1)**2 + 1)**(3/2))
- 1.32420503922957*math.sin(phi1)*math.cos(phi1)/(math.sqrt(1 + 0.00277446931380131*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**2)*(0.0526732314729343*math.sin(phi
1)**2 + 1)**(3/2)))/math.sqrt((0.0533469650208306*math.sin(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(3/2) - 0.0533469650208306*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)*
*(3/2) + 0.00842987111076221*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(5/2))**2 + 1.02574523451281/(0.0526732314729343*math.sin(phi1)**2 + 1)) + 0.53930891890412
2*(-0.039788672413643*math.sqrt(632.0196/(0.0526732314729343*math.sin(phi1)**2 + 1) + 1.75351898592098*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**3) + 1 + (-0.0526
732314729343*math.sin(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1) + 0.0526732314729343*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1) - 0.00554893862760263*math.sin(phi1)**2*math
.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**2)/(1 + 0.00277446931380131*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**2))/math.sqrt((0.0533469650208306*
math.sin(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(3/2) - 0.0533469650208306*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(3/2) + 0.00842987111076221*math.sin(phi1)**2*mat
h.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(5/2))**2 + 1.02574523451281/(0.0526732314729343*math.sin(phi1)**2 + 1)) + 0.0383663379417315*(0.039788672413643*math.sqrt(632.0196/(0.052673
2314729343*math.sin(phi1)**2 + 1) + 1.75351898592098*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**3) - 1 - (-0.0526732314729343*math.sin(phi1)**2/(0.0526732314729343
*math.sin(phi1)**2 + 1) + 0.0526732314729343*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1) - 0.00554893862760263*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**
2 + 1)**2)/(1 + 0.00277446931380131*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**2))/math.sqrt((0.0533469650208306*math.sin(phi1)**2/(0.0526732314729343*math.sin(phi
1)**2 + 1)**(3/2) - 0.0533469650208306*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)**2 + 1)**(3/2) + 0.00842987111076221*math.sin(phi1)**2*math.cos(phi1)**2/(0.0526732314729343*math.sin(phi1)*
*2 + 1)**(5/2))**2 + 1.02574523451281/(0.0526732314729343*math.sin(phi1)**2 + 1))
结果:-0.0013788634210874233
页: [1]
查看完整版本: python求解复杂方程