zhenfupu.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
load 'sjjl';
step=h;
y1=v(2,:);
y2=v(3,:);
%step=0.001
sf=1/step;
%t=0:step:4;
%y1=sin(2*pi*10*t);
%y2=sin(2*pi*5*t);
ng=length(y1);
nfft=2^nextpow2(ng);
ff=(0:sf/nfft:sf/2-sf/nfft)*2*pi;
nn=1:length(ff);
%nfft=1024;
yy1=fft(y1,nfft);
yy2=fft(y2,nfft);
absyy1(nn)=(abs(real(yy1(nn)).^2+imag(yy1(nn))).^2).^0.5;
absyy2(nn)=(abs(real(yy2(nn)).^2+imag(yy2(nn))).^2).^0.5;
figure(1);
plot(ff(nn),absyy1(nn));
figure(2);
plot(ff(nn),absyy2(nn));
figure(3);
plot(ff(nn),real(yy1(nn)));
figure(4);
plot(ff(nn),imag(yy1(nn)));
figure(5);
plot(ff(nn),real(yy2(nn)));
figure(6);
plot(ff(nn),imag(yy2(nn)));
figure(7);
plot(real(yy1(nn)),imag(yy1(nn)));
figure(8);
plot(real(yy2(nn)),imag(yy2(nn)));