计算信号的Wigner-Ville分布 附M文件 出问题
clear allclose 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分布
=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 编辑 ] 还有200＋HZ处多出了峰值，请高手帮我看看 另外 信号s 写成两列（调频调幅一列 正弦一列） 帮助里说是求两信号的“cross wv” 怎么理解 谢谢 我怀疑是信号的问题 大家帮帮忙啊 问题已解决 谢谢大家
页:
[1]