资源描述：
《数字信号处理上机报告一》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、数字信号处理第一次上机实验报告实验一：设给定模拟信号乙⑴之一网“，/的单位是ms。(1)利用MATLAB绘制岀其时域波形和频谱图(傅里叶变换)，估计英等效带宽(忽略谱分量降低到峰值的3%以下的频谱)。⑵用两个不同的采样频率对给定的乙⑴进行采样。①以£=5000样木/秒采样％⑴得到x.(h)o画出西0)及其频谱②以/；=1000样本/秒采样£(/)得到西(/?)0画出西(町及其频谱x】(严)。比较两种采样率下的信号频谱，并解释。实验一MATLAB程序：(1)N=10;Fs=5;Ts=l/Fs;n=[-N:Ts:N];xn=exp(-abs(n));w=-4*pi:0・01:4*p
2、i;X=xn*exp(-j*(n'*w));subplot(211)plot(n,xn);title(*x_a(t)时域波形*);xlabel(ft/ms*);ylabel(*x_a(t)');axis([-10,10,0,1]);subplot(212);plot(w/pizabs(X));title(1x_a(t)频谱图1);xlabel(1omega/pi');ylabel('X_a(eA(jomega))*);ind=find(X〉=0・03*max(X))*0.01;eband=(max(ind)-min(ind));fprintf(*等效带宽为%fKHZ*
3、zeband);运行结果：aFigure1I=
4、回l・^J等效带宽为12.110000KHZ(2).N=10;omega二-3*pi:0.01:3*pi;%Fs=5000Fs二5;Ts二1/Fs;n=-N:Ts:N;xn=exp(-abs(n));X二xn*exp(-j*(n‘Iomega));subplot(221);stem(n,xn);gridon;axis([T0,10,0,1.25]);titleC吋域波形(f_s二5000)');xlabelCn');ylabel('x_l(n)');subplot(222);plot(omega/pi,abs(X));titleC
5、频谱图(f_s二5000)');xlabelComega/pi');ylabel('X_1(f)');%Fs=1000Fs=l;Ts=l/Es;n二-N:Ts:N;xn=exp(~abs(n));X二xn*exp(-j*(n，Iomega));subplot(223);stem(n,xn);gridon;axis([-10,10,0,1.25]);title(，时域波形(f_s二1000)');xlabel('n);ylabel('x2(n)');subplot(224);plot(omega/pi,abs(X));title('频谱图(f_s=1000)');xlabel
6、('omega/pi');ylabel('X_2(f)');运行结杲：❻Figure1["回实验二给定一指数型衰减信号曲)=严cos(2龙加)，采样率丁为采样周期。为方便起见，重写成复指数形式x(t)=e-atej2^o采样后的信号为x(nT)=严％対砂，加窗后长度为厶的形式为:xL(hT)=x(hT),m=0丄…，厶一1这3个信号兀（小兀（M）,xL（nT）的幅度谱平方分别为:模拟信号:]/+（2兀（/_九）『采样信号:xlabelCf/Hz');ylabel('
7、X(f)「2');legendC模拟信号幅度谱平方
8、X(f)「2','f_s=lHz时,釆样信号幅度谱平方
9、
10、TX(f)「2','f_s=2Hz时,釆样信号幅度谱平方ITX(f)「2');运行结果：(2)f=0:0.01:3;alpha=0.2;f0=0.5;Ll=5;L2=10;L3=20;T1=0.5Xa=l./(alpha*2+(2*pi*(f-fO)).2);Xsl=Tl*(l-2*exp(-alpha*Tl*Ll)*cos(2*pi*(f-fO)*T1*L1)+exp(-2*alpha*Tl*Ll))./(1-2*exp(-alpha*Tl)*cos(2*pi*(f-fO)*T1)+exp(-2*alpha*Tl));Xs2=Tl*(l-2*exp(-alpha*Tl*L2)
11、*cos(2*pi*(f-f0)*Tl*L2)+exp(-2*alpha*Tl*L2))./(l-2*exp(~alpha*Tl)*cos(2*pi*(f-f0)*Tl)+exp(~2*alpha*Tl));Xs3=Tl*(12*exp(-alpha*Tl*L3)*cos(2*pi*(f-fO)*T1*L3)+exp(-2*alpha*Tl*L3))./(1-2*exp(-alpha*Tl)*cos(2*pi*(f-fO)*T1)+exp(-2*alpha*Tl));plot(f,Xa,'b
