信电学院MATLAB第二次试验.doc

实验四 连续时间信号的频域分析 1、（1）编程实现信号的幅度频谱 M文件： syms t w; Gt1=sym(Heaviside(t+1)-Heaviside(t-1)); Gt2=sym(Heaviside(2*t+1)-Heaviside(2*t-1)); Fw1=fourier(Gt1,t,w) Fw2=fourier(Gt2,t,w) FP1=abs(Fw1); FP2=abs(Fw2); subplot(2,1,1);ezplot(FP1,[-10*pi,10*pi]);grid; axis([-10*pi,10*pi,0,2.2]);xlabel(w);ylabel(Fw1); %title(Gt1的幅频特性曲线); subplot(2,1,2);ezplot(FP2,[-20*pi,20*pi]);grid; axis([-20*pi,20*pi,0,1.2]);xlabel(w);ylabel(Fw2); %title(Gt2的幅频特性曲线); 运行结果： Fw1 = 2/w*sin(w) Fw2 = 2/w*sin(1/2*w) 2、(2)利用ifourier()函数，求频谱函数的傅氏反变换。 M文件： syms t w; Fw=sym((-w^2+5*i*w-8)/(-w^2+6*i*w+5)); ft=ifourier(Fw,w,t) 运行结果： ft = dirac(t)+(-3*exp(-t)+2*exp(-5*t))*heaviside(t) 实验五 系统的零极点及频域响应特性 1、（1）由系统函数H（s），求其零极点，并画出零极点图。 M文件： num=[1 0 1]; den=[1 2 5]; [z,p]=tf2zp(num,den) zplane(z,p) 运行结果： z = 0 + 1.0000i 0 - 1.0000i p = -1.0000 + 2.0000i -1.0000 - 2.0000i 2、(1) 已知系统函数H(s)，求其频率特性。 M文件： num=[2 0]; den=[1 2^0.5 1]; w=logspace(-2,2); freqs(num,den,w) 运行结果： 实验六 LTI系统的响应 1、（1）已知系统微分方程和激励信号e(t)，用MATLAB绘出系统单位冲激响应和系 统零状态响应的波形。 M文件： a=[1 4 4]; b=[1 3]; subplot(2,1,1);impulse(b,a,4) p=0.01; t=0:p:4; x=exp(-t); subplot(2,1,2);lsim(b,a,x,t) 运行结果： 2、（2）用MATLAB求出差分方程所描述的离散系统，在0~20时间范围内的单位函数 相应、阶跃响应和系统零状态响应的数值解，并绘出其波形。 M文件： a=[1 -0.7 0.1]; b=[7 -2 0]; figure(1); impz(b,a,-1:20) hk=impz(b,a,0:20) k=0:20; x=ones(size(k)); y=filter(b,a,x) figure(2); subplot(2,1,1);dstep(b,a,20) subplot(2,1,2);stem(k,y) 运行结果： hk = 7.0000 2.9000 1.3300 0.6410 0.3157 0.1569 0.0783 0.0391 0.0195 0.0098 0.0049 0.0024 0.0012 0.0006 0.0003 0.0002 0.0001 0.0000 0.0000 0.0000 0.0000 y = Columns 1 through 10 7.0000 9.9000 11.2300 11.8710 12.1867 12.3436 12.4218 12.4609 12.4805 12.4902 Columns 11 through 20 12.4951 12.4976 12.4988 12.4994 12.4997 12.4998 12.4999 12.5000 12.5000 12.5000 Column 21 12.5000