《控制系统仿真与CAD》
实验课程报告
                         
一、实验教学目标与基本要求
上机实验是本课程重要的实践教学环节。实验的目的不仅仅是验证理论知识，更重要的是通过上机加强学生的实验手段与实践技能，掌握应用MATLAB/Simulink 求解控制问题的方法，培养学生分析问题、解决问题、应用知识的能力和创新精神，全面提高学生的综合素质。
通过对MATLAB/Simulink进行求解，基本掌握常见控制问题的求解方法与命令调用，更深入地认识和了解MATLAB语言的强大的计算功能与其在控制领域的应用优势。
上机实验最终以书面报告的形式提交，作为期末成绩的考核内容。
二、题目及解答
第一部分：MATLAB 必备基础知识、控制系统模型与转换、线性控制系统的计算机辅助分析
1.
>>f=inline('[-x(2)-x(3);x(1)+a*x(2);b+(x(1)-c)*x(3)]','t','x','flag','a','b','c');[t,x]=ode45(f,[0,100],[0;0;0],[],0.2,0.2,5.7);plot3(x(:,1),x(:,2),x(:,3)),grid,figure,plot(x(:,1),x(:,2)),grid
2.
>>y=@(x)x(1)^2-2*x(1)+x(2);ff=optimset;ff.LargeScale='off';ff.TolFun=1e-30;ff.TolX=1e-15;ff.TolCon=1e-20;x0=[1;1;1];xm=[0;0;0];xM=[];A=[];B=[];Aeq=[];Beq=[];[x,f,c,d]=fmincon(y,x0,A,B,Aeq,Beq,xm,xM,@wzhfc1,ff)
Warning: Options LargeScale = 'off' and Algorithm =
'trust-region-reflective' conflict.
Ignoring Algorithm and running active-set algorithm. To run
网络对生活的影响
trust-region-reflective, set
LargeScale = 'on'. To run active-set without this warning, use
Algorithm = 'active-set'.
> In fmincon at 456
Local minimum possible. Constraints satisfied.
fmincon stopped because the size of the current search direction is less than
twice the selected value of the step size tolerance and constraints are
satisfied to within the selected value of the constraint tolerance.
<stopping criteria details>
Active inequalities (to within options.TolCon = 1e-20):
  lower      upper    ineqlin  ineqnonlin
    2                               
x =
    1.0000
        0
    1.0000
f =
  -1.0000
c =
    4
192.168.0.1修改WIFI密码d =
iterations: 5
funcCount: 20
lssteplength: 1
stepsize: 3.9638e-26
金时厚
怎么取消快捷支付
algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
firstorderopt: 7.4506e-09
constrviolation: 0
message: [1x766 char]
3.
(a) >> s=tf('s');G=(s^3+4*s+2)/(s^3*(s^2+2)*((s^2+1)^3+2*s+5))
G =
                      s^3 + 4 s + 2
  ------------------------------------------------------
  s^11 + 5 s^9 + 9 s^7 + 2 s^6 + 12 s^5 + 4 s^4 + 12 s^3
Continuous-time transfer function.
(b)
>> z=tf('z',0.1);H=(z^2+0.568)/((z-1)*(z^2-0.2*z+0.99))
H =
          z^2 + 0.568
  -----------------------------
  z^3 - 1.2 z^2 + 1.19 z - 0.99
Sample time: 0.1 seconds
Discrete-time transfer function.
4.
>> A=[0 1 0;0 0 1;-15 -4 -13];B=[0 0 2]';C=[1 0 0];D=0;G=ss(A,B,C,D),Gs=tf(G),Gz=zpk(G)
G =
a =
        x1  x2  x3
  x1    0    1    0
梦见乌龟是什么意思  x2    0    0    1
  x3  -15  -4  -13
洛克王国海豹旗手 b =
      u1
  x1  0
  x2  0
  x3  2
c =
      x1  x2  x3
  y1  1  0  0
d =
      u1
  y1  0