求助--小波脊线提取

yao1983

请高手帮我看一下程序:
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);

结果好象不对!
不知道什么原因

yao1983

没有高手自己顶个

zxn880625

回复 1 # yao1983 的帖子

你的意思是什么不对，求出的结果不对，不是预期值，还是？