zonghemoxing.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
a=pi/2;
h=0.001;
t=5;
w=2*pi*100;
bt=rem(t,h);
bbt=t-bt;
n1=bbt/h;
    k=0;
for i=1:n1
    t(:,i)=(i-1)*h;


    oo=w*t(:,i)/60;
   
    if oo<=2*k*pi+0.5*pi-a&oo>=2*k*pi

        ff(:,i)=1;
    elseif oo<=2*k*pi+0.5*pi+a&oo>2*k*pi+0.5*pi-a

        ff(:,i)=0.5*(1+cos((0.5*(oo-0.5*pi+a)/a)*pi));
    elseif oo<=2*k*pi+3*pi/2-a&oo>2*k*pi+0.5*pi+a

        ff(:,i)=0;
    elseif oo<=2*k*pi+3*pi/2+a&oo>2*k*pi+3*pi/2-a

        ff(:,i)=0.5*(1-cos((0.5*(oo-3*pi/2+a)/a)*pi));
    %elseif oo>2*k*pi+3*pi/2+a&oo<=2*pi+2*k*pi
    elseif 2*k*pi+3*pi/2+a<oo<=2*pi+2*k*pi
         ff(:,i)=1;
    end
     if oo>2*(k+1)*pi
        k=k+1;
    end
end
save('kaiguan','t','ff');
plot(t(:),ff(:));
ylim([-0.02 1.02]);
xlabel('时间(s)');
ylabel('幅值');
title('方波模型开关函数');