1、一、设计题目及要求太阳光源跟踪系统利用伺服系统控制太阳电池帆板的移动，使其跟踪并始终垂直于太阳光线，最大程度地接受太阳能。太阳光源跟踪系统由感光器与检测线路和电机的功率放大器（可以简化视为一个增益放大环节） ，太阳帆板（作为直流力矩电机的负载，可以近似看作常值转动惯量加到电机轴上） ，电机位置传感器（其输出与电机转角成正比的电压信号）和直流力矩电机组成。太阳光源跟踪系统如题图 a 所示。计算机控制系统方块图如题图 b 所示。试用连续域-离散化设计方法设计数字控制器，满足如下指标要求：（1）超调量 15%（2）上升时间 ；r0.ts（3）调节时间 。s（4）静态速度误差系数 。v5K(a)(b)
2、图 太阳光源跟踪计算机控制系统 E(z)TD(z)TZOH 2615.9(0)sT 设计要求：（1）计算未加控制器时的性能指标，并绘出仿真曲线；（2）设计连续域控制器 D(s)，写出设计步骤，验算加控制器后的性能指标，并绘出仿真曲线；（3）选用两种合适的离散化方法，将 D(s)离散为 D(z)。并绘制采样周期 T 分别为0.01s，0.05s ，0.1s 时，计算机控制系统的单位阶跃响应仿真曲线，记录时域指标，计算 。比较两种离散化方法的性能，并说明连续域-离散化设计与采样周rs%,t和 vK期 T 的关系。比较离散化前后系统的阶跃响应曲线，分析离散化后系统性能变化的原因。（4）最终选定你认为
3、最合适的一种离散化方法和采样周期。说明：所有的仿真都需有程序清单或 simulink 模型2、设计过程1、在未加控制器时对系统性能进行检测 2.813s0ss2）（ ）（）（ G编程如下：a=1 20 1231.82;b=0 0 1231.82;step(b,a,0:0.01:2)输出曲线如下图所示：Step ResponseTime (sec)Amplitude0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 200.20.40.60.811.21.4System: sysTime (sec): 0.0921Amplitude: 1.39System: sysTime
4、(sec): 0.0554Amplitude: 0.998System: sysTime (sec): 0.446Amplitude: 1.01由输出曲线可知 %=39%， rt=0.0554s， t=0.446s，计算得 Kv=61,可知超调量无法满足设计要求，即需要加控制器 D(s)进行控制，它的主要作用是降低超调量。2、设计符合要求的连续域控制器拟定超调量为 8.6%，调节时间为 0.2s，设计出【,555s29s555s1ss2ee ）（）（）（GG（ ）为 加入控制器后有的闭环系统。编程检验闭环控制系统性能：a=1 29 555;b=0 0 555;step(b,a,0:0.01:2
5、)输出曲线如下图所示：0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 200.20.40.60.811.21.4System: sysTime (sec): 0.17Amplitude: 1.09System: sysTime (sec): 0.12Amplitude: 0.999System: sysTime (sec): 0.29Amplitude: 1Step ResponseTime (sec)Amplitude由输出曲线可知 =9%，上升时间 0.12s,及调节时间 0.29s 均远远满足要求，并通过计算得 vK=19。故经控制器校正后的连续系统，其性能是符合
6、要求的。3、将 G（s）离散化为 G（z）（1）T=0.01s 时，编程如下：num=1231.82;den=1, 20,0; c,d=c2dm(num,den,0.01, zoh)输出结果为：c =0 0.0577 0.0540d =1.0000 -1.8187 0.8187则 212 z87.0z.85471.0z87.154)z( G（2）T=0.05s 时，编程如下：num=1231.82;den=1, 20,0; c,d=c2dm(num,den,0.05, zoh)输出结果为：c =0 1.1329 0.8137d =1.0000 -1.3679 0.3679则 212 z3679
7、.0z.829367.0z9.18)z( G（3）T=0.1s 时，编程如下：num=1231.82;den=1, 20,0; c,d=c2dm(num,den,0.1, zoh)输出结果为：c =0 3.4963 1.8292d =1.0000 -1.1353 0.1353则 212 z35.0z.9864315.0z3.29864)z( G4、对 D（s）进行离散化（1）使用阶跃响应不变法离散 D（s ）T=0.01s 时，离散过程如下：num=0.4505, 0.4505*20;den=1, 29; c,d=c2dm(num,den,0.01, zoh)输出结果为：c = 0.4505
8、-0.3723d = 1.0000 -0.7483则所以仿真分析其阶跃响应特性：num=0,0.0268,0.0029,-0.0207;den=1,-2.5402,2.1825,-0.6333;r=ones(1,101);k=0:100;c=filter(num,den,r);plot(k,c,r-);输出的单位阶跃响应仿真曲线如下：0 10 20 30 40 50 60 70 80 90 10000.20.40.60.811.21.4X: 15Y: 1.122过输出曲线可以观察到 =12.2%，上升时间及调节时间均明显满足要求。T=0.05s 时，离散过程如下：num=0.4505, 0.4
9、505*20;den=1, 29; c,d=c2dm(num,den,0.05, zoh)输出结果为：c = 0.4505 -0.2127d =1.0000 -0.2346所以在 T=0.05s 时， 212 z3679.0z.829367.0z9.18)z( G所以仿真分析其阶跃响应特性：num=0,0.5104,0.1256,-0.17307;den=1,-0.8575,0.8144,-0.25938;r=ones(1,101);k=0:100;c=filter(num,den,r);plot(k,c,r-);0 10 20 30 40 50 60 70 80 90 10000.20.40
10、.60.811.21.4X: 2Y: 1.074通过输出曲线可以观察到 =7.4%，能满足设计要求，但上升时间及调节时间均明显不满足要求， 也不满足要求。.)()(1)1(lim11zGzDzTKzvT=0.1s 时，离散过程如下：num=0.4505, 0.4505*20;den=1, 29; c,d=c2dm(num,den,0.1, zoh)输出结果为：c =0.4505 -0.1569d = 1.0000 -0.0550所以在 T=0.1s 时， 212 z35.0z.9864315.0z3.29864)z( G仿真分析其阶跃响应特性：num=0, 0.575,0.27549,-0.2
11、870;den=1,-0.3847,0.4732,-0.29444;r=ones(1,101);k=0:100;c=filter(num,den,r);plot(k,c,r-);0 10 20 30 40 50 60 70 80 90 10000.20.40.60.811.21.4X: 2Y: 1.072有图可以看出是不满足要求的。（2）使用零极点匹配法离散 D（s ）T=0.01s 时，离散过程如下：num=0.4505, 0.4505*20;den=1, 29; sysc=tf(num,den);sysd=c2d(sysc,0.01,matched)输出结果为：Transfer funct
12、ion:0.4315 z - 0.3533-z - 0.7483Sampling time: 0.01即 11z483.701z533.30315.40483.70z533.30z315.40z）（DT=0.01s 时， 212 878547.187.54)( G所以， 3213216317.01826.25421.2101908.0003.00249.0)()(1)()z(z zzzzzzzDzGzDG）（仿真分析其阶跃响应特性：num=0,0.0249,0.003,-0.01908;den=1,-2.5421,2.1826,-0.6317;r=ones(1,101);k=0:100;c=f
13、ilter(num,den,r);plot(k,c,r-);输出曲线如下图所示：0 10 20 30 40 50 60 70 80 90 10000.20.40.60.811.21.4X: 16Y: 1.129通过输出曲线可以观察到 =12.9%，上升时间及调节时间均明显满足要求。T=0.05s 时，离散过程如下：num=0.4505, 0.4505*20;den=1, 29; sysc=tf(num,den);sysd=c2d(sysc,0.05,matched)Transfer function:0.3762 z - 0.1384-z - 0.2346Sampling time: 0.05
14、即 11z346.201z384.10762.30346.20z.13840z.37620z）（DT=0.05s 时， 212 3679829367.9.18)( G所以， 3213211989.08381.01763.111126.01493.04262.0)()(1)()z(z zzzzzzzDzGzDG）（仿真分析其阶跃响应特性：num=0,0.4262,0.1493,-0.1126;den=1,-1.1763,0.8381,-0.1989;r=ones(1,101);k=0:100;c=filter(num,den,r);plot(k,c,r-);输出曲线如下图所示：0 10 20 3
15、0 40 50 60 70 80 90 10000.20.40.60.811.21.4X: 3Y: 1.372通过输出曲线可以观察到 =37.2%不能满足设计要求，上升时间及调节时间均明显满足要求。T=0.1s 时，离散过程如下：num=0.4505, 0.4505*20;den=1, 29; sysc=tf(num,den);sysd=c2d(sysc,0.1,matched)Transfer function:0.3395 z - 0.04595-z - 0.05502Sampling time: 0.111z5502.001z4595.00395.305502.00z.045950z.3
16、3950z）（D切有 T=0.1s 时， 212 z35.z.986435.z13.29864)( G所以， 3213210915.06582.00033.0108405.04604.0187.1)()(1)()z(z zzzzzzzDzGzDG）（仿真分析其阶跃响应特性：num=0,1.187,0.4604,-0.08405;den=1,-0.0033,0.6582,-0.0915;r=ones(1,101);k=0:100;c=filter(num,den,r);plot(k,c,r-);输出曲线如下图所示：0 10 20 30 40 50 60 70 80 90 10000.20.40.60.811.21.41.61.8有输出曲线可以观察到各性能指标均无法满足要求。4 选择：由上面的分析可以看出，采用阶跃响应不变法离散 D（s ） ，且 T=0.01 时效果最好。
