lipingheng.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
load 'sjjl';
x=v(2:3,:);
nn=length(x(2,:));
y=zeros(2,nn);
f=zeros(2,nn);
f(1,:)=0;
f(2,:)=F0(2,:);
h=0.01;
dt=h/2^20;
for i=1:nn
    %H is change along with time
    kt=K-door1(i)*dk1;
    H=[-c/2/m 1/m;0.25*c^2/m-kt -0.5*c/m];
    y(:,i)=H*x(:,i)+f(:,i);
end
xx(1,:)=x(1,:);
xx(2,:)=(x(2,:)-0.5*c*x(1,:))/m;
yy(1,:)=y(1,:);
yy(2,:)=(y(2,:)-0.5*c*y(1,:))/m;
wy=xx(1,:);
sd=xx(2,:);
jsd=yy(2,:);
for i=1:nn
    kt=K-door1(i)*dk1;
bph(i)=f(2,i)-m*jsd(i)-c*sd(i)-kt*wy(i);
end
figure(1)
plot(v(1,:),xx(1,:),'k',v(1,:),xx(2,:),'-.k',v(1,:),yy(2,:),'--k');
ylim([-0.5,0.7]);
h=legend('位移','速度','加速度',1);
xlabel('时间(s)');
ylabel('幅值');
title('步长为1位移，速度和加速度响应图');
set(h,'box','on');
figure(2);
plot(v(1,1:2:100),bph(1:2:100),'k');
xlabel('时间(s)');
ylabel('不平衡幅值');
title('步长为1系统不平衡时域图');