杨2009 发表于 2009-1-8 20:35

帮我看下用matlab求解非线性方程时哪里出了问题

程序如下:其中s1,s2,s3,s4,s5,s6,s7,s8是方程;k2 y1 y2 m1 m2 x1 x2为相应的7个变量,但是求解时,有如下显示,没求出解,需要高手指导,谢谢!!!
syms k2 y1 y2 m1 m2 x1 x2;
>> s1=k2+8.5847142026029796983461767280000;
>> s3=y1-(-2.-2.040*k2)/(-.2996e-1-.200e-3*k2);
>> s4=y2-(1+.20e-1*k2)/(-.2996e-1-.200e-3*k2);
s5=m1-(-.6e-1*(-2+.10e-1*y1)*(1+.10e-1*y2)+1.960+.1920e-1*y2-.200e-3*y1+.40e-3*x1*(1+.10e-1*y2))/(.6e-1*(1+.10e-1*y2)^2+.40e-1+.400e-3*y2);
s6=m2-((-.6e-1-.120e-2*k2)*(1+.10e-1*y2)-.20e-1-.400e-3*k2+.80e-3*x2*(1+.10e-1*y2))/(.6e-1*(1+.10e-1*y2)^2+.40e-1+.400e-3*y2);
s7=x1-2.5*(-1020000.+6.*k2*m1*y2+3.*m2*y2*y1+800.*m2*m1*y2+3.*m2*y2^2*m1-20400.*k2+40000.*m1-1040000.*m2+300.*y1+40800.*x2-10200.*m2*y2+300.*m1*y2+800.*k2*m1+6.*k2*y1+50000.*m2*m1+400.*m2*y1-400.*x2*m1-4.*x2*y1)/(5.*m2*y2+1498.+10.*k2+500.*m2);
s8=x2-.625*(30000.+1200.*k2+80000.*m2+600.*m2*y2+12.*k2^2+1600.*k2*m2+12.*k2*m2*y2+50000.*m2^2+800.*m2^2*y2+3.*m2^2*y2^2)/(1499.+10.*k2+500.*m2+5.*m2*y2);
>> solve(s1,s3,s4,s5,s6,s7,s8,' k2 y1 y2 m1 m2 x1 x2')

ans =
    k2:
    m1:
    m2:
    x1:
    x2:
    y1:
    y2:
页: [1]
查看完整版本: 帮我看下用matlab求解非线性方程时哪里出了问题