刚学matlab ，用它求解一个非线性方程组，结果提示如下错误： fsolve stopped because it exceeded the functi on evaluation limit,function f=group(t) x(1)=t(1); %x(1)=Ft3% x(2)=t(2); %x(2)=FAa% x(3)=t(3); %x(3)=FBa% x(4)=t(4); %x(4)=FCa% x(5)=t(5); %x(5)=FDa% x(6)=t(6); %x(6)=FAr% x(7)=t(7); %x(7)=FBr% x(8)=t(8); %x(8)=FCr% x(9)=t(9); %x(9)=FDr% x(10)=t(10); %x(10)=aA% x(11)=t(11); %x(11)=aB% x(12)=t(12); %x(12)=aC% x(13)=t(13); %x(13)=aD% x(14)=t(14); %x(14)=Fa3% x(15)=t(15); %x(15)=Fr3% a1=0.058; %a1=LA% a2=0.058; %a2=LB% a3=0.111; %a3=LC% a4=0.111; %a4=LD% a5=0.036; %a5=LA'% a6=0.036; %a6=LB'% a7=0.027; %a7=L1% a8=0.027; %a8=L3% b1=0.033; %b1=R% b2=0.102; %b2=R1% b3=0.126; %b3=R3% c1=40.352; %c1=TH% d1=0.15; %d1=Jy% d2=0.311; %d2=Jz% d3=0.67; %d3=Ya% d4=0.67; %d4=Yc% d5=74.351; %d5=WH% d6=20*pi/180; %d6=β% d7=20*pi/180; %d7=α% d8=80*pi/180; %d8=δ3% d9=116*pi/180; %d9=δ1% c4=c1/b2*cos(d7); %c4=Ft1% c2=c4*tan(d7)*sin(d9); %c2=Fa1% c3=c4*tan(d7)*cos(d9); %c3=Fr1% f=zeros(15,1); f(1)=x(7)*sin(x(11))*a2-x(6)*sin(x(10))*a1-x(9)*sin(x(13))*a4-x(8)*sin(x(12))*a3; f(2)=-x(8)*cos(x(12))*a3-x(9)*cos(x(10))*a4-x(6)*cos(x(10))*a1*cos(d6)-x(7)*cos(x(11))*a2*cos(d6); f(3)=-c1-x(6)*cos(x(10))*a1*sin(d6)+x(7)*cos(x(11))*a2*sin(d6); f(4)=-x(6)*cos(x(10))-x(8)*cos(x(12))-x(7)*cos(x(11))-x(9)*cos(x(13)); f(5)=-x(6)*sin(x(10))*cos(d6)-x(7)*sin(x(11))*cos(d6)-x(9)*sin(x(13))-x(8)*sin(x(12))+x(3)*sin(d6)-x(2)*sin(d6); f(6)=x(4)-x(5)+x(3)*cos(d6)-x(2)*cos(d6)+x(6)*sin(x(10))*sin(d6)-x(7)*sin(x(11))*sin(d6); f(7)=c2*b2+x(14)*b3+c3*a7+x(15)*a8+x(2)*b1*sin(x(10))-x(3)*b1*sin(x(11))+x(6)*a5*sin(x(10))-x(7)*a6*sin(x(11))-sin(d6)*((1-cos(d6))*d2+d1*cos(d6))*d5^2; f(8)=-x(6)*cos(x(10))*a5+x(7)*cos(x(11))*a6-x(2)*cos(x(10))*b1+x(3)*cos(x(11))*b1-c4*a7+x(1)*a8; f(9)=-c4*b2+x(1)*b3; f(10)=x(2)-x(6)/(2*d3); f(11)=x(3)-x(6)/(2*d3)-c2-x(14); f(12)=x(4)-x(8)/(2*d4); f(13)=x(5)-x(8)/(2*d4)-(x(3)-x(2))*cos(d6)-(x(7)*sin(x(11))+x(6)*sin(x(10)))*sin(d6); f(14)=x(14)-x(1)*tan(d7)*sin(d8); f(15)=x(15)-x(1)*tan(d7)*cos(d8);
更多回帖
