哈尔滨工业大学信号与系统第三时间机器实验.docVIP

实验五 连续系统的频域分析及连续信号的采样与重构 二、利用MATLAB分析系统频响特性 实验内容 FH图5-4 二阶低通滤波器电路图 图5-4所示的电路为最平坦幅度型二阶低通滤波器。试用MATLAB程序画出系统响应的幅度响应及相频响应，并与理论分析的结果进行比较。的截止频率 F H 图5-4 二阶低通滤波器电路图 b=[0 0 1]; a=[2 2.82 2]; [h,w]=freqs(b,a,100); h1=abs(h); h2=angle(h); subplot(211); plot(w,h1); grid xlabel(角频率(W)); ylabel(幅度); title(H(jw)的幅频特性); subplot(212); plot(w,h2*180/pi); grid xlabel(角频率(w)); ylabel(相位(度)); title(H(jw)的相频特性); 三、连续信号的采样与重构 实验内容 设，由于不是严格的带限信号，但其带宽可根据一定的精度要求做一近似。试根据以下三种情况用MATLAB实现由采样信号重构并求出两者误差，分析三种情况下的结果。 （1），，； 运行程序： wm=5000*pi; wc=1*wm; Ts=pi/wm; ws=2*pi/Ts; n=-100:100; nTs=n*Ts; f=exp(-1000.*abs(nTs)); Dt=0.005;t=-15:Dt:15; fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t)))); error=abs(fa-exp(-1000.*abs(t))); t1=-15:0.5:15; f1=exp(-1000.*abs(t1)); subplot(311); stem(t1,f1); ylabel(f(kTs)); title(exp(-1000.*abs(t))的采样信号); subplot(312); plot(t,fa) ylabel(fa(t)); title(由exp(-1000.*abs(t))的采样信号重构)); grid; subplot(313); plot(t,error); ylabel(error(t)); title(采样信号与原信号的误差error(t)); （2），，； 运行程序： wm=10000*pi; wc=1.1*wm; Ts=pi/wm; ws=2*pi/Ts; n=-100:100; nTs=n*Ts; f=exp(-1000.*abs(nTs)); Dt=0.005;t=-15:Dt:15; fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t)))); error=abs(fa-exp(-1000.*abs(t))); t1=-15:0.5:15; f1=exp(-1000.*abs(t1)); subplot(311); stem(t1,f1); ylabel(f(kTs)); title(exp(-1000.*abs(t))的采样信号); subplot(312); plot(t,fa) ylabel(fa(t)); title(由exp(-1000.*abs(t))的过采样信号重构)); grid; subplot(313); plot(t,error); ylabel(error(t)); title(过采样信号与原信号的误差error(t)); （3），，； 运行程序： wm=2500*pi; wc=0.9*wm; Ts=pi/wm; ws=2*pi/Ts; n=-100:100; nTs=n*Ts; f=exp(-1000.*abs(nTs)); Dt=0.005;t=-15:Dt:15; fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t)))); error=abs(fa-exp(-1000.*abs(t))); t1=-15:0.5:15; f1=exp(-1000.*abs(t1)); subplot(311); stem(t1,f1); ylabel(f(kTs)); title(exp(-1000.*abs(t))的采样信号); subplot(312); plot(t,fa) ylabel(fa(t)); title(由exp(-1000.*abs(t))的采样信号重构)); grid; subplot(313); plot(t,error); ylabel(error(t)); title(采