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

《有限长序列频谱DTFT的性质.doc》由会员分享，可在线阅读，更多相关《有限长序列频谱DTFT的性质.doc（4页珍藏版）》请在优知文库上搜索。

1、-实验报告一、实验目的和要求通过演示实验，建立对典型信号及其频谱的直观认识，理解DTFT的物理意义、主要性质。二、实验容和步骤2-1用MATLAB，计算得到五种共9个序列： 2-1-1单位阶跃序列 *(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从正弦信号*(t)=sin(2pft+delta)抽样得到的正弦序列*(n)=sin(2pfnT+delta)。如，信号频率f=1Hz，初始相位delta=0，抽样间隔T

2、=0.1秒，序列长length=10。2-1-5含两个频率分量的复合函数序列*(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 观察同种序列取不同参数时的频谱，

3、发现它们的差异，给予解释。实验名称：_：_*：_ 2.3 用MATLAB，对下面系统进展分析因果系统, 求（1） 求出, 画出零极点示意图（2） 画出和三、主要仪器设备MATLAB编程。四、操作方法和实验步骤参见“二、实验容和步骤五、实验数据记录和处理列出MATLAB程序清单，加注释。2-12-1-1N=15;*=ones(1,N);%生成行矩阵*n=5:N-1+5;%横坐标figure(1);*label(*(n);ylabel(n);stem(*n,2*);%纵坐标变为2倍a*is on;a*is(-1 14 0 4);set(gca,YTick,0:1:4);set(gca,*Tick,

4、-1:1:14);title(单位阶跃序列);2-1-2N=10;%可修改参数a=0.5;%可修改参数*n=0:N-1;%横坐标figure(1);*label(*(n);ylabel(n);stem(*n,a.(*n);%画指数函数a*is on;a*is(-3 N+3 0 a0);set(gca,*Tick,-1:1:N+2);title(a=0.5 length=10);2-1-3N=20;%可修改参数a=0.4;%可修改参数b=0.6;%可修改参数*n=0:N-1;%横坐标figure(1);*label(*(n);ylabel(n);stem(*n,(a+j*b).(*n);%画复指

5、数函数a*is on;a*is(-3 N+3 -0.5 (a+j*b)0);set(gca,*Tick,-1:1:N+2);title(a=0.4 b=0.6 length=20);2-1-4f=1;%可修改参数delta=0;%可修改参数T=0.1;%可修改参数N=10;%可修改参数*n=0:T:T*(N-1);%横坐标figure(1);stem(*n,sin(2*pi*f*n+delta);%画采样后的结果hold on;%停留t=-0.5:0.001:1.5;%画原始函数plot(t,sin(2*pi*f*t);a*is on;a*is(-0.5 1.5 -1.2 1.2);grid

6、on;title(Sample);2-1-5f1=1;%可修改参数f2=3;%可修改参数delta=0.5;%可修改参数phi=0;%可修改参数T=0.1;%可修改参数N=10;%可修改参数*n=0:T:T*(N-1);%横坐标figure(1);stem(*n,sin(2*pi*f1*n)+delta*sin(2*pi*f2*n+phi);%画采样后的结果hold on;%停留t=-0.5:0.001:1.5;%画原始函数plot(t,sin(2*pi*f1*t)+delta*sin(2*pi*f2*t+phi);a*is on;a*is(-0.5 1.5 -1.2 1.2);grid on

7、;title(2-1-5(1);2-2-1(5)f1=1;%可修改参数f2=3;%可修改参数delta=0.5;%可修改参数phi=180;%可修改参数T=0.1;%可修改参数N=10;%可修改参数*n=0:T:T*(N-1);%横坐标figure(1);yn=sin(2*pi*f1*n)+delta*sin(2*pi*f2*n+phi);yn0=real(yn);%实部yn1=imag(yn);%虚部yn2=abs(yn);%模yn3=angle(yn);%角subplot(2,2,1);stem(*n,yn0,filled);title(real);subplot(2,2,2);stem(

8、*n,yn1,filled);title(imag);subplot(2,2,3);stem(*n,yn2,filled);title(abs);subplot(2,2,4);stem(*n,yn3,filled);title(angle);六、实验结果与分析观察实验结果数据及图形的特征，做必要的记录，做出解释。包括：6-1 各种序列的图形时域和频谱频域各有何特征，给予解释。2-12-1-12-1-22-1-32-1-42-1-52-2-1yn=sin(2*pi*f1*n)+delta*sin(2*pi*f2*n+phi);其中f1=1,f2=3,delta=0.5,phi=180实部、虚部、模、角. z.