1、实验一序列频谱DFT的性质专业:_信息与通信工程姓名:_ _学号_日期:_地点_实验报告课程名称: 数字信号处理 指导老师: 成绩:_实验名称: FIR序列、频谱、DFT的性质 实验类型:_演示_同组学生姓名: 一、实验目的和要求设计通过演示实验,建立对典型信号及其频谱的直观认识,理解DFT的物理意义、主要性质。二、实验内容和步骤2-1用MATLAB,计算得到五种共9个序列:2-1-1实指数序列例如,a=0.5, length=10 a=0.9, length=10 a=0.9, length=202-1-2复指数序列例如,a=0.5, b=0.8, length=102-1-3从正弦信号x(
2、t)=sin(2ft+delta)抽样得到的正弦序列x(n)=sin(2fnT+delta)。如,信号频率f=1Hz,初始相位delta=0,抽样间隔T=0.1秒,序列长length=10。2-1-4从余弦信号x(t)=cos(2ft + delta)抽样得到的余弦序列x(n)=cos(2fnT + delta)。如,信号频率f=1Hz,初相位delta=0,抽样间隔T=0.1秒,序列长length=10。2-1-5含两个频率分量的复合函数序列x(n)=sin(2f1nT)+deltasin(2f2nT+phi)。如,频率f1(Hz)频率f2(Hz)相对振幅delta初相位phi (度)抽样间
3、隔T(秒)序列长length130.500.110130.5900.110130.51800.1102-2 用MATLAB,对上述各个序列,重复下列过程。2-2-1画出一个序列的实部、虚部、模、相角;观察并记录实部、虚部、模、相角的实验名称:FIR序列、频谱、DFT的性质 姓名:_ 特征。2-2-2 计算该序列的幅度谱、频谱实部、频谱虚部;观察和并记录它们的特征,给予解释。2-2-3 观察同种序列取不同参数时的频谱,发现它们的差异,给予解释。三、主要仪器设备MATLAB编程。四、操作方法和实验步骤(参见“二、实验内容和步骤”)五、实验数据记录和处理列出MATLAB程序清单,加注释。2-1-1a
4、 (a=0.5, length=10)程序n=0:9;xn=(0.5).n).*(0=n&n=9);xw=dftmtx(10)*xn; %用DFT求频谱f=n/10.*(0=n&n=5)+(10-n)/10.*(6=n&n=9); %求出对应频率figure(1); %画出序列的实部、虚部、模、相角subplot(2,2,1);stem(n,real(xn);xlabel(n);ylabel(real(xn);subplot(2,2,2);stem(n,imag(xn);xlabel(n);ylabel(imag(xn);subplot(2,2,3);stem(n,abs(xn);xlabel
5、(n);ylabel(abs(xn);subplot(2,2,4);stem(n,angle(xn);xlabel(n);ylabel(angle(xn);figure(2); %画出序列的幅度谱、频谱实部、频谱虚部subplot(3,1,1);stem(f,abs(xw); xlabel(f/Hz);ylabel(abs(xw);实验名称:FIR序列、频谱、DFT的性质 姓名:_ 3subplot(3,1,2);stem(f,real(xw);xlabel(f/Hz);ylabel(real(xw);subplot(3,1,3);stem(f,imag(xw);xlabel(f/Hz);yl
6、abel(imag(xw);2-1-1b(a=0.9, length=10)程序n=0:9;xn=(0.9).n).*(0=n&n=9);xw=dftmtx(10)*xn; %用DFT求频谱f=n/10.*(0=n&n=5)+(10-n)/10.*(6=n&n=9); %求出对应频率figure(1); %画出序列的实部、虚部、模、相角subplot(2,2,1);stem(n,real(xn);xlabel(n);ylabel(real(xn);subplot(2,2,2);stem(n,imag(xn);xlabel(n);ylabel(imag(xn);subplot(2,2,3);st
7、em(n,abs(xn);xlabel(n);ylabel(abs(xn);subplot(2,2,4);stem(n,angle(xn);xlabel(n);ylabel(angle(xn);figure(2); %画出序列的幅度谱、频谱实部、频谱虚部subplot(3,1,1);stem(f,abs(xw); xlabel(f/Hz);ylabel(abs(xw);subplot(3,1,2);stem(f,real(xw);xlabel(f/Hz);ylabel(real(xw);subplot(3,1,3);stem(f,imag(xw);xlabel(f/Hz);ylabel(ima
8、g(xw);2-1-1cn=0:19;xn=(0.9).n).*(0=n&n=19);实验名称:FIR序列、频谱、DFT的性质 姓名:_邵振江_学号_3080102350_P. 4xw=dftmtx(20)*xn; %用DFT求频谱f=n/20.*(0=n&n=10)+(20-n)/20.*(11=n&n=0&n=9);f=n/10.*(0=n&n=5)+(10-n)/10.*(6=n&n=0&n=9);xw=dftmtx(10)*xn; %用DFT求频谱f=n.*(0=n&n=5)+(10-n).*(6=n&n=0&n=9);xw=dftmtx(10)*xn; %用DFT求频谱f=n.*(0=n&n=5)+(10-n).*(6=n&n=9); %求出对应频率figure(1); %画出序列的实部、虚部、模、相角subplot(2,2,1);stem(n,real(xn);xlabel(n);ylabel(real(xn);subplot(2,2,2);stem(n,imag(xn);xlabel(n);ylabel(imag(xn);subplot(2,2,3);stem(n,abs(xn);xlabel(n);ylabel(abs(xn);subplot(2,2,4);stem(n,angle(xn);xlabel(n);y
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1