DSP实验-电信技巧.doc

实验一 连续时间系统的时域和频域分析 实验内容： 已知描述系统的微分方程和激励信号如下 ， 要求：（1）从理论上求解系统的冲激响应和零状态响应，并根据求解结果用MATLAB绘制其时域波形； （2）分别用MATLAB的impulse()函数和lsim()函数绘制系统的冲激响应和零状态响应，验证（1）中的结果。 程序： a=[1 4 4]; b=[0 1 3]; t=0:0.1:5; subplot(2,2,1) impulse(b,a); subplot(2,2,2); u=exp(-t); lsim(b,a,u,t); subplot(2,2,3) y1=exp(-2*t)+t.*exp(-2*t); plot(t,y1); xlabel(t); ylabel(y1); subplot(2,2,4); y2=(1+t).*(exp(-2*t)-exp(-3*t)); plot(t,y2); xlabel(t); ylabel(y2); 下图是用RLC元件构成的二阶低通滤波器。设，，，试用MATLAB的freqs()函数绘出该系统的幅频特性曲线和相频特性曲线。（求模：abs()函数，求相角：angle()函数） 程序： a=[0.08 0.4 1]; b=[0 0 1]; w=0:0.1:50; h=freqs(b,a,w); H=abs(h); subplot(2,1,1); plot(H); xlabel(w); ylabel(|H|); P=angle(h); subplot(2,1,2); plot(P); xlabel(w); ylabel(P); 已知某二阶系统的零极点分别为，，（二重零点），试用MATLAB绘出该系统在0 ~ 1kHz频率范围内的幅频特性曲线（要求用频率响应的几何求解方法实现），并说明该系统的作用。 程序： w=0:2*pi*1000; a1=100+i*w; a2=200+i*w; b1=0+i*w; b2=0+i*w; A1=abs(a1); A2=abs(a2); B1=abs(b1); B2=abs(b2); H=(B1.*B2)./(A1.*A2); plot(w,H); xlabel(w); ylabel(|H|); 实验二 离散时间系统的时域和频域分析 实验内容： 1. 已知某系统的系统函数为，，要求：（1）从理论上求解系统的单位冲激响应和零状态响应，并根据求解结果用MATLAB绘制其时域波形；（2）试分别用MATLAB的impz()函数和filter()函数绘制系统的单位冲激响应和零状态响应。 程序： N=20; n=0:N-1; h=-7/4*(-0.6).^n+11/4*(0.2).^n; yzs=-11/16*(0.2).^n-21/32*(-0.6).^n+75/32; subplot(221) stem(n,h,.);title(理论计算的单位冲激响应); subplot(222) stem(n,yzs,.);title(理论计算的零状态响应); b=[1 2]; a=[1 0.4 -0.12]; x=ones(1,N); yzs1=filter(b,a,x); subplot(224) stem(n,yzs1,.);title(系统的零状态响应); subplot(223) y=impz(b,a,20); stem(n,y,.);title(系统的单位冲激响应); 已知某系统的系统函数为 （1）绘制其零极点图 （2）用freqz()函数绘出该系统的幅频特性曲线和相频特性曲线，并说明该系统的作用。 程序： b=[1]; a=[1 -0.5]; subplot(221) zplane(b,a);title(系统的零极点分布); n=(0:50)*pi/50; [h,w]=freqz(b,a,n); subplot(223) plot(w/pi,abs(h));title(幅频特性曲线); subplot(224) plot(w/pi,angle(h));title(相频特性曲线); 3. 已知某系统的差分方程为，其中，，，设输入信号为，绘出该系统的幅频特性曲线和相频特性曲线，以及系统的输入信号和响应，并说明该系统的性能。 程序： b=[1 -1.1 0.6]; a=[0.6 -1.1 1]; n=(0:50)*pi/50; [h,w]=freqz(b,a,n); subplot(221) plot(w/pi,a