《控制系统仿真与CAD》.docVIP

实验一 Matlab使用方法和程序设计 二、 实验内容 6、基本程序设计 （1）编写命令文件：计算1+2+…+n2000时的最大n值； n=0;s=0; while s2000 n=n+1; s=s+n; end n=n-1 （2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。 （2.1）function [sum]=mi(n) sum=0; for m=0:n sum=sum+2^m; end （2.2）function [sum]=miwhile(n) sum=0;m=0; while m=n sum=sum+2^m; m=m+1; end （3）如果想对一个变量x自动赋值。当从键盘输入y或Y时（表示是），x自动赋为1； 当从键盘输入n或N时（表示否），x自动赋为0；输入其他字符时终止程序。 s=input(please input a key:,s) while(s==n||s=N||s==y||s=Y) if(s==y||s=Y) clc; x=1; disp(x=); disp(x); else(s==n||s=N) clc; x=0; disp(x=); disp(x); end s=input(please input a key,s) clc; end 实验二 控制系统分析（一） 二、 实验内容 1、时域分析 （1）根据下面传递函数模型：绘制其单位阶跃响应曲线并从图上读取最大超调量，绘制系统的单位脉冲响应、零输入响应曲线。 num=5*[1 ,5 ,6]; den=[1, 6 ,10, 8]; sys=tf(num,den); step(sys) figure(2) sys1=ss(sys); x0=[1;0;0]; initial(sys1,x0) （2）典型二阶系统传递函数为： 当ζ=0.7,ωn取2、4、6、8、10、12的单位阶跃响应。 zet=0.7;wn=[2,4,6,8,10,12]; hold on for i=1:length(wn) sys=tf(wn(i)^2,[1,2*zet*wn(i),wn(i)^2]); step(sys) end （3）典型二阶系统传递函数为： 当ωn =6,ζ取0.2、0.4、0.6、0.8、1.0、1.5、2.0的单位阶跃响应。 wn=6;zet=[0.2,0.4,0.6,0.8,1.0,1.5,2.0]; hold on for i=1:length(zet) sys=tf(wn^2,[1,2*zet(i)*wn,wn^2]); step(sys) end 2、频域分析 （1）典型二阶系统传递函数为： 当ζ=0.7,ωn取2、4、6、8、10、12的伯德图。 zet=0.7;wn=[2,4,6,8,10,12]; hold on for i=1:length(wn) sys=tf(wn(i)^2,[1,2*zet*wn(i),wn(i)^2]); bode(sys) end （2）典型二阶系统传递函数为： 当ωn =6,ζ取0.2、0.4、0.6、0.8、1.0、1.5、2.0的伯德图。 wn=6;zet=[0.2,0.4,0.6,0.8,1.0,1.5,2.0]; hold on for i=1:length(zet) sys=tf(wn^2,[1,2*zet(i)*wn,wn^2]); bode(sys) end 3、根轨迹分析 根据下面负反馈系统的开环传递函数，绘制系统根轨迹，并分析系统稳定的K值范围。 sys=tf(1,[1, 3, 2 ,0]) rlocus(sys) [x,y]=ginput(1) p=x+y*i k=rlocfind(sys,p) 实验三 控制系统分析（二） 二、实验内容 1、系统稳定性分析 （1）代数法稳定性判据：（用求分母多项式的根和routh函数两种方法来判定） 已知系统的开环传递函数为： 试对系统闭环判别其稳定性 k=100; p=[0,-1,-20]; z=-2; sys=zpk(z,p,k) sys=tf(sys) G_c=feedback(sys,1);%转换为闭环传递函数 roots(G_c.den{1})%求闭环特征方程的根 [rtab,msg]=routh(G_c.den{1}) （2）根轨迹法判断系统稳定性： 已知一个单位负反馈系统开环传递函数为：