基于MATLAB实现的对机械振动信号用三维能量谱进行分析+使用说明文档.rar

N=1024; fs=64; T=1/fs; tt=(0:N-1)*T; s1=sin(2*pi*2*tt)+5*sin(2*pi*5*tt); imf=emd(s1); [m,n]=size(imf); for i=1:m subplot(m,1,i); plot(imf(i,:)); end figure(3) dt=1/fs; h = nspab(imf',500,0,500,dt); subplot(211) surf(h(1:200,50:end-20)) shading interp xlabel('时间(点数)','FontSize',12); Ylabel('频率(Hz)','FontSize',12); zlabel('幅值','FontSize',12); title('a','FontSize',12) view([-75,25]) 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('幅值','FontSize',12); title('c','FontSize',12)