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

沁*丿、象实验报告 课程名称： 数字信号处理 指导老师： 英 成绩： 实验名称： 有限长序列、频谱、 DTFT的性质 实验类型: 演示 同组学生: — 、实验目的和要求 通过演示实验，建立对典型信号及其频谱的直观认识，理解 DTFT的物理意义、主要性质。 、实验容和步骤 2-1用MATLAB，计算得至U五种共 9个序列: 2-1-1单位阶跃序列 x( n)=2u( n-5) 2-1-2实指数序列x(n) 0 n len gth 1 otherwise 例如，a=0.5, length=1O a=0.9, le ngth=10 a=0.9, le ngth=20 2-1-3复指数序列x(n) (a jb)n 0 n len gth otherwise 1 例如, a=0.4, b=0.6, le ngth=20 2-1-4从正弦信号x(t)=sin(2 ft+delta)抽样得到的正弦序列 x(n)=sin(2 fnT+delta)。女口，信号频率 f=1Hz，初始相位delta=0，抽样间隔 T=0.1秒，序列长length=10 。 2-1-5含两个频率分量的复合函数序列 x(n)=sin(2 f1nT)+delta Xsin(2 f2nT+phi)。如, 频率f1 频率f2 相对振幅 初相位phi 抽样间隔T 序列长 (Hz) (Hz) delta (度) 渺) len gth 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观察同种序列取不同参数时的频谱，发现它们的差异，给予解释。 实验名称： : 2.3用MATLAB，对下面系统进行分析 i 因果系统 y(n) y(n 1) y(n 2) x(n)，求 6 求出H(z)，画出零极点示意图 画出 H(ejw)和 H(ejw) 三、 主要仪器设备 MATLAB 编程。 四、 操作方法和实验步骤 (参见“二、实验容和步骤”) 五、 实验数据记录和处理 列出MATLAB程序清单，加注释。 2-1 2-1-1 N=15; x=on es(1,N); %生成行矩阵 xn=5:N-1+5; %黄坐标 figure(1); xlabel(x (n)); ylabel( n); stem(x n,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-2 N=10; 涮修改参数 a=0.5; %可修改参数 xn=0:N-1; %横坐标 figure(1); xlabel(x(n)); ylabel( n); stem(x n, a.A(x n)); %画指数函数axis on; axis([-3 N+3 0 aA0]); set(gca,XTick,[-1:1:N+2]); title(a=0.5 length=10); 2-1-3 N=20; 涮修改参数 a=0.4; 涮修改参数 b=0.6; 涮修改参数 xn=0:N-1; %黄坐标 figure(1); xlabel(x (n)); ylabel( n); stem(x n,(a+j*b).A(x n)); %画复指数函数 axis on; axis([-3 N+3 -0.5 (a+j*b)A0]); set(gca,XTick,[-1:1:N+2]); title(a=0.4 b=0.6 len gth=20); 2-1-4 f=1; 9可修改参数 delta=0; %可修改参数 T=0.1; %可修改参数 N=10; %可修改参数 xn=0:T:T*(N-1); %黄坐标 figure(1); stem(x n,sin (2*pi*f*x n+delta)); %画采样后的结果 hold on; %亭留 t=-0.5:0.001:1.5; %画原始函数 plot(t,si n( 2*pi*f*t)); axis on; axis([-0.5 1.5 -1.2 1.2]); grid on; title(Sample); 2-1-5 f1=1; %可修改参数 f2=3; %可修改参数 delta=0.5; %可修改参数 phi=0; %可修改参数 T=0.1; %可修改参数 N=10; %可修改参数 xn=0:T:T*(N-1); %黄坐标 figure(1); s