《数字信号处理实验一.docVIP

河北联合大学 信息工程学院 通信工程系 数字信号处理 实验名称：系统响应及系统稳定性 班 级： 11电科1 学 号： 201124420115 姓 名： 王晓云 指导老师：崔东艳 实验名称：系统响应及系统稳定性 学生姓名：王晓云 班级：11电科1 预习报告 实验目的： 掌握求系统响应的方法。 掌握时域离散系统的时域特性。 分析、观察及检验系统的稳定性。 实验要求： 简述在时域求系统响应的方法。 简述通过实验判断系统稳定性的方法。分析上面第三个实验的稳定输出的波形。 对各实验所得结果进行简单分析和解释。 简要回答思考题。 打印程序清单和要求的各信号波形。 实验报告 实验目的: 掌握求系统响应的方法。 掌握时域离散系统的时域特性。 分析、观察及检验系统的稳定性 实验要求： 简述在时域求系统响应的方法。 简述通过实验判断系统稳定性的方法。分析上面第三个实验的稳定输出的波形。 对各实验所得结果进行简单分析和解释。 简要回答思考题。 打印程序清单和要求的各信号波形。 实验方法（步骤及结果）： 编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter函数或conv函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。 给定一个低通滤波器的差分方程为y(n)=0.05x(n)+0.05x(n-1)+0.9y(n-1) 输入信号 x1(n)=R8(n),x2(n)=u(n) 分别求出x1(n)=R8(n)和x2(n)=u(n)的系统响应，并画出其波形。 求出系统的单位脉冲响应，画出其波形。 程序：close all;clear all %======内容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性====== A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量B和A x1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号x1(n)=R8(n) x2n=ones(1,128); %产生信号x2(n)=u(n) hn=impz(B,A,58); %求系统单位脉冲响应h(n) subplot(2,2,1);stem(hn,.); %调用函数tstem绘图 title((a) 系统单位脉冲响应h(n));box on y1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n) subplot(2,2,2);stem(y1n,.); title((b) 系统对R8(n)的响应y1(n));box on y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n) subplot(2,2,4);;stem(y2n,.); title((c) 系统对u(n)的响应y2(n));box on 结果： 给定系统的单位脉冲响应为 h1(n)=R10(n),h2(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3),用线性卷积法求x1（n）=R8（n）分别对系统h1(n)和h2(n)的输出响应，并画出波形。 程序：%===内容2：调用conv函数计算卷积============================ n=0:19 x1n=[1 1 1 1 1 1 1 1 ]; %产生信号x1(n)=R8(n) h1n=[ones(1,10) zeros(1,10)]; h2n=[1 2.5 2.5 1 zeros(1,10)]; y21n=conv(h1n,x1n); y22n=conv(h2n,x1n); figure(2) subplot(2,2,1);;stem(n,h1n,.); %调用函数tstem绘图 title((d) 系统单位脉冲响应h1(n));box on subplot(2,2,2);stem([0:length(y21n)-1],y21n,.); title((e) h1(n)与R8(n)的卷积y21(n));box on subplot(2,2,3);stem([0:length(h2n)-1],h2n,.); %调用函数tstem绘图 title((f) 系统单位脉冲响应h2(n));box on subplot(2,2,4);stem([0: