基于DSP的数字信号实验一.docVIP

一、实验目的 1、掌握用MATLAB产生基本离散时间信号的方法； 2、掌握离散线性系统的差分方程表示法； 3、加深对单位脉冲响应和离散卷积的理解； 4、掌握离散系统的频率响应分析方法； 5、加深对零极点分布的理解。 二、实验原理 信号的分析处理及MATLAB软件的应用 三、实验内容 1、用MATLAB产生单位脉冲序列、单位阶跃序列、矩形序列、实指数序列(底数a＝0.5)和正弦序列(ω0＝π/8)，并用stem函数分别做出各序列的图形(-4≤n≤16)。 2、求以下差分方程所描述系统的单位脉冲响应h(n)，0≤n≤30 3、计算序列x(n)和y(n)的离散卷积，并作图表示卷积结果： ； 4、求以下因果系统的频率响应及零极点分布图： 四、实验结果与分析 (1)、单位脉冲序列 程序： n=-4:1:16; x=[0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; stem(n,x); xlabel(n); ylabel(x(n)); grid; set(gca,xtick,n); 图像： （2）、单位阶跃序列 程序： n=-4:1:16; x=[0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; stem(n,x); xlabel(n); ylabel(x(n)); grid; set(gca,xtick,n); 图像： （3）、矩形序列 程序： n=-4:1:16; x=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; stem(n,x); xlabel(n); ylabel(x(n)); grid; set(gca,xtick,n); 图像： （4）、实指数序列(底数a＝0.5) 程序： n=-4:1:16; x=0.5.^n; stem(n,x); xlabel(n); ylabel(x(n)); grid; set(gca,xtick,n); 图像： （5）、正弦序列(ω0＝π/8) 程序： n=-4:1:16; x=sin(pi/8*n); figure; stem(n,x); xlabel(n); ylabel(x); title(sin(pi/8*n)); 图像： 2、程序： N=31; a=[1 0.1 -0.01]; b=[1 -0.2]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(b,a,x); stem(k,y); xlabel(n); ylabel(y(n)); 图像： 3、程序： x=[1 2 3 6 10]; y=[4 1 -3 -2]; z=conv(x,y); n=0:(length(x)+length(y)-2); stem(n,z); xlabel(n); ylabel(z(n)); 图像： 4、程序： a=[1 0 0.3333]; b=[0.3333 0.6667 0.3333]; zplane(b,a); k=256; w=0:pi/k:2*pi; H=freqz(b,a,w); figure; subplot(2,1,1); plot(w,abs(H)); %plot(w,20*log10(abs(H)));%对数谱 grid; xlabel(\omega); ylabel(幅度); title(幅度响应); subplot(2,1,2); plot(w,angle(H)); grid; xlabel(\omega); ylabel(弧度); title(相位响应); 图像： 五、实验小结 实验进行的还算比较顺利，因为有范文，而其中的原理虽然能够看明白、弄清楚，但是如果全部自己写，估摸着够呛，须得多加练习。