资源描述:
《DSP实验1:有限长序列、频谱、DFT地性质》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实用标准文案课程名称:数字信号处理实验名称:有限长序列、频谱、DFT的性质一、实验目的和要求设计通过演示实验,建立对典型信号及其频谱的直观认识,理解DFT的物理意义、主要性质。二、实验内容和步骤2-1用MATLAB,计算得到五种共9个序列:2-1-1实指数序列例如,a=0.5,length=10a=0.9,length=10a=0.9,length=202-1-2复指数序列例如,a=0.5,b=0.8,length=102-1-3从正弦信号x(t)=sin(2pft+delta)抽样得到的正弦序列x(n)=
2、sin(2pfnT+delta)。如,信号频率f=1Hz,初始相位delta=0,抽样间隔T=0.1秒,序列长length=10。2-1-4从余弦信号x(t)=cos(2pft+delta)抽样得到的余弦序列x(n)=cos(2pfnT+delta)。如,信号频率f=1Hz,初相位delta=0,抽样间隔T=0.1秒,序列长length=10。2-1-5含两个频率分量的复合函数序列x(n)=sin(2pf1nT)+delta×sin(2pf2nT+phi)。如,频率f1(Hz)频率f2(Hz)相对振幅del
3、ta初相位phi(度)抽样间隔T(秒)序列长length130.500.110130.5900.110130.51800.1102-2用MATLAB,对上述各个序列,重复下列过程。2-2-1画出一个序列的实部、虚部、模、相角;观察并记录实部、虚部、模、相角的特征。2-2-2计算该序列的幅度谱、频谱实部、频谱虚部;观察和并记录它们的特征,给予解释。2-2-3观察同种序列取不同参数时的频谱,发现它们的差异,给予解释。文档大全实用标准文案三、主要仪器设备MATLAB编程。四、操作方法和实验步骤(参见“二、实验内容
4、和步骤”)五、实验数据记录和处理1.实指数序列1-1.a=0.5,length=10clc;clf;clear;%清除n=0:9;%设置自变量区间xn=0.5.^n;%计算相应的x(n)k=0:9;%设置DFT采样长度xw=dftmtx(10)*xn';%DFT变换figure(1);%画出原序列的实部、虚部、模、相角subplot(2,2,1);stem(n,real(xn),'filled');xlabel('x');ylabel('real(xn)');title('xn实部');subplot(2,
5、2,2);stem(n,imag(xn),'filled');xlabel('x');ylabel('imag(xn)');title('xn虚部');subplot(2,2,3);stem(n,abs(xn),'filled');xlabel('x');ylabel('abs(xn)');title('xn模');subplot(2,2,4);stem(n,angle(xn),'filled');xlabel('x');ylabel('angle(xn)');title('xn相角');figure(2)
6、;%画出频谱的幅度谱、频谱实部、频谱虚部subplot(3,1,1);stem(k,abs(xw),'filled');xlabel('k');ylabel('abs(xw)');title('幅度谱');subplot(3,1,2);stem(k,real(xw),'filled');xlabel('k');ylabel('real(xw)');title('频谱实部');subplot(3,1,3);stem(k,imag(xw),'filled');xlabel('k');ylabel('imag(x
7、w)');title('频谱虚部');1-2.a=0.9,length=10clc;clf;clear;%清除n=0:9;%设置自变量区间xn=(0.9).^n;%计算相应的x(n)k=0:9;%设置DFT采样长度xw=dftmtx(10)*xn';%DFT变换figure(1);%画出原序列的实部、虚部、模、相角subplot(2,2,1);stem(n,real(xn),'filled');xlabel('n');ylabel('real(xn)');title('xn实部');subplot(2,2,
8、2);stem(n,imag(xn),'filled');xlabel('n');ylabel('imag(xn)');title('xn虚部');文档大全实用标准文案subplot(2,2,3);stem(n,abs(xn),'filled');xlabel('n');ylabel('abs(xn)');title('xn模');subplot(2,2,4);stem(n,angle(xn),'filled')