gonglp3.m

内燃机转子仿真

%function mpyy1=gonglp1()
clear all;
clc;
close all;
format long;
load 'nb';
y1=xdx1(2000:20000)-mean(xdx1(2000:20000));
step=h;

sf=1/step;

%y1=sin(2*pi*10*t);
%y2=sin(2*pi*5*t);
ng=length(y1);
%nfft=2^nextpow2(ng)*32;
nfft=2^10;
%w1=xcorr(y1,z1,200);
%w2=xcorr(y2,z2,200);
%x=-nn*step:step:nn*step;
yy1=fft(y1,nfft);
pyy1=yy1.*conj(yy1)/nfft;

ff=(0:sf/nfft:sf/2-sf/nfft);%%%%%%%%%%%%%%%%%%%%%%%%
nn=1:length(ff);
figure
plot(ff(nn),pyy1(nn),'k');
title('m1功率谱');
xlabel('频率(Hz)');
ylabel('幅值');
%xlim([-1,6]);
xlim('auto');
ylim('auto');