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

实验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 all subplot(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(理想采样信号序列的幅度谱; 0204060-2000 200理想采样信号序列 0204060 0500 1000理想采样信号序列的幅度谱 0204060-50 5理想采样信号序列的相位谱0204060 -1 1理想采样信号序列 0204060 1 2理想采样信号序列的幅度谱 上机实验答案: 分析理想采样信号序列的特性 产生在不同采样频率时的理想采样信号序列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(相位; -50 050 -5 05 7 X a 序列 n X a -50 050012 8X a 的幅度谱k 幅度 -50 050-5 5 X a 的相位谱k 相位 -50 050 -50518 X b 序列 n 相位 -50 05005 1018X b 的幅度谱k 幅度 -50 050-5 5 X b 的相位谱 k 相位 -50 050 -505x 1026 X c 序列 n X c -50 050 05 10x 10 26X c 的幅度谱k 幅度 -50 050 -5 05 X c 的相位谱 k 相位 由图可以看出: 当采样频率为1000Hz时,采样序列在折叠频率附近处,无 明显混叠。 当采样频率进一步降低时,主瓣宽度逐渐变宽,频率混叠现 象也逐渐严重,存在较明显的失真现象。 原因:信号理想