新人求教：课题基于HHT的汽轮机碰磨振动故障诊断

小飞侠1992

RT
老师吧代码给我们了

x=ZT3PM1;
N=length(x);
fs=7147.2;
figure(1);%波形图
plot(x);
y=fft(x);%进行fft变换
mag=abs(y).*2/N;%求幅值
m=length(y);
f=(0:m/2-1).*fs/m;%进行对应的频率转换
figure(2);  %频谱图
plot(f,mag(1:m/2));
xlabel('频率(Hz)');
ylabel('幅值');
title('仿真信号频谱图');

imf=emd(x);
[m,n]=size(imf);
for i=1:m
    figure(3)  % EMD分解IMF分量图
    subplot(m,1,i);
    plot(imf(i,:));   
end
dt=1/fs;

h = nspab(imf',500,0,500,dt);
figure(4)  % 三维HHT图
subplot(211)
surf(h(3:200,50:end-30))
shading interp
xlabel('时间(点数)','FontSize',12);
Ylabel('频率(Hz)','FontSize',12);
zlabel('幅值(mV)','FontSize',12);
title('a','FontSize',12)
view([-75,25])

figure(5)  % 二维HHT图和边际谱图
yt=subplot(223)
imagesc(h(1:200,:))
Ylabel('频率(Hz)','FontSize',12);
xlabel('时间(点数)','FontSize',12);
set(yt,'ydir','nor')
title('b','FontSize',12)
ms=mspc(h)
subplot(224)
plot((1:length(ms)),ms)
xlabel('频率(Hz)','FontSize',12);
ylabel('幅值(mV)','FontSize',12);
title('c','FontSize',12)




请大神帮我看看代码有没有问题
9.066357

9.066357

11.151146

11.151146

11.056383

5.678573

5.678573

9.090047

10.440422

10.440422

7.147402

7.147402

9.824462

12.004015

12.004015

9.800771

9.800771

12.288304

15.533943

16.197287

16.197287

16.126213

16.126213

20.058886

22.901781

22.901781

21.669859

21.669859

27.474104

31.312014

31.312014

25.815748

23.896793

23.896793

29.890566

29.890566

27.924229

22.570110

22.333202

22.333202

23.920485

23.920485

19.561378

15.770851

15.770851

14.728456

12.999028

10.393041

10.227205

9.942916

7.952889

4.162362

3.143658

2.361861

-0.812705

-1.096995

小飞侠1992

因为新人，不能发图，有、运行出来有5张图，需要解释，但是我不太能解释清楚，请问有大哥大姐帮忙稍微教下吗？ 可以加我Q 573631136

小飞侠1992

师哥师姐帮帮忙。。。

小飞侠1992

师哥师姐帮帮忙。。。

airazor

我水平有限。但是你应该先弄一个很简单的模拟信号，然后调试各个程序，这样才能看出有没有问题。一上来就拿实际采集的信号分析，很难知道程序是否正确。

airazor

如果还有问题，这几天赶快回帖，我给你详细说。