离散信号和离散系统分析.docVIP

实验一 离散信号与离散系统分析 一、实验目的 1、掌握用MATLAB产生基本离散时间信号的方法； 2、掌握简单的绘图命令，学会用matlab创建函数； 3、加深对单位脉冲响应和离散卷积的理解； 4、掌握离散系统的频率响应分析方法； 5、加深对零极点分布的理解； 6、掌握离散线性系统的差分方程表示法。 二、实验原理 1、stem(n, x)：绘制离散序列x(n)。 2、y=conv(h, x): 计算h和x的线性卷积(离散卷积)，结果放在y中。 3、y=filter(b, a, x): 对输入序列x进行数字滤波，b和a分别是滤波器分子(对应输入信号x)和分母(对应输出信号y)的系数向量。 4、H=freqz(b, a, w): 计算数字滤波器的频率响应，w为频率向量，通常通过对2π取N等份得到。[H w]=freqz(b, a): 对0～π进行512点采样，采样频率输出在向量w中。如果省略左边的H或[H w]，则直接绘出滤波器的对数频谱和相频响应。 5、zplane(b, a)：绘制零极点分布图，b和a分别是滤波器分子(对应输入信号x)和分母(对应输出信号y)的系数向量。 三、实验内容 1、用MATLAB产生单位脉冲序列、单位阶跃序列、矩形序列(N＝4)、实指数序列(底数a＝0.5)和正弦序列(ω0＝π/8)，并用stem函数分别做出各序列的图形(-4≤n≤16)。 2、求以下差分方程所描述系统的单位脉冲响应h(n)，0≤n≤30 3、计算序列x(n)和y(n)的离散卷积，并作图表示卷积结果： ； 4、求以下因果系统的频率响应及零极点分布图： 四、实验结果与分析 第一题 代码： x =-4:16; y1 =[zeros(1,4) 1 zeros(1,16)];%-4到16的单位脉冲序列 y2 =[zeros(1,4) ones(1,17)];%单位阶跃序列 y3 =[zeros(1,4) ones(1,4) zeros(1,13)];%矩形脉冲 y4 = 0.5.^(x);%x实指数序列 y5 =sin(pi/8. *x);%正弦序列 figure; stem(x,y1); xlabel(n); ylabel(y1); title(单位脉冲序列); grid; figure; stem(x,y2); xlabel(n); ylabel(y2); title(单位阶跃序列); grid; figure; stem(x,y3); xlabel(n); ylabel(y3); title(矩形脉冲); grid; figure; stem(x,y4); xlabel(n); ylabel(y4); title(实指数序列); grid; figure; stem(x,y5); xlabel(n); ylabel(y5); title(正弦序列); grid; 图形： 分析： 产生单位脉冲序列和单位阶跃序列的脉冲式要学会用矩阵输入来简化带码；一些运算符号(如*和^)前面要加上“.”， 第二题： 代码： 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)); 图形： 分析：使用数字滤波函数filter时，要注意滤波器分子(对应输入信号x)和分母(对应输出信号y)的系数向量的输入顺序，注意不要颠倒。 第三题： 代码： 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)); 图形： 分析：使用线性卷积(离散卷积) 函数conv时，绘制卷积结果时要记得加上绘制的长度。 第四题： 代码： a=[0.3333 0.6667 0.3333]; b=[1 0 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(相位响应); 图形： 分析：使用系统的频率响应及零极点分布freqz和zplane函数时，要注意滤波器分子(对应输入信号x)和分母(对应输出信号y)的系数向量的输