| 查看: 394 | 回复: 0 | ||
[求助]
solve函数解方程得出的解是错误的，是软件算错了吗？
|
|
syms r_2 d_1 d_2 d_3 f_SYS eq=sym('- (d_3 - f_SYS)/f_SYS - ((f_SYS*(4*f_SYS*d_1^2*d_2^2 - 2*d_1*d_2*r_2^3 + 4*f_SYS*d_1*d_2*r_2^2 + d_1*r_2^4 - d_2*r_2^4))^(1/2) - 2*d_1*d_2*f_SYS)/(f_SYS*r_2^2)=0'); solve(eq,'r_2'); 运行后，得出r_2的两个解： ((d_1*d_2*f_SYS*(4*d_3*f_SYS^2 - 8*d_3^2*f_SYS + 4*d_3^3 + d_1*d_2*f_SYS - 4*d_1*d_3*f_SYS + 4*d_2*d_3*f_SYS))^(1/2) - d_1*d_2*f_SYS)/(d_2*f_SYS - d_1*f_SYS - 2*d_3*f_SYS + d_3^2 + f_SYS^2) -((d_1*d_2*f_SYS*(4*d_3*f_SYS^2 - 8*d_3^2*f_SYS + 4*d_3^3 + d_1*d_2*f_SYS - 4*d_1*d_3*f_SYS + 4*d_2*d_3*f_SYS))^(1/2) + d_1*d_2*f_SYS)/(d_2*f_SYS - d_1*f_SYS - 2*d_3*f_SYS + d_3^2 + f_SYS^2) 第一个解是正确的；但第二个解是错误的，代入原方程不成立！ 奇怪的是，如果把方程中d_1 d_2 d_3 f_SYS赋值，比如：d_1=-1700;d_2=1700;d_3=-1700;f_SYS=-5000; 则方程eq='((17000000*r_2^4 - 28900000000*r_2^3 - 289000000000000*r_2^2 + 835210000000000000000)^(1/2) - 28900000000)/(5000*r_2^2) + 33/50=0' 这时方程仅计算出一个解： -2.2907e+03 不知道出现上述情况是什么原因，大家帮忙看看。 |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
matlab 用solve解决不了的参数方程 怎么办？求大牛们指点指点哈
已经有4人回复
- solve解超越方程少了一个解
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人