zhmxbc.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
a1=pi/4;
a2=pi;
a3=5*pi/4;
h=0.0001;
tt=3;
w=1000;
bt=rem(tt,h);
bbt=tt-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<a1&oo>=2*k*pi

        ff(:,i)=0;
    elseif (oo-2*k*pi)<a2&oo-2*k*pi>=a1

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

        ff(:,i)=0;
    elseif oo-2*k*pi<2*pi&oo-2*k*pi>=a3

        ff(:,i)=0.5*(1-cos((0.5*(oo-3*pi/2+a1)/a1)*pi));

       end
     if oo-2*k*pi>=2*pi
        k=k+1;
    end
 
end
plot(t(:),ff(:));