离散时间系统及离散卷积.pdf

实验一、离散时间系统及离散卷积

1、单位脉冲响应

源程序：

functionpr1()%定义函数pr1

a=[1,-1,0.9];%定义差分方程y(n)-y(n-1)+0.9y(n-2)=x(n)

b=1;

x=impseq(0,-20,120);%调用impseq函数（matlab软件的函数库）

n=[-20:120];%定义n的范围，从-20到120

h=filter(b,a,x);%调用函数给纵坐标赋值

figure(1)%绘图figure1(冲激响应)

stem(n,h);%在图中绘出冲激

title(单位冲激响应（耿海锋）);%定义标题为:冲激响应（耿海锋）

xlabel(n);%绘图横座标为n

ylabel(h(n));%绘图纵座标为h(n)

figure(2)%绘图figure2

[z,p,g]=tf2zp(b,a);%绘出零极点图

zplane(z,p)

function[x,n]=impseq(n0,n1,n2)%声明impseq函数

n=[n1:n2];

x=[(n-n0)==0];

结果：

Figure1:

Figure2:

2、离散系统的幅频、相频的分析

源程序：

functionpr2()

b=[0.0181,0.0543,0.0543,0.0181];

a=[1.000,-1.76,1.1829,-0.2781];

m=0:length(b)-1;%m的范围，从0到3

l=0:length(a)-1;%l的范围，从0到3

K=5000;

k=1:K;

w=pi*k/K;%角频率w

H=(b*exp(-j*m*w))./(a*exp(-j*l*w));%对系统函数的定义

figure(1)

magH=abs(H);%magH为幅度

angH=angle(H);%angH为相位

plot(w/pi,magH-耿海锋);%绘制w(pi)-magH-耿海锋的图形

figure(2)

axis([0,1,0,1]);%限制横纵座标从0到1

xlabel(w(pi));%x座标为w(pi)

ylabel(|H|);%y座标为angle(H)-耿海锋

title(幅度，相位响应（耿海锋）);%图的标题为:幅度，相位响应

（耿海锋）

plot(w/pi,angH);%绘制w(pi)-angH的图形

grid;%为座标添加名称

xlabel(w(pi));%x座标为w(pi)

ylabel(angle(H));%y座标为angle(H)

结果：

Figure1

Figure2

3、卷积计算

源程序：

functionpr3()

n=-5:50;%声明n的范围，从-5到50

u1=stepseq(0,-5,50);%调用stepseq函数声用明u1=u(n)

u2=stepseq(10,-5,50);%调用stepseq函数声用明u2=u(n-10)

%输入x(n)和冲激响应h(n