lnly 发表于 2007-9-12 10:47

波形图分析

下面是一个画转子非线性振动波形图的程序,但不知道怎么回事画出的x方向的图好象不对?这是说明系统失稳了吗?
请大家指教! 谢谢!
function xdot=zhuanzi1(t,x,flag,omega,e)
A0=1.2e-4;N=80;C0=2e-4;I0=4;miu0=4*pi*(10.^-7);m=0.434;g=10;
Kxi=96.5097;Kyi=96.5097;Kxx=1.9302e+6;Kyy=1.9311e+6;Kp=20.5;Kd=0.006;
i0=m*g*C0^2/miu0*A0*N^2*I0;Kiy1=-2*miu0*A0*N^2*i0/C0^3;Kyy1=3*miu0*A0*N^2*I0*i0/C0^4;
C=-Kyi*Kd;
Kx1=Kxi*Kd*sqrt(C*g)/m*g;Kx2=C*(Kxi*Kp+Kxx)/m*g;Ky1=Kyi*Kd*sqrt(C*g)/m*g;Ky2=C*Kiy1*Kd*sqrt(C*g)/m*g;
Ky3=C*(Kyi*Kp+Kyy)/m*g;Ky4=C^2*(Kyy1+Kiy1*Kp)/m*g;tao=3230.447*t;
xdot=[x(3);
    x(4);
   -Kx1*(omega^-1)*x(3)-Kx2*(omega^-2)*x(1)+e*cos(tao);
    -Ky1*(omega^-1)*x(4)-Ky2* (omega^-1)*x(2)*x(4)-Ky3*(omega^-2)*x(2)-Ky4*(omega^-2)*x(2)^2+e*sin(tao)];
%%%%%%%%%%%%%%%%%%%%%%%%%
clear;
clc;
syms omega
t_final=0.6;x0=;
omega=14.4471;e=0.1;
=ode45('zhuanzi1',,x0,[],omega,e);
figure;plot(t,x(:,1))% x向波形图


咕噜噜 发表于 2007-9-12 10:50

回复 #1 lnly 的帖子

把图贴出来吧，我看看，我刚重装系统还没来的及装matlab那，你把图贴出来也方便大家帮你分析

lnly 发表于 2007-9-12 11:05

不好意思刚才图没传上,附件中是画的图.




无水1324 发表于 2007-9-12 11:10

回复 #3 lnly 的帖子

附件还是没有看到

lnly 发表于 2007-9-12 11:41

不知道怎么回事图传不上去了！不好意思！

无水1324 发表于 2007-9-12 11:43

回复 #5 lnly 的帖子

看一下图片的格式和大小。论坛有一定的限制

octopussheng 发表于 2007-9-12 11:48

回复 #6 无水1324 的帖子

我正在给她试算，算一次时间还蛮长的！

octopussheng 发表于 2007-9-12 15:51

不知道程序有没有问题，算了5个小时还没出来结果！晕死！

咕噜噜 发表于 2007-9-12 19:24

回复 #8 octopussheng 的帖子

所以楼主还是把自己的图贴出来吧

无水1324 发表于 2007-9-12 19:35

syms omega

这句话应该是多余的，

lnly 发表于 2007-9-13 08:53

波形图1

咕噜噜 发表于 2007-9-13 09:19

回复 #11 lnly 的帖子

楼主是不是不会上传附件啊？

octopussheng 发表于 2007-9-14 08:09

回复 #12 咕噜噜 的帖子

算了很久都没有结果，我已经放弃了！呵呵！楼主把图重新贴一下吧！
页: [1]
查看完整版本: 波形图分析