数字信号处理实验报告一 离散信号及其MATLAB实验.docVIP

南昌航空大学实验报告 2012 年 04 月 06 日 课程名称： 数字信号处理 实验名称： 离散信号及其MATLAB实现 班级： 090423班 学号： 姓名： 张 丽 实验一 离散信号及其MATLAB实验 实验目的 熟悉MATLAB的主要操作命令； 学会离散信号的表示方法及其基本运算； 掌握简单的绘图命令； 用Matlab编程并学会创建函数。 实验内容 序列的加、减、乘、除和乘方运算。输入A=[1 2 3 4 ]，B=[3 4 5 6]，起点n=0，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B，并用stem语句画出A，B，C，D，E，F，G。 用MATLAB实现下列序列： x(n)= 0≤n≤15 x(n)= 0≤n≤15 x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π) 0≤n≤15 编写函数如stepshift（n0，n1，n2）实现u(n)的移动序列u(n-n0)，n1≤n≤n2，给出该函数的图形。设n1=0，n2=9，n0=2 x(n)=[1，-1，3，5]，起点n=0，产生并绘出下列序列的样本： x1(n)=2x(n+2)-x(n-1)-2x(n) 三、实验程序及实验图形 实验（1） 程序： n=[0:3]; A=[1 2 3 4 ]; %定义序列A，即一行四列的矩阵 B=[3 4 5 6]; %定义序列B C=A+B;D=A-B;E=A.*B;F=A./B;G=A.^B; figure(1); subplot(2,2,1) %将图形界面分成2行2列，其中第一个显示下列图形 stem(n,A,r*); %绘制二维离散数据的火柴杆图，用红线和*号标出 xlabel(n); %x轴标签为n ylabel(A); %y轴标签为A grid on; %绘制网格 subplot(2,2,2) stem(n,B);xlabel(n);ylabel(B);grid on; subplot(2,2,3) stem(n,C);xlabel(n);ylabel(C);grid on;title(序列的运算); subplot(2,2,4); stem(n,D);xlabel(n);ylabel(D);grid on; figure(2); subplot(3,1,1) stem(n,E);xlabel(n);ylabel(E);grid on;title(序列的运算); subplot(3,1,2) stem(n,F);xlabel(n);ylabel(F);grid on; subplot(3,1,3) stem(n,G);xlabel(n);ylabel(G);grid on; 运行结果 图1 A、B、C、D的波形 图2 E、F、G的波形 实验结果分析： 由图1和图2可知：序列的加、减、乘、除是在n上的对应点的加、减、乘、除，.*、./、.^也是矩阵的对应点相乘、除、乘方，对离散序列进行运算可以看作是对两个行向量的运算。 实验（2） 程序： clear; close all; %初始化 n=[0:15]; x1=(0.8).^n; x2=exp((0.2+3j).*n); x3=3*cos(0.175*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi); % 定义x1,x2.x3 figure(1); subplot(3,1,1); stem(n,x1); xlabel(n);ylabel(x1); grid on; title(序列的表示); %显示图形的名字 subplot(3,1,2); stem(n,x2); %画火柴图 xlabel(n);ylabel(x2); grid on; subplot(3,1,3);stem(n,x3);xlabel(n);ylabel(x3);grid on; 2、运行结果： 图3 序列的基本运算 3、实验结果分析： 从图3上可以看出：离散序列可以看作是连续函数在一些等间隔点上的采样，即在连续函数的基础上对图形进行离散化，使得原来可以取定义域内任意值的自变量t转换成只能取定义域内的整数的自变量n，如此进行抽样，但该离散序列图形的包络和与之对应的连续函数的波形是一致的。 实验（3） 程序 % 方法一 下标为0--9 clear;close all; n1=0;n2=9;n0=2; N=n2-n1+1; x=ones(1,N); %生成一个1行N列的单位矩阵 [n,x1]=stepshift(n0,n1,n2); %移位 figure(1); subplot(2,1,1);