典型离散时间序列上机报告选编.docx

曾子叶典型离散时间序列上机报告 曾子叶 第  PAGE \* MERGEFORMAT 14 页 共  NUMPAGES \* MERGEFORMAT 16 页 实验目的 掌握各种常用的序列，理解其数学表达式和波形表示。 掌握在计算机中生成及绘制数字信号波形的方法。 掌握序列的相加、相乘、移位、反褶等基本运算及计算机实现与作用。 掌握线性卷积软件实现的方法。 掌握计算机的使用方法和常用系统软件及应用软件的使用。 通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。 实验原理 序列的基本概念 离散时间信号在数学上可用时间序列{x(n)}来表示，其中x(n)代表序列的第n个数字，n代表时间的序列，n的取值范围为R的整数，n取其它值x(n)没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)x(t)进行等间隔采样，采样间隔为T，得到{x(nt)}一个有序的数字序列就是离散时间信号，简称序列。 常用序列 常用序列有：单位脉冲序列（单位抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。 序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 主要实验仪器及材料 微型计算机、Matlab、TC编程环境。 实验内容 知识准备 认真复习以上基础理论，理解本实验所用到的实验原理。 离散时间信号（序列）的产生 利用MATLAB编程产生下列有限长序列： 单位脉冲（抽样）序列 单位阶跃序列 矩形序列 正弦序列 实指数序列 复指数序列 Sinc函数 序列的运算 利用MATLAB编程实现离散时间信号的代数运算，移位运算，翻折运算，卷积运算。 音频实验 对给定两个音频进行叠加。 上机调试并打印或记录实验结果。 完成实验报告。 实验程序及结果 离散时间信号 %单位抽样序列 k=-30:30; delta=[zeros(1,30),1,zeros(1,30)]; figure(name,单位抽样序列); stem(k,delta,filled); %单位阶跃序列 uk=[zeros(1,30),1,ones(1,30)]; figure(name,单位阶跃序列); stem(k,uk,filled); %矩形序列 rk=[zeros(1,30),ones(1,20),zeros(1,11)]; figure(name,矩形序列); stem(k,rk,filled); %正弦序列 k=0:39; fk=sin(pi/6*k); figure(name,正弦序列); stem(k,fk,filled); %实指数序列 c=1; a=5/4; k=0:20; x=c*(a.^k); uk=[zeros(1,30),1,ones(1,30)]; figure(name,实指数序列); stem(k,x,filled); %复指数序列 r=0.8; w=pi/4; k=0:20; f=(r*exp(1i*w)).^k; xr=real(f); xi=imag(f); xa=abs(f); xw=angle(f); figure(name,复指数序列); subplot(2,2,1),stem(k,xr,filled),title(实部); subplot(2,2,2),stem(k,xi,filled),title(虚部); subplot(2,2,3),stem(k,xa,filled),title(模); subplot(2,2,4),stem(k,xw,filled),title(相角); %sinc函数 x=-30:0.1:30; y=sin(x)./x; figure(name,sinc函数); plot(x,y); 序列的运算 % %序列移位 k=-30:30; delta=[zeros(1,35),1,zeros(1,25)]; figure(name,序列移位); stem(k,delta,filled); % %序列的翻转 k=-30:30; uk=[ones(1,31),zeros(1,30)]; figure(name,序列的翻转); stem(k,uk,filled); % %序列相加 k=0:39; f1=0.3*sin(pi/6*k); f2=0.2*cos(pi/4*k); figure(name,序列相加); fa=f1+f2; subplot(2,2,1),stem(k,f1,filled),title(序列f1); subplot(2,2,2),stem(k,f2,filled),title(序列f2); subplot(2,2,3),stem(k,fa,fille