24小时热门版块排行榜

返回列表

寒竹冷霜

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 82.5
帖子: 25
在线: 98.2小时
虫号: 1823743
注册: 2012-05-18
专业: 动力学与控制

[求助] 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）

回复此楼