| 查看: 1576 | 回复: 2 | ||
[求助]
MATLAB求解多项式形式的非线性方程组（2个方程） 已有1人参与
|
|
求助在MATLAB中求解两个高次多项式组成的非线性方程组，用fsolve命令求解 方程的参数和形式如下，下面的P1和P2为两个多项式，需要求解的变量是P1和P2中的as（1）和as（2） function G=coupled_FR2(as,s) C=38.9636; omega10= 325.0643; omega20= 865.1030; r1=136.9815; F=1.2875*10^4; x11=1; x12= 1; x21=0.6390; x22= -1.5649; P1=(4*x11*x22*C*as(2)*omega20-4*C*as(2)*omega20*x12*x21)^2/(-x11^3*r1*x21*as(1)^3+... x11*r1*x21^3*as(1)^3)^2+((3*s*omega10-omega20)*(8*omega20*x11*x22-8*omega20*x12*x21)+3*r1*x21*x12^3*as(2)^3-... 3*r1*x11*x22^3*as(2)^3+6*r1*x21*x11^2*as(1)^2*x12*as(2)-6*r1*x11*x21^2*as(1)^2*x22*as(2))^2/(-x11^3*r1*x21*as(1)^3+x11*r1*x21^3*as(1)^3)^2 -1; P2=(-x11^2*x22*as(1)^2*omega10*x21+as(1)^2*omega10*x12*x21^2*x11+... 3*x11*x22^2*x12*as(2)^2*omega20-3*x22*x12^2*as(2)^2*omega20*x21)^2*C^2/(as(1)^2*F^2*x22^2*x11^2*x21^2)+... (1/16)*(3*x22*r1*x11^4*as(1)^4*x21-12*as(2)^2*r1*x22*x11^2*as(1)^2*x12^2*x21-... 8*x11^2*(s*omega10-omega10)*omega10*x21*as(1)*x22+24*x11*x22^2*x12*as(2)*(omega20-3*omega10)*omega20+... 8*x11*(s*omega10-omega10)*omega10*x21^2*as(1)*x12-3*r1*x12*x21^4*as(1)^4*x11-... 72*x11*x22^2*x12*as(2)*(s*omega10-omega10)*omega20+9*x11*x22^4*x12*as(2)^4*r1+... 12*as(2)^2*r1*x12*x21^2*as(1)^2*x22^2*x11-24*x22*x12^2*as(2)*(omega20-3*omega10)*omega20*x21-... 9*x22*x12^4*as(2)^4*r1*x21+72*x22*x12^2*as(2)*(s*omega10-omega10)*omega20*x21)^2/(as(1)^2*F^2*x22^2*x11^2*x21^2) -1; G=[P1 P2]; 下面利用fsolve命令求解 因为上述P1和P2中还有一个变量是s 求解命令为 clc clear all s=1.1; % define parameter first as = fsolve(@(as) coupled_FR2(as,s),[13;10^(-2)],optimset('display','iter','MaxFunEvals',20000)) 但是MATLAB最后结果为 No solution found. fsolve stopped because the last step was ineffective. However, the vector of function values is not near zero, as measured by the default value of the function tolerance. <stopping criteria details> as = 14.7209 -1.0705 而且这里我只需要求出as（1）和as（2）的正实数解就可以， 请教论坛中的高手给点意见，怎么能求出这里的as（1）和as（2） |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有124人回复
-
物理学I论文润色/翻译怎么收费?
已经有184人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有0人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2015-12-10 17:44:40
3楼2015-12-10 18:27:10
200