资源描述:
《数字信号处理实验作业完全版》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实验1:理想采样信号的序列,幅度谱,相位谱,以及改变参数后的图像。源程序:clc;n=0:50;A=444.128;a=50*sqrt(2.0)*pi;T=0.001;w0=50*sqrt(2.0)*pi;x=A*exp(-a*n*T).*sin(wO*n*T);closeallsubplot(3,2,l);stem(x,J・');title(理想采样信号序列');k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).A(n'*k);magX=abs(X);subplot(3,2⑵;stem(magX,‘・'理想采样信号序列的幅
2、度谱');angX=angle(X);subplot(3,2,3);stem(angX);titleCH想采样信号序列的相位谱)n=0:50;A=l;a=0.4?w0=2.0734;T=l;x=A*exp(-a*n*T).*sin(wO*n*T);subplot(3,2?4);stem(x/・,);title。理想采样信号序列1;k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).A(n'*k);magX=abs(X);subplot(3,2,5);stem(magX,'・,)title。理想采样信号序列的幅度谱J;理想采样信
3、号序列的幅度谱理想采样信号序列1理想采样信号序列的相位谱■1o2040600204060上机实验答案:分析理想采样信号序列的特性产生在不同采样频率时的理想采样信号序列Xa(n),并记录各自的幅频特性,观察频谱“混淆”现象是否明显存在,说明原因。源程序:A=444.128;a=50*pi*sqrt(2.0);W0=50*pi*sqrt(2.0);n=-50:l:50;Tl=l/1000;Xa=A*(exp(a*n*Tl)).*(sin(WO*n*Tl));subplot(3,3,l);plot(n,Xa);title(Xa序列);xlabel('n');ylabel
4、('Xa');k=-25:25;X1=Xa*(exp(-j*pi/12.5)).A(n'*k);subplot(3,3,2);stem(k,abs(Xl),'.);titleCXa的幅度谱');xlabelCk^ylabelC・幅度');subplot(3,3,3);stem(k,angle(X1),'.');title('Xa的相位谱);xlabel(k);ylabelC相位);T2=1/300;Xb=A*(exp(a*n*T2)).*(sin(W0*n*T2));subplot(3,3,4);plot(n,Xb);title('Xb序列');xlabel('n
5、*);ylabel('位);k=・25:25;X2=Xb*(exp(-j*pi/12.5)).A(n,*k);subplot(3,3,5);stem(k,abs(X2)/.*);titleCXb的幅度谱);xlabelCk*);ylabel('・幅度);subplot(3,3,6);stem(k,angle(X2)/.');titleCXb的相位谱);xlabel(k);ylabel(相位);T3=l/200;Xc=A*(exp(a*n*T3)).*(sin(W0*n*T3));subplot(3,3,7);plot(n,Xc);title(Xc序列);xlabe
6、l('n,);ylabel('Xc');k=-25:25;X3=Xc*(exp(-j*pi/12.5)).A(nr*k);subplot(3,3,8);stem(k,abs(X3)/.*);titleCXc的幅度谱);xlabel('k');ylabel('度);subplot(3,3,9);stem(k,angle(X3)title('Xc的相位谱);xlabel(k);ylabel('相位');0-5-50050n5心列2x榔勺幅度谱0-50050Xa的相位谱50-5-50050nxb知序列翌1SkXb的相位谱50x10先序列-5-50050nk由图可以看出:
7、当采样频率为1000Hz时,采样序列在折叠频率附近处,无明显混叠。当采样频率进一步降低时,主瓣宽度逐渐变宽,频率混叠现象也逐渐严重,存在较明显的失真现象。原因:信号理想采样后的频谱是原信号频谱的周期严拓,其严拓周期等于采样角频率。根据时域采样定理,如果原信号是带限信号,且采样频率高于原信号最高频率的2倍才不会发生频率混叠。否则频谱会发生明显混叠如采样频率为200和300Hz时。1、离散信号、系统和系统响应的分析1)观察信号Xb(n)和Hb(n)的时域和频谱特性,比较系统输出Y(n)与Hb(n)的时域及幅频特性,绘出图形。clcn=l:50;hb=zeros(l,5
8、0);hb