求助关于平滑伪WVD

nanwenya

新建 Microsoft Word 文档 (3).doc (52 KB, 下载次数: 9)

2009-3-30 10:12 上传

点击文件名下载附件

请问为什么我的时频图纵坐标不对呢？请大家帮我看下吧 谢拉
sig=xlsread('D:\WRI\both\1\WTIx_both1.xls');
fs=20;
%时域波形
figure(1);
plot(sig);
xlabel('时间t');
ylabel('幅值A');
title('WTIx_both1时域');
[tfr,t,f]=tfrspwv(sig);  
figure(2);
contour(t,f,tfr);
xlabel('t');
ylabel('f');
title('WTIx_both1平滑伪Wigner-Ville分布');

lingyunzhi

给你一个例子，你可以看一下，这是CWD分布，要是用平滑伪可以换成tfrspwv(sig);  
clc
clear
N=1024;
fs=1024;
t=1/fs*(0:N/4-1);
x1=sin(2*pi*150*t);
t=1/fs*(N/4:N/2-1);
x2=sin(2*pi*350*t);
t=1/fs*(N/2:3*N/4-1);
x3=sin(2*pi*50*t);
t=1/fs*(3*N/4:N-1);
x4=sin(2*pi*20*t);

x=[x1 x2 x3 x4];
t=1/fs*(0:(N-1));
figure(1)
plot(t,x);xlabel('时间 t/s');ylabel('幅值 A/um');
%计算Choi-Williams分布
sig=x';
sig1=hilbert(sig);
[tfr,t,f]=tfrcw(sig1);
figure(2)
contour(t/fs,f(1:N)*fs,abs(tfr));
xlabel('Choi-Williams分布时频图   时间 t/s');
ylabel('频率 f/Hz');
grid on