计算机控制系统上机报告..doc

计算机控制系统大作业 完成时间：2012年6月13日星期三 上机作业1 3-9： 实验题目 如下图所示计算机控制系统，其中，，。 实验要求： 将连续部分变换为等效的脉冲传递函数G(z)，然后分别列写D(z)和G(z)对应的状态方程和输出方程，最终列出闭环系统的状态方程和输出方程。 直接写出连续部分的连续状态方程，然后将该连续状态方程离散化，给出相应的离散状态方程和输出方程，再结合（1）中计算出的D(z) 对应的状态方程和输出方程，最终给出闭环系统的状态方程和输出方程。 分别对上述两种方法得到的闭环系统状态空间描述做Matlab仿真，并分别绘制闭环系统在单位阶跃信号输入的情况下的响应曲线，以及状态变化曲线。 对两种方法得到的仿真结果进行对比分析，给出仿真结论（即：两种数学描述的等价性）。 用Matlab提供的simulink仿真工具对该闭环系统进行仿真，给出单位阶跃响应曲线。（选作） 理论分析： （1）被控对象离散化： = 依串行法写状态方程： 控制器离散化: 状态方程为 闭环系统方程: 1 0] （2）由写出连续状态方程如下: = 最后求得的状态方程为： 实验程序源代码： %% N = 20; r = ones(N,1); %% selected state variables 1 state1 = zeros(3,N); x1 = zeros(3,1); % initial value of state variable y1 = zeros(N,1); F1 = [1 -1.25 1.25;0.4 -0.25 1.25;0 -0.3 -0.2]; G1 = [1.25 1.25 0.3]; C1 = [0 1 0]; for i = 1 : N y1(i) = C1 * x1; state1(:,i) = x1; x1 = F1 * x1 + G1 * r(i); end %% selected state variables 2 state2 = zeros(3,N); x2 = zeros(3,1); % initial value of state variable y2 = zeros(N,1); F2 = [0.75 0.05 0.005;-5 0 0.1;-15 -3 -0.2]; G2 = [0.005 0.1 0.3]; C2 = [50 10 0]; for i = 1 : N y2(i) = C2 * x2; state2(:,i) = x2; x2 = F2 * x2 + G2 * r(i); end %% figure(1); plot([1:N],y1,ko-); xlabel(iteration index k); ylabel(ouput y(k)); figure(2); plot([1:N],state1,-o); legend(x_1(k),x_2(k),x_3(k)); xlabel(iteration index k); ylabel(state variables x(k)); %% figure(3); plot([1:N],y2,k*-); xlabel(iteration index k); ylabel(ouput y(k)); figure(4); plot([1:N],state2,-*); legend(x_1(k),x_2(k),x_3(k)); xlabel(iteration index k); ylabel(state variables x(k)); 仿真结果： 用第一种方法绘制闭环系统在单位阶跃信号输入的情况下的响应曲线，以及状态变化曲线。 用第二种方法绘制闭环系统在单位阶跃信号输入的情况下的响应曲线，以及状态变化曲线。 结论： 由仿真得出对于同一个系统，由于状态变量的选择不同，可以得到不同的状态方程和输出方程，但最后的输出响应相同。所以状态变量的选取不影响系统的输出响应。 上机作业2 5.24 实验题目： 太阳光源跟踪系统利用伺服系统控制太阳电池帆板的移动，使其跟踪并始终垂直于太阳光线，最大程度地接收太阳能。太阳光源跟踪系统由感光器与检测线路和电机功率放大器（简化为增益放大环节），太阳帆板（电机负载），电机位置传感器和直流力矩电机组成。 太阳光源跟踪系统如下图（a）所示，计算机控制系统方框图如图（b）所示，其中。设计数字控制器，满足如下指标要求： （1