图1-22是仓库大门自动控制系统原理图。试说明系统自动控制大门开闭的工作原理并画出系统方框图。
可以先去尝试着分析这个系统是怎么运作的,可以看到,右下角是一个对称的电路(电位器桥式测量电路),当我们合上开门开关时,从开门开关处采集的电压值与从大门位置采集的电压值共同进入放大器,由于电阻分压原理,二者不等,产生偏差电压,经放大器放大后,驱动伺服电动机带动绞盘转动,提起大门,直到与大门相连的电位器电刷与开门开关位置处的电位器电刷水平。
图1-24为水温控制系统原理示意图。冷水在热交换器中由通入的蒸汽加热,从而得到一定温度的热水。冷水流量变化用流量计测量。试绘制系统方块图,并说明为了保持热水温度为期望值,系统是如何工作的?系统的被控对象和控制装置各是什么?
先分析工作原理,首先,水温控制系统,自然是要求水的温度可控,明确被控量与给定量都应该与温度有关。逐个分析:
之后,依据反馈系统的基本框图,绘制方框图:
结论:系统的被控对象是热交换器,被控量是热水温度,控制装置是温度控制器
此类系统可用线性微分方程式描述,一般形式为:
a 0 d n d t n c ( t ) + a 1 d n − 1 d t n − 1 c ( t ) + . . . + a n − 1 d d t c ( t ) + a n c ( t ) = b 0 d m d t m r ( t ) + b 1 d m − 1 d t m − 1 r ( t ) + . . . + b m − 1 d d t r ( t ) + b m r ( t ) a_0\frac{{\rm d}^n}{{\rm d}t^n}c(t)+a_1\frac{{\rm d}^{n-1}}{{\rm d}t^{n-1}}c(t)+...+a_{n-1}\frac{{\rm d}}{{\rm d}t}c(t)+a_nc(t)=b_0\frac{{\rm d}^m}{{\rm d}t^m}r(t)+b_1\frac{{\rm d}^{m-1}}{{\rm d}t^{m-1}}r(t)+...+b_{m-1}\frac{{\rm d}}{{\rm d}t}r(t)+b_mr(t) a0dtndnc(t)+a1dtn−1dn−1c(t)+...+an−1dtdc(t)+anc(t)=b0dtmdmr(t)+b1dtm−1dm−1r(t)+...+bm−1dtdr(t)+bmr(t)
线性定常系统还可以根据输入量变化规律分为恒指控制系统,随动系统和程序控制系统(见课本12页)
离散系统是指系统的某处或多处的信号为脉冲序列或数码形式,因而信号在时间上是离散的。连续信号经过采样开关的采样就可以转换成离散信号。一般,在离散系统中既有连续的模拟信号,也有离散的数字信号。
因此,离散系统要用差分方程描述,线性差分方程的一般形式为:
a 0 c ( k + n ) + a 1 c ( k + n − 1 ) + . . . + a n − 1 c ( k + 1 ) + a n c ( k ) = b 0 r ( k + m ) + b 1 r ( k + m − 1 ) + . . . + b m − 1 r ( k + 1 ) + b m r ( k ) a_0c(k+n)+a_1c(k+n-1)+...+a_{n-1}c(k+1)+a_nc(k)=b_0r(k+m)+b_1r(k+m-1)+...+b_{m-1}r(k+1)+b_mr(k) a0c(k+n)+a1c(k+n−1)+...+an−1c(k+1)+anc(k)=b0r(k+m)+b1r(k+m−1)+...+bm−1r(k+1)+bmr(k)
线性系统系统特性:
满足叠加性和齐次性。
如果输入 r 1 ( t ) ⟹ r_1(t)\implies r1(t)⟹输出 y 1 ( t ) y_1(t) y1(t),输出 r 2 ( t ) ⟹ r_2(t)\implies r2(t)⟹输出 y 2 ( t ) y_2(t) y2(t)
则输入 a r 1 ( t ) + b r 2 ( t ) ⟹ ar_1(t)+br_2(t)\implies ar1(t)+br2(t)⟹输出 a y 1 ( t ) + b y 2 ( t ) ay_1(t)+by_2(t) ay1(t)+by2(t)
系统中只要有一个元部件的输入-输出特性是非线性的,这类系统就称为非线性控制系统。这时,要用非线性微分(或差分)方程描述其特性。
非线性方程的特点:系数与变量有关,或者方程中含有变量及其导数的高次幂或乘积项。
下列各式是描述系统的微分方程,其中 c ( t ) c(t) c(t)为输出量, r ( t ) r(t) r(t)为输入量,试判断那些是线性定常或时变系统,哪些是非线性系统?
(1) c ( t ) = 5 + r 2 ( t ) + t d 2 r ( t ) d t 2 c(t)=5+r^2(t)+t\frac{{\rm d}^2r(t)}{{\rm d}t^2\\} c(t)=5+r2(t)+tdt2d2r(t)
(2) d 3 c ( t ) d t 3 + 3 d 2 c ( t ) d t 2 + 6 d c ( t ) d t + 8 c ( t ) = r ( t ) \frac{{\rm d}^3c(t)}{{\rm d}t^3}+3\frac{{\rm d}^2c(t)}{{\rm d}t^2}+6\frac{{\rm d}c(t)}{{\rm d}t}+8c(t)=r(t) dt3d3c(t)+3dt2d2c(t)+6dtdc(t)+8c(t)=r(t)
(3) t d c ( t ) d t + c ( t ) = r ( t ) + 3 d r ( t ) d t t\frac{{\rm d}c(t)}{{\rm d}t}+c(t)=r(t)+3\frac{{\rm d}r(t)}{{\rm d}t} tdtdc(t)+c(t)=r(t)+3dtdr(t)
(4) c ( t ) = r ( t ) c o s w t + 5 c(t)=r(t)coswt+5 c(t)=r(t)coswt+5
(5) c ( t ) = 3 r ( t ) + 6 d r ( t ) d t + 5 ∫ − ∞ t r ( τ ) d τ c(t)=3r(t)+6\frac{{\rm d}r(t)}{{\rm d}t}+5\int_{-\infty}^tr(\tau){\rm d}\tau c(t)=3r(t)+6dtdr(t)+5∫−∞tr(τ)dτ
(6) c ( t ) = r 2 ( t ) c(t)=r^2(t) c(t)=r2(t)
(7) c ( t ) = { 0 , t<6 r ( t ) , t>=6 c(t)=\begin{cases}0, & \text{t<6} \\ r(t), & \text{t>=6} \\ \end{cases} c(t)={0,r(t),t<6t>=6
(4)还可以用叠加定理验证:
令 c 1 ( t ) = r 1 ( t ) c o s w t + 5 , c 2 ( t ) = r 2 ( t ) c o s w t + 5 c_1(t)=r_1(t)coswt+5,c_2(t)=r_2(t)coswt+5 c1(t)=r1(t)coswt+5,c2(t)=r2(t)coswt+5
若 r 3 ( t ) = a r 1 ( t ) + b r 2 ( t ) , 可 推 导 c 3 ( t ) = a r 1 ( t ) c o s w t + b r 2 ( t ) c o s w t + 5 r_3(t)=ar_1(t)+br_2(t),可推导c_3(t)=ar_1(t)coswt+br_2(t)coswt+5 r3(t)=ar1(t)+br2(t),可推导c3(t)=ar1(t)coswt+br2(t)coswt+5
但 a c 1 ( t ) + b c 2 ( t ) = a r 1 ( t ) c o s w t + b r 2 ( t ) c o s w t + 5 ( a + b ) ac_1(t)+bc_2(t)=ar_1(t)coswt+br_2(t)coswt+5(a+b) ac1(t)+bc2(t)=ar1(t)coswt+br2(t)coswt+5(a+b)
二者明显不等
L = J w L=Jw L=Jw
M ⃗ = d L ⃗ d t \vec{M}=\frac{{\rm d}\vec{L}}{{\rm d}t} M=dtdL
单根:如果分母多项式有n个单根,分别是 p 1 , p 2 , . . . , p n p_1,p_2,...,p_n p1,p2,...,pn
对于 F ( s ) = N ( s ) D ( s ) = a 0 s m + a 1 s m − 1 + . . . + a m b 0 s n + b 1 s n − 1 + . . . + b n F(s)=\frac{N(s)}{D(s)}=\frac{a_0s^m+a_1s^{m-1}+...+a_m}{b_0s^n+b_1s^{n-1}+...+b_n} F(s)=D(s)N(s)=b0sn+b1sn−1+...+bna0sm+a1sm−1+...+am
可以展开为: F ( s ) = K 1 s − p 1 + K 2 s − p 2 + 。 。 。 + K n s − p n F(s)=\frac{K_1}{s-p_1}+\frac{K_2}{s-p_2}+。。。+\frac{K_n}{s-p_n} F(s)=s−p1K1+s−p2K2+。。。+s−pnKn
得 K i = [ ( s − p i ) F ( s ) ] s = p i K_i=[(s-p_i)F(s)]_{s=p_i} Ki=[(s−pi)F(s)]s=pi
共轭复根:略
重根:略
例: G ( s ) = b 1 s + b 2 a 0 s 2 + a 1 s + a 2 G(s)=\frac{b_1s+b_2}{a_0s^2+a_1s+a^2} G(s)=a0s2+a1s+a2b1s+b2
由传递函数可得s的代数方程 ( a 0 s 2 + a 1 s + a 2 ) C ( s ) = ( b 1 s + b 2 ) R ( s ) (a_0s^2+a_1s+a_2)C(s)=(b_1s+b_2)R(s) (a0s2+a1s+a2)C(s)=(b1s+b2)R(s),在零初始条件下,用微分算符 d d t \frac{{\rm d}}{{\rm d}t} dtd置换s,得到相应的微分方程
a 0 d 2 d t 2 c ( t ) + a 1 d d t c ( t ) + a 2 c ( t ) = b 1 d d t r ( t ) + b 2 r ( t ) a_0\frac{{\rm d}^2}{{\rm d}t^2}c(t)+a_1\frac{{\rm d}}{{\rm d}t}c(t)+a_2c(t)=b_1\frac{{\rm d}}{{\rm d}t}r(t)+b_2r(t) a0dt2d2c(t)+a1dtdc(t)+a2c(t)=b1dtdr(t)+b2r(t)
传递函数解题:对于初态不为0的情况,先用交叉相乘法还原为微分方程,再用拉氏变换法求解;对于初态为0的情况,直接求解
试证明图2-49(a)的电网络与图2-49(b)的机械系统有相同的数学模型
首先,写出左图的传递函数。
回顾知识点:复频域中的电路定律与电路模型(运算电路)
转换为复频域模型:直接列写传递函数方程
其次,再来看右图,首先,我们要明确阻尼器的运动机理。
阻尼器:产生摩擦力,大小与速度成正比,方向与运动方向相反。
观察等式,对于系统上部而言,上部产生了向下的位移,我们假设是由一个向下的力造成的,则第一个等式的左半部分,便是与这个力相平衡的力,由于力具有相互作用性,等式左端的力会给系统中间部分产生一个向下的反作用力,使得中间部分向下移动,此时,下边的阻尼器为保持平衡,产生向上的力与之平衡。同理,下面的阻尼器也会产生反作用力作用与最下面的弹簧,也就是第二个等式。
我们将二者的传递函数抽离出来:
G a ( s ) = R 1 R 2 C 1 C 2 s 2 + ( R 1 C 1 + R 2 C 2 ) s + 1 R 1 R 2 C 1 C 2 s 2 + ( R 1 C 1 + R 2 C 2 + R 1 C 2 ) s + 1 G_a(s)=\frac{R_1R_2C_1C_2s^2+(R_1C_1+R_2C_2)s+1}{R_1R_2C_1C_2s^2+(R_1C_1+R_2C_2+R_1C_2)s+1} Ga(s)=R1R2C1C2s2+(R1C1+R2C2+R1C2)s+1R1R2C1C2s2+(R1C1+R2C2)s+1
G b ( s ) = f 1 f 2 K 1 K 2 s 2 + ( f 1 K 1 + f 2 K 2 ) s + 1 f 1 f 2 K 1 K 2 s 2 + ( f 1 K 1 + f 2 K 2 + f 1 K 2 ) s + 1 G_b(s)=\frac{\cfrac{f_1f_2}{K_1K_2}s^2+(\cfrac{f_1}{K_1}+\cfrac{f_2}{K_2})s+1}{\cfrac{f_1f_2}{K_1K_2}s^2+(\cfrac{f_1}{K_1}+\cfrac{f_2}{K_2}+\cfrac{f_1}{K_2})s+1} Gb(s)=K1K2f1f2s2+(K1f1+K2f2+K2f1)s+1K1K2f1f2s2+(K1f1+K2f2)s+1
可以看到,二者的微分方程模型阶次相同。而且如果选择适当的参数之后: R 1 = f 1 , R 2 = f 2 , C 1 = 1 K 1 = 1 , C 2 = 1 K 2 = 1 R_1=f_1,R_2=f_2,C_1=\cfrac{1}{K_1}=1,C_2=\cfrac{1}{K_2}=1 R1=f1,R2=f2,C1=K11=1,C2=K21=1。当两系统的输入按照相同的规律变化时,两系统的输出相应曲线也相同。
满足上述条件的系统被称为相似系统(即本题中的相同数学模型),需要注意的是仅仅满足阶次相同是不够的,只有当令微分方程对应的系数相等,能够求解出适当的系统参数时,才可以称作相似系统。