数字信号处理实验报告95152.docVIP

数字信号处理课程设计报告 通信11级16班 学号 姓名：薛焕福 实验一 离散时间系统及离散卷积 一、实验目的 1、熟悉MATLAB软件的使用方法。 2、熟悉系统函数的零极点分布、单位脉冲响应和系统频率响应等概念。 3、利用MATLAB绘制系统函数的零极点分布图、系统频率响应和单位脉冲响应。 4、熟悉离散卷积的概念，并利用MATLAB计算离散卷积。 二、实验原理 本实验利用MATLAB的库函数，实验冲击响应零极点、幅度相位谱的绘制。熟悉MATLAB的使用。并利用离散卷积的概念编写一个计算离散卷积的程序。 三、实验内容及要求 1.离散时间系统的单位冲击响应 （1）选择一个离散时间系统； （2）用笔进行差分方程的递推计算； （3）编制差分方程的递推计算程序； （4）在计算机上实现递推运算； （5）将程序计算结果与笔算的计算结果进行比较，验证程序运行的正确性； 2.离散系统的幅频、相频的分析方法 （1）离散系统的幅频、相频的分析方法 （2）用笔计算几个特殊的幅频、相频的值，画出示意曲线图； （3）编制离散系统的幅频、相频的分析程序； （4）在计算机上进行离散系统的幅频、相频特性计算，并画出曲线； （5）通过比较，验证程序的正确性； 3. 离散卷积的计算 （1）选择两个有限长序列，用笔计算其线性卷积； （2）编制有限长序列线性卷积程序； （3）利用计算程序对（1）选择的有限长序列进行卷积运算； （4）比较结果验证程序的正确性。 四、实验过程 1、离散时间系统的单位脉冲响应程序: a=[1,-1,0.91]; % 系统的系数 b=[1,2,3]; n0=0; n=-20:120; x=[(n-n0)==0]; % 利用filter函数进行滤波，从而求得系统的单位冲击响应 h=filter(b,a,x); figure(1) stem(n,h); % 利用stem绘制曲线图 title(冲击响应); xlabel(n); ylabel(h(n)); h figure(2) [z,p,g]=tf2zp(b,a); % 求传递函数 zplane(z,p) % 求零极点 （2）实验结果： 2、离散系统的幅频、相频的分析程序： b=[0.6,1.2,0.4,0.9]; % 分子的系数 a=[1.0,-1.7,1.1,-0.5] ; % 分母的系数 m=0:length(b)-1; l=0:length(a)-1; K=500; k=1:K; w=pi*k/K; H=(b*exp(-j*m*w))./(a*exp(-j*l*w)); % 求传递函数 magH=abs(H); % 求传递函数的幅值 angH=angle(H); % 求传递函数的相位 figure(1) subplot(2,1,1); plot(w/pi,magH); % 定义坐标轴的区间 grid; axis([0,1,0,1]); xlabel(w(pi)); ylabel(|H|); title(·幅度，相位响应); subplot(2,1,2); plot(w/pi,angH); grid; xlabel(w(pi)); ylabel(angle(H)); （2）实验结果： 3、离散卷积的计算程序： function pr13() n=-5:50; u1=stepseq(0,-5,50); u2=stepseq(10,-5,50); x=u1-u2; h=((0.9).^n).*u1; figure(1) subplot(3,1,1); stem(n,x); axis([-5,50,0,2]); title(输入序列); xlabel(n); ylabel(x(n)); subplot(3,1,2); stem(n,h); axis([-5,50,0,2]); title(冲击响应序列); xlabel(n); ylabel(h(n)); [y,ny]=conv_m(x,n,h,n); subplot(3,1,3); stem(ny,y); axis([-5,50,0,8]); title(输出响应); xlabel(n); ylabel(y(n)); function [x,n]=stepseq(n0,n1,n2) n=n1:n2; x=[(n-n0)=0]; function [y,ny]=conv_m(x,nx,h,