有限长序列、频谱、DF.doc

实验报告 课程名称： 数字信号处理 指导老师： 成绩：__________________ 实验名称： 有限长序列、频谱、DFT的性质 实验类型：___演示_______同组学生姓名： —— 一、实验目的和要求 设计通过演示实验，建立对典型信号及其频谱的直观认识，理解DFT的物理意义、主要性质。 二、实验内容和步骤 2-1用MATLAB，计算得到五种共9个序列： 2-1-1实指数序列 例如，a=0.5, length=10 a=0.9, length=10 a=0.9, length=20 2-1-2复指数序列例如，a=0.5, b=0.8, length=10 2-1-3从正弦信号x(t)=sin(2(ft+delta)抽样得到的正弦序列x(n)=sin(2(fnT+delta)。如，信号频率f=1Hz，delta=0，抽样间隔T=0.1秒，序列长length=10。 2-1-4从余弦信号x(t)=cos(2(ft + delta)抽样得到的余弦序列x(n)=cos(2(fnT + delta)。如，信号频率f=1Hz，delta=0，抽样间隔T=0.1秒，序列长length=10。 2-1-5含两个频率分量的复合函数序列x(n)=sin(2(f1nT)+delta×sin(2(f2nT+phi)。如， 频率f1 (Hz) 频率f2 (Hz) 相对振幅 delta 初相位phi (度) 抽样间隔T (秒) 序列长 length 1 3 0.5 0 0.1 10 1 3 0.5 90 0.1 10 1 3 0.5 180 0.1 10 2-2 用MATLAB，对上述各个序列，重复下列过程。 2-2-1画出一个序列的实部、虚部、模、相角；观察并记录实部、虚部、模、相角的特征。 2-2-2 计算该序列的幅度谱、频谱实部、频谱虚部；观察和并记录它们的特征，给予解释。 2-2-3 观察同种序列取不同参数时的频谱，发现它们的差异，给予解释。 三、主要仪器设备 MATLAB编程。 四、操作方法和实验步骤 （参见“二、实验内容和步骤”）  五、实验数据记录和处理 （一）实指数序列 （1） a=0.5, length=10 %program 2.1.1a clear;clf;clc;%清除缓存 n=0:9;%设置区间 xn=((0.5).^n).*(0=nn=9); xw=dftmtx(10)*xn; %用DFT求频谱 f=n/10.*(0=nn=5)+(10-n)/10.*(6=nn=9); %求出对应频率 figure(1); %画出序列的实部、虚部、模、相角 subplot(2,2,1);stem(n,real(xn)); xlabel(n);ylabel(real(xn));title(序列的实部); subplot(2,2,2);stem(n,imag(xn)); xlabel(n);ylabel(imag(xn));title(序列的虚部); subplot(2,2,3);stem(n,abs(xn)); xlabel(n);ylabel(abs(xn));title(序列的模); subplot(2,2,4);stem(n,angle(xn)); xlabel(n);ylabel(angle(xn));title(序列的相角); figure(2); %画出序列的幅度谱、频谱实部、频谱虚部 subplot(3,1,1);stem(f,abs(xw)); xlabel(f/Hz);ylabel(abs(xw));title(序列的幅度谱); subplot(3,1,2);stem(f,real(xw)); xlabel(f/Hz);ylabel(real(xw));title(频谱实部); subplot(3,1,3);stem(f,imag(xw)); xlabel(f/Hz);ylabel(imag(xw));title(频谱的虚部); (2) a=0.9, length=10 %program 2.1.1B clear;clf;clc;%清除缓存 clear n=0:9; xn=((0.9).^n).*(0=nn=9); xw=dftmtx(10)*xn; %用DFT求频谱 f=n/10.*(0=nn=5)+(10-n)/10.*(6=nn=9); %求出对应频率 figure(1); %画出序列的实部、虚部、模、相角 subplot(2,2,1);stem(n,real(xn)); xlabel(n);ylabel(real(xn));title(序列的实部); subplot(2,2,2);st