一维连续小波变换实验报告.docxVIP

连续小波变换实验报告实验目的：通过matlab编程实现一维连续小波变换，更好地理解连续小波变换的算法和作用，以及小波系数矩阵的含义。同时通过小波系数矩阵对原始信号进行频谱分析，并了解小波小波系数在尺度和位移两个分量上的意义。实验原理：一维连续小波变换公式：当小波函数为实函数时在给定尺度下，对待分析信号和小波函数按照，进行采样，其中为采样间隔，则小波变换可近似如下：=对给定的值，依次求出不同值下的一组小波系数，由于数据采样间隔?t为0.03（常量），所以可以把这个系数忽略,并通过公式下面对小波变换矩阵进行归一化处理。实验结果：实验程序及注释（1）主程序load(data.mat);n=length(dat);amax=70; %尺度a的长度a=zeros(1,amax);wfab=zeros(amax,n); %小波系数矩阵，均以零矩阵形式赋初值mexhab=zeros(1,n); %某尺度下小波系数for s=1:amax %s表示尺度 for k=1:nmexhab(k)=mexh(k/s); endfor t=1:n % t 表示位移wfab(s,t)=(sum(mexhab.*dat))/sqrt(s); %将积分用求和代替mexhab=[mexh(-1*t/s),mexhab(1:n-1)]; %mexhab修改第一项并右移endendwfab_abs=abs(wfab);figure(3); %画三维图colormap(pink(255)); surfc(wfab_abs); xlabel(TRANSLATION) ylabel(SCALE)zlabel(AMPLITUDE) for index=1:amax %小波系数矩阵归一化处理 max_coef=max(wfab_abs(index,:)); min_coef=min(wfab_abs(index,:)); ext=max_coef-min_coef; wfab_abs(index,:)=255*(wfab_abs(index,:)-min_coef)/ext;endfigure(1);plot(dat); %画原始数据图title(ORIGINAL DATA);xlabel(TIME)ylabel(AMPLITUDE) figure(2); image(wfab_abs); %画尺寸-位移图colormap(pink(255)); title(COEFFS ABSOLUTE);xlabel(TIME) ylabel(SCALE)（2）墨西哥帽小波函数function Y=mexh(x) %单独用.M文件定义此函数if abs(x)=5Y=((pi^(-1/4))*(2/sqrt(3)))*(1-x*x)*exp(-(x*x)/2);elseY=0;end;