zfp1.m

内燃机转子仿真

%function mpyy1=gonglp1()
clear all;
clc;
close all;
format long;
load 'bneiranji4';
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:Ts-step;
ni=length(t)-1;
% w=pw(6)  ;
% pimx2(w);
%load 'sjjl';
%subplot(2,2,1);
%plot(t,xyz1(2,:));
%title('m1振动时域波形图');
%xlabel('时间(s)');
%ylabel('振幅');
%xlim([0,Ts]);
 %subplot(2,2,2);
%plot(t,xyz1(6,:));
%title('m2振动时域波形图');
%xlabel('时间(s)');
%ylabel('振幅');
%xlim([0,Ts]);

y1=xyz4(2,:);
y2=xyz4(6,:);
y3=xyz4(9,:);
sf=1/step;
%y1=sin(2*pi*10*t);
%y2=sin(2*pi*5*t);
ng=length(y1);
nfft=2^nextpow2(ng);
%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;
pyy1=yy1.*conj(yy1)/nfft;
yy2=fft(y2,nfft);
pyy2=yy2.*conj(yy2)/nfft;
yy3=fft(y3,nfft);
%pyy3=(yy3.*conj(yy3)).^0.5;
pyy3=yy3.*conj(yy3)/nfft;
%ff=2*pi*(0:sf/nfft:sf/2-sf/nfft);%%%%%%%%%%%%%%%%%%%%%%%%
ff=(0:sf/nfft:sf/2-sf/nfft);
nn=1:length(ff);
subplot(3,1,1);
plot(ff,pyy1(nn),'k');
title('曲轴1扭振幅值谱');
xlabel('频率(Hz)');
ylabel('|FFT|');
xlim([0,200]);
%xlim('auto');
ylim('auto');
subplot(3,1,2);
plot(ff,pyy2(nn),'k');
title('曲轴5扭振幅值谱');
xlabel('频率(Hz)');
ylabel('|FFT|');
xlim([0,200]);
%xlim('auto');
ylim('auto');
subplot(3,1,3);
plot(ff,pyy3(nn),'k');
title('曲轴9扭振幅值谱');
xlabel('频率(Hz)');
ylabel('|FFT|');
xlim([0,200]);
%xlim('auto');
ylim('auto');

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