声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1302|回复: 1

[综合] 求助关于平滑伪WVD

[复制链接]
发表于 2009-3-30 10:12 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
新建 Microsoft Word 文档 (3).doc (52 KB, 下载次数: 9) 请问为什么我的时频图纵坐标不对呢?请大家帮我看下吧 谢拉
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分布');
回复
分享到:

使用道具 举报

发表于 2009-6-15 16:40 | 显示全部楼层
给你一个例子,你可以看一下,这是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
回复 支持 1 反对 0

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 16:24 , Processed in 0.078754 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表