求解非线性齐次方程组
如图所示，也就是转子动力学--钟一谔书里面的那个例子
下面是我编的，解的答案是错误的
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:47 编辑
检查了楼主的程序 语法上应该是没有错误的 尝试运行了一下 得到结果 但是似乎与楼主提供的答案对不上（楼主提供的答案显示不太清楚啊）
我也不知道怎么回事了 期待高手解决
附上一张图 显示的A矩阵更清楚一些 本帖最后由 ME! 于 2013-6-21 11:44 编辑
补充：里面有几个参数为复数，公式打不出来
下面这张图片清楚一点貌似没有那么简单，要非开写成8个方程，实部和实部相等，虚部和虚部相等
本帖最后由 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]