小波分析实验报告.doc

文档介绍：小波分析实验报告姓名: 班级: 学号: 成绩: 教师签名:实验一小波函数的Fourier变换和Fourier逆变换实验目的 用Matlab实现函数的Fourier变换和Fourier逆变换实验内容用Matlab实现卜列函数的Fourier变换和Fourier逆变换工Morlct小波 0(兀)=e2€叫690>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(-l))A2/2)»f=ifourier(F)f=exp((iA2*wOA2)/2-(t-i*wO)A2/2)rMarr小波 {/(x)=(1-x2)e2»symsx;»f=(l-xA2)*exp(-xA2/2);»F=fourier(f)(2A(l/2)*piA(1/2)*wA2)/exp(wA2/2)»f=ifourier(F)f=(2A(l/2)*((2A(l/2)*piA(l/2))/exp(xA2/2)-(2A(l/2)*piA(l/2)*xA2)/exp(xA2/2)))/(2*piA(1/2))DOG小波»symsx;f=exp(-xA2/2)-exp(-xA2/8)/2;»F=fouricr(f)(2A(l/2)*piA(l/2))/exp(wA2/2)-(2A(l/2)*piA(l/2))/exp(2*wA2)»f=ifourier(F)f=((2*pi)/exp(xA2/2)-pi/exp(xA2/8))/(2*pi)1工 1 1Gauss函数族ga(x)=—e4aQ=l,—,—416»symsx;»gl=l/(2*piA(l/2))*exp(-xA2/4);»Fl=fourier(gl)Fl=(5081767996463981*piA(l/2))/(9007199254740992*exp(wA2))»gl=ifourier(Fl)gl=5081767996463981/(18014398509481984*exp(xA2/4))»symsx;»g2=1/(piA(1/2))*exp(-xA2);»F2=fourier(g2)F2=(5081767996463981*piA(1/2))/(9007199254740992*exp(wA2/4))»g2=ifburier(F2)g2=5081767996463981/(9007199254740992*exp(xA2))a=l/16时»symsx;»g3=2/(piA(l/2))*exp(-4*xA2);»F3=fourier(g3)F3=(5081767996463981*piA(l/2))/(9007199254740992*exp(wA2/l6))»g3=ifourier(F3)g3=508176799646398l/(4503599627370496*cxp(4*xA2)): y/(x)=e2e,(Oi'xco{}>5及其频域形式»x=-6::6;»w0=6;»y=exp(-l/2.*-wO*i.*x);»plot(x,y); ^(x)=(1-x2)e?及其频域形式»x=-6::6;»y=(1-).*exp(-);»plot(x,y);x2[— 讥x)=e'-— 及其频域形式2»x=-6::6;»y=exp(-)-(exp(-))/2;»plot();i =—e4aa=1,—,—及其频域2yj7ra 416vvXP5OO53vvgl也、(2沃p>(1/2))<?xp(,x・>2、4);vvg2Hl、(pi>(l、2))*exp(—x>2)八vvg3==2、(pi>(l、2))*exp(-4*x・>2rvvplot(x9gl><aQ2><aQ3)-5・wx51iF<Jlafsnsin(2a)+sin(4a)+sin(lom)wjfs'jtBffl^®vvfHsin(2*pi<5+sin(4*pi&+sin(10*p*t);VVP105O实验三:,并显示各层低频高频图形,加以比较。具体格式如上例单尺度:»loadnoissin;»s=noissin(1:1000);»[cAl,cDI]=dwt(s;db4,);»Al=upcoef('a',cA1,'db4',l);>>DI二upcoef('