《数字信号处理实验作业完全版.docxVIP

实验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(w0*n*T);close allsubplot(3,2,1);stem(x,’.’);title(理想采样信号序列);k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n*k);magX=abs(X);subplot(3,2,2);stem(magX,’.’);title(理想采样信号序列的幅度谱);angX=angle(X);subplot(3,2,3);stem(angX);title(理想采样信号序列的相位谱)n=0:50;A=1;a=0.4,w0=2.0734;T=1;x=A*exp(-a*n*T).*sin(w0*n*T);subplot(3,2,4);stem(x,’.’);title(理想采样信号序列);k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n*k);magX=abs(X);subplot(3,2,5);stem(magX,’.’)title(理想采样信号序列的幅度谱);上机实验答案：分析理想采样信号序列的特性产生在不同采样频率时的理想采样信号序列Xa(n)，并记录各自的幅频特性，观察频谱“混淆”现象是否明显存在，说明原因。源程序：A=444.128;a=50*pi*sqrt(2.0);W0=50*pi*sqrt(2.0);n=-50:1:50;T1=1/1000;Xa=A*(exp(a*n*T1)).*(sin(W0*n*T1));subplot(3,3,1);plot(n,Xa);title(Xa序列);xlabel(n);ylabel(Xa);k=-25:25;X1=Xa*(exp(-j*pi/12.5)).^(n*k);subplot(3,3,2); stem(k,abs(X1),.);title(Xa的幅度谱);xlabel(k);ylabel(·幅度);subplot(3,3,3);stem(k,angle(X1),.); title(Xa的相位谱);xlabel(k);ylabel(相位);T2=1/300;Xb=A*(exp(a*n*T2)).*(sin(W0*n*T2)); subplot(3,3,4);plot(n,Xb);title(Xb序列);xlabel(n);ylabel(相位);k=-25:25;X2=Xb*(exp(-j*pi/12.5)).^(n*k);subplot(3,3,5); stem(k,abs(X2),.);title(Xb的幅度谱);xlabel(k);ylabel(·幅度);subplot(3,3,6);stem(k,angle(X2),.); title( Xb的相位谱);xlabel(k);ylabel(相位);T3=1/200;Xc=A*(exp(a*n*T3)).*(sin(W0*n*T3)); subplot(3,3,7);plot(n,Xc);title(Xc序列);xlabel(n);ylabel(Xc);k=-25:25;X3=Xc*(exp(-j*pi/12.5)).^(n*k);subplot(3,3,8); stem(k,abs(X3),.);title(Xc的幅度谱);xlabel(k);ylabel(幅度);subplot(3,3,9);stem(k,angle(X3),.); title(Xc的相位谱);xlabel(k);ylabel(相位);由图可以看出：当采样频率为1000Hz时，采样序列在折叠频率附近处，无明显混叠。当采样频率进一步降低时，主瓣宽度逐渐变宽，频率混叠现象也逐渐严重，存在较明显的失真现象。原因：信号理想采样后的频谱是原信号频谱的周期严拓，其严拓周期等于采样角频率。根据时域采样定理，如果原信号是带限信号，且采样频率高于原信号最高频率的2倍才不会发生频率混叠。否则频谱会发生明显混叠如采样频率为200和300Hz时。离散信号、系统和系统响应的分析观察信号Xb(n)和Hb(n)的时域和频谱特性，比较系统输出Y(n)与Hb(n)的时域及幅频特性，绘出图形。clcn=1:50;hb=zeros(1,50);hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;subplot(3,3,1);stem(hb);title(系统hb[n]);m=1:50;T=0.001;A=444.128;a=sqrt(2.0)*pi;w0=50*sqrt(2.0)*pix=A*