资源描述：
《计算机控制系统上机报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、计算机控制系统大作业完成时间：2012年6月13日星期三上机作业13-9：实验题目如下图所示计算机控制系统，其中，，。实验要求：（1）将连续部分变换为等效的脉冲传递函数G(z)，然后分别列写D(z)和G(z)对应的状态方程和输出方程，最终列出闭环系统的状态方程和输出方程。（2）直接写出连续部分的连续状态方程，然后将该连续状态方程离散化，给出相应的离散状态方程和输出方程，再结合（1）中计算出的D(z)对应的状态方程和输出方程，最终给出闭环系统的状态方程和输出方程。（1）分别对上述两种方法得到的闭环系统状态空间描述做Matlab仿真，并分别绘制闭环系统在单位阶跃信号输
2、入的情况下的响应曲线，以及状态变化曲线。（2）对两种方法得到的仿真结果进行对比分析，给出仿真结论（即：两种数学描述的等价性）。（3）用Matlab提供的simulink仿真工具对该闭环系统进行仿真，给出单位阶跃响应曲线。（选作）理论分析：（1）被控对象离散化：=依串行法写状态方程：控制器离散化:状态方程为闭环系统方程:10]（2）由写出连续状态方程如下:=最后求得的状态方程为：实验程序源代码：%%N=20;r=ones(N,1);%%selectedstatevariables1state1=zeros(3,N);x1=zeros(3,1);%initialval
3、ueofstatevariabley1=zeros(N,1);F1=[1-1.251.25;0.4-0.251.25;0-0.3-0.2];G1=[1.251.250.3]';C1=[010];fori=1:Ny1(i)=C1*x1;state1(:,i)=x1;x1=F1*x1+G1*r(i);end%%selectedstatevariables2state2=zeros(3,N);x2=zeros(3,1);%initialvalueofstatevariabley2=zeros(N,1);F2=[0.750.050.005;-500.1;-15-3-0.2
4、];G2=[0.0050.10.3]';C2=[50100];fori=1:Ny2(i)=C2*x2;state2(:,i)=x2;x2=F2*x2+G2*r(i);end%%figure(1);plot([1:N],y1,'ko-');xlabel('iterationindexk');ylabel('ouputy(k)');figure(2);plot([1:N],state1,'-o');legend('x_1(k)','x_2(k)','x_3(k)');xlabel('iterationindexk');ylabel('statevariablesx(k
5、)');%%figure(3);plot([1:N],y2,'k*-');xlabel('iterationindexk');ylabel('ouputy(k)');figure(4);plot([1:N],state2,'-*');legend('x_1(k)','x_2(k)','x_3(k)');xlabel('iterationindexk');ylabel('statevariablesx(k)');仿真结果：用第一种方法绘制闭环系统在单位阶跃信号输入的情况下的响应曲线，以及状态变化曲线。用第二种方法绘制闭环系统在单位阶跃信号输入的情况下的响应曲线，以及
6、状态变化曲线。结论：由仿真得出对于同一个系统，由于状态变量的选择不同，可以得到不同的状态方程和输出方程，但最后的输出响应相同。所以状态变量的选取不影响系统的输出响应。上机作业25.24实验题目：太阳光源跟踪系统利用伺服系统控制太阳电池帆板的移动，使其跟踪并始终垂直于太阳光线，最大程度地接收太阳能。太阳光源跟踪系统由感光器与检测线路和电机功率放大器（简化为增益放大环节），太阳帆板（电机负载），电机位置传感器和直流力矩电机组成。太阳光源跟踪系统如下图（a）所示，计算机控制系统方框图如图（b）所示，其中。设计数字控制器，满足如下指标要求：（1）超调量，上升时间，调节时间
7、。（2）稳态速度误差系数，采样周期。（a）（b）实验要求：1）画出期望极点的允许范围；2）给出设计的具体推导过程，给出设计出的控制器脉冲传递函数D(z)（提示：注意系统稳定性）；绘出设计好的闭环系统阶跃响应曲线和误差曲线（建议使用simulink）。程序代码：clc;clear;closeall;%%%绘制根轨迹%num=[2*615.91];%den=[129.68690];%[n,d]=c2dm(num,den,0.1,'zoh');%[z,p,k]=tf2zp(n,d)%rlocus(n,d);Gs=tf([2*615.91],[1,29.6869,0]);
8、Gz=c2
