《数字信号处理》实验及代码.doc

实验一 用Matlab进行数字信号处理 一、实验目的 1．掌握 MATLAB 基本的操作； 2．学习典型的离散信号的Matlab实现方法。 3．学习离散时间序列的基本运算：相加、相乘、移位等； 二、实验内容 1. 练习把y=sin(x), z=cos(x),u=2sin(x),v=sin(x)/cos(x)在[0，2pi]区间内的4个子图分别用不同的颜色、点型和线型绘制在同一个窗口中，并加上纵坐标、标题、图例和网格线。 2. 利用MATLAB编程产生和绘制下列有限长序列： 　　 (1)单位脉冲序列 　 　(2)单位阶跃序列 　　 (3)矩形序列 　　 (4)正弦序列 (5)实指数序列 3．上机调试并打印或记录实验结果。 4．完成实验报告。 三、实验结果 1. 实验程序如下： x=0:pi/10:2*pi; y=sin(x); plot(x,y,r*:); grid on %绘制网格线 hold on z=cos(x); plot(x,z,y+-); hold on u=2*sin(x); plot(x,u,bx-); hold on v=sin(x)/cos(x); plot(x,v,ko-); hold on xlabel(x 轴); %x轴注释 ylabel(y 轴); %y轴注释 legend({y=sin(x),z=cos(x), u=2sin(x) , v=sin(x)/cos(x) }); %图形注解 2. 程序如下： n=-20:20;n0=0;n1=8; w0=pi/5; w1=pi/3; x1=[(n-n0)==0]; x2=[(n-n0)=0]; x3=[(n-n0)=0 (n-n1)=0]; x4=sin(w0*n+w1); x5=0.9.^n.*x2; subplot(511);stem(n,x1);axis([ -20 20 0 2]);ylabel(\sigma(n)); subplot(512);stem(n,x2);axis([ -20 20 0 2]);ylabel(u(n)); subplot(513);stem(n,x3);axis([ -20 20 0 2]);ylabel(R8(n)); subplot(514);stem(n,x4);axis([ -20 20 -2 2]);ylabel( sin(w0n+w1) ); subplot(515);stem(n,x5);axis([ -20 20 0 2]);ylabel(0.9nu(n)); xlabel(n); 实验二 离散信号与系统 一、实验目的 1．掌握卷积定理、熟悉离散信号和系统的时域特性； 2．学习Matlab进行卷积计算方法； 3．学习Matlab求解差分方程 二、实验内容 1．离散信号的基本运算：对序列x(n)={2,3,4,1,2,5} ,n=0,1,2,3,4,5,的移位、乘法、加法、翻转及尺度变换。编写程序、并画出图形。 2．用matlab求解系统的单位脉冲响应、单位阶跃响应。编写程序、并画出图形。要求分别用filter、conv、impz三种函数完成。频率响应和零极点分布。 y(n)-0.4y(n-1)+0.5y(n-2)=x(n)-0.5x(n-1) 3．上机调试并打印或记录实验结果。 4．完成实验报告。 三、实验结果 1. 离散信号的基本运算：对序列x(n)={2,3,4,1,2,5} ,n=0,1,2,3,4,5,的移位、乘法、加法、翻转及尺度变换。编写程序、并画出图形。 程序如下： n=-10:10;k=2;N=length(n); x=zeros(size(n)); x1=zeros(size(n)); y=zeros(size(n)); y1=zeros(size(n)); y2=zeros(size(n)); y3=zeros(size(n)); y4=zeros(size(n)); x(11:16)=[2 3 4 1 2 5];x1(11:16)=[1 2 3 4 5 6]; y(1+k:N)=x(1:N-k); y1(1:N-k)=x(1+k:N); y2=x.*x1; y3=x+x1; y4=fliplr(x); y5=x(1:2:N);n1=fix(n(1:2:N)/2); subplot(421);stem(n,x);xlabel(n);ylabel(x(n)); subplot(422);stem(n,x1);xlabel(n);ylabel(x1(n)); subplot(423);stem(n,y);xlabel(n);ylabel(x(n-2)); subplot(424);stem(n,y1);xlabel(n);ylabe