信号与系统实验报告材料6.docVIP

实用 文档 信号与系统实验报告 实验六：离散信号与系统的时域分析 一、实验目的 1．学会用 M AT L A B表示常用离散信号的方法； 2．学会用 M AT L A B实现离散信号卷积的方法； 3． 学会用 M AT L A B求解离散系统的单位响应； 4． 学会用 M AT L A B求解离散系统的零状态响应； 二、实验原理 1．离散信号的 M AT L A B表示 表示离散时间信号 f(k) 需要两个行向量，一个是表示序号k=[ ]，一个是表示相应 函数值 f=[ ]，画图命令是 stem 。 2. 离散信号的卷积和 两个有限长序列 f1，f2 卷积可调用 MATLAB 函数 conv，调用格式是 f=conv(f1,f2), f 是卷积结果，但不显示时间序号，可自编一个函数dconv给出 f和 k，并画图。 function [f,k]=dconv(f1,f2,k1,k2) 3． 离散系统的单位响应 MATLAB提供画系统单位响应函数impz，调用格式是 impz(b,a) 式中b和 a是表示离散系统的行向量； impz(b,a,n) 式中 b和 a是表示离散系统的行向量，时间范围是0~n； impz(b,a,n1,n2) 时间范围是n1~n2 ；y=impz(b,a,n1,n2) 由y给出数值序列； 4． 离散系统的零状态响应 MATLAB 提供求离散系统零状态响应数值解函数 filter，调用格式为 filter(b,a,x),式中 b 和 a是表示离散系统的向量，x 是输入序列非零样值点行向量，输出向量序号同x一样。 三 实验验证（截图） 2. 3. 4. 5. 6. 7. 8. 解答： 代码 a=[2,-2,1]; b=[1,3,2]; impz(b,a) impz(b,a,60) impz(b,a,-10:40) 图形如下： 已知y(k)+y(k-1)+0.25*y(k-2)=f(k),输入f(t)=e(k),画输出波形，范围0~15。 解答： 代码： a=[1 1 0.25]; b=[1]; t=0:20; x=heaviside(t); y=filter(b,a,x) subplot(2,1,1) stem(t,x) title(输入序列) subplot(2,1,2) stem(t,y) title(响应序列) 图形如下：