爱吱声

标题: 拜读晨大《控制之道》及问题请教 [打印本页]

作者: 可梦之    时间: 2024-3-17 14:18
标题: 拜读晨大《控制之道》及问题请教
刚读完第一章《反馈、动态与稳定性》。虽然是讲控制的,但是对我现在研究电路也很有帮助。/ r+ y7 P4 |, B2 o0 S4 ^
4 Q6 r/ t5 j( r3 G' R9 M
这章讲明白了好几个点。虽然我之前也接触过,数学公式、程序验证也对,但是总觉得隔着一层,就是不知道为什么这样。晨大的书透到了本质,方便理解。
: b+ F; a5 G/ c/ w: ?* I1. 极点的磨平效应,零点的锐化作用。之前理解极点还容易点,零点一直不知道有什么作用。
6 j! l9 F$ F& _2 E2. 极点和时间常数互为倒数。之前我也用时间常数表达式,但是不知道这个概念。
( T2 F0 b4 B% D$ V; h9 z" v3. 欠阻尼,临界阻尼,过阻尼。我现在就遇到了临界阻尼的情况。8 [- t# Z& k% A# k( q% B
4. 极点实部和虚部的作用。我们的极点都是实数,那应该不会出现我之前担心的震荡(欠阻尼)现象。之前也看过类似的分析,但是这次看书才联系起来,解决了我的困惑。
' w( Y. h& n4 `8 C! I$ M' J% @: H2 Q) A' s1 e3 p% w
但还有一些问题没有想明白,想请教一下晨大
: @3 a+ b' z* P0 m( S7 o6 `1. 极点数多于零点数:没太搞明白解释。输出不能超前于输入,这个没问题。但是磨平作用不仅和极点的多少有关,和极点的大小也有关吧。零点的锐化也类似。那么有没有可能我有几个很大的极点,更多但是更小的零点,使得输出还是之后输入呢?3 `+ a, s8 y% Y0 N3 c
4 X+ _: X( t6 `* m& ~3 e3 b
我们电路分析中也有tranfer function,里面也有零点书和极点数。有的书说的是小于,有的地方说的是小于等于。之前也没当一回儿事。现在好奇了:什么时候会等于,电路输入输入的transfer function和闭环、开环控制有哪些想通之处?
' [8 l$ X) z4 k/ h7 U* k% q" p2 F4 z& x& p4 u% i
2. 我们对系统降阶之后,其对应的波形在0点附近会先下降出现负值再回来,类似临界阻尼中峰值会超过1. 当然我们输入不是step,而是ramp函数。这个肯定是降阶造成的误差,原系统中不可能为负数。这种现象有什么名称吗?. n/ B; D8 K2 x+ g2 |

8 X) v' r- y; g) h) J& z3. 图1-8中,0.35上面的实线是“稳态值”,0.35下面的虚线没有标注,我猜是“设定值”。那上升时间图片标的是到设定值的时间,下面文字写的是穿过“稳态值”的时间。这个有点困惑,但不影响理解意思。" C% o8 A/ d7 |' ?! H4 \7 }! X
& }; d6 s" l* Y% ^- i& Y
感谢感谢
作者: 数值分析    时间: 2024-3-17 20:26
都已经拿到书开始学习了?羡慕了,我的书还不知道在哪儿呢。。。
作者: 晨枫    时间: 2024-3-17 22:52
本帖最后由 晨枫 于 2024-3-17 09:19 编辑 0 b) K% Q2 R( m7 P" ]0 x: Y

( b. y: i$ z! ]8 E- d$ r9 [+ e能帮梦兄改善对零极点的理解,是我的荣幸。后续章节还请多提宝贵意见。
: ?% p7 Z4 c" Y2 {5 x# S9 v$ T/ p" e7 u0 y* n) \
从频域的角度,可以把每一个极点理解为多一个低通滤波环节,每一个零点为高通滤波环节。每一个频域环节有两个特征:幅度和相位。磨平还是锐化是从幅度来说的,但加一个极点,就在相位上滞后一个90度,加一个零点就超前90度,不受幅度影响。这个90度可以从正弦微分就变余弦,余弦积分就变正弦来理解。极点好比积分,零点好比微分。
3 }% W8 j0 l. h! Q+ D) b- @. R+ {' v8 d5 P& ?4 X6 w8 x
实数极点只是开环时不震荡,闭环了还是可以震荡的。在根轨迹上,实数极点永远在实轴上,不震荡;增益增加到一定程度后,即将分叉,那就是临界阻尼点了。
4 i( X7 {) n5 G: Y
4 t1 g" _5 P% N7 \这样好理解一点吗?
6 Z) y1 y( J8 f2 x% @6 k0 ]& Q
1 L: e7 i( @! [/ s- p极点数大于等于零点数的叫proper functino(正定函数),严格大于、不能等于的叫strictly proper function(严格正定函数)。物理世界里想的出来的都是严格正定函数的传递函数关系,这正是从输入不可能超前于输出来理解的。正定和非正定只有数学上的意义,就和时域的纯超前一样。很多定理只对正定或者严格正定函数适用,在实际上并不构成限制。2 W5 e: F! P$ b1 i. M! J; p
3 Y. Y6 G: I. a/ {, i) W6 s
传递函数只是定义动态系统的输入输出关系,干什么用都可以。传递函数本身的动态行为与开环控制相同,开环的输入-输出行为直接由传递函数决定。闭环把输出反馈回来,修正传递函数的行为。开环是一厢情愿的,想要达到某一结果,但是否达到既看不到也不关心;闭环是见招拆招的,一面实施修正动作,一面观察结果,作为进一步修正的基础。稳定的闭环最终达到想要的结果,不稳定的闭环弄巧成拙,越描越黑。
- p+ w4 b6 w" d5 W/ b, F7 l# h) C( D0 y, U1 S, K  c0 J4 q
我们对系统降阶之后,其对应的波形在0点附近会先下降出现负值再回来,类似临界阻尼中峰值会超过1.

9 K2 c4 Q. x/ u- K7 L8 f1 F. r- {2 i7 Q& ^- ^$ K) d, Y" U4 [1 ?0 D+ }
这个要想一下,能给一个波形图看看吗?
2 }% i& C- x. ^  E
9 z: o" u! p8 ?上升时间(rise time)的定义并不统一,用稳态值、设定值、设定值加减5%、10%的都有,我自己也没有统一起来,不好意思。
, j* y2 V# e' m- a
3 ]5 n& O# Z/ Z! O在实用中,这个只是定性地用一用,所以这些定义差异并不打紧。稳定时间(settling time)更加有用。
作者: 可梦之    时间: 2024-3-18 09:36
本帖最后由 可梦之 于 2024-3-18 10:11 编辑 1 q6 o( p5 `2 e3 A# ?$ }

0 P5 K) y$ [4 Y- v. }万分感谢!
; W8 ^- W$ u/ T+ R3 E, z* Q$ O8 A/ ]' Y$ }1 k
正定函数这个,书上证明系统是passive的时候用到了,不过理论证明太抽象,认怂跳过去了。你这样解释就清楚明白好多了。7 R- \* s0 c: c0 t6 ]5 i
: Y) p# E  b0 p/ t# H8 i
函数图片是这样的2 D; C* W1 V: d; D, |1 c
8 O  p2 {7 S* Z

3 N5 L* i; `1 l  c/ V2 E+ Y* }3 m7 Z+ y' m% L1 F- J5 _
原系统是单调的,降阶到二维的结果(一般不会这么低,为了突出这个现象)。poles是:[-0.2580, -6.5426], residue是 [0.3019   -1.1122]. zeros:    -2.0830(手算出来的,可能不对)。我们用的ramp输入。
- Y  O+ \* z# \7 ~2 c' f$ `) z
0 s3 \1 H8 M' E6 d2 h* g1 G后面还没有全看完,但有一个小问题。很多图片是彩色的,很漂亮。但是对我这种色弱的残障人士,波形图颜色和图示颜色匹配有些困难,尤其是颜色多的时候。之后出书是否可以考虑除了颜色,还可以用不同线段(实线,虚线。。。)或者符号(+.*)区分表示?
5 l6 k/ C) K1 \; ~0 w9 n. C" g" _+ \- Y3 p

作者: 晨枫    时间: 2024-3-18 09:54
可梦之 发表于 2024-3-17 19:36* A: N. n6 @8 X8 q
万分感谢!
  Z7 y6 y7 o( f' x* d' R/ C$ Y* H" y" q9 z" Z2 k6 [2 }
正定函数这个,书上证明系统是passive的时候用到了,不过理论证明太抽象,认怂跳过去了。你这 ...

+ s# k7 s' s! s, d. v* \7 _$ v; Y0 L
梦兄太客气。4 e. z* X6 ^! e7 `, J2 v
: E5 b0 \8 ]3 `- O
这个好像是non-minimum phase的样子,你是不是降阶后产生一个右半平面的零点?应该就是这东西在作怪。
5 @' e$ I* E* v! e2 ]8 ^6 T+ Y& D! B; H. o
用各种line style的办法已经用上了,但有时线太多,还是区分不过来。这个问题我以后会注意,尽量用粗细线、点划线等来区分。多谢提醒。
作者: 可梦之    时间: 2024-3-18 10:09
晨枫 发表于 2024-3-18 09:54% k0 G( l. a' N$ x) l
梦兄太客气。2 r0 r+ l6 Q  B( o% F, G
6 d  `) `/ r1 I3 N/ G
这个好像是non-minimum phase的样子,你是不是降阶后产生一个右半平面的零点?应该就是这 ...

: H& |' Q) f( n, V8 q7 U又学习到了,去了解了一下minimum phase。刚才贴图搞错了没显示出来。
# A, r1 c3 r2 H0 T" i1 y$ s2 p# b) f/ U3 I
我们用的pole-residue的格式。这个transfer function是 0.3019/(s+0.2580) - 1.1122/(s+6.5426) = -(s-2.0830)/(s+0.2580)*(s+6.5426). 的确零点在右平面,刚才算错了。牛!!!




欢迎光临 爱吱声 (http://aswetalk.net/bbs/) Powered by Discuz! X3.2