xinmoxing.m

内燃机转子仿真

clear;
clc;
close all;
ts=1;
h=0.0001
ddt=rem(ts,h);
dt=ts-ddt;
nn=dt/h;
T=0.1;
tt=(0:h:T);
ddi=rem(T,h);
di=T-ddi;
np=di/h;
w=600;
for ii=1:nn
    t(:,ii)=(ii-1)*h;
    i=rem(ii,np);
    if i>=0&i<=0
    ff7(:,ii)=1;
    elseif i>=0&i<=400
    ff7(ii)=0.5*(1+cos(2*pi*100*(tt(i)-tt(1))/8));
elseif i>=400&i<=600
    ff7(:,ii)=0;
elseif i>=600&i<=1000
    ff7(ii)=0.5*(1-cos(2*pi*100*(tt(i)-tt(600))/8));
elseif i>=1000&i<=1000
    ff7(:,ii)=1;
    end
end
save('zijidekaiguan7','t','ff7');
plot(t,ff7,'k');
ylim([-0.1 1.1]);
xlim([0 0.6]);
xlabel('时间(s)');
ylabel('幅值');
title('横向裂纹开关函数模型');