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

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

数字信号处理第一次上机实验报告

实验一：

设给定模拟信号()1000t a x t e -=，的单位是ms 。

(1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的进行采样。 ○1

。 ○2

。

比较两种采样率下的信号频谱，并解释。 实验一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)');

t ()a x t ()()15000s a f x t x n =以样本秒采样得到。()()11j x n X e ω画出及其频谱()()11000s a f x t x n =以样本秒采样得到。()()

11j x n X e ω画出及其频谱

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）.

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)');

运行结果：

实验二：

给定一指数型衰减信号，采样率，为采样周期。为方便起见，重写成复指数形式。

采样后的信号为，加窗后长度为的形式为：

这3个信号，，的幅度谱平方分别为： 模拟信号：

采样信号：

()()0cos 2at x t e f t π-=1

s f T

=T ()02j f t at x t e e π-=()02j f nT anT x nT e e π-=L ()(),0,1,,1L x nT x nT n L ==-L ()x t ()x nT ()L x nT ()()()

2

2

2

01

2X f a f f π=

+-()()()2

201

ˆ12cos 2aT

aT

X

f e

f f T e

π--=--+

加窗（取有限个采样点）信号：

且满足如下关系： 实验内容

（1） 在同一张图上画出：模型号幅度谱平方；

（2） 在同一张图上画出：模型号幅度谱平方；

；改变值，

结果又如何？

（1） f=0:0.01:3; alpha=0.2; f0=0.5; L=10; T1=1; T2=0.5;

Xa=1./(alpha^2+(2*pi*(f-f0)).^2); Xs1=T1*(1-2*exp(-alpha*T1*L)*cos(2*pi*(f-f0)*T1*L)+exp(-2*alpha*T1*L))./(1-2*exp(-alpha*T1)*cos(2*pi*(f-f0)*T1)+exp(-2*alph a*T1));

()()()()()22

02012cos 2ˆ12cos 2aTL aTL L

aT

aT

e f f TL e X f e

f f T e

ππ------+=--+()()()()ˆˆˆlim ,lim s L

L f X f X f TX f X f →∞

→∞

==100.2sec ,0.5Hz,1Hz 2Hz =10s s a f f f L -====取采样频率分别取和，。()2

X f ()()2

ˆ1Hz 2Hz 0Hz 3Hz s s f f TX

f f ==≤≤和时，采样信号幅度谱平方()2

X f ()()2

ˆ2Hz 0Hz 3Hz s f TX

f f =≤≤时，采样信号幅度谱平方L