实验一数字信号处理实验报告.doc

实 验 报 告 课 程： 数字信号处理实验 专业班级： 中兴通信131班 学生姓名： 赵 庆 学 号： 6103113160 实验一?系统响应及系统稳定性 一、实验目的 1.掌握常见信号和系统的表示方法 2.掌握求系统响应的方法 3.掌握时域离散系统的时域特性 4.分析观察系统稳定性 二、实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在平域可以用系统函数描述系统特性。已知输入信号，可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域。在计算机适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱filter函数。也可以用MATLAB语言的工具箱函数conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统响应。系统的稳态性有系统的差分方程的系数决定。实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的。系统的稳态输出是指当n为无穷大时，系统的输出。如果系统稳定，信号加入系统后系统输出的开始一段称为暂态响应，随n的增大，幅度趋于稳定，达到稳态输出。 三、实验内容和步骤 1.编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter函数或conv函数求解系统输出响应的子程序； 2.给定离散时域系统的差分方程，求系统单位脉冲响应和系统响应； （1）产生输入信号X1（n)=R8(n)以及其系统响应： 程序如下： %program1 n=-10:30 x=[n=0n=7] subplot(211)??????????????????? stem(n,x);title(x1) axis([-10 30 0.5 1.5]) b=[0.05 0.05]; a=[1 -0.9]; y=filter(b,a,x); subplot(212) stem(n,y);title(x1（n）的系统响应y1（n）) （2）产生输入信号X2(n)=u(n)以及其系统响应： 程序如下： %program2 n=-100:100 x=[n=0] subplot(211) stem(n,x);title(x2) b=[0.05 0.05]; a=[1 -0.9]; y=filter(b,a,x); subplot(212) stem(n,y);title(x2（n）的系统响应y2（n）) （3）求出系统的单位脉冲响应h（n）并画出波形： 程序如下： %program3 clear; clc n=-50:50; b=[0.05 0.05]; a=[1 -0.9]; h=impz(b,a,n); subplot(111); stem(n,h);title(单位脉冲响应h（n）) ???? ?3.给定离散时域系统的单位脉冲响应，求系统的输出响应，其中x1（n）=R8（n）； （1）用卷积法求出x1对h1的系统响应： 程序如下： %program4 n=-10:10 x=[n=0n=7] subplot(311) stem(n,x);title(x1) n=-20:20 h=[n=0n=9] subplot(312) stem(n,h);title(h1) y=conv(h,x) subplot(313) stem(y);title(y1=x1*h1) ?????? ? ? 用卷积法求出x1对h2的系统响应： 程序如下： %program5 n=-10:10 x=[n=0n=7] subplot(311) stem(n,x);title(x1) n=-10:10 z=[n==0]; b=[1 2.5 2.5 1]; a=[1]; h=filter(b,a,z) subplot(312) stem(n,h);title(h2) y=conv(h,x) subplot(313) stem(y);title(y2=x1*h2 ? 分析观察给定系统的稳定性给定差分方程： Y(n)=1.8237y(n-1)-0.9801y(n-2)+b0x(n)-b0(n-2) ; 令b0=1/100.49 程序如下： %program4 clear; clc s=1/100.49; n=-50:50; b=[s 0 -s]; a=[1 -1.8237 0.9801]; h=impz(b,a,n); subplot(111); stem(n,h);title(系统的单位脉冲响