实验一：离散时间序列卷积和MATLAB实现

实验一：离散时间序列卷积和MATLAB实现

（一）实验目的

1、掌握离散卷积计算方法；

2、学会线性常系数差分方程的迭代解法； 3、学会针对具体系统设计程序；

（二）实验内容

1、用MATLAB画出如下表达式的脉冲序列

f(n)?8.0?(n)?3.4?(n?1)?1.8?(n?2)?5.6?(n?3)?2.9?(n?4)?0.7?(n?5)

例如：x(1)=8.0; n=0:5;

fn=[8.0 3.4 1.8 5.6 2.9 0.7]; stem(n,fn); xlabel('n'); ylabel('fn');

876543210fn00.511.522.5n33.544.55

2、 计算卷积

用MATLAB计算序列{1 2 3 3 2 1 1}和序列{-1 2 0 1}的离散卷积。 首先用手工计算，然后用MATLAB函数conv(x,y)编程验证。 MATLAB部分参考程序如下： a=[1 2 3 3 2 1 1]; c=conv(a,b); M=length(c)-1; stem(n,c);

xlabel('n'); ylabel('幅度');

subplot(3,1,1); a=[1 2 3 3 2 1 1]; stem(a); subplot(3,1,2)

b=[-1 2 0 1]; stem(b); subplot(3,1,3); c=conv(a,b); M=length(c)-1;

幅度3210210-11050-512345n67891011.522.533.541234567stem(c);

xlabel('n'); ylabel('幅度');

3、参考教材P14页，使用通用卷积convu函数，计算如下卷积

已知h(n)=R5(n+2)，x(n)=3R4(n-3)，计算y(n)=x(n)*h(n)编程实现并在一幅图内分别画出h(n)，x(n)，y(n)，标好横坐标与纵坐标。

subplot(3,1,1); h=ones(1,5);nh=-2;2; stem(h);

xlabel('n'); ylabel('幅度'); subplot(3,1,2);

幅度10.50x=3*ones(1,4);nx=3:6; stem(x);

xlabel('n'); ylabel('幅度'); subplot(3,1,3); [y,ny]=convu(h,nh,x,nx); stem(y);

xlabel('n'); ylabel('幅度');

11.522.53n3.544.554幅度2011.522.5n33.5420幅度1001234n5678

4、编写教材P31第15题的MATLAB程序。 a=0.5;ys=0;

xn=[1 2 3 0 4 2 1];

B=[1 0.5 3 0 0.25];A=[1 0.5]; xi=filtic(B,A,ys); yn=filter(B,A,xn,xi); n=0:length(yn)-1; stem(n,yn,'.')

title('(a)');xlabel('n');ylabel('y(n)')

5、编写教材P32第18题的MATLAB程序。

（三）心得体会及总结（要求每个人的总结都不一样）

(a)161412108y(n)6420-20123n456

本文来源：https://www.bwwdw.com/article/k8td.html