数字信号处理第三次实验报告.docVIP

DSP第三次实验实验报告 实验内容： 给定连续时间信号。 以足够小的时间间隔，在足够长的时间内画出信号时域图形。 用公式计算信号的频谱。以足够小的频率间隔，在足够大的频率范围内，画出其频谱图，估计信号的带宽。 以抽样频率3000Hz对抽样，得到离散时间信号，画出其图形，标明坐标轴。 用DFT计算的频谱，画出频谱图形，标明坐标轴。 由的频谱得到原信号的频谱的估计，计算均方误差在模拟频域上考察对原信号频谱的逼近程度。 将经理想内插后得到原信号的估计，从连续时间域上考察信号的恢复程度，计算均方误差。 抽样频率为00 Hz，重做3。 对比和分析，验证时域抽样定理。 t1=-0.1;%时间的起始位置 t2=0.1;%时间的终止位置 t=t1:0.00001:t2;%间隔为0.00001 x=exp(-1000*abs(t)); plot(t,x); xlabel(t/s); ylabel(x); title(x(t)=exp(-1000|t|)); 1.2运行结果 2.1编程原理和公式： 求的： 2.2程序脚本和注释： clear all; omega1=-3000*pi;%频率范围取-3000pi~3000pi omega2=3000*pi; omega=omega1:0.5*pi:omega2;%频率间隔取0.5pi N=length(omega); for(k=1:N) X(k)=2000/(1000000+omega(k)*omega(k)); end plot(omega,X); xlabel(omega（rad/s）); ylabel(X); title( X(omega)=2000/(1000000+omega(k)*omega(k))); 2.3程序运行结果： 信号的3DB带宽为W=208pi 即 104HZ 20DB带宽为W=956pi, 即 477HZ 3.1编程原理和公式 以3000HZ的抽样频率对x(t)进行抽样得x(n)： 3.2程序脚本和注释 clear all; f=3000; T=1/f; n1=-300;%n从-300变化到300 n2=300; n=n1:n2; xn=exp(-1000*abs(n*T)); stem(n,xn); xlabel(n); ylabel(x(n)); title(x(n)=exp(-1000|nT|)); 3.3运行结果 该题的图形与第一题的图形相似，第一题的图画的是t的连续信号，该图画的是时间离散的信号。 4.1编程原理和思路 用第二次实验中编写的dtft子函数计算x(n)的dtft,为了方便观察图形，利用dtft的周期性，画出-4pi-4pi四个周期内的图像。 4.2程序脚本和注释 clear all; f=3000; T=1/f; n1=-300; n2=300; n=n1:n2; M=128%每个周期内取128个点 xn=exp(-1000*abs(n*T));%得到离散时间信号x(n) Xww=dtft(xn,M);%计算x(n)的DTFT（在一个周期内0-2pi) Xww_abs=abs(Xww); ww=0:2*pi/M:2*pi*(M-1)/M; for(i=1:M)%利用周期性画出x(n)的DTFT在-4pi-4pi内的图像 Xw(i)=Xww_abs(i); Xw(M+i)=Xww_abs(i); Xw(M*2+i)=Xww_abs(i); Xw(M*3+i)=Xww_abs(i); end L=length(Xw); w=-4*pi:2*pi/M:2*pi*(L-1)/M-4*pi; plot(w,Xw); xlabel(w/rad); ylabel(Xw); title(x(n)的DTFT在-4pi-4pi内的图像)); 4.3运行结果 X(n)的DTFT在-4pi-4pi四个周期内的图像 5．1编程原理和公式 由可以得到离散时间信号x(n)的傅里叶变化 显然是以为周期的函数 通过理想低通滤波器即得到恢复出的信号的频谱。 利用公式（）计算均方误差。编程时，为第二小题中计算所得频谱，为由变换所得的频谱。 5.2程序脚本和注释 clear all; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %原信号频谱 omegaA1=-3000*pi;%频率范围取-3000pi~3000pi omegaA2=3000*pi; jg=0.5*pi; omegaA=omegaA1:jg:omegaA2;%频率间隔取0.5pi N=length(omegaA); for(k=1