新人求教：课题基于HHT的汽轮机碰磨振动故障诊断
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);
=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
因为新人，不能发图，有、运行出来有5张图，需要解释，但是我不太能解释清楚，请问有大哥大姐帮忙稍微教下吗？ 可以加我Q 573631136 师哥师姐帮帮忙。。。{:3_50:} 师哥师姐帮帮忙。。。{:3_50:} 我水平有限。但是你应该先弄一个很简单的模拟信号，然后调试各个程序，这样才能看出有没有问题。一上来就拿实际采集的信号分析，很难知道程序是否正确。 如果还有问题，这几天赶快回帖，我给你详细说。
页:
[1]