小波上机实验选读.doc

小波分析上机实验报告 实验报告一 实验目的 运用傅里叶变换知识对常用的基本函数做基本变换。 加深对因果滤波器的理解，并会判断因果滤波器的类型。 运用卷积公式对基本信号做滤波处理并作出分析，以加深理解 熟悉Matlab中相关函数的用法 实验原理 1．运用傅里叶正、反变换的基本公式： 及其性质，对所要处理信号做相应的傅里叶变换和逆变换。 2．运用卷积的定义式： 对所求信号做滤波处理。 实验步骤与内容 实验题目： Butterworth滤波器，其冲击响应函数为 1. 求 2. 判断是否因果；是低通、高通、带通还是带阻？ 3. 对于信号,画出图形 4. 画出滤波后图形，比较滤波前后图形，你会发现什么，这里取 5. 取采用不同的变量值(初始设定A=α=10) 画出原信号图形与滤波后图形，比较滤波效果。 实验步骤及分析过程： 1.求 由傅里叶变换的定义式可得： （1） 2. 判断是否因果；是低通、高通、带通还是带阻？ 该滤波器的幅频特性为：，假定，绘制该滤波器的幅频特性曲线如下： 图1.1滤波器的幅频特性曲线 （1）观察滤波器响应函数可知，只有在输入信号到达后，该滤波器才会有输出响应，此外实际应用的滤波器均是因果滤波器，非因果不可用；所以，题中滤波器是因果滤波器。 （2）由图1可知，该滤波器为低通滤波器。 3. 对于信号，画出图形 编写matlab程序（见附录1），绘制信号的图形如下： 图1.2 信号f(t) 曲线 4. 画出滤波后图形，比较滤波前后图形，你会发现什么，这里取。 根据卷积运算的滤波特性，编写matlab程序，取时间离散点数n=10000；可绘制该参数下的滤波信号如下图所示： 图1.3 f(t) 滤波后信号曲线 比较图1.2和图1.3中，可以看出：经滤波处理后，信号f（t） 取采用不同的变量值(初始设定A==10) 画出原信号图形与滤波后图形，比较滤波效果。 原始波形如下： （1）根据题意，绘制信号曲线；并取A==1、3、5、1、15 、20(初始设定A==10)采用给定滤波器，进行滤波处理，结果如下： 由1问可知：该滤波器的幅频特性为：转折频率；所以当A=分别取1、3、5、10、15、20时，滤波器对于该频率成分的信号起到抑制作用。 通过观察不同滤波曲线可知：可以看出随着A、a值逐渐增大，滤波后信号毛刺（高频波动信号）增多，即对高频信号的抑制效果变差，同时也可以看出滤波器输出信号中的低频成分也呈增大趋势。 总之，滤波器在A、a值为1时对高频的抑制效果最好，但必须指出在这种情况下低频信号也受到一定的削弱，效果并不一定最好，需要根据实际使用需求设定参数。 程序附录 %信号图绘制 t=linspace(0,pi,100000); f=exp(-t/3).*(sin(2*t)+2*sin(4*t)+0.4*sin(2*t).*sin(40*t)); plot(t,f);grid on xlabel(时间/t); ylabel(幅值/f(t)); %滤波处理后信号 A=10;a=10; h=A*exp(-a*t); fy=conv(f1,h); plot(fy);grid on xlabel(时间/t); ylabel(幅值/f(t)); %信号图绘制 t=linspace(0,pi,100000); f1=exp(-t).*(sin(5*t)+sin(3*t)+sin(t)+sin(40*t)); plot(t,f1);grid on xlabel(时间/t); ylabel(幅值/f(t)); %滤波处理后信号 A=1;a=1; h=A*exp(-a*t); fy=conv(f1,h); plot(fy);grid on xlabel(时间/t); ylabel(幅值/f(t)); title(A=a=10时信号波形); A=3;a=3; h=A*exp(-a*t); fy=conv(f1,h); plot(fy);grid on xlabel(时间/t); ylabel(幅值/f(t)); title(A=a=3时信号波形); A=5;a=5; h=A*exp(-a*t); fy=conv(f1,h); plot(fy);grid on xlabel(时间/t); ylabel(幅值/f(t)); title(A=a=5时信号波形); A=10;a=10; h=A*exp(-a*t); fy=conv(f1,h); plot(fy);grid on xlabel(时间/t); ylabel(幅值/f(t)); title(A=a=10时信号波形); A=15;a=15; h=A*exp(-a*t); fy=conv(f1,h); plot(fy);grid on