期末复习学习信号处理原理实验.pdfVIP

本次实验为函数x(t)进行离散Harr小波变换，对其进行逆变换并求其PSNR；

并对x(t)的FT的逆变换求PSNR。代码如下：

f.m

%x(t)assignedintheHW9

functionx=f(t)

x=(cos(10*pi*t)+2*sin(30*pi*t)).*(0=tt5)+

cos(40*pi*t).*(5=tt10)+(cos(60*pi*t)+0.6*sin(90*

pi*t)).*(10=tt15)+sin(100*pi*t).*(15=tt20);

test.m

x=0:0.01:20;

y=f(x);

%DWT

figure(1);

plot(y);

[ca,cd]=dwt(y,haar);

figure(2);

plot(ca);

figure(3);

plot(cd);

Idwt=idwt(ca,cd,haar);

figure(4);

plot(Idwt);

psnr_dwt=psnr(Idwt,[y,y(length(y))])

%FFT

x2=0:0.01:20;

y2=f(x2);

psnr_fft=psnr(ifft(fft(f(x2))),y2)

其中原函数图像如下：

本次实验为函数x（t）进行离散harr小波变换，对其进行逆变换并求其psnr;并

对x（t）的ft的逆变换求psnr。代码如下：

f.m

在HW9函数中分配的％x（t）

x=f（t）

x=（cos（10*pi*t）+2*sin（30*pi*t））。*（10=t＆t15）

+sin（100*pi*t）。*（15=t＆t20）;

test.m

X=0：0.01：20;y=f（x）;％DWT图

（1）;绘图（y）;[ca，cd]=dwt（y，

haar）;图（2）;情节（CA）;图（3）;图

（CD）;IDWT=IDWT（CA，CD，

HAAR）;图（4）;情节（IDWT）;

psnr_dwt=psnr（idwt，[y，y

（length（y））]）％fftx2=0：

0.01：20;y2=f（x2）;psnr_fft=

psnr（ifft

（fft（f（x2）））），y2）其中原函数图

像如下：

进行离散Harr小波变换后，近似分量如下图：

细节分量如下图：

Harr小波变换后，近似分量如下图：

细节分量如下图：