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

《数字信号处理matlab实验报告2》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验报告课程名称：数字信号处理实验专业班级：姓名：学号：实验名称数字信号处理实验地点实验时间实验成绩一、实验目的及任务•学会运川MATLAB求解离散时间系统的；状态响应；•学会运用MATLAB求解离散时间系统的单位取样响应；•学会运用MATLAB求解离散时间系统的卷积和。二、实验内容与步骤3.2.1离散时间系统的响应离散时间LTI系统可用线性常系数差分方程来描述，即NM=x(n-j)(3-1)/=o7=0其中，％•(Z=0,1，…，AO和(7=0,1,•••,M)为实常数。MATLAB屮函数filter可对式(13-1)的差分方程在指

2、定时间范围内的输入序列所产生的响应进行求解。函数filtei•的语句格式为y=filter(b，a，x)其中，x为输入的离散序列；y力输出的离散序列；y的长度与x的长度一样；b与a分别力差分方程右端与左端的系数向量。3.2.2离散时间系统的单位取样响应系统的单位取样响应定义为系统在激励K系统的零状态响应，用似＞)表示。MATLAB求解单位取样响应可利用函数filter,并将激励设为前而所定义的impDT函数。MATLAB另一种求单位取样响应的方法是利用控制系统工具箱提供的函数impz来实现。impz函数的常用语句格式为impz(b，a

3、，N)其中，参数N通常为正整数，代表计算单位取样响应的样值个数。3.2.3离散时间信号的卷积和运算由于系统的零状态响应是激励与系统的单位取样响应的卷积，因此卷积运算在离散时间信号处理领域被广泛应用。离散时间信号的卷积定义为OOy(/?)=x(n)*h(n)=x(m)h(n-m)(3-2)m=-<»可见，离散时间信号的卷积运算是求和运算，因而常称为“卷积和”。MATLAB求离散时间信号卷积和的命令为conv,其语句格式为y=conv(x，h)其中，x与h表示离散吋间信号值的向量；y为卷积结果。用MATLAB进行卷积和运算吋，无法实现无限

4、的累加，只能计算时限信号的卷积。对于给定函数的卷积和，我们应计算卷积结果的起始点及其长度。两个时限序列的卷积和长度一般等于两个序列长度的和减1。上机练习：1.试用MATLAB命令求解以下离散时间系统的单位取样响应。(1)3y(n)+4y(n-1)+y(n-2)=x(n)+x(n-1)程序代码：a=l341J;b=flH；n=0:35;x=impDT(n);h=filter(b,a,x);stem(n,h,'fill'),gridonxlabel(’nf)，titlef单位取样响座h(n)MATLAB仿真•)(2)—y(/i)+6y(n

5、-1)+1Qy(n-2)=x(n)程序代码：a=[5/2610];b:[lj;n=0:30;x=impDT(n);h=filter(b，a，x);stem(n,h,'fill'),gridonxlabelCn’)，title(f单位取样响应h⑻MATLAB仿真•)2.己知某系统的单位取样响应为/7（M=试用MATLAB求当激励8信号为5）时，系统的零状态响应。程序代码：nx=-l:8;nh=-2:15;x=uDT(nx)-uDT(nx-5);h二0.8.八nh.*(uDT(nh)-uDT(nh-10));y=conv(x，h);nyl

6、=nx(l)+nh(l);ny=nyl+(0:(length(nx)+length(nh)-2));subplot(311)stem(nx，x，’fiir)，gridonxlabel(’n’)，title(’x(n)’)axis([-41603])subplot(312)stem(nh，h’，'fiir)，gridonxlabelfn')，title(’h(n)’)axis(卜41603])subplot(313)stem(ny，y，'fiir)，gridonxlabel('n')，title(’y(n)=x(n)*h(n)’)axi

7、s([-41605])实验心得: