网站服务器怎么建设,美容养生连锁东莞网站建设,广州模板建站软件,深圳建站公司专业公司实验目的#xff1a; 在理解连续小波变换原理的基础上#xff0c;通过编程实现对一维信号进行连续小波变换#xff0c;#xff08;实验中采用的是墨西哥帽小波#xff09;#xff0c;从而对连续小波变换增加了理性和感性的认识#xff0c;并能提高编程能力#xff0c;为…实验目的 在理解连续小波变换原理的基础上通过编程实现对一维信号进行连续小波变换实验中采用的是墨西哥帽小波从而对连续小波变换增加了理性和感性的认识并能提高编程能力为今后的学习和工作奠定基础。 实验工具 计算机matlab6.5 程序附录 (1) 墨西哥帽小波函数,按照(***)式编程 % mexh.m function Ymexh(x) if abs(x)8 Yexp(-x*x/2)*(1-x^2); else Y0; End (2) 实验程序按照(**)式编程详细过程请参考“本实验采取的一些小技巧” % clc;clear; load(data.mat); lenlength(dat); lna70; % 尺度a的长度 azeros(1,lna); wfabzeros(lna,len); %小波系数矩阵 mexhabzeros(1,len); % 离散化小波系数矩阵 for s1:lna %s 表示尺度 for k1:len mexhab(k)mexh(k/s); end for t1:len % t 表示位移 wfab(s,t)(sum(mexhab.*dat))/sqrt(s); %将积分用求和代替 mexhab[mexh(-1*t/s),mexhab(1:len-1)]; %mexhab修改第一项并右移 end end figure(1); plot(dat); title(原始数据图); figure(2); %小波系数谱 image(wfab); colormap(pink(128)); title(小波系数图); %surf(wfab); %title(小波系数谱网格图); %pwfabwfab.*wfab; %%瞬态功率谱 %figure(3); %subplot(1,2,1); %surf(pwfab); %title(瞬态功率谱网格图); %subplot(1,2,2); %contour(pwfab); %title(瞬态功率谱等值线); (3)test函数。 %test 函数 clc;clear; for i1:200 dat(i)sin(2*pi*i*0.05); %正弦波函数 end lenlength(dat); lna40; wfabzeros(lna,len); mexhabzeros(1,len); for s1:lna %s 表示尺度 for k1:len mexhab(k)mexh(k/s); end for t1:len % t 表示位移 wfab(s,t)(sum(mexhab.*dat))/sqrt(s); %将积分用求和代替 mexhab[mexh(-1*t/s),mexhab(1:len-1)]; %mexhab修改第一项并右移 end end figure(1); plot(dat); title(orignal dat); figure(2); %小波系数谱 image(wfab); colormap(pink(128)); title(正弦波的小波系数图); 4用fft实现cwt %按照圆周卷积定理原周卷积和线性卷积的关系LMN-1 %按照圆周卷积的定义相关和线性卷积的关系原始算法和线性卷积的关系 %注意画图理解 clc;clear; t1cputime; load(data.mat); lenlength(dat); lna70; % a尺度的长度 azeros(1,lna); % a 表示尺度 bzeros(1,len); % b 表示位移 wfabzeros(lna,len); %小波系数矩阵 mexhabzeros(1,2*len-1); data[zeros(1,len-1),dat]; Ydatafft( data ,4*len); for s1:lna for k1:2*len-1 mexhab(k)mexh((k-len)/s); end tempifft( Ydata.*fft( mexhab,4*len ) ,4*len); wfab(s,:)real(temp(2*len-1:3*len-2))/sqrt(s); %为什么要取实部而不是取模我也不是很清楚可是有种感觉 end figure(1); plot(dat); title(原始数据图); figure(2); %小波系数谱 image(wfab); colormap(pink(128)); title(小波系数谱 ); cputime-t1 4fft快速计算cwt %按照圆周卷积的定义 %注意画图理解 clc;clear; t1cputime; load(data.mat); lenlength(dat); lna70; % a尺度的长度 a5; data[dat,zeros(1,len)]; Ydatafft(dat,2*len); for s1:lna mexhabzeros(1,2*len); k[-a*s:1:a*s]; mexhab(klen)mexh2(k./s); tempifft( Ydata.*fft( mexhab,2*len ) ,2*len); wfab(s,:)real(temp(len1:2*len))/sqrt(s); %要取实部而不是取模呵呵 end figure(1); plot(dat); title(原始数据图); figure(2); %小波系数谱 image(wfab); colormap(pink(128)); title(小波系数谱 ); cputime-t1 5保存为mexh2.m function Ymexh2(x) Yexp(-x.*x/2).*(1-x.^2); Torstan 2005.09.16转载于:https://www.cnblogs.com/Torstan/archive/2011/08/31/2161427.html
