用希爾伯特黃變換(hht)求時頻譜和邊際譜(編輯修改稿)
2025-06-12 06:41 本頁面
【文章內(nèi)容簡介】 oman39。)set(gca,39。fontsize39。,)ylabel([39。imf39。,int2str(i)])endsubplot(m+1,1,m+1)。set(gcf,39。color39。,39。w39。)plot(t,c(m,:),39。k39。)set(gca,39。fontname39。,39。times New Roman39。)set(gca,39。fontsize39。,)ylabel([39。r39。,int2str(m1)])%畫出每個IMF分量及剩余分量residual的幅頻曲線figure(2)subplot(m+1,1,1)set(gcf,39。color39。,39。w39。)[f,z]=fftfenxi(t,z)。plot(f,z,39。k39。)set(gca,39。fontname39。,39。times New Roman39。)set(gca,39。fontsize39。,)ylabel([39。initial signal39。,int2str(m1),39。Amplitude39。])for i=1:m1subplot(m+1,1,i+1)。set(gcf,39。color39。,39。w39。)[f,z]=fftfenxi(t,c(i,:))。plot(f,z,39。k39。)set(gca,39。fontname39。,39。times New Roman39。)set(gca,39。fontsize39。,)ylabel([39。imf39。,int2str(i),39。Amplitude39。])endsubplot(m+1,1,m+1)。set(gcf,39。color39。,39。w39。)[f,z]=fftfenxi(t,c(m,:))。plot(f,z,39。k39。)set(gca,39。fontname39。,39。times New Roman39。)set(gca,39。fontsize39。,)ylabel([39。r39。,int2str(m1),39。Amplitude39。])hx=hilbert(z)。xr=real(hx)。xi=imag(hx)。%計算瞬時振幅sz=sqrt(xr.^2+xi.^2)。%計算瞬時相位sx=angle(hx)。%計算瞬時頻率dt=diff(t)。dx=diff(sx)。sp=dx./dt。figure(6)plot(t(1:N1),sp)title(39。瞬時頻率39。)%計算HHT時頻譜和邊際譜[A,fa,tt]=hhspectrum(c)。[E,tt1]=toimage(A,fa,tt,length(tt))。figure(3)disp_hhs(E,tt1) %二維圖顯示HHT時頻譜，E是求得的HHT譜pausefigure(4)for i=1:size(c,1)faa=fa(i,:)。[FA,TT1
公司管理相關(guān)推薦
鄂ICP備17016276號-1