matlab小波脊线,小波脊线提取，模极大值法。运行的结果不太对，代码有些地方我也没完全看懂...

本帖最后由 1393107100 于 2019-5-1 11:10 编辑

clear,clc

close all

%%%%%%　小波变换　%%%%%%%%%%%%%%%%

fs=1024;

t=1/fs:1/fs:1;

f1=100;f2=200;f3=300;

s=cos(2*pi*f1*t.*(t>=0&t<0.3))+2*cos(2*pi*f2*t.*(t>=0.3&t<0.8))+3*sin(2*pi*f3*t.*(t>=0.8&t<=1));

%s=cos(2*pi*f1*t);

figure(1);

plot(t,s,'b')

title('原始单频信号')

Y = awgn(s,6);%加高斯白噪声

subplot(412);

plot(t,Y,'r')

title('加信噪比为6的白噪声')

%%%%%%%%%%%%%%%%%小波时频图绘制 %%%%%%%%%%%%%%%%%%%%%

wavename='cmor4-4';    %%选用带宽参数和中心频率均为4 的复morlet小波

totalscal=256;          %尺度序列的长度，即scal的长度

fc=centfrq(wavename);%% 小波的中心频率

%%%%%%%%% 绘制时间域小波波形 %%%%%%%%%%%%%%%%%%%

cparam=2*fc*totalscal;%为得到合适的尺度所求出的参数

a=tot