发布时间 : 星期一 文章2016-《数字信号处理》实验讲义(学生)更新完毕开始阅读
图1-3 系统单位取样响应
1.2.3 离散时间信号的卷积和运算
由于系统的零状态响应是激励与系统的单位取样响应的卷积,因此卷积运算在离散时间信号处理领域被广泛应用。离散时间信号的卷积定义为
?y(n)?x(n)*h(n)?m????x(m)h(n?m) (1-2)
可见,离散时间信号的卷积运算是求和运算,因而常称为“卷积和”。
MATLAB求离散时间信号卷积和的命令为conv,其语句格式为
y=conv(x,h)
其中,x与h表示离散时间信号值的向量;y为卷积结果。用MATLAB进行卷积和运算时,无法实现无限的累加,只能计算时限信号的卷积。
例如,利用MALAB的conv命令求两个长为4的矩形序列的卷积和,即
g(n)?[u(n)?u(n?4)]*[u(n)?u(n?4)],其结果应是长为7(4+4-1=7)的三角
序列。用向量[1 1 1 1]表示矩形序列,MATLAB源程序为
>>x1=[1 1 1 1]; >>x2=[1 1 1 1]; >>g=conv(x1,x2) g=
1 2 3 4 3 2 1 如果要绘出图形来,则利用stem命令,即
4
>>n=1:7;
>>stem(n,g,'fill'),grid on,xlabel('n') 程序运行结果如图1-4所示。
图1-4 卷积结果图
对于给定函数的卷积和,我们应计算卷积结果的起始点及其长度。两个时限序列的卷积和长度一般等于两个序列长度的和减1。
【实例1-3】 已知某系统的单位取样响应为h?n??0.8?u?n??u?n?8??,试
n用MATLAB求当激励信号为x(n)?u(n)?u(n?4)时,系统的零状态响应。
解:MATLAB中可通过卷积求解零状态响应,即x(n)*h(n)。由题意可知,描述h(n)向量的长度至少为8,描述x(n)向量的长度至少为4,因此为了图形完整美观,我们将h(n)向量和x(n)向量加上一些附加的零值。MATLAB源程序为
>>nx=-1:5; %x(n)向量显示范围(添加了附加的零值) >>nh=-2:10; %h(n)向量显示范围(添加了附加的零值) >>x=uDT(nx)-uDT(nx-4);
>>h=0.8.^nh.*(uDT(nh)-uDT(nh-8)); >>y=conv(x,h);
>>ny1=nx(1)+nh(1); %卷积结果起始点
>>%卷积结果长度为两序列长度之和减1,即0到(length(nx)+length(nh)-2) >>%因此卷积结果的时间范围是将上述长度加上起始点的偏移值 >>ny=ny1+(0:(length(nx)+length(nh)-2)); >>subplot(311)
5
>>stem(nx,x,'fill'),grid on >>xlabel('n'),title('x(n)') >>axis([-4 16 0 3]) >>subplot(312)
>>stem(nh,h','fill'),grid on >>xlabel('n'),title('h(n)') >>axis([-4 16 0 3]) >>subplot(313)
>>stem(ny,y,'fill'),grid on
>>xlabel('n'),title('y(n)=x(n)*h(n)') >>axis([-4 16 0 3]) 程序运行结果如图1-5所示。
图1-5 利用卷积和法求解系统的零状态响应
6
1.3 编程练习
1.
试用MATLAB命令求解以下离散时间系统的单位取样响应,并判断系统的稳定性。
(1)3y(n)?4y(n?1)?y(n?2)?x(n)?x(n?1) (2)2.
5y(n)?6y(n?1)?10y(n?2)?x(n) 278n已知某系统的单位取样响应为h?n??()?u?n??u?n?10??,试用MATLAB求当激励信号为x(n)?u(n)?u(n?5)时,系统的零状态响应。
3. 计算两个从n = 0开始的序列
x(n)?u(n)?any(n)?Aeu(n)的卷积和,其中和
A=40,a = 0.5。取50个样值点,作出序列x(n)、y(n)及卷积和f(n)的图像。
第2章 z变换及离散时间LTI系统的z域分析
2.1 实验目的
? 学会运用MATLAB求离散时间信号的z变换和z反变换; ? 学会运用MATLAB分析离散时间系统的系统函数的零极点;
? 学会运用MATLAB分析系统函数的零极点分布与其时域特性的关系; ? 学会运用MATLAB进行离散时间系统的频率特性分析。
2.2 实验原理及实例分析
2.2.1 z正反变换
序列x?n?的z变换定义为
X?z??Z?x?n???n????x?n?z??n (2-1)
其中,符号Z表示取z变换,z是复变量。相应地,单边z变换定义为
7