有限长序列、频谱、DTFT的性质.docx

实验报告课程名称：数字信号处理指导老师：刘英成绩：__________________实验名称：有限长序列、频谱、DTFT的性质实验类型：演示同组学生姓名：——一、实验目的和要求通过演示实验，建立对典型信号及其频谱的直观认识，理解DTFT的物理意义、主要性质。二、实验内容和步骤2-1用MATLAB，计算得到五种共9个序列： 2-1-1单位阶跃序列 x(n)=2u(n-5)2-1-2实指数序列例如，a=0.5, length=10a=0.9, length=10a=0.9, length=202-1-3复指数序列例如，a=0.4, b=0.6, length=202-1-4从正弦信号x(t)=sin(2ft+delta)抽样得到的正弦序列x(n)=sin(2fnT+delta)。如，信号频率f=1Hz，初始相位delta=0，抽样间隔T=0.1秒，序列长length=10。2-1-5含两个频率分量的复合函数序列x(n)=sin(2f1nT)+delta×sin(2f2nT+phi)。如，频率f1(Hz)频率f2(Hz)相对振幅delta初相位phi (度)抽样间隔T(秒)序列长length130.500.110130.5900.110130.51800.1102-2用MATLAB，对上述各个序列，重复下列过程。2-2-1画出一个序列的实部、虚部、模、相角；观察并记录实部、虚部、模、相角的特征。2-2-2计算该序列的幅度谱、频谱实部、频谱虚部；观察和并记录它们的特征，给予解释。2-2-3观察同种序列取不同参数时的频谱，发现它们的差异，给予解释。实验名称：_______________________________姓名：______________学号：__________________ 2.3 用MATLAB，对下面系统进行分析因果系统, 求求出, 画出零极点示意图画出和三、主要仪器设备MATLAB编程。四、操作方法和实验步骤（参见“二、实验内容和步骤”）五、实验数据记录和处理列出MATLAB程序清单，加注释。2-12-1-1N=15;x=ones(1,N);%生成行矩阵xn=5:N-1+5;%横坐标figure(1);xlabel(x(n));ylabel(n);stem(xn,2*x);%纵坐标变为2倍axis on;axis([-1 14 0 4]);set(gca,YTick,[0:1:4]);set(gca,XTick,[-1:1:14]);title(单位阶跃序列);2-1-2N=10;%可修改参数a=0.5;%可修改参数xn=0:N-1;%横坐标figure(1);xlabel(x(n));ylabel(n);stem(xn,a.^(xn));%画指数函数axis on;axis([-3 N+3 0 a^0]);set(gca,XTick,[-1:1:N+2]);title(a=0.5 length=10);2-1-3N=20;%可修改参数a=0.4;%可修改参数b=0.6;%可修改参数xn=0:N-1;%横坐标figure(1);xlabel(x(n));ylabel(n);stem(xn,(a+j*b).^(xn));%画复指数函数axis on;axis([-3 N+3 -0.5 (a+j*b)^0]);set(gca,XTick,[-1:1:N+2]);title(a=0.4 b=0.6 length=20);2-1-4f=1;%可修改参数delta=0;%可修改参数T=0.1;%可修改参数N=10;%可修改参数xn=0:T:T*(N-1);%横坐标figure(1);stem(xn,sin(2*pi*f*xn+delta));%画采样后的结果hold on;%停留t=-0.5:0.001:1.5;%画原始函数plot(t,sin(2*pi*f*t));axis on;axis([-0.5 1.5 -1.2 1.2]);grid on;title(Sample);2-1-5f1=1;%可修改参数f2=3;%可修改参数delta=0.5;%可修改参数phi=0;%可修改参数T=0.1;%可修改参数N=10;%可修改参数xn=0:T:T*(N-1);%横坐标figure(1);stem(xn,sin(2*pi*f1*xn)+delta*sin(2*pi*f2*xn+phi));%画采样后的结果hold on;%停留t=-0.5:0.001:1.5;%画原始函数plot(t,sin(2*pi*f1*t)+delta*sin(2*pi*f2*t+phi));axis on;axis([-0.5 1.5 -1.2 1.2]);grid on;title(2-1-5(1));2-2-1(5)f1=1;%可修改参数f