现代数字信号信号处理matlab编程范例.docVIP

PAGE PAGE 39 信号处理仿真实验 目录 TOC \o 1-3 \h \z \u 一、 离散傅立叶变换 2 1、补零DFT 2 2、高密度频谱和高分辨率频谱 4 3、用DFT进行谱分析 6 二、 数字滤波器的设计 11 4、IIR滤波器的设计 11 5、FIR滤波器的设计 14 三、 短时傅立叶变换分析信号特性 17 6、短时傅立叶变换分析信号特性 17 四、 功率谱估计与应用 18 7、采用周期图法估计信号的功率谱 18 8、采用平均周期图法估计信号的功率谱 20 9、采用分段平均周期图法估计信号的功率谱 21 五、 同态信号处理 23 10、同态信号处理 24 六、 非平稳信号分析与处理 26 七、 Kalman滤波在机动目标跟踪中的应用 31 八、 机载R-D成像 55 离散傅立叶变换 1、补零DFT ，求N0分别取10，20时的X(k)。 解：设N0=10，则 设N0=20，则 用Matlab进行仿真，得图1-1。 Matlab源程序： %X(k)与X(ejw)的关系 clear; n=0:4;x=[ones(1,5)]; k=0:999;w=(pi/500)*k; X=x*(exp(-j*pi/500)).^(n*k);%计算离散时间傅立叶变换X(ejw) Xe=abs(X); subplot(321);stem(n,x);ylabel(x(n)); subplot(322);plot(w/pi,Xe);ylabel(|X(ejw)|); N=15;x=[ones(1,5),zeros(1,N-5)]; n=0:N-1; k=0:9; X=x*(exp(-j*pi/5)).^(n*k);%N=10点离散傅立叶变换 magX=abs(X); subplot(323);stem(n,x);ylabel(x(n)); subplot(324);stem(k,magX);axis([0,10,0,5]);ylabel(|X(k)|); N=20;x=[ones(1,5),zeros(1,N-5)]; n=0:1:N-1; k=0:N-1; X=x*(exp(-j*pi/10)).^(n*k);;%N=20点离散傅立叶变换 magX=abs(X); subplot(325);stem(n,x);ylabel(x(n)); subplot(326);stem(k,magX);axis([0,20,0,5]);ylabel(|X(k)|); 图1-1 与的关系 结论： 填零是给原始序列的运算，给给原始序列的离散时间傅立叶变换提供间隔较密的样本。 为画出，只需要5点的用内插公式即可得到，但实际上是用10或20点的来填充的值。 填零运算提供了较密的频谱，而没有增加任何新的信息，因此它不能提供高分辨率的频谱。 为得到高分辨率的频谱，需从实验或观察中取得更多的数据。 2、高密度频谱和高分辨率频谱 为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 当时，确定并画出x(n)的离散傅立叶变化。 当时，确定并画出x(n)的离散傅立叶变换。 当时，确定并画出x(n)的离散傅立叶变换。 用Matlab进行仿真，得图1-2。 Matlab源程序： %高密度频谱与高分辨率频谱 clear; figure(1); N=40;M=10; n=0:N-1; x=2*cos(0.35*pi*n)+cos(0.5*pi*n); n1=0:M-1;x1=x(1:M); subplot(231);stem(n1,x1);title(没有足够采样点的信号); k1=0:M-1;w1=2*pi*k1/M; X1=x1*(exp(-j*2*pi/M)).^(n1*k1);%计算M点DFT Y1=abs(X1); subplot(234);stem(w1/pi,Y1);title(信号的频谱); n2=0:N-1;x2=[x1,zeros(1,N-M)]; subplot(232);stem(n2,x2);title(填零信号); k2=0:N-1;w2=2*pi*k2/N; X2=x2*(exp(-j*2*pi/N)).^(n2*k2);%计算N点DFT Y2=abs(X2); subplot(235);stem(w2/pi,Y2);title(高密度频谱); subplot(233);stem(n,x);title(有足够采样点的信号); k=0:N-1;w=2*pi*k/N;