fuliyebian.m

内燃机转子仿真

%ti xing bo fu li ye zhan kai

clear;
clc;
close all;
format long;
AA=42975;
zhuansu=1000;
d=60.0/zhuansu/360*3.0; 
w=2*pi*zhuansu/60;
a=63
k=(1:2:91);
n1=length(k);
t=(0:0.001:6.5);
tx1=42975*(0:0.)
n2=length(t);
xx=zeros(n1,n2);
x=zeros(1,n2);
y=zeros(1,n2);
ctt=(1:2:91);

for k=1:1:n1
    xx(k,:)=xx(k,:)+sin((2*k-1)*(0.5*pi-a*pi/180))*sin((2*k-1)*(1.0005755-6.5*0.5/2/pi)*t)/(2*k-1)^2;
%%1.0005755-6.5*0.5/2/pi= 0.48332193495134
%k=(1:2:91)*0.48332193495134
end
x(:)=sum(xx);
y(:)=(4*AA/(0.5*pi-a*pi/180)/pi)*x(:);
plot(t,y);
ctt1=ctt';
for i=1:n2
    it=(i-1)*0.001;
    t(1,i)=it;
    fffx(1,i)=(4*AA/(0.5*pi-a*pi/180)/pi)*sum(sin(ctt1*(0.5*pi-a*pi/180)).*sin(ctt1*0.48332193495134*it)./ctt1.^2);
end
figure(2);
plot(t,fffx);