西电数字信号处理上机实验报告.docx

. . . . 数字信号处理上机实验报告 14020710021 张吉凯 第一次上机 实验一： 设给定模拟信号  1000 t xa t e  ， t 的单位是 ms 。 利用 MATLAB 绘制出其时域波形和频谱图（傅里叶变换） ，估计其等效带宽 （忽略谱分量降低到峰值的 3% 以下的频谱）。 用两个不同的采样频率对给定的 xa t 进行采样。 1○以f s 5000样本 秒采样 xa t 得到x1 n 。 画出x n 及其频谱 X e j 。 1 1 2○以f s 1000样本 秒采样 xa t 得到x1 n 。 画出x n 及其频谱 X e j 。 1 1 比较两种采样率下的信号频谱，并解释。 （ 1 ） MATLAB 程序： N=10; Fs=5; Ts=1/Fs; n=[-N:Ts:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t) 时域波形 '); xlabel('t/ms');ylabel('x_a(t)'); axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t) 频谱图'); xlabel('\omega/\pi');ylabel('X_a(e^(j\omega))'); ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf(' 等效带宽为 %fKHZ\n',eband); 运行结果： 等效带宽为 12.110000KHZ （ 2 ） MATLAB 程序： N=10; omega=-3*pi:0.01:3*pi; %Fs=5000 Fs=5; Ts=1/Fs; n=-N:Ts:N; xn=exp(-abs(n)); X=xn*exp(-j*(n'*omega)); subplot(2,2,1);stem(n,xn);grid on;axis([-10, 10, 0, 1.25]); title(' 时域波形 (f_s=5000)'); xlabel('n');ylabel('x_1(n)'); subplot(2,2,2);plot(omega/pi,abs(X)); title(' 频 谱 图 (f_s=5000)'); xlabel('\omega/\pi');ylabel('X_1(f)'); grid on; %Fs=1000 Fs=1; Ts=1/Fs; n=-N:Ts:N; xn=exp(-abs(n)); X=xn*exp(-j*(n'*omega)); subplot(2,2,3);stem(n,xn);grid on;axis([-10, 10, 0, 1.25]); title(' 时域波形 (f_s=1000)'); xlabel('n');ylabel('x_2(n)'); grid on; subplot(2,2,4); plot(omega/pi,abs(X)); title(' 频谱图 (f_s=1000)'); xlabel('\omega/\pi'); ylabel('X_2(f)'); grid on; 运行结果： 实验二： 给定一指数型衰减信号  x t e  at cos 2  f 0t ，采样率 fs  1 ，T 为采样周期。 为 T 方便起见，重写成复指数形式 x t e at ej 2 f 0 t 。 采样后的信号为 x nT e anT ej 2  f0 nT ，加窗后长度为 L 的形式为： xL nT x nT , n 0,1,L , L 1 这 3 个信号 x t ， x nT ， xL nT 的幅度谱平方分别为： 模拟信号： 2 X f a 2 2 1 2 f f0 采样信号： 2X? f 2  1 2e aT  cos 2 1 f f0 T e  2aT 加窗（取有限个采样点）信号： 2X? f 2 1 2e aTL cos 2 f f0 TL e 2 aTL L 1 2e aT cos 2 f f0 T e 2aT 且满足如下关系： XLlim ? X L L f X? f , lim f TX? f X f s 实验内容 1取a 0.2sec 1 , f0 0.5Hz, 采样频率分别取 fs 1Hz和f s 2Hz ， L =10。 ） 在同一张图上画出：模型号幅度谱平方 2 X f ； fs 1Hz和f s 2Hz 时，采样信号幅度谱平方 TX?