求助--小波脊线提取
请高手帮我看一下程序:clc
clear
%设置a的初值；
a(1)=0.02;
Te=0.001;
t=0:0.001:2;
s=cos(2*pi*30*(t.^2));
s=hilbert(s);
hudu=4*pi;
%用morlet小波进行分析；
kk=length(t);
k=0;
aa=0;
bb=0;
for j=1:1000
wa=cwt(s,a(1),'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);
break;
end
end
ar(k)=abs(aa);%赋值；
a(1)=ar(k);%继续求a;
end
plot(ar);
结果好象不对!
不知道什么原因 没有高手自己顶个 回复 1 # yao1983 的帖子
你的意思是什么不对，求出的结果不对，不是预期值，还是？
页:
[1]