ME! 发表于 2013-6-20 19:15

求解非线性齐次方程组


如图所示,也就是转子动力学--钟一谔书里面的那个例子
下面是我编的,解的答案是错误的
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=;

x=inv(A)*b

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

本帖最后由 米斯兰达 于 2013-6-20 22:47 编辑

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

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

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

ME! 发表于 2013-6-21 11:15

本帖最后由 ME! 于 2013-6-21 11:44 编辑

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

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

ME! 发表于 2013-6-21 15:14

本帖最后由 ME! 于 2013-6-21 15:39 编辑

米斯兰达 发表于 2013-6-20 22:46 static/image/common/back.gif
检查了楼主的程序 语法上应该是没有错误的 尝试运行了一下 得到结果 但是似乎与楼主提供的答案对不上(楼主 ...
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=;
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    %x1theta_y1 x2theta_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=;
x2=inv(B2)*b2       %y1theta_x1 y2theta_x2
是不是应该这样写,但是得出的数据x1,x2还只是数据要转换成Z,和fai
但是我不知道怎么变换成指数形式,请教下就是这个公式的变化欧拉公式e^jx=cosx+jsinx
页: [1]
查看完整版本: 求解非线性齐次方程组