《DSP实验1：有限长序列、频谱、DFT的性质》由会员分享，可在线阅读，更多相关《DSP实验1：有限长序列、频谱、DFT的性质（13页珍藏版）》请在金锄头文库上搜索。
1、课程名称: 数字信号解决 实验名称: 有限长序列、频谱、F的性质 一、实验目的和规定设计通过演示实验，建立对典型信号及其频谱的直观结识，理解DFT的物理意义、重要性质。二、实验内容和环节2-1用MTLAB,计算得到五种共9个序列:21-实指数序列例如,a=0.5, legth=1 =09， length1 a9, lengh2-12复指数序列例如，a=.5,=.8，th121-3从正弦信号（t)=s(2pfdela)抽样得到的正弦序列x(n)=s(2pfnT+la)。如,信号频率f=1Hz，初始相位ta=0,抽样间隔1秒，序列长lengh=10。2-14从余弦信号(t)=os(2pf et)抽
2、样得到的余弦序列x（n）=c（2pT+dela)。如,信号频率f1z,初相位delt，抽样间隔T=.1秒,序列长lengt10。2-1-5含两个频率分量的复合函数序列x(n)=n(2p1n)+delsin(pf2nphi)。如，频率f1(z)频率f2(H)相对振幅dla初相位ph (度)抽样间隔T(秒)序列长lngt130.50.100.51010510.110 用ML，对上述各个序列,反复下列过程。2-2-1画出一种序列的实部、虚部、模、相角；观测并记录实部、虚部、模、相角的特性。-2-2 计算该序列的幅度谱、频谱实部、频谱虚部;观测和并记录它们的特性，予以解释。2-23 观测同种序列取不同
3、参数时的频谱,发现它们的差别，予以解释。 三、重要仪器设备MA编程。四、操作措施和实验环节(参见“二、实验内容和环节”）五、实验数据记录和解决1. 实指数序列-. a=., ength1clc;lf;lar;%清除n=0:9；%设立自变量区间x05.;%计算相应的x(n)k=0:9;%设立FT采样长度xw=mt（0）*xn;%DF变换figure(1）;%画出原序列的实部、虚部、模、相角subplt(2,2,1)；stem(n,rel（xn)，filed)；label(x);ylabel（a(xn）);tile(xn实部）;supot(,2,2）;stem(n，imag(n),filled);
4、xlbel（）；yael(mg(xn);tl（x虚部);ubplot(2,2,3);stem（n,as（xn），fill);xlabel(）;ylabe(a（xn)）；ie(模）;ubot（，,4）;te(n，ngl（），lle）;xabel();ylbe(angle(x）；tile(相角);fgre（)；画出频谱的幅度谱、频谱实部、频谱虚部subplo(3,1，)；sem（k,a（xw),fild);xlbel(k)；lel(abs(xw）)；ti(幅度谱)；sublot(,1,);s（k，rea(x),fld)；xlabl(k）;lal(rea(xw);tite（频谱实部）;subpot(
5、3,1,3）;e(,imag(),filled);labl(k);lab(imag（w)）;itl（频谱虚部）;-2 =0.9， lenth=10cc;clf；clar;%清除n=0:9;设立自变量区间xn=()n;%计算相应的x()=0：9;%设立DFT采样长度xwdtmx(10）*xn;%D变换fgre();画出原序列的实部、虚部、模、相角suplot(2,2，1);t(,real(x),filed)；();yael(rel(xn);title(x实部);ublo（,2,2)；st（,mg(),illed);lb(n）;ylbel（i(xn);titl（xn虚部);sublt(2,3)；s
6、em（n,as(n）,filed）;lbel（n）；yabel(bs(xn)；titl（n模);sl(2，2,）;sem(n，ngl(xn)，flled);xlabe();ybel（anle(）;title(xn相角)；igue（2);画出频谱的幅度谱、频谱实部、频谱虚部suplot(3,1，1);tem(k,bs(x）,flld);xlabl()；ylal（bs(xw)；titl（幅度谱);suplt(，1,2);em(,real(xw)，filed)；xae(k)；ybel（rea（xw）);il(频谱实部);upot（3，1,)；em(k,iag(xw）,ld);xlabe（);ylab
7、(imag（w)；tite(频谱虚部); a=.9， ength=20clc；cf;cear;%清除n=0:19;设立自变量区间=(0.9).n；%计算相应的x(）k=0:;设立T的采样长度xwdtmx（20）*n;%DFT变换igue（);%画出原序列的实部、虚部、模、相角subplt(,1);stm(n,real(xn)，file)；xlabel（n);ylbe（eal(x)；til(xn实部);subplt(2,2,）;stem（n,mag(xn),filed);xabel()；ael（ag(xn）)；title(n虚部);subplot（2，2,3)；ste（n,as（xn）,fld)
8、；lael（n)；ylabel(bs()；titl(n模);subpot(2，2,4）;stem（n,nle(x）,filled);xlbl（n）;lbel(angle(x);ile（n相角）;igure（2)；%画出频谱的幅度谱、频谱实部、频谱虚部suplo(3,1)；t（k,ab(xw）,file)；xlabel（k);ylabl(ab(x）)；itl（幅度谱);ubplot（,1,）;stem（，al(w),filed);xlae(k)；ylbel(rea(xw）);e（频谱实部);sulot（3,3);stem（k,ima(),illed）;xle(k);labl(ia（w);tt（频
9、谱虚部);. 复指数序列0.5, b8， ngth=0clc;f；lar；%清除n=0：9;%设立自变量区间=(05+1.8)n; 计算相应的x()k=:;设立DF的采样长度xw=dftmt（10)*;%DF变换figure(1）;%画出原序列的实部、虚部、模、相角sublo(,2，1)；tem(n，ral（xn），filled）；lae(n）;lael（ra(xn);it(xn实部);ublot(2,，2)；stem(n,iag（n),fled);xbl(n)；ybl(ima();titl（n虚部);ubp(2，2，3);stem(,abs(xn)，file);lab();lbe(as（x）
10、);ttle(n模）;subpt（2,2,)；tm(，ng（xn),filed）;xlabl(n)；ylel（nle(n);il（x相角);fiu（2)；%画出频谱的幅度谱、频谱实部、频谱虚部subpo（，,1);tm(k,abs（xw)，lled)；xb(k）;ylbel(as（xw);tite（幅度谱)；ubplt(,);m(k,al(xw)，filld)；xlbel(k）;yael(rl(w）);e(频谱实部);sbplot(,1,3);stem(k,mg(xw),lled);lbl（k);yabl(imag（xw）)；ite(频谱虚部);. 从正弦信号x(t)si(2pft+eta)抽
11、样得到的正弦序列(n)=sn(2pfdelt)。如，信号频率=1H,初始相位dlta=0,抽样间隔T=0.1秒,序列长et=0。cl;cf;er;%清除0:;设立自变量区间xn=si(2*n*0.）；计算相应的x(n）k=0:9;设立DFT的采样长度xwttx(10)*n;%FT变换figur(1)；%画出原序列的实部、虚部、模、相角ubot(2,2,）;stem(n,rel(xn),illed);xlbel(n);ylbel(real(x）;te(n实部）;blo（2,2,2)；(,mg(n)，filled）;xlael（);ylbl（ag(x)）；(x虚部）;suplot(2，3）;sem(n,abs（n)，filed）；xlabel(n)；lbel（ab(x）);tile(xn模）;st(2,2，4)；ste(n,ngle（n)，fil);xlael(n);ylabe（angle(n）)；ite(n相角);figre（）； 画出频谱的幅度谱、频谱实部、频谱虚部u（3,1,);em(k,s（xw),flld)；xlbe(k);ylabel(as（）)；tile(幅度谱);ubplot(3,2);sem(k,real(x)，filled)；xle(k);ylab（real（)）;tile（频谱实部)；bplot（，,3）;stm(k,img(x),e
