太阳光源跟踪系统设计题

1、一、设计题目及要求太阳光源跟踪系统利用伺服系统控制太阳电池帆板的移动，使其跟踪并始终垂直于太阳光线，最大程度地接受太阳能。太阳光源跟踪系统由感光器与检测线路和电机的功率放大器（可以简化视为一个增益放大环节），太阳帆板（作为直流力矩电机的负载，可以近似看作常值转动惯量加到电机轴上），电机位置传感器（其输出与电机转角成正比的电压信号）和直流力矩电机组成。太阳光源跟踪系统如题图a所示。计算机控制系统方块图如题图b所示。试用连续域-离散化设计方法设计数字控制器，满足如下指标要求：（1）超调量（2）上升时间；（3）调节时间。（4）静态速度误差系数。(a)E(z)TD(z)TZOH (b)图 太阳光源跟踪

2、计算机控制系统 设计要求：（1）计算未加控制器时的性能指标，并绘出仿真曲线；（2）设计连续域控制器D(s)，写出设计步骤，验算加控制器后的性能指标，并绘出仿真曲线；（3）选用两种合适的离散化方法，将D(s)离散为D(z)。并绘制采样周期T分别为0.01s，0.05s，0.1s时，计算机控制系统的单位阶跃响应仿真曲线，记录时域指标，计算。比较两种离散化方法的性能，并说明连续域-离散化设计与采样周期T的关系。比较离散化前后系统的阶跃响应曲线，分析离散化后系统性能变化的原因。（4）最终选定你认为最合适的一种离散化方法和采样周期。说明：所有的仿真都需有程序清单或simulink模型2、 设计过程1、

3、在未加控制器时对系统性能进行检测编程如下：a=1 20 1231.82;b=0 0 1231.82;step(b,a,0:0.01:2)输出曲线如下图所示：由输出曲线可知%=39%，=0.0554s，=0.446s，计算得Kv=61,可知超调量无法满足设计要求，即需要加控制器D(s)进行控制，它的主要作用是降低超调量。2、 设计符合要求的连续域控制器拟定超调量为8.6%，调节时间为0.2s，设计出【555s29s555s1ss2ee+=+=F）（）（）（GG,F（）为加入控制器后有的闭环系统。编程检验闭环控制系统性能：a=1 29 555;b=0 0 555;step(b,a,0:0.01:2

4、)输出曲线如下图所示：由输出曲线可知=9%，上升时间0.12s,及调节时间0.29s均远远满足要求，并通过计算得=19。故经控制器校正后的连续系统，其性能是符合要求的。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则（2）T=0.05s时，编程如下：num=1231.82;den=1, 20,0; c,d=c2dm(num,den,0.05, 'zo

5、h')输出结果为：c =0 1.1329 0.8137d =1.0000 -1.3679 0.3679则（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则4、 对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')输

6、出结果为：c = 0.4505 -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-'); 输出的单位阶跃响应仿真曲线如下：过输出曲线可以观察到=12.2%，上升时间及调节时间均明显满足要求。T=0.05s时，离散过程如下：num=0.4505, 0.4505*20;den=1, 29; c,d=c2dm(num,den,0.05,

7、 'zoh')输出结果为：c = 0.4505 -0.2127d =1.0000 -0.2346所以在T=0.05s时，所以仿真分析其阶跃响应特性：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-');通过输出曲线可以观察到=7.4%，能满足设计要求，但上升时间及调节时间均明显不满足要求，.)()(1)1(lim11+-=®zGzDzTKzv也不满足要求。T=0.1s时，离散过

8、程如下：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时，仿真分析其阶跃响应特性：num=0, 0.575,0.27549,-0.2870;den=1,-0.3847,0.4732,-0.29444;r=ones(1,101);k=0:100;c=filter(num,den,r);plot(k,c,'r-'); 有图可以看出是不满足要求的。（2） 使用零极点匹配法离散D（s）T=0.01

9、s时，离散过程如下：num=0.4505, 0.4505*20;den=1, 29; sysc=tf(num,den);sysd=c2d(sysc,0.01,'matched')输出结果为：Transfer function:0.4315 z - 0.3533- z - 0.7483 Sampling time: 0.01即11z483.701z533.30315.40483.70z533.30z315.40z-=-=）（DT=0.01s时，所以， 3213216317.01826.25421.2101908.0003.00249.0)()(1)()z(z-+-+=+=Fzzz

10、zzzzDzGzDG）（仿真分析其阶跃响应特性：num=0,0.0249,0.003,-0.01908;den=1,-2.5421,2.1826,-0.6317;r=ones(1,101);k=0:100;c=filter(num,den,r);plot(k,c,'r-');输出曲线如下图所示： 通过输出曲线可以观察到=12.9%，上升时间及调节时间均明显满足要求。T=0.05s时，离散过程如下：num=0.4505, 0.4505*20;den=1, 29; sysc=tf(num,den);sysd=c2d(sysc,0.05,'matched')Trans

11、fer function:0.3762 z - 0.1384- z - 0.2346 Sampling time: 0.05>>即11z346.201z384.10762.30346.20z.13840z.37620z-=-=）（DT=0.05s时，所以， 3213211989.08381.01763.111126.01493.04262.0)()(1)()z(z-+-+=+=FzzzzzzzDzGzDG）（仿真分析其阶跃响应特性：num=0,0.4262,0.1493,-0.1126;den=1,-1.1763,0.8381,-0.1989;r=ones(1,101);k=0:1

12、00;c=filter(num,den,r);plot(k,c,'r-');输出曲线如下图所示： 通过输出曲线可以观察到=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.05502 Sampling time: 0.111z5502.001z4595.00395.305502.00z.045950z.33950z-=-=）（D切有T=0.1s时，所以， 3213210915.06582.00033.0108405.04604.0187.1)()(1)()z(z-+-+=+=FzzzzzzzDzGzDG）（仿真分析其阶跃响应特性：num