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

实验报告课程名称： 数字信号处理 指导老师： 刘英 成绩：_实验名称： 有限长序列、频谱、DTFT的性质 实验类型： 演示 同组学生姓名： 一、实验目的和要求通过演示实验，建立对典型信号及其频谱的直观认识，理解DTFT的物理意义、主要性质。二、实验内容和步骤2-1用MATLAB，计算得到五种共9个序列： 2-1-1单位阶跃序列 x(n)=2u(n-5)2-1-2实指数序列例如，a=0.5, length=10 a=0.9, length=10 a=0.9, length=202-1-3复指数序列例如，a=0.4, b=0.6, length=202-1-4从正弦信号x(t)=sin(2pft+delta)抽样得到的正弦序列x(n)=sin(2pfnT+delta)。如，信号频率f=1Hz，初始相位delta=0，抽样间隔T=0.1秒，序列长length=10。2-1-5含两个频率分量的复合函数序列x(n)=sin(2pf1nT)+deltasin(2pf2nT+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，对下面系统进行分析因果系统, 求（1） 求出, 画出零极点示意图（2） 画出和三、主要仪器设备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 a0);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;%可修改参数f2=3;%可修改参数delta=0.5;%可修改参数phi=180;%可修改参数T=0.1;%可修改参数N=10;%可修改参数xn=0:T:T*(N-1);%横坐标figure(1);yn=sin(2*pi*f1*xn)+delta*sin(2*pi*f2*xn+phi);yn0=real(yn);%实部yn1=imag(yn);%虚部yn2=abs(yn);%模yn3=angle(yn);%角subplot(2,2,1);stem(xn,yn0,filled);title(real);subplot(2,2,2);stem(xn,yn1,filled);title(imag);subplot(2,2,3);stem(xn,yn2,filled);title(abs);subplot(2,2,4);stem(xn,yn3,filled);title(angle);六、实验结果与分析观察实验结果（数据及图形）的特征，做必要的记录，做出解释。包括：6-1 各种序列的图形（时域）和频谱（频域）各有何