sjjlf.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
load 'sjjl';
%v is input data
x=v(3,:);
sf=1024;
np=100000;
t=0:1/sf:(np-1)/sf;
nt=length(x);
s=1.5*std(x);
m=0;
y=zeros(1,np);
w=nt-np;
for k=2:w
    a=abs(x(k-1)-s);
    b=abs(x(k)-s);
    c=abs(x(k+1)-s);
    if b<a&b<c
        y(1:np)=y(1:np)+x(k:k+np-1);
        m=m+1;
    end
end
y=y./m;
save('outsjjl','y');
plot(t,y);
xlabel('time(s)');
ylabel('amplitude');
grid on;