实验范例(2009春).pptVIP

ts=0;te=5;dt=0.01; 实验一示例1 sys=tf([2,4],[1,4,3]);%sys是LTI系统模型 t=ts:dt:te; y=impulse(sys,t);%用impulse()函数来仿真分析系统冲击响应，如果是仿真分析阶跃响应，则用step()函数 plot(t,y); xlabel(t(sec)); ylabel(y(t)); grid on; ts=0;te=5;dt=0.01; 实验一示例2 sys=tf([6],[1,5,6]); t=ts:dt:te; f=10*sin(2*pi*t); y=lsim(sys,f,t);%用lsim()函数来仿真分析系统零状态响应 plot(t,y); xlabel(t(sec));ylabel(y(t)); title(零状态响应); grid on; f=sym(‘exp(-2*t)*heaviside(t)’); 实验二示例 F=fourier(f); FT=subs(F,[-2*pi:0.01*pi:2*pi]);%用subs()函数来得到傅立叶变换的数值解从而绘出幅频和相频特性曲线 subplot(2,1,1); plot([-2*pi:0.01*pi:2*pi],abs(FT)); title(f(t)的幅频特性曲线); xlabel(w);ylabel(幅度); subplot(2,1,2); plot([-2*pi:0.01*pi:2*pi],angle(FT)); title(f(t)的相频曲线);xlabel(w);ylabel(相位(弧度)); b=[1]; 实验四示例1 a=[1,-1]; k=0:10; subplot(1,2,1); zplane(b,a);%用zplane()函数绘制系统零极点图 grid on; h=impz(b,a,k);%用impz()函数来绘制单位函数响应。 subplot(1,2,2); stem(h);%离散序列绘图用stem()函数 grid on; b=[1]; 实验四示例2 a=[1,-0.5]; subplot(2,2,1); zplane(b,a); h=freqz(b,a);%用freqz()函数来分析系统的频率响应特性 subplot(2,2,2); plot(h); subplot(2,2,3); plot(abs(h)); subplot(2,2,4); plot(angle(h)); * *