| 查看: 234 | 回复: 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 ] |
回复此楼» 猜你喜欢
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有176人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有38人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有1人回复
-
基金申请
已经有44人回复
2楼2008-09-16 14:32:44