| 查看: 399 | 回复: 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 不知道出现上述情况是什么原因，大家帮忙看看。 |
回复此楼» 猜你喜欢
遇见不省心的家人很难过
已经有17人回复
-
退学或坚持读
已经有25人回复
-
博士延得我，科研能力直往上蹿
已经有4人回复
-
免疫学博士有名额，速联系
已经有14人回复
-
面上基金申报没有其他的参与者成吗
已经有4人回复
-
多组分精馏求助
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
matlab 用solve解决不了的参数方程 怎么办？求大牛们指点指点哈
已经有4人回复
- solve解超越方程少了一个解
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人