《有限长序列、频谱、DTFT的性质》由会员分享，可在线阅读，更多相关《有限长序列、频谱、DTFT的性质（9页珍藏版）》请在金锄头文库上搜索。
1、课程名称：数字信号处理指导老师： 刘英 成绩:实验名称：有限长序列、频谱、DTFT的性质 实验类型:一同组学生姓名：一、实验目的和要求通过演示实验，建立对典型信号及其频谱的直观认识，理解DTFT的物理意义、主要性质。二、实验内容和步骤2-1用MATLAB，计算得到五种共9个序列：2-1-1 单位阶跃序列 x(n)=2u(n-5)2-1-2 实指数序列 x(n) = |an0 n 叽h-1例如，a=0.5 length=100otherwisea=0.9, length=10a=0.9, length=202-1-3 复指数序列 x(n) = | (a + 力)0 - n -length-1 例
2、如，a=0.4, b=0.6, length=200otherwise2-1-4从正弦信号x(/)=sin(2劝1+delta)抽样得到的正弦序列x(n)=sin(2物T+delta)。如，信号频率f=1Hz，初始相位delta=0，抽样间隔T=0.1秒，序列长length=10。2-1-5含两个频率分量的复合函数序列x(n)=sin(2硏nT)+deltax sin(2利nT+phi)。如,频率f(Hz)频率f2(Hz)相对振幅delta初相位phi(度)抽样间隔T (秒)序列长 length130.500.110130.5900.110130.51800.1102-2用MATLAB，对上述
3、各个序列，重复下列过程。2-2-1 画出一个序列的实部、虚部、模、相角；观察并记录实部、虚部、模、相角的特征。2-2-2 计算该序列的幅度谱、频谱实部、频谱虚部；观察和并记录它们的特征，给予解释2-2-3 观察同种序列取不同参数时的频谱，发现它们的差异，给予解释。实验名称：姓名：学号： 2.3用MATLAB，对下面系统进行分析51因果系统 y(n) + y(n -1) + y(n - 2) = x(n),求66(1) 求出H(z)，画出零极点示意图(2) 画出 H(ejw)和 ZH(ejw)三、主要仪器设备MATLAB 编程。四、操作方法和实验步骤(参见“二、实验内容和步骤”)五、实验数据记录
4、和处理列出 MATLAB 程序清单，加注释2-12-1-1%生成行矩阵%横坐标%纵坐标变为2倍N=15; x=ones(1,N); xn=5:N-1+5;figure(1); xlabel(x(n); ylabel(n);stem(xn,2*x);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);st em(xn,a/(xn);%画
5、指数函数axis on;axis(-3 N+3 0 aA0);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).A(xn);%画复指数函数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 length=20);f=1；%可修改参数delta=0;%可修改参数T=0.1;%可修改参数N=10;%可修改参数xn=0:T:T*(N-1);% 横
6、坐标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-4%可修改参数 %可修改参数 %可修改参数 %可修改参数 %可修改参数 %可修改参数 %横坐标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*
7、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);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
8、=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 各种序列的图形（时域）和频谱（频域）各有何特征，给予解释。2-12-1-12-1-2单位阶跃序列a = 0 5 length=Wa = 0.9 length=202-1-5yn=sin(2*pi*f1*xn)+delta*sin(2*pi*f2*xn+phi); 其中 f1=1,f2=3,delta=0.5,phi=180 实部、虚部、模、角
