求解非线性齐次方程组（钟一谔--转子动力学）

ME!

下面是我编的程序，结构是错的
a=0.4;
O=250;
w=O;
m=102;
EI=6.136*1e4;
Jd=6.377;
Jp=2*6.377;
e1=0.005*1e-2;
e2=0.003*1e-2;
v1=pi/6;
v2=pi/3;
A =[ (21*EI)/(2*a^3) - m*w^2,                  (3*EI)/(2*a^2),               (9*EI)/(2*a^3),                -(3*EI)/(2*a^2);
          (3*EI)/(2*a^2),           (13*EI)/(2*a) - Jd*w^2 + Jp*O*w,          (3*EI)/(2*a^2),                      -EI/(2*a);
         (9*EI)/(2*a^3),                  (3*EI)/(2*a^2),                (15*EI)/(2*a^3) - m*w^2,                -(9*EI)/(2*a^2);
         -(3*EI)/(2*a^2),                       -EI/(2*a),                    -(9*EI)/(2*a^2),         (7*EI)/(2*a) - Jd*w^2 + Jp*O*w]
f1=m*e1*O^2*exp(i*v1);
f2=m*e2*O^2*exp(i*v2);
b=[f1;0;f2;0];

x=inv(A)*b

玉林

1.这怎么是非线性方程组？
2.你给出的那个解对么？z=7.718*1e-3*exp(-250+0.745i)=1.5143e-111 +1.3967e-111i，如此之小

ME!

玉林 发表于 2013-6-21 10:32

                               
登录/注册后可看大图

1.这怎么是非线性方程组？
2.你给出的那个解对么？z=7.718*1e-3*exp(-250+0.745i)=1.5143e-111 +1.3967e-1 ...

答案应该是对的z=7.718*1e-3*exp（i*(-250t+0.745)）

ME!

玉林 发表于 2013-6-21 10:32

                               
登录/注册后可看大图

1.这怎么是非线性方程组？
2.你给出的那个解对么？z=7.718*1e-3*exp(-250+0.745i)=1.5143e-111 +1.3967e-1 ...

不知道怎么有两个帖子，可能是我按错了

http://forum.chinavib.com/thread-127256-1-1.html

你看看是不是应该这样http://forum.chinavib.com/thread-127256-1-1.html