nbxunhuan.asv

内燃机转子仿真

clear all;
clc;
close all;
format long;
for i=1:30
w=(i-1)*3;
w=(i-1)*3;
for i=1:30
w=2*pi*(i-1)*3;
pim2x(w);
x=(xdx(1,10000:40000)-mean(xdx(1,10000:40000)))*180/pi;
step=h;
sf=1/step;
ng=length(x);
nfft=2^11;
yy=fft(x,nfft);
pyy=yy.*conj(yy)/nfft;
ff=(0:sf/nfft:sf/2-sf/nfft);
nn=length(ff);
for j=1:nn
if ff(j-1)<=(i-1)*3
if ff(j)>=(i-1)*3
pm=max[pyy(j-1);pyy(j-1)];
end
end
pmm(:,i)=pm;
end
plot