matlab求解系统响应filter,数字信号处理matlab实验报告2.doc

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

实验报告课程名称数字信号处理实验专业班级姓名学号实验名称数字信号处理实验地点实验时间实验成绩一、实验目的及任务学会运用MATLAB求解离散时间系统的零状态响应；学会运用MATLAB求解离散时间系统的单位取样响应；学会运用MATLAB求解离散时间系统的卷积和。二、实验内容与步骤321离散时间系统的响应离散时间LTI系统可用线性常系数差分方程来描述，即(31)MJNIJNXBINYA00其中，(，1，，N)和(，1，，M)为实常数。IA0JBMATLAB中函数FILTER可对式(131)的差分方程在指定时间范围内的输入序列所产生的响应进行求解。函数FILTER的语句格式为YFILTERB,A,X其中，X为输入的离散序列；Y为输出的离散序列；Y的长度与X的长度一样；B与A分别为差分方程右端与左端的系数向量。322离散时间系统的单位取样响应系统的单位取样响应定义为系统在激励下系统的零状态响应，用表示。NNHMATLAB求解单位取样响应可利用函数FILTER，并将激励设为前面所定义的IMPDT函数。MATLAB另一种求单位取样响应的方法是利用控制系统工具箱提供的函数IMPZ来实现。IMPZ函数的常用语句格式为IMPZB,A,N其中，参数N通常为正整数，代表计算单位取样响应的样值个数。323离散时间信号的卷积和运算由于系统的零状态响应是激励与系统的单位取样响应的卷积，因此卷积运算在离散时间信号处理领域被广泛应用。离散时间信号的卷积定义为(32)MNHXNHXY可见，离散时间信号的卷积运算是求和运算，因而常称为“卷积和”。MATLAB求离散时间信号卷积和的命令为CONV，其语句格式为YCONVX,H其中，X与H表示离散时间信号值的向量；Y为卷积结果。用MATLAB进行卷积和运算时，无法实现无限的累加，只能计算时限信号的卷积。对于给定函数的卷积和，我们应计算卷积结果的起始点及其长度。两个时限序列的卷积和长度一般等于两个序列长度的和减1。上机练习1试用MATLAB命令求解以下离散时间系统的单位取样响应。(1)12143NXNYNY程序代码A341B11N035XIMPDTNHFILTERB,A,XSTEMN,H, FILL ,GRIDONXLABEL N ,TITLE 单位取样响应HNMATLAB仿真 (2)21065NXYNY程序代码A5/2610B1N030XIMPDTNHFILTERB,A,XSTEMN,H, FILL ,GRIDONXLABEL N ,TITLE 单位取样响应HNMATLAB仿真 2已知某系统的单位取样响应为，试用MATLAB求当激1087NUNHN励信号为时，系统的零状态响应。5UNX程序代码NX18NH215XUDTNXUDTNX5H08NHUDTNHUDTNH10YCONVX,HNY1NX1NH1NYNY10LENGTHNXLENGTHNH2SUBPLOT311STEMNX,X, FILL ,GRIDONXLABEL N ,TITLE XN AXIS41603SUBPLOT312STEMNH,H , FILL ,GRIDONXLABEL N ,TITLE HN AXIS41603SUBPLOT313STEMNY,Y, FILL ,GRIDONXLABEL N ,TITLE YNXNHN AXIS41605实验心得