资源描述：
《离散时间系统及离散卷积》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、实验一、离散时间系统及离散卷积1、单位脉冲响应源程序：functionpr1()%定义函数pr1a=[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到120h=filter(b,a,x);%调用函数给纵坐标赋值figure(1)%绘图figure1(冲激响应)stem(n,h);%在图中绘出冲激title('单位冲激响应（耿海锋）');%定义标题为:'冲激响应（
2、耿海锋）'xlabel('n');%绘图横座标为nylabel('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:1、离散系统的幅频、相频的分析源程序：functionpr2()b=[0.0181,0.0543,0.0543,0.0181];a=[1.000,-1.76,1.182
3、9,-0.2781];m=0:length(b)-1;%m的范围，从0到3l=0:length(a)-1;%l的范围，从0到3K=5000;k=1:K;w=pi*k/K;%角频率wH=(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到1xlabel('w(pi
4、)');%x座标为w(pi)ylabel('
5、H
6、');%y座标为angle(H)-耿海锋title('幅度，相位响应（耿海锋）');%图的标题为:'幅度，相位响应（耿海锋）'plot(w/pi,angH);%绘制w(pi)-angH的图形grid;%为座标添加名称xlabel('w(pi)');%x座标为w(pi)ylabel('angle(H)');%y座标为angle(H)结果：Figure1Figure21、卷积计算源程序：functionpr3()n=-5:50;%声明n的范围，从-5到50u1=stepseq(0,-5,50);%调
7、用stepseq函数声用明u1=u(n)u2=stepseq(10,-5,50);%调用stepseq函数声用明u2=u(n-10)%输入x(n)和冲激响应h(n)x=u1-u2;%x(n)=u(n)-u(n-10)h=((0.9).^n).*u1;%h(n)=0.9^n*u(n)figure(1)subplot(3,1,1);%绘制第一个子图stem(n,x);%绘制图中的冲激axis([-5,50,0,2]);%限定横纵座标的范围title('输入序列-52101702-耿海锋');%规定标题为:'输入序列-52101702-耿海锋'xla
8、bel('n');%横轴为nylabel('x(n)');%纵轴为x(n)subplot(3,1,2);%绘制第二个子图stem(n,h);%绘制图中的冲激axis([-5,50,0,2]);%限定横纵座标的范围title('冲激响应序列-52101702-耿海锋');%规定标题为:'冲激响应序列-52101702-耿海锋'xlabel('n');%横轴为nylabel('h(n)');%纵轴为h(n)%输出响应[y,ny]=conv_m(x,n,h,n);%调用conv_m函数subplot(3,1,3);%绘制第三个子图stem(ny,y)
9、;axis([-5,50,0,8]);title('输出响应-52101702-耿海锋');%规定标题为:'输出响应-52101702-耿海锋'xlabel('n');ylabel('y(n)');%纵轴为y(n)%stepseq.m子程序%实现当n>=n0时x(n)的值为1function[x,n]=stepseq(n0,n1,n2)n=n1:n2;x=[(n-n0)>=0];%con_m的子程序%实现卷积的计算function[y,ny]=conv_m(x,nx,h,nh)nyb=nx(1)+nh(1);nye=nx(length(x))+
10、nh(length(h));ny=[nyb:nye];y=conv(x,h);结果：实验二、离散傅立叶变换与快速傅立叶变换1、离散傅立叶变换（DFT）
