用希爾伯特黃變換(hht)求時頻譜和邊際譜(已修改)
2025-05-28 06:41 本頁面
【正文】 . . . .用希爾伯特黃變換（HHT）求時頻譜和邊際譜？HHT就是先將信號進行經(jīng)驗模態(tài)分解（EMD分解），然后將分解后的每個IMF分量進行Hilbert變換，得到信號的時頻屬性的一種時頻分析方法。 ( KB)201025 18:30 ( KB)201025 18:30 ( KB)201025 18:30EMD分解的流程圖如下： ( KB)201025 18:49。給定頻率分別為10Hz和35Hz的兩個正弦信號相疊加的復合信號，采樣頻率fs=2048Hz的信號，表達式如下：y=5sin(2*pi*10t)+5*sin(2*pi*35t)(1)為了對比，先用fft對求上述信號的幅頻和相頻曲線。復制內(nèi)容到剪貼板代碼:function fftfenxiclear。clc。N=2048。%fft默認計算的信號是從0開始的t=linspace(1,2,N)。deta=t(2)t(1)。1/detax=5*sin(2*pi*10*t)+5*sin(2*pi*35*t)。% N1=256。N2=512。w1=*2*pi。w2=*2*pi。w3=*2*pi。% x=(t=200amp。t=200+N1*deta).*sin(w1*t)+(t200+N1*detaamp。t=200+N2*deta).*sin(w2*t)+(t200+N2*detaamp。t=200).*sin(w3*t)。y = x。m=0:N1。f=1./(N*deta)*m。%可以查看課本就是這樣定義橫坐標頻率范圍的%下面計算的Y就是x(t)的傅里葉變換數(shù)值%Y=exp(i*4*pi*f).*fft(y)%將計算出來的頻譜乘以exp(i*4*pi*f)得到頻移后[2,2]之間的頻譜值Y=fft(y)。z=sqrt(Y.*conj(Y))。plot(f(1:100),z(1:100))。title(39。幅頻曲線39。)xiangwei=angle(Y)。figure(2)plot(f,xiangwei)title(39。相頻曲線39。)figure(3)plot(t,y,39。r39。)%axis([2,2,0,])title(39。原始信號39。) ( KB)201025 18:42 ( KB)201025 18:42 ( KB)201025 18:42（2）用Hilbert變換直接求該信號的瞬時頻率復制內(nèi)容到剪貼板代碼:clear。clc。clf。%假設待分析的函數(shù)是z=t^3N=2048。%fft默認計算的信號是從0開始的t=linspace(1,2,N)。deta=t(2)t(1)。fs=1/deta。x=5*sin(2*pi*10*t)+5*sin(2*pi*35*t)。z=x。hx=hilbert(z)。xr=real(h
公司管理相關推薦
鄂ICP備17016276號-1