AMD.rar_amd_故障诊断_机械_目标识别_识别

N=10000; fs=12000; t=0:1/fs:(N-1)/fs; f=(0:1/N:1/2-1/N)*fs; imf=emd(aaa); plot_emd_result(imf); Y1=fft(imf(1,1:10000)); Y2=fft(imf(2,1:10000)); Y3=fft(imf(3,1:10000)); Y4=fft(imf(4,1:10000)); Y5=fft(imf(5,1:10000)); Y6=fft(imf(6,1:10000)); Y7=fft(imf(7,1:10000)); Y8=fft(imf(8,1:10000)); Y9=fft(imf(9,1:10000)); Y10=fft(imf(10,1:10000)); Y11=fft(imf(11,1:10000)); Y10=fft(imf(10,1:10000)); Y11=fft(imf(11,1:10000)); Y12=fft(imf(12,1:10000)); Y13=fft(imf(13,1:10000)); mag1=abs(Y1); mag2=abs(Y2); mag3=abs(Y3); mag4=abs(Y4); mag5=abs(Y5); mag6=abs(Y6); mag7=abs(Y7); mag8=abs(Y8); mag9=abs(Y9); mag10=abs(Y10); mag11=abs(Y11); mag12=abs(Y12); mag13=abs(Y13); figure(2) subplot(911) plot(f,mag1(1:N/2)*2/N) subplot(912) plot(f,mag2(1:N/2)*2/N) subplot(913) plot(f,mag3(1:N/2)*2/N) subplot(914) plot(f,mag4(1:N/2)*2/N) subplot(915) plot(f,mag5(1:N/2)*2/N) subplot(916) plot(f,mag6(1:N/2)*2/N) subplot(917) plot(f,mag7(1:N/2)*2/N) subplot(918) plot(f,mag8(1:N/2)*2/N) subplot(919) plot(f,mag9(1:N/2)*2/N) figure(3) subplot(411) plot(f,mag10(1:N/2)*2/N) subplot(412) plot(f,mag11(1:N/2)*2/N) subplot(413) plot(f,mag12(1:N/2)*2/N) subplot(414) %plot(f,mag13(1:N/2)*2/N)