python连续小波变换_2-1 连续小波变换程序

2-1是对分形数据进行墨西哥帽子小波变换，实现连续小波变换。程序文件列表，见表1。

表1    程序文件列表

文件

作用

Singularity_Detection.m

对已知信号进行墨西哥帽子变换

test.m

对分形数据vonkoch进行变换

源程序-Singularity_Detection.m

%某个尺度的连续小波变换的M函数

%delta小波变换的尺度

%N小波函数的长度

%s原始信号

%g原始信号某个尺度下的小波变换系数

function g= Siguarity_Detection(delta, N, s);

%原始信号长度

n= length(s);

%构造墨西哥帽子小波函数

for index_x= 1: N;

x= index_x-(N-1)/2;

phi_x(index_x)= ((pi^ (-1/4))*(2/sqrt(3)))*(1-x.*x/(delta^2))*exp(-(x.*x)/(2*delta^2));

end;

%对信号做卷积

phi_x= phi_x/ norm(phi_x);%能量归一化

g= conv(s,phi_x);%卷积

g= wkeep(g, n);%保持信号长度

源程序-test

%多个尺度连续小波变换的实现

clc;clear

%下载信号

load vonkoch

vonkoch= vonkoch(1: 510);

%尺度1-32的连续小波变换

S_Min&#