文档介绍：小波分析实验报告姓名: 班级: 学号: 成绩: 教师签名: 实验一 小波函数的Fourier变换和Fourier逆变换实验目的 用Matlab实现函数的Fourier变换和Fourier逆变换实验内容用Matlab实现下列函数的Fourier变换和Fourier逆变换_x2Morlet小波 肖⑴=e2elf,,l'xco0>5»symsxiwO;»f=exp(-xA2/2)*exp(i*wO*x);»F=fourier(f,x)(2A(l/2)*piA(l/2))/exp((x+i*wO*sqrt(-1))A2/2)»f=ifourier(F)f=cxp((i人2*w0人2)/2-(t-i*wO)人2/2)Marr小波 =(1-x2)e2»symsx;»f=(l-xA2)*exp(-xA2/2);»F=fourier(f)F=(2A(l/2)*piA(l/2)*wA2)/exp(wA2/2)»f=ifourier(F)f==e(2A(1/2)*((2A(1/2)*piA(l/2))/exp(xA2/2)-(2A(1/2)*piA(l/2)*xA2)/exp(xA2/2)))/(2*piA(1/2))DOG小波»symsx;f=exp(-xA2/2)-exp(-xA2/8)/2;»F=tourier(f)(2A(l/2)*piA(l/2))/exp(wA2/2)-(2A(l/2)*piA(l/2))/exp(2*wA2)»f=ifourier(F)((2*pi)/exp(xA2/2)-pi/exp(xA2/8))/(2*pi)] X1 jjGauss函数族ga(x)=— 《4"Q=l,—,—2y]7ra 416»symsx;»gl=l/(2*piA(l/2))*exp(-xA2/4);»Fl=fourier(gl)Fl=(5081767996463981*pi人(1⑵)/(9007199254740992*exp(wA2))»gl=ifourier(Fl)gl=5081767996463981/(18014398509481984*exp(xA2/4))»symsx;»g2=l/(piA(l/2))*exp(-xA2);»F2=fourier(g2)F2=(5081767996463981*piA(l/2))/(9007199254740992*exp(wA2/4))»g2=ifourier(F2)g2=5081767996463981/(9007199254740992*exp(xA2))a=l/16时»symsx;»g3=2/(piA(1/2))*exp(-4*xA2);»F3=fourier(g3)F3=(5081767996463981*piA(l/2))/(9007199254740992*exp(wA2/16))»g3=ifourier(F3)g3=508176799646398l/(4503599627370496*exp(4*xA2)): 肖0)=幺2 ©)n5及其频域形式»x=-6::6;»w0=6;»y=exp(-1/2.*+wO*i.*x);»plot(x,y); =(1-x2)e?及其频域形式»x=-6::6;»y=(l-).*exp(-);»plot(x,y);A | . y/(x)=e~-―幺飞及其频域形式»x=-6::6;»y=exp(-)-(exp(-))/2;»plot(x,y);4-Gauss函数族同-坐标系中MX躺八 归,誌及其频域形式»x=-5::5;»gl=l/(2*piA(l/2))*exp(-);»g2=1/(piA(1/2))*exp(-);»g3=2/(piA(l/2))*exp(-4*);»plot(x,gl,x,g2,x,g3);(t)=sin(2加)+sin(4加)+sin(10加),并画出图形»(=1::5;»f^sin(2*pi*t)+sin(4*pi*t)+sin(10*pi*t);»plot(t,f);5实验三:一•,并显示各层低频高频图形,加以比较。具体格式如上例单尺度:»loadnoissin;»s=noissin(l:1000);»[cAI,cDI]=dwt(s,,db4,);»A1=upcoef('a\cA1,'db4\l);»Dl=upcoefCd\cD1;db4;l);»subp