复小波变换提取小波脊线

复小波变换提取小波脊线，可较好识别频率

clc;clear;close all
t1=0:1/40:7.5-1/40;
t2=7.5:1/40:25.6-1/40;
t=[t1t2];
x1=sin（1.6*2*pi*t1）+sin（3*2*pi*t1+0.04*2*pi*t1.^2+4*sin（0.1*2*pi*t1））;
x2=sin（1.6*2*pi*（t2-7.5）+0.055*2*pi*（t2-7.5）.^2）+sin（3*2*pi*t2+0.04*2*pi*t2.^2+4*sin（0.1*2*pi*t2））;
x=[x1x2];
N=length（x）;
y=awgn（x11）;
fs=40;fc=3;totalscal=200;
FreqBins=linspace（1/fs6/fstotalscal）;
fw=FreqBins*fs;
scales=fc./FreqBins;

coes_y=cwt（yscales‘cmor 4-3‘）;
imagesc（tscalesabs（coes_y））;
figure（2）
imagesc（tfwabs（coes_y））;
a=zeros（1N）;
b=zeros（1N）;
b_p=zeros（1N）;


[a（1）b（1）]=max（abs（coes_y（81:1211）））;
temp1=80;b_p（1）=b（1）+temp1;
for i=2:N
    [a（i）b（i）]=max（abs（coes_y（temp1+b（i-1）-10:temp1+b（i-1）+10i）））;
    temp1=b_p（i-1）-10-1;
    b_p（i）=b（i）+temp1;
end
figure（3）
plot（tfw（b_p））
c=zeros（1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      84917  2010-08-12 11:04  原始信号小波量图.png

     文件      38786  2010-08-12 11:04  提取出的小波脊线.png

     文件       1846  2010-08-10 09:19  wtsimulation.m

----------- ---------  ---------- -----  ----

               125549                    3