specgram

xueyebaixuan · 发表于 2009-4-24 15:03

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

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

x

specgram(xx,1199,500,hanning(256),0);
为什么算出来的横轴--时是轴不是从0开始的,而是从-0.2开始的呢?

ChaChing · 发表于 2009-4-24 21:46

时域解析度不够造成这麽明显!
LZ的例子明显仅有4次

xueyebaixuan · 发表于 2009-4-26 19:12

那长为1200的序列划成几组比较适当呢?

ChaChing · 发表于 2009-4-26 22:50

看LZ的需求而定吧! 或许可试试使用overlaps

xueyebaixuan · 发表于 2009-4-27 12:43

哦,谢谢!
还有我想请教一下.
毕设老师让求相干的时频表示.
function [S,F,T]=stft(xx,yy);
[Sxx,Fxx,Txx]=specgram(xx,1199,500,hanning(200),0);
Syy=specgram(yy,1199,500,hanning(200),0);
Sxx=abs(Sxx).^2;%信号xx的自谱
Syy=abs(Syy).^2;%信号yy的自谱
Sxy=Syy.*conj(Sxx)%%%%%%请问这个式子求的是两个信号的互谱吗?
S=(abs(Sxy).^2)./(Sxx.*Syy); %相干函数
F=Fxx;
T=Txx;
imagesc(T,F/2/pi,20*log(S));axis xy; colormap(jet);
由此求出来的结果有几点问题:
(1)S有大于1的数出现.
(2)试了好几组信号,却发现求出来的图像都是在某一特定的时域处有最大值,不知道是不是程序有什么问题?

[ 本帖最后由 xueyebaixuan 于 2009-4-27 12:45 编辑 ]

ChaChing · 发表于 2009-4-27 14:33

我并非讯号处理专家! 数学公式也模糊了, 怕说错误导, 待高人路过确定! 但直觉下列有问题, LZ自己再审视吧
Sxx=abs(Sxx).^2; Syy=abs(Syy).^2; % 已经为实数
Sxy=Syy.*conj(Sxx); 何来conj?
是否应为
Pxx=abs(Sxx).^2; Pyy=abs(Syy).^2; Pxy=Syy.*conj(Sxx); S=(abs(Pxy).^2)./(Pxx.*Pyy);

xueyebaixuan · 发表于 2009-4-27 15:17

哦,是的.我也刚刚发现了这个问题.可是改完之后,出来的数据全是1.:'( 我要晕死了~

xueyebaixuan · 发表于 2009-4-27 15:25

哦,我刚才也发现这个问题,可是改完后得到的相干全变成1了.有人说相干要分段平滑之类的.可是短时傅里叶已经是分段的了,不能再分段了.唉~

账号		自动登录	找回密码
密码			我要加入

[其他] specgram

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

回复沙发 ChaChing 的帖子

回复板凳 xueyebaixuan 的帖子

回复地板 ChaChing 的帖子

回复 5楼 xueyebaixuan 的帖子

回复 6楼 ChaChing 的帖子

回复 6楼 ChaChing 的帖子

[其他] specgram

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

回复 沙发 ChaChing 的帖子

回复 板凳 xueyebaixuan 的帖子

回复 地板 ChaChing 的帖子

回复 5楼 xueyebaixuan 的帖子

回复 6楼 ChaChing 的帖子

回复 6楼 ChaChing 的帖子

回复沙发 ChaChing 的帖子

回复板凳 xueyebaixuan 的帖子

回复地板 ChaChing 的帖子