声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1064|回复: 1

[小波] 请帮忙看一个wavelab工具箱awt函数的问题,谢谢!

[复制链接]
发表于 2012-8-1 19:28 | 显示全部楼层 |阅读模式

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

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

x
我目前正在研究如何用小波得到非平稳信号的瞬时频率。我大致读了“Wavelet Tour of signal processing”这本书的4.3和4.4节。为了测试解析小波,我做了一个简单的例子。我用了WAVELAB工具箱提供的AWT函数。然而,测试结果和预期相距甚远。下面是我的matlab代码:

////////////////////////////MATLAB PROGRAM/////////////////////////////////////////////////////////
t=1/800:1/800:1;
s=sin(2*pi*50*t)+sin(2*pi*100*t);  %the signal has a duration of 1 sec, with 800 sample rate, and with two frequency elements, 50 and 100%
awt = AWT(s);   %analytic wavelet transform with default voice number of 12%
%set and figure out the scalogram of signal s%
delta = 1/15;
unit = (1-3*delta)/3;
h3 = [1.5*delta delta 1-2.5*delta 2*unit];
h2 = [delta 2*(unit+delta) 1-2*delta 2*unit];
h1 = [delta 4*unit+3*delta 1-2*delta unit];
h0 = [delta delta 1-2*delta 1-2*delta];
sz     = size(awt);
logawt = awt';
linawt = log2linChap4(flipud(logawt),12);
figure(1);clf
set(gcf, 'NumberTitle','off', 'Name','Window 1')
axes('position',h3);
ImageAWTnew(linawt','Overall','gray','lin');
//////////////////////////////////////////////////////////////////////////////////////////////////////////////

以上的代码是我参照WAVELAB提供的“Wavelet Tour of signal processing”这本书的参考代码写的,这个参考代码就是图4.16的matlab仿真,参考代码叫做wt04fig16.m可以在wavelab850\\Books\WaveTour\WTCh04 目录中找到。 下面贴出结果:
tt.jpg
从图中可以看出,小波系数的大部分能量集中在80hz和160hz左右,而不是预想的50hz和100hz。请问谁能告诉我这个代码有什么问题吗?还是AWT这个函数有问题呢?
谢谢大家!
回复
分享到:

使用道具 举报

 楼主| 发表于 2012-8-1 22:21 | 显示全部楼层
up ~ up~ up~
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 05:23 , Processed in 0.082627 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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