nbxunhuan.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
ww=[150;300;500;750;1000;1200;1400];
II=[2.5;5.6;9.2;13;16.5;19.2;22];
www=(150:5:1400);
II1=spline(ww,II,www);
for iu=1:251
w=2*pi*www(iu);
[xdx1,h]=pim2x(w,iu,II1);
x=(xdx(1,20000:40000)-mean(xdx(1,20000:40000)))*180/pi;

niubx(iu,:)=x;
save('niubai','niubx','h');
end

for ip=1:251
    x=(niubx(ip,:)-mean(niubx(ip,:)))*180/pi;
step=h;
sf=1/step;
ng=length(x);
nfft=2^11;
yy=fft(x,nfft);
pyy=yy.*conj(yy)/nfft;
yyy(ip,:)=pyy;
end
yyy=yyy';
myyy=max(yyy);
plot(www,myyy);