声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3773|回复: 9

[综合讨论] 求小波脊线的程序

[复制链接]
发表于 2007-4-7 09:51 | 显示全部楼层 |阅读模式

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

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

x
紧急求助,由于之前没怎么接触小波,现在急需得到小波脊线,哪位大侠有程序可以帮帮忙吗?就是想用小波变换得到信号的时频图,请大家多多指点。谢谢。
回复
分享到:

使用道具 举报

发表于 2007-4-7 09:56 | 显示全部楼层

小波脊线的提取方法??

clc
clear
%设置a的初值;
a(1)=0.003;
Te=0.001;k=1;
N=1000;
t=0:0.001:1;
s=exp(-t.^2).*exp(i*50.*t+10*t.^2);
hudu=50;
wucha=0.001;
count=0;
%用morlet小波进行分析;
kk=length(t);
k=0;
flag=1;
while flag>0
flag=0;
wa=cwt(s,a,'morl');
fai=angle(wa);

for i=2:kk
a(i)=hudu/((fai(i)-fai(i-1))/Te); %求Db;
if abs(a(i)-a(i-1))/(abs(a(i-1)))<0.001 %判断是否满足条件;
k=k+1;
aa=a(i);
bb=i*Te;
break;
else
end
end
ar(k)=aa;%赋值;
bbb(k)=bb;
a=0;aa=0;bb=0;
a(1)=ar(k);%继续求a;
flag=flag+1;
count=count+1;
if count>100%求所有的a;
break;
end
end
plot(ar,bbb);
不知道是不是和你的一个概念???

另外有关小波的资料:

http://forum.vibunion.com/thread-9141-1-1.html

评分

1

查看全部评分

 楼主| 发表于 2007-4-7 21:35 | 显示全部楼层
我所指的小波脊线就是对信号做小波变换后,取小波系数的模极大值形成小波脊线,我也用上述程序做过分析,但是改用我所分析的信号后,程序就无法运行,总是出错了。还请高手指点!谢谢。
发表于 2007-12-27 18:13 | 显示全部楼层

回复 #3 xinquan 的帖子

这个程序的确存在问题。也一直没有解决...
发表于 2007-12-27 21:02 | 显示全部楼层

回复 #4 fudy10 的帖子

什么问题?
发表于 2007-12-28 10:09 | 显示全部楼层

回复 #5 花如月 的帖子

将原程序中的s=exp(-t.^2).*exp(i*50.*t+10*t.^2);
替换别的信号就会出现问题。比如
s=exp(-2*pi*3t).*exp(i*50.*t)。
还有,我觉得这个程序反映的是尺度与时间的关系。并非小波系数与时间的关系。
不知对否?请大家多多讨论啊

评分

1

查看全部评分

发表于 2007-12-28 12:52 | 显示全部楼层

回复 #6 fudy10 的帖子

相信你自己的判断吧,2楼在给的时候本来就加了“?”号

因此希望大家批判地吸收
发表于 2008-4-9 19:33 | 显示全部楼层
再顶,请问有没有高手提供一下取小波系数的模极大值形成小波脊线的程序,谢谢!:@)
发表于 2010-9-2 22:10 | 显示全部楼层
程序没写完整,看上面程序是用迭代算法求的小波脊线.
发表于 2011-10-20 16:17 | 显示全部楼层
回复 3 # xinquan 的帖子

小波脊线求出来了吗?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-26 00:46 , Processed in 0.063402 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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