求助：局域波时频法
有哪位高手懂得局域波时频分析方法，求有关的 matlab 源程序！！<BR>小弟先谢过了！！！！！！！！！！！！！ 你指的局域波，是大连理工大学的一个博士所提出的 “局域波”吗？如果是，那么我感觉就是 Huang等人提出的经验模态分解，或者 HHT(Hilbert-Huang Transform)。 本帖最后由 wdhd 于 2016-3-11 11:17 编辑基于傅里叶变换的功率谱图上的任一条谱线代表一个同频率的周期信号, 该信号分布在整个时间域内, 称之为全域波. 而非平稳信号的主要特征是其时变性, 其频率是瞬变的, 仅仅是在某一局部时间内才存在, 称之为局域波。
%用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;
=meshgrid(fnew,tnew);
subplot(223),mesh(F,T,abs(TF))
subplot(224),contour(F,T,abs(TF))
不知道是不
页:
[1]