实验离散时间信号的分析实验.docVIP

数字信号处理——实验一 武汉工程大学 电气信息学院 通信工程 红烧大白兔 实验目的 认识常用的各种信号，理解其数学表达式和波形的表示。 掌握在计算机中生成及绘制数字信号波形的方法。 掌握序列的简单运算及计算机实现与作用。 理解离散时间傅里叶变换、Z变换及它们的性质和信号的频域特性。 实验设备 计算机，MATLAB语言环境。 实验基础理论 序列的相关概念 离散时间信号用一个称为样本的数字序列来表示。一般用{x[n]}表示，其中自变量n的取值范围是负无穷到正无穷之间的整数。为了表示方便，序列通常直接用x[n]表示。 离散时间信号可以是一个有限长序列，也可以是一个无限长序列。 满足x[n+kN]=x[n]（对于所有n）的序列称为周期为N的周期序列，其中N取任意正整数；k取任意整数。 常见序列 单位取样信号、单位阶跃信号、矩形序列、斜变序列、单边指数序列、正弦序列、复指数序列等。 序列的基本运算 加法、乘法、倒置、移位、尺度变换、卷积等。 离散傅里叶变换的相关概念 5、Z变换的相关概念 实验内容与步骤 （1）信号的平滑 R=51; D=0.8*(rand(R,1)-0.5); m=0:R-1; s=2*m.*(0.9.^m); x=s+D; subplot(2,1,1); plot(m,D,r-,m,s,g--,m,x,b-.); xlabel(Time index n);ylabel(Amplitude); legend(D[n],s[n],x[n]); x1=[0 0 x];x2=[0 x 0];x3=[x 0 0]; y=(x1 + x2 + x3)/3; subplot(2,1,2); plot(m,y(2: R+1),r-,m,s,g--); legend(y[n],s[n]); xlabel(Time index n);ylabel(Amplitude); （2）调制信号的产生 n=0:0.01:10;m=10;A=8; wl=pi;wh=10*pi; xl=cos(wl*n); xh=cos(wh*n); subplot(3,1,1); plot(n,xl); subplot(3,1,2); plot(n,xh); y=A.*(1+m.*xl).*xh; subplot(3,1,3); plot(n,y); 卷积的计算 h=[3 2 1 -2 1 0 -4 0 3]; x=[1 -2 3 -4 3 2 1]; y=conv(h,x); n=0 : 14; subplot(2,1,1); stem(n,y); xlabel(Time index n);ylabel(Amlitude); title(Output Obtained by Convolution);grid; （4）信号离散傅里叶的正反变换 x=[0,1,2,3,4,5,6,7]; N=length(x); n=0:(N-1); k=0:(N-1); Xk=xn*exp(-j*2*pi/N).^(n*k); subplot(1,2,1); stem(n,abs(x)); title(IDFT|X(k)|); subplot(1,2,2); stem(k,abs(Xk)); title(|X(k)|); （5）圆周移位 M=6; a=[0 1 2 3 4 5 6 7 8 9]; b=circshift(a,[0 M]); L=length(a)-1; n=0:L; subplot(2,1,1); stem(n,a);axis([0,L,min(a),max(a)]); title(Original Sequence); subplot(2,1,2); stem(n,b);axis([0,L,min(a),max(a)]); title([Sequence Obtained by Cirularly Shifting by,num2str(M),Samples]); （6）时域移位 x=[0 2 4 6 8 10 12 14 16]; N=length(x)-1;n=0:N; y=circshift(x,5); XF=fft(x); YF=fft(y); subplot(2,2,1); stem(n,abs(XF));grid title(Magnitude of DFT of Original Sequence); subplot(2,2,2); stem(n,abs(YF));grid title(Magnitude of DFT of Cicularly Shifted Sequence); subplot(2,2,3); stem(n,angle(XF));grid