zfp2.m

内燃机转子仿真

%function mpyy1=gonglp1()
clear all;
clc;
close all;
format long;
%load 'neiranji1';
load 'sjjl';
pw=[0 2^0.5 2^0.5+(5^0.5-2^0.5)/2 5^0.5 4 6];
step=h;
t=0:step:tf-step;
ni=length(t)-1;
% w=pw(6)  ;
% pimx2(w);

subplot(2,2,1);
plot(t,v(2,:));
title('m1振动时域波形图');
xlabel('时间(s)');
ylabel('振幅');
xlim([0,12]);
 subplot(2,2,2);
plot(t,v(3,:));
title('m2振动时域波形图');
xlabel('时间(s)');
ylabel('振幅');
xlim([0,12]);

y1=v(2,:);
y2=v(3,:);
sf=1/step;
%y1=sin(2*pi*10*t);
%y2=sin(2*pi*5*t);
ng=length(y1);
nfft=2^nextpow2(ng)*32;
%nfft=512;
%w1=xcorr(y1,z1,200);
%w2=xcorr(y2,z2,200);
%x=-20000*step:step:20000*step;
yy1=fft(y1,nfft);
pyy1=(yy1.*conj(yy1)).^0.5;
yy2=fft(y2,nfft);
pyy2=(yy2.*conj(yy2)).^0.5;
ff=2*pi*(0:sf/nfft:sf/2-sf/nfft);%%%%%%%%%%%%%%%%%%%%%%%%
nn=1:length(ff);
subplot(2,2,3);
plot(ff(nn),pyy1(nn));
title('m1幅值谱');
xlabel('频率(Hz)');
ylabel('|FFT|');
xlim([-1,12]);
%xlim('auto');
ylim('auto');
subplot(2,2,4);
plot(ff(nn),pyy2(nn));
title('m2幅值谱');
xlabel('频率(Hz)');
ylabel('|FFT|');
xlim([-1,12]);
%xlim('auto');
ylim('auto');

%figure(3);
%plot(x,w1);
%figure(4);
%plot(x,w2);
%mpyy1=max(pyy1);