| 查看: 454 | 回复: 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 不知道出现上述情况是什么原因，大家帮忙看看。 |
回复此楼» 猜你喜欢
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有8人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
关于Filecode分析方法
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
找到一些相关的精华帖子，希望有用哦~
matlab 用solve解决不了的参数方程 怎么办？求大牛们指点指点哈
已经有4人回复
- solve解超越方程少了一个解
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
