| 查看: 261 | 回复: 1 | |||
| 当前主题已经存档。 | |||
[交流]
[求助]懂MATLAB解非线性方程组得虫虫近来看一下~
|
|||
|
3个方程，3个未知数为x, y ,z K1，K2，K3，a 分别为常数；. 编如下程序求解: clear clc syms x y z a k1 k2 k3; expr1 = '(4*x+y-3*z)^4*(2*x-y-z)^2-k1*(1/(a+1)-x)*(a/(a+1)-x-y+z)*(1+4*x-2*z)^4'; expr2 = '(4*x+y-3*z)*y-k2*(2*x-y-z)*(a/(a+1)-x-y+z)'; expr3 = 'z*(a/(a+1)-x-y+z)*(1+4*x-2*z)^2-k3*(2*x-y-z)*(4*x+y-3*z)^3'; fun1 = subs(expr1,{a,k1,k2,k3},{1,1,1,1}); fun2 = subs(expr2,{a,k1,k2,k3},{1,1,1,1}); fun3 = subs(expr3,{a,k1,k2,k3},{1,1,1,1}); s=solve(fun1,fun2,fun3); result=[s.x,s.y,s.z]; 当a,k1,k2,k3取值如：{1，1，1，1，}可以算出解；当当a,k1,k2,k3取值如{3，21000，12，18255}时，怎么结果如下： Warning: Explicit solution could not be found. > In D:\Program\matlab\toolbox\symbolic\solve.m at line 133 In D:\Program\matlab\toolbox\symbolic\@sym\solve.m at line 49 In D:\Program\matlab\work\zyh.m at line 10 ??? Access to an object's fields is only permitted within its methods. Error in ==> D:\Program\matlab\work\zyh.m On line 11 ==> result=[s.x,s.y,s.z] 是不是意味着没有解呢？而且解一个要30分钟（CPU 是T7300 ），太慢了。 可是按理应该有解啊？用fsolve 程序是不是可以求解（不太会用 ），而且更快？？能帮忙解决一下吗？20个金币相谢！急！ [ Last edited by kele888888 on 2008-9-16 at 13:28 ] |
回复此楼» 猜你喜欢
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
物理学I论文润色/翻译怎么收费?
已经有194人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有2人回复
2楼2008-09-16 14:32:44