压缩机动力学分析

m=0.7;%动子质量

c=0.05;%阻尼系数

k=44;%弹簧弹性系数

f=40;%电流频率

A=pi/2;%电流相位

x1=dsolve('m*D2x+c*Dx+k*x=100*sin(2*pi*f*t+A)','x(0)=0','Dx(0)=0','t'); %微分方程

t=0:0.001:0.2;

y1=100*sin(2*pi*f*t+A);

plot(t,y1);%电磁力图

x1_1=eval(x1);

plot(t,x1_1);%位移图

x2=diff(x1);

x2_1=eval(x2);

plot(t,x2_1);%速度图

x3=diff(x2);

x3_1=eval(x3);

plot(t,x3_1);%加速度图

电磁力：

位移：

速度：

加速度：

稳定状态下的相位关系

在t=[200*pi，200*pi+0.2]

红色是加速度；

紫红色是电磁力；

绿色是速度；

蓝色是位移；

由于数量级上的差距，所以乘上了系数。实际的速度和位移是：

m=0.7;

c=0.005;

k=44;

f=40;

A=pi/2;

x1=dsolve('m*D2x+c*Dx+k*x=100*sin(2*pi*f*t+A)','x(0)=0','Dx(0)=0','t'); t=200*pi:0.001:200*pi+0.2;

x1_1=10000*eval(x1);

plot(t,x1_1,'b');

hold on;

x2=diff(x1);

x2_1=100*eval(x2);

plot(t,x2_1,'g');

hold on;

y1=100*sin(2*pi*f*t+A); plot(t,y1,'m');

hold on;

x3=diff(x2);

x3_1=eval(x3);

plot(t,x3_1,'r');

hold on;