数字信号处理课程设计实验报告..docVIP

基础实验 实验一 离散时间系统及离散卷积 实验原理 利用Matlab软件计算出系统函数的零极点分布、单位脉冲响应和系统频率响应等的图像并于笔算结果进行比较，找出异同。编译合适程序能计算取值范围不同的离散卷积。 实验目的 （1）熟悉MATLAB软件的使用方法。 （2）熟悉系统函数的零极点分布、单位脉冲响应和系统频率响应等概念。 （3）利用MATLAB绘制系统函数的零极点分布图、系统频率响应和单位脉冲响应。 （4）熟悉离散卷积的概念，并利用MATLAB计算离散卷积。 三、实验源程序及实验结果 1.离散时间系统的单位脉冲响应 function pr1() %定义函数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函数 n=[-20:120]; %定义n从-20 到120 h=filter(b,a,x); %输入x为滤波前序列，b为分子，a为分母 figure(1) %绘图figure 1 (冲激响应) stem(n,h); %在图中绘出冲激 title(冲激响应); %定义标题为:冲激响应 xlabel(n); %绘图横座标为n ylabel(h(n)); %绘图纵座标为h(n) figure(2) %绘图figure 2 [z,p,g]=tf2zp(b,a); %绘出零极点图 zplane(z,p) function [x,n]=impseq(n0,n1,n2) %声明impseq函数 n=[n1:n2]; x=[(n-n0)==0]; 离散系统的幅频、相频的分析方法 function pr2() 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=500;k=1:K;w=pi*k/K; %角频率w H=(b*exp(-j*m*w))./(a*exp(-j*l*w));%对系统函数的定义 magH=abs(H); %magH为幅度 angH=angle(H); %angH为相位 figure(1) subplot(2,1,1); %在同一窗口的上半部分绘图 plot(w/pi,magH); %绘制w(pi)-magH的图形 grid; axis([0,1,0,1]); %限制横纵座标从0到1 xlabel(w(pi)); %x座标为 w(pi) ylabel(|H|); %y座标为 angle(H) title(幅度，相位响应); %图的标题为:幅度，相位响应 subplot(2,1,2); %在同一窗口的下半部分绘图 plot(w/pi,angH); %绘制w(pi)-angH的图形 grid; %为座标添加名称 xlabel(w(pi)); %x座标为 w(pi) ylabel(angle(H)); %y座标为 angle(H) 离散卷积的计算 function pr3() 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=u1-u2; %x(