希望能够得到大侠帮助，看看我的程序（小波脊线），不胜感激！！
能运行，结果不对，自己看不出来问题，急啊！！！a=0.3; % 设置a的初值；
Te=1;
k=0;
t=0:1000-1
for ttt=1:length(t)
fs(ttt)=0;
end
s=exp(i*2*pi*10/1000*t+i*2*pi*10*(t/1000).^2);
hudu=5;
for k1=1:lenth(t)-1
x(k1)=0;
if (k1==1)
y(1)=0.3;
else
y(k1)=0;
end
end
%用morlet小波进行分析；
kk=length(t);
count=1;
for count=1:lenth(t)-1
wa=cwt(s,a,'morl');
fai=angle(wa);
for k2=2:kk
y(k1)=hudu/((fai(k1)-fai(k1-1))/Te); %求Db；
if ( abs(y(k1)-y(k1-1))/abs(y(k1-1))<0.001) %判断是否满足条件；
k=k+1;
x(k)=y(k1-1); %赋值；
a=x(k); %继续求a;
count=count+1;
break;
end
end
end
for k4=1:999
fs(k4)=hudu/x(k4)/(2*pi)*1000;
end
[ 本帖最后由 zhangnan3509 于 2007-7-4 14:12 编辑 ]
页:
[1]