直流电机三种控制方法示例及比较

1、直流电机的三种控制方法讨论比较了三种跟踪定位点命令及减小负载扰动敏感性的技术。前馈控制积分反馈控制LQR最优二次型算法控制一、问题描述在电枢控制的直流电动机中，外加电压Va 控制电机转轴的转角速度。如图1 所示图 1 带负载电机工作示意图图 1 中显示了两种减小角速度对负载变化（改变电机负载的反向转矩）敏感性的方法。图 2 所示的是一个简单的直流电机模型，转矩 Td 作为电机的负载扰动。在该扰动下，必须使转速的变化减到最小。图2带负载电机结构模型模型参数入下所示：R = 2.0;L = 0.5;Km = 0.1; Kb = 0.1;Kf = 0.2;J = 0.02;% Ohms% Henry

2、s% torque and back emf constants% Nms 粘滞摩擦系数% kg.m2/s2转动惯量转矩和反电势常数首先构造一个直流电机的状态空间模型，有两个输入（Va， Td）和一个输出（ w）：h1= tf(Km,L R);% armature 电枢传递函数h2= tf(1,J Kf);% eqn of motion转动负载的传递函数dcm = ss(h2) * h1 , 1;% w = h2 * (h1*Va + Td)角速度函数dcm = feedback(dcm,Kb,1,1);% close back emf loop现在绘制阶跃输入电压为Va 的角速度响应曲线，如

3、图3 所示。右键点击图形，选择Characteristics可查看具体响应参数，上升时间，调节时间，峰值时间等。stepplot(dcm(1);图3带负载电机在输入单位阶跃电压信号Va 时的响应曲线二、前馈控制设计图 4 所示，是一个简单的前馈控制结构，通过控制器调节使角速度w 达到w_ref的指定值。图 4 前馈控制示意图前馈增益Kff应该设置为电压到角速度的直流增益的倒数Kff = 1/dcgain(dcm(1)。为了估计消除负载干扰的前馈控制器参数，模拟一个时间从5 秒到 10 秒范围的干扰转矩，干扰转矩为Td=-0.1 牛米，设定参考角速度为w_ref 1。t = 0:0.1:15;T

4、d = -0.1 * (t>5 & t<10);%装在扰动信号 Tdu = ones(size(t) ; Td;%输入信号为指定角速度w_ref=1 and 扰动信号 Tdcl_ff = dcm * diag(Kff,1);%在系统中增加前馈增益set(cl_ff,'InputName','w_ref','Td','OutputName','w');h = lsimplot(cl_ff,u,t);title('Setpoint tracking and disturbance rejec

5、tion')legend('cl_ff')图 5 前馈控制系统对干扰信号的响应曲线显然，前馈控制处理负载扰动不佳。三、 反馈控制设计图 6 所示为反馈控制结构示意图。图 6 反馈控制结构示意图为了实施零状态误差，我们在系统中使用积分控制C(s) = K/s ，式中的K 是待确定的。为了确定增益K，可以在电压到角速度的开环传递函数中使用根轨迹的方法，用 1/s 乘以开环传递函数。利用反馈控制系统开环传递函数在根轨迹图中确定积分环节的增益K。图 7 所示为加入积分的反馈控制系统的根轨迹图。h = rlocusplot(tf(1,1 0) * dcm(1);setoption

6、s(h,'FreqUnits','rad/sec');set(gca,'Xlim',-15 5,'Ylim',-15 15);图 7加入积分的反馈控制系统的根轨迹图按下曲线可读取增益值及相关信息。这里合理的K 值选择为K=5。注意单入单出系统设计工具提供了一个综合的图形用户界面来完成这样的设计。（相关内容请参考sisotool帮助）在同一个测试系统中，将这个新的反馈控制设计和前面的前馈设计进行比较。对应的时域响应曲线如图8所示。K=5;C = tf(K,1 0);%积分控制器的传递函数C=K/scl_rloc = feedback

7、(dcm * append(C,1),1,1,1);h = lsimplot(cl_ff, cl_rloc,u,t);set(cl_rloc,'InputName','w_ref','Td','OutputName','w');title('Setpoint tracking and disturbance rejection')legend('feedforward','feedback w/ rlocus','Location','Nort

8、hWest')图 8 前馈控制与反馈控制响应曲线图从图中看到绿色线是积分反馈控制的响应曲线，蓝色是前馈控制曲线。用根轨迹设计的积分反馈控制器对负载扰动的抑制效果比前馈控制要好很多。三、 LQR（线性二次型最优控制 算法）控制器设计为了进一步提高性能，如图9 所示，在原反馈控制结构中设计一个线性二次型最优控制器。图 9 LQR 控制系统结构示意图除了积分误差，LQR方案还用状态向量x=(i, w)来综合处理驱动电压Va。电压处理后表示为：Va = K1 * w + K2 * w/s + K3 * i式中的电流i 为电枢电流。为了有更好的干扰抑制效果，用了一个代价函数对大的积分误差进行处理

9、：C(20 q(t )2w(t )20.01Va (t) 2 )dt0式中 q(s)w( s) / s 。由这个代价函数计算出最优的LQR增益：dc_aug = 1 ; tf(1,1 0) * dcm(1);%将输入 w/s 增加到直流电机模型中K_lqr = lqry(dc_aug,1 0;0 20,0.01);%计算最后LQR增益为仿真实现建立闭环模型：P = augstate(dcm);%模型输入为 :Va, Td输出为 :w,xC = K_lqr * append(tf(1,1 0),1,1);%将积分控制器1/s 包含到 LQR控制器中OL = P * append(C,1);% open loop开环传递函数为OL三中设计的闭环频率响应曲线图如图 10 所示，左边图为设定角速度的频率响应曲线，右边图中为干扰信号 Td 的频率响应曲线。在曲线上点击可确定各曲线所代表的系统或者查看相关数据。bodeplot(cl_ff,cl_rloc,cl_lqr);图 10三种控制器波特图四、设计比较最后我们在仿真系统中比较这三种设计：h = lsimplot(cl_ff,cl_rloc,cl_lqr,u,t);title('Setpoint tracking and disturbance rejection')lege