数字信号处理办法实验一离散时间信号分析.pptVIP

测控技术与仪器研究室 指导教师：林晓钢 xglin@ 测控技术与仪器研究室 实验一 离散时间信号分析 2010年11月 测控技术与仪器研究室 实验目的 1．掌握各种常用的序列，理解其数学表达式和波形表示。 2．掌握在计算机中生成及绘制数字信号波形的方法。 3．掌握序列的相加、相乘、移位、反转等基本运算及计算机实现与作用。 4．掌握线性卷积软件实现的方法。 5．掌握计算机的使用方法和相关应用软件的使用。 6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。 测控技术与仪器研究室 实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列{x(n)}来表示，其中x(n)代表序列的第n个数字，n代表时间的序列。注意： x(n)只 在 n 为 整 数 时 才 有 意 义, n 不 是 整 数 时 无 定 义, 但 不 能 认 为 是 0。 离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号xa(t)进行等间隔采样，采样间隔为T，得到一个有序的数字序列{xa(nT)}就是离散时间信号，简称序列。 测控技术与仪器研究室 实验原理 2．常用序列 单位脉冲序列（单位抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。 测控技术与仪器研究室 方法1 n1 = -5:5; x1=[(n1-0)==0]; %显然,用逻辑式是比较高明的方法 stem(n1,x1); %作用是绘制以n1为横坐标，x1为纵坐标的离散杆图 title(单位脉冲序列); %绘图标题 axis([-5,5,0,1]); %横坐标从-5到5，纵坐标从0到1 方法2 x=[1 zeros(1,N-1)]; % 产生长度为N的单位脉冲序列 实验原理 2．常用序列 单位脉冲序列（单位抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。 测控技术与仪器研究室 clear; n1 = -5:5; x1=[(n1-0)=0]; % 用逻辑式是比较高明的方法 stem(n1,x1); title(‘阶跃序列); axis([-5,5,0,1.1]); 实验原理 2．常用序列 测控技术与仪器研究室 clear; n1 = -5:5; x1=[(n1-0)=0]; x2=[(n1-3=0)]; x3=x1-x2; stem(n1,x3);title(矩形序列) axis([-5,5,0,1.1]) 实验原理 2．常用序列 ，式中ω0为数字频率 测控技术与仪器研究室 clear; n1 = 0:40; x1=exp((0.1+(pi/6)*i)*n1); subplot(2,1,1); stem(n1,real(x1));title(‘复指数序列’); ylabel(‘实部’); subplot(2,1,2); stem(n1,imag(x1)); ylabel(‘虚部’); 实验原理 2．常用序列 测控技术与仪器研究室 clear; n = 0:40; x1=3*cos(0.1*pi*n+pi/3); x2=2*sin(0.2*pi*n); subplot(2,1,1); stem(n,x1);title(余弦序列); subplot(2,1,2); stem(n,x2);title(正弦序列); 实验原理 2．序列的基本运算 2.1 序列的加法 ：z(n)=x (n)+y(n) 加法用于序列的合成。它把两个序列中位置序号相同的样本相加，形成新的样本序列。 测控技术与仪器研究室 function [y,n] = seqadd(x1,n1,x2,n2) % [y,n] = seqadd(x1,n1,x2,n2) % ----------------------------- % 实现y(n) = x1(n)+x2(n) % y = 在包含n1和n2的n点上求序列和, % x1 = 在位置向量n1上的第一序列 % x2 = 在位置向量n2上的第二序列(n2可与 n1不同) n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的长度 y1 = zeros(1,length(n)); y2 = y1; % 初始化 y1(find((n=min(n1))(n=max(n1))))=x1; % 具有y的长度的x1 y2(find((n=min(n2))(n=max(n2))))=x2; % 具有y的长度的x2