数字信号处理上机报告-一.pdfVIP

数字信号处理第一次上机实验报告 实验一 ：  1000t t 设给定模拟信号x t e ， 的单位是ms。 a (1)利用MATLAB绘制出其时域波形和频谱图〔傅里叶变换〕，估计其等效带宽 〔忽略谱分量降低到峰值的3%以下的频谱〕。 (2)用两个不同的采样频率对给定的 进行采样。 x t a 1    ○以f 5000样本 秒采样x t 得到x n 。 s a 1     画出x n 及其频谱X ej 。 1 1 2    ○以f 1000样本 秒采样x t 得到x n 。 s a 1     画出x n 及其频谱X ej 。 1 1 比较两种采样率下的信号频谱，并解释。 实验一MATLAB程序： 〔1〕 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); 运行结果： 等效带宽为 〔2〕 . 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(221);stem(n,xn);grid on;axis([-10, 10, 0, 1.25]); title('时域波形(f_s=5000)'); xlabel('n');ylabel('x_1(n)'); subplot(222);plot(omega/pi,abs(X)); title('频谱图(f_s=5000)'); xlabel('\omega/\pi');ylabel('X_1(f)'); %Fs=1000 Fs=1; Ts=1/Fs; n=-N:Ts:N; xn=exp(-abs(n)); X=xn*exp(-j*(n'*omega)); subplot(223);stem(n,xn);grid on;axis([-10, 10, 0, 1.25]); title('时域波形(f_s=1000)'); xlabel('n');ylabel('x_2(n)'); subplot(224); plot(omega/pi,abs(X)); title('频谱图(f_s=1000)'); xlabel('\omega/\pi');ylabel('X_2(f)'); 运行结果： 实验二：  at   f  1  T 给定一指数型衰减信号x t