zzcjingque.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
w=60;
h=0.01;
ts=12;
dt=rem(ts,h);
ddt=ts-dt;
nn=ddt/h;
xy=zeros(3,nn);
xy(:,1)=[0;0;0]
for i=2:nn
    t=(i-1)*h;
    x1=(10/(6-3*w^2)-10/(15-3*w^2))*cos(w*t)-10/(6-3*w^2)*cos(2^0.5*t)+10/(15-3*w^2)*cos(5^0.5*t);
    x2=(10/(6-3*w^2)+20/(15-3*w^2))*cos(w*t)-10/(6-3*w^2)*cos(2^0.5*t)-20/(15-3*w^2)*cos(5^0.5*t);
    xy(:,i)=[t;x1;x2];
end
figure(1);
plot(xy(1,:),xy(2,:));
xlabel('时间(s)');
ylabel('x1');
title('质量1的响应');
figure(2);
plot(xy(1,:),xy(3,:));
xlabel('时间(s)');
ylabel('x2');
title('质量2的响应');
save('jingque1','xy');