求助：局域波时频法

lz1736

有哪位高手懂得局域波时频分析方法，求有关的 matlab 源程序！！<BR>小弟先谢过了！！！！！！！！！！！！！

峰之巅

你指的局域波，是大连理工大学的一个博士所提出的 “局域波”吗？如果是，那么我感觉就是 Huang等人提出的经验模态分解，或者 HHT(Hilbert-Huang Transform)。

zhyuer

　　基于傅里叶变换的功率谱图上的任一条谱线代表一个同频率的周期信号, 该信号分布在整个时间域内, 称之为全域波. 而非平稳信号的主要特征是其时变性, 其频率是瞬变的, 仅仅是在某一局部时间内才存在, 称之为局域波。

　　%用WD变换解析信号x=exp(j*pi*k*t.^2)

　　clear,close all

　　k=4;T=4;

　　fc=k*T;fs=4*fc;%采样频率大于4倍的信号宽度

　　Ts=1/fs;N=T/Ts;

　　x=zeros(1,N);

　　t=0:N-1;

　　x=exp(j*k*pi*(t*Ts).^2);

　　subplot(221),plot(t*Ts,real(x));

　　X=fftshift(fft(x));

　　subplot(222),plot((t-N/2)*fs/N,abs(X))

　　R=zeros(N);

　　for n=0:N-1

　　M=min(n,N-1-n);

　　for k=0:M

　　R(n+1,k+1)=x(n+k+1)*conj(x(n-k+1));

　　end

　　for k=N-1:-1:N-M

　　R(n+1,k+1)=conj(R(n+1,N-k+1));

　　end

　　end

　　TF=zeros(N);

　　for n=0:N-1

　　temp=fftshift(fft(R(n+1,:)));

　　TF(n+1,:)=temp;

　　end

　　fnew=(t-N/2)*fs/2/N;

　　tnew=(0:N-1)*Ts;

　　[F,T]=meshgrid(fnew,tnew);

　　subplot(223),mesh(F,T,abs(TF))

　　subplot(224),contour(F,T,abs(TF))

　　不知道是不