计算信号的Wigner-Ville分布 附M文件 出问题

ddy41

clear all
close all
%采样频率
fs=1920;
%信号产生
for i=0:511
t=i/fs;
f=(1+0.2*sin(2*pi*7.5*t))*cos(2*pi*30*t+0.5*sin(2*pi*15*t))+sin(2*pi*120*t);
j=i+1;
s(j)=f;
end
figure;
plot(s);
grid on;
%计算信号的Wigner-Ville分布
[tfr,t,f]=tfrwv(s');
%绘图
%%绘制等高线图
figure;
contour(t,(0:255)/512*fs,abs(tfr(1:256,:)));
xlabel('时间 t');
ylabel('频率 f');
%%三维图
figure;
mesh(t,(0:255)/512*fs,abs(tfr(1:256,:)));
xlabel('时间 t');
ylabel('频率 f');
zlabel('幅值 A');

结果与信号有出入，120HZ处无峰值（感觉整个就不对 原信号是一个调幅调频信号＋正弦信号） 不知什么原因 请高手帮忙看看！

[ 本帖最后由 ddy41 于 2008-6-2 18:57 编辑 ]

ddy41

还有200＋HZ处多出了峰值，请高手帮我看看

ddy41

另外 信号s 写成两列（调频调幅一列 正弦一列） 帮助里说是求两信号的“cross wv” 怎么理解 谢谢

ddy41

我怀疑是信号的问题 大家帮帮忙啊

ddy41

问题已解决 谢谢大家