求解非线性齐次方程组

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

米斯兰达

检查了楼主的程序 语法上应该是没有错误的 尝试运行了一下 得到结果 但是似乎与楼主提供的答案对不上（楼主提供的答案显示不太清楚啊）

我也不知道怎么回事了 期待高手解决

附上一张图 显示的A矩阵更清楚一些

ME!

补充：  里面有几个参数为复数，公式打不出来

下面这张图片清楚一点貌似没有那么简单，要非开写成8个方程，实部和实部相等，虚部和虚部相等

ME!

米斯兰达 发表于 2013-6-20 22:46

                               
登录/注册后可看大图

检查了楼主的程序 语法上应该是没有错误的 尝试运行了一下 得到结果 但是似乎与楼主提供的答案对不上（楼主 ...

a=0.4;
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;
O=250;
w=O;

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*(cos(v1));
f2=m*e2*O^2*(cos(v2));
b1=[f1;0;f2;0];
B1=[A(1,1) A(1,2) A(1,3) A(1,4);
    A(2,1) A(2,2) A(2,3) A(2,4);
    A(3,1) A(3,2) A(3,3) A(3,4);
    A(4,1) A(4,2) A(4,3) A(4,4)];

x1=inv(B1)*b1    %x1  theta_y1 x2  theta_y2



B2=[A(1,1) -A(1,2) A(1,3) -A(1,4);
    A(2,1) -A(2,2) A(2,3) -A(2,4);
    A(3,1) -A(3,2) A(3,3) -A(3,4);
    A(4,1) -A(4,2) A(4,3) -A(4,4)];
f3=m*e1*O^2*(sin(v1));
f4=m*e2*O^2*(sin(v2));
b2=[f3;0;f4;0];
x2=inv(B2)*b2       %y1  theta_x1 y2  theta_x2
是不是应该这样写，但是得出的数据x1,x2还只是数据要转换成Z，和fai
但是我不知道怎么变换成指数形式，请教下就是这个公式的变化欧拉公式e^jx=cosx+jsinx