《数字信号处理报告—实习报告》.docVIP

数字信号处理报告学 校： 中国地质大学（武汉）姓 名：班 级：学 号：指导老师：目录 1-1：离散卷积的运算 1-2：离散傅里叶变换及其应用 1-3：ＩＩＲ滤波器的设计 1-4：ＦＩＲ低通滤波器的设计 1-5：用ＭＡＴＬＡＢ设计一信号处理系统 1-6:频率抽选 1-7：实习心得体会1-1离算卷积的应用 1：基本要求h(n)，输入序列为x(n) 1、h(n)=(0.8)n，0≤n≤4； x(n)=u(n)-u(n-4) 2、h(n)=(0.8)nu(n), x(n)=u(n)-u(n-4) 3、h(n)=(0.8)nu(n), x(n)=u(n) 求以上三种情况下系统的输出y(n)，显示输入和输出波形。 2：实验目的及算法原理 实验目的： （1）：掌握离散卷积计算机实现。 （2）：进一步对离散信号卷积算法的理解。 算法原理：把冲激响应h(n)与输入序列x(n)分别输入到程序中，然后调用离散卷积函数y=conv(x.,h)即可得到所要求的结果。 离散卷积定义为 当序列为有限长时则 3：程序代码 x1=[1 1 1 1 ];nx1=0:3; h1=[1 0.8 0.64 0.8^3 0.8^4];nh1=0:4; y1=conv(x1,h1); subplot(3,3,1);stem(nx1,x1);title(序列x1); xlabel(n);ylabel(x1(n)); subplot(3,3,2);stem(nh1,h1);title(序列h1); xlabel(n);ylabel(h1(n)); subplot(3,3,3);stem(y1);title(序列y1); xlabel(n);ylabel(y1(n)); x2=[1 1 1 1];nx2=0:3; nh2=0:1:20; h2=(0.8).^nh2; y2=conv(x2,h2); subplot(3,3,4);stem(nx2,x2);title(序列x2); xlabel(n);ylabel(x2(n)); subplot(3,3,5);stem(nh2,h2);title(序列h2); xlabel(n);ylabel(h2(n)); subplot(3,3,6);stem(y2);title(序列y2); xlabel(n);ylabel(y2(n)) ; nx3=0:1:20; x3=1.^nx3; nh3=0:1:20; h3=(0.8).^nh3; y3=conv(x3,h3); subplot(3,3,7);stem(nx3,x3);title(序列x3); xlabel(n);ylabel(x3(n)); subplot(3,3,8);stem(nh3,h3);title(序列h3); xlabel(n);ylabel(h3(n)); subplot(3,3,9);stem(y3);title(序列y3); xlabel(n);ylabel(y3(n)); 4：结果分析1-2离散傅里叶变换及应用 1：基本要求x(n)=cos(0.48πn)+cos(0.52πn) 分析下列三种情况下的幅频特性。 采集数据长度N=16，分析16点的频谱，并画出幅频特性。采集数据长度N=16，并补零到64点，分析其频谱，并画出幅频特性。 采集数据长度N=64，分析64点的频谱，并画出幅频特性。 观察三幅不同的幅频特性图，分析和比较它们的特点及形成原因。 2：实习目的及相关算法 实习目的： （1）了解DFT及FFT的性质和特点。 （2）利用FFT算法计算信号的频谱。 相关算法： 读入离散序列x(n) 采集长度为N=16的数据，调用matlab中的函数fft（x,16）与fft（x,64）对其作离散傅里叶变换得到16点、64点的频谱 采集数据长度为N=64，，调用matlab中的函数fft（x,46）对其作离散傅里叶变换得到64点的频谱。 3：相关实现代码1-3：IIR滤波器的设计 1：基本要求 IIR滤波器。 2：实习目的及相关算法 实习目的： （1）掌握不同IIR滤波器的性质、特点。 （2）通过实验学习如何设计各种常用的IIR滤波器，以便在实际工作中能根据具体情况使用IIR滤波器。 相关算法： 输入通带截止频率Wp，阻带截止频率Ws，通带衰减Rp，阻带衰减Rs，通过这些数值调用[N Wn]=buttord(Wp,Ws,Rp,Rs) 函数计算巴特沃斯数字滤波器的阶数N和截止频率wn，再根据阶数N通过函数[b,a]=butter(N,Wn)，即可得到所要的巴特沃斯滤波器。设计一个正弦波信号，再调用函数A=filter(b,a,I)