Chương 2 : Điều khiển tối ưu

Chia sẻ: Nguyễn Thị Hà Hạnh | Ngày: | Loại File: PDF | Số trang:132

Thêm vào BST

Báo xấu

153
lượt xem 31
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Một hệ điều khiển được thiết kế ở chế độ làm việc tốt nhất là hệ luôn ở trạng thái tối ưu theo một tiêu chuẩn chất lượng nào đó ( đạt được giá trị cực trị ). Trạng thái tối ưu có đạt được hay không tùy thuộc vào yêu cầu chất lượng đặt ra, vào sự hiểu biết về đối tượng và các tác động lên đối tượng

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Chương 2 : Điều khiển tối ưu

Chương 2: ði u khi n t i ưu Chương 2 ðI U KHI N T I ƯU 2.1 CH T LƯ NG T I ƯU 2.1.1 ð c ñi m c a bài toán t i ưu 1. Khái ni m M t h ñi u khi n ñư c thi t k ch ñ làm vi c t t nh t là h luôn tr ng thái t i ưu theo m t tiêu chu n ch t lư ng nào ñó ( ñ t ñư c giá tr c c tr ) . Tr ng thái t i ưu có ñ t ñư c hay không tùy thu c vào yêu c u ch t lư ng ñ t ra , vào s hi u bi t v ñ i tư ng và các tác ñ ng lên ñ i tư ng , vào ñi u ki n làm vi c c a h ñi u khi n … M t s ký hi u s d ng trong chương 2. Hình 2.1: Sơ ñ h th ng ñi u khi n . H th ng ñi u khi n như hình trên bao g m các ph n t ch y u : ñ i tư ng ñi u khi n ( ðTðK ) , cơ c u ñi u khi n ( CCðK ) và vòng h i ti p ( K ) . V i các ký hi u : r : tín hi u ñ u vào, m c tiêu ñi u khi n, ñáp ng mong mu n c a h th ng. u : tín hi u ñi u khi n, lu t ñi u khi n. x : tín hi u ñ u ra, ñáp ng ra c a h th ng. ε = r – x : sai l ch c a h th ng. f : tín hi u nhi u Ch tiêu ch t lư ng J c a m t h th ng có th ñư c ñánh giá theo sai l ch c a ñ i lư ng ñư c ñi u khi n x so v i tr ñáp ng mong mu n r , lư ng quá ñi u khi n ( tr s c c ñ i xmax so v i tr s xác l p x ( ∞ ) tính theo ph n trăm ) , th i gian quá ñ … hay theo m t ch tiêu h n h p trong ñi u ki n làm vi c nh t ñ nh như h n ch v công su t , t c ñ , gia t c … Do ñó vi c ch n m t lu t ñi u khi n và cơ c u ñi u khi n ñ ñ t ñư c ch ñ làm vi c t i ưu J ñ t c c tr còn tùy thu c vào lư ng thông tin ban ñ u mà ta có ñư c. ñây chúng ta có th th y ñư c s khác bi t v k t qu nh n ñư c ch t lư ng t i ưu khi lư ng thông tin ban ñ u thay ñ i ( Hình 2.2 ) . Trang 125
Chương 2: ði u khi n t i ưu Hình 2.2 : T i ưu c c b và t i ưu toàn c c . Khi tín hi u ñi u khi n u gi i h n trong mi n [u1,u2] , ta có ñư c giá tr t i ưu c c ñ i J1∗ c a ch tiêu ch t lư ng J ng v i tín hi u ñi u khi n u1 . ∗ Khi tín hi u ñi u khi n u không b ràng bu c b i ñi u ki n u1 ≤ u ≤ u2 , ta có ñư c giá tr t i ưu J 2 > J1∗ ng v i u2 . Như v y giá tr t i ưu th c s ∗ ∗ bây gi là J 2 . ∗ T ng quát hơn , khi ta xét bài toán trong m t mi n [um , un ] nào ñó và tìm ñư c giá tr t i ưu J i∗ thì ñó là giá tr t i ưu c c b . Nhưng khi bài toán không có ñi u ki n ràng bu c ñ i v i u thì giá tr t i ưu là J ∗ = extremum( J i∗ ) v i J i∗ là các giá tr t i ưu c c b , giá tr J ∗ chính là giá tr t i ưu toàn c c . ði u ki n t n t i c c tr : • ð o hàm b c m t c a J theo u ph i b ng 0 : ∂J =0 ∂u • Xét giá tr ñ o hàm b c hai c a J theo u t i ñi m c c tr : ∂2J > 0 : ñi m c c tr là c c ti u ∂u 2 ∂2J < 0 : ñi m c c tr là c c ñ i ∂u 2 2. ði u ki n thành l p bài toán t i ưu Trang 126
Chương 2: ði u khi n t i ưu ð thành l p bài toán t i ưu thì yêu c u ñ u tiên là h th ng ph i có ñ c tính phi tuy n có c c tr . Bư c quan tr ng trong vi c thành l p m t h t i ưu là xác ñ nh ch tiêu ch t lư ng J . Nhi m v cơ b n ñây là b o ñ m c c tr c a ch tiêu ch t lư ng J . Ví d như khi xây d ng h t i ưu tác ñ ng nhanh thì yêu c u ñ i v i h là nhanh chóng chuy n t tr ng thái này sang tr ng thái khác v i th i gian quá ñ nh nh t , nghĩa là c c ti u hóa th i gian quá ñ . Hay khi tính toán ñ ng cơ tên l a thì ch tiêu ch t lư ng là vư t ñư c kho ng cách l n nh t v i lư ng nhiên li u ñã cho . Ch tiêu ch t lư ng J ph thu c vào tín hi u ra x(t) , tín hi u ñi u khi n u(t) và th i gian t . Bài toán ñi u khi n t i ưu là xác ñ nh tín hi u ñi u khi n u(t) làm cho ch tiêu ch t lư ng J ñ t c c tr v i nh ng ñi u ki n h n ch nh t ñ nh c a u và x . Ch tiêu ch t lư ng J thư ng có d ng sau : T J = ∫ L[ x(t ), u (t ), t ]dt 0 Trong ñó L là m t phi m hàm ñ i v i tín hi u x , tín hi u ñi u khi n u và th i gian t . L y ví d v bài toán ñi u khi n ñ ng cơ ñi n m t chi u kích t ñ c l p Φ kt = const v i tín hi u ñi u khi n u là dòng ñi n ph n ng iu và tín hi u ra x là góc quay ϕ c a tr c ñ ng cơ . Hình 2.3 : ð ng cơ ñi n m t chi u kích t ñ c l p . Ta có phương trình cân b ng moment c a ñ ng cơ : dω (1) k M iu − M c = M q dt dϕ (2) ω= dt Trang 127
Chương 2: ði u khi n t i ưu trong ñó k M = CM Φ = const ; Mq là moment quán tính ; ω là t c ñ góc ;ϕ là góc quay . Gi s b qua ph t i trên tr c ñ ng cơ ( M c = 0 ) thì : d 2ϕ (3) kM iu = M q dt 2 N u xét theo th i gian tương ñ i b ng cách ñ t : τ = t kM / M q thì (3) có d ng : d 2ϕ (4) = iu dτ 2 T ñó ta có : d 2x (5) =u dτ 2 V y phương trình tr ng thái c a ñ ng cơ ñi n là m t phương trình vi phân c p hai v i tín hi u ñi u khi n u là dòng ñi n ph n ng iư, tín hi u ra là góc quay ϕ . • Bài toán t i ưu tác ñ ng nhanh ( th i gian t i thi u ) : Tìm lu t ñi u khi n u(t) v i ñi u ki n h n ch u ≤ 1 ñ ñ ng cơ quay t v trí ban ñ u có góc quay và t c ñ ñ u b ng 0 ñ n v trí cu i cùng có góc quay b ng ϕ0 và t c ñ b ng 0 v i m t kho ng th i gian ng n nh t . Vì c n th i gian ng n nh t nên ch tiêu ch t lư ng J s là : T J = ∫ L[ x(t ), u (t ), t ]dt = T 0 Rõ ràng t phương trình trên ta ph i có L[ x(t ), u (t ), t ] = 1 . Như v y , ñ i v i bài toán t i ưu tác ñ ng nhanh thì ch tiêu ch t lư ng J có d ng : T J = ∫ 1dt = T 0 • Bài toán năng su t t i ưu : Năng su t ñây ñư c xác ñ nh b i góc quay l n nh t c a ñ ng cơ trong th i gian T nh t ñ nh . Khi ñó ch tiêu ch t lư ng J có d ng : T T J = ∫ L[ x(t ), u (t ), t ]dt = ϕT − ϕ0 = ∫ ϕ (t )dt & 0 0 Do ñó L[ x(t ), u (t ), t ] = ϕ (t ) = x(t ) và ta s có ch tiêu ch t lư ng J ñ i v i & & bài toán năng su t t i ưu như sau : T J = ∫ x ( t )dt & 0 Trang 128
Chương 2: ði u khi n t i ưu • Bài toán năng lư ng t i thi u : T n hao năng lư ng trong h th ng : T Q = ∫ U u iu dt 0 D a vào phương trình cân b ng ñi n áp : U u = iu Ru + keω dω và phương trình cân b ng moment : k M iu − M c = M q dt Ta tính ñư c : T T ke M c Q = ∫ U u iu dt = (ϕT − ϕ0 ) + ∫ Ru iu2 dt kM 0 0 ð có ñư c tiêu hao năng lư ng t i thi u , ta ch c n tìm c c ti u c a J : T T J = ∫ L[ x(t ), u (t ), t ]dt = ∫ iu2 dt 0 0 Mà dòng ñi n ph n ng iu ñây chính là tín hi u ñi u khi n u . Vì v y ch tiêu ch t lư ng J ñ i v i bài toán năng lư ng t i thi u có d ng : T J = ∫ u 2 (t )dt 0 3. T i ưu hoá tĩnh và ñ ng Chúng ta c n phân bi t hai d ng bài toán t i ưu hoá tĩnh và t i ưu hóa ñ ng . T i ưu hóa tĩnh là bài toán không ph thu c vào th i gian . Còn ñ i v i t i ưu hóa ñ ng thì th i gian cũng là m t bi n mà chúng ta c n ph i xem xét ñ n. 2.1.2 Xây d ng bài toán t i ưu 1. T i ưu hóa không có ñi u ki n ràng bu c M t hàm ch tiêu ch t lư ng vô hư ng L ( u ) ñư c cho trư c là m t hàm c a m t vector ñi u khi n hay m t vector quy t ñ nh u ∈ R m . Chúng ta c n ch n giá tr c a u sao cho L(u) ñ t giá tr nh nh t . ð gi i bài toán t i ưu , ta vi t chu i Taylor m r ng cho ñ bi n thiên c a L(u) như sau : 1 dL = LT du + du T Luu du + O(3) (2.1) u 2 V i O(3) là s h ng th 3. Grad c a L theo u là m t vector m c t : Trang 129
Chương 2: ði u khi n t i ưu  ∂L / ∂u1    ∂L  ∂L / ∂u 2  ∆ (2.2) = Lu ∂u   M   ∂L / ∂u m  và ñ o hàm c p 2 c a L theo u là m t ma tr n m x m ( còn g i là ma tr n Hessian ) :  ∂2L  2 ∆∂ L =  (2.3) Luu ∂u 2  ∂u i ∂ u j    Luu ñư c g i là ma tr n u n . M t ñi m c c tr ho c ñi m d ng xu t hi n khi s bi n thiên dL v i thành ph n th nh t ti n v 0 v i m i bi n thiên du trong quá trình ñi u khi n . Vì v y , ñ có ñi m c c tr thì : Lu = 0 (2.4) Gi s ñang t i ñi m c c tr , có Lu = 0 như (2.4) . ð ñi m c c tr tr thành ñi m c c ti u , chúng ta c n có : 1 dL = du T Luu du + O(3) (2.5) 2 là xác ñ nh dương v i m i s bi n thiên du . ði u này ñư c ñ m b o n u ma tr n u n Luu là xác ñ nh dương : Luu > 0 (2.6) N u Luu là xác ñ nh âm thì ñi m c c tr chính là ñi m c c ñ i ; còn n u Luu là không xác ñ nh thì ñi m c c tr chính là ñi m yên ng a . N u Luu là bán xác ñ nh thì chúng ta s xét ñ n thành ph n b c cao hơn trong (2.1) ñ xác ñ nh ñư c lo i c a ñi m c c tr . Nh c l i : Luu là xác ñ nh dương ( ho c âm ) n u như các giá tr riêng c a nó là dương ( ho c âm ) , không xác ñ nh n u các giá tr riêng c a nó v a có dương v a có âm nhưng khác 0 , và s là bán xác ñ nh n u t n t i giá tr riêng b ng 0 . Vì th n u Luu = 0 , thì thành ph n th hai s không hoàn toàn ch ra ñư c lo i c a ñi m c c tr . 2. T i ưu hóa v i các ñi u ki n ràng bu c Cho hàm ch tiêu ch t lư ng vô hư ng L( x, u ) , v i vector ñi u khi n u ∈ R m và vector tr ng thái x ∈ R n . Bài toán ñưa ra là ch n u sao cho hàm Trang 130
Chương 2: ði u khi n t i ưu ch tiêu ch t lư ng L(x,u) ñ t giá tr nh nh t và th a mãn ñ ng th i các phương trình ñi u ki n ràng bu c . f ( x, u ) = 0 (2.7) Vector tr ng thái x ñư c xác ñ nh t m t giá tr u cho trư c b ng m i quan h (2.7) , vì th f là m t h g m n phương trình vô hư ng , f ∈ R n . ð tìm ñi u ki n c n và ñ c a giá tr c c ti u , ñ ng th i th a mãn f ( x, u ) = 0 , ta c n làm chính xác như trong ph n trư c . ð u tiên ta khai tri n dL dư i d ng chu i Taylor , sau ñó xác ñ nh s h ng th nh t và th hai là Lu & Luu. Th a s Lagrange và hàm Hamilton . T i ñi m c c tr , dL v i giá tr th nh t b ng 0 v i m i s bi n thiên c a du khi df b ng 0 . Như v y chúng ta c n có: (2.8) dL = LT du + LT dx = 0 u x df = f u du + f x dx = 0 (2.9) T (2.7) ta xác ñ nh ñư c x t giá tr u ñã có, ñ bi n thiên dx ñư c xác ñ nh b i (2.9) t giá tr bi n thiên du v i ñi u ki n ma tr n Jacobi là không kỳ d f x ≠ 0 . Như v y , ma tr n Jacobi fx không kỳ d và : dx = − f x−1 f u du (2.10) Thay dx vào (2.8) ta ñư c : dL = ( LT − LT f x−1 f u )du (2.11) u x ð o hàm riêng c a L theo u ch a h ng s f ñư c cho b i phương trình : ∂L ( ) T = LT − LT f x−1 f u = Lu − f uT f x−T L x (2.12) u x ∂u df =0 v i f x−T = ( f x−1 ) . Lưu ý r ng : T ∂L (2.13) = Lu ∂u dx =0 ð thành ph n th nh t c a dL b ng không v i giá tr du tùy ý khi df = 0 , ta c n có : (2.14) Lu − f uT f x−T L x = 0 ðây là ñi u ki n c n ñ có giá tr c c ti u . Trư c khi ñi tìm ñi u ki n ñ , chúng ta hãy xem xét thêm m t vài phương pháp ñ có ñư c (2.14) . Vi t (2.8) và (2.9) dư i d ng: dL   LT LT   dx    = 0 (2.15) x u  df  =   f x f u  du  Trang 131
Chương 2: ði u khi n t i ưu H phương trình tuy n tính này xác ñ nh m t ñi m d ng , và ph i có m t [ ]T k t qu dx T du T . ði u này ch x y ra n u ma tr n h s (n + 1) × (n + m ) có h ng nh hơn n+1 . Có nghĩa là các hàng c a ma tr n tuy n tính v i nhau ñ t n t i m t vector λ có n s h ng như sau: T L LT  T [ ] 1 λ . x =0 (2.16) u f x fu   Hay: (2.17) LT + λT f x = 0 x (2.18) LT + λT f u = 0 u Gi i (2.17) ta ñư c λ : λT = − LT f x−1 (2.19) x và thay vào (2.18) ñ có ñư c (2.14) . Vector λ ∈ R n ñư c g i là th a s Lagrange , và nó s là công c h u ích cho chúng ta sau này . ð hi u thêm ý nghĩa c a th a s Lagrange ta xét du = 0 , t (2.8) và (2.9) ta kh dx ñ ñư c : dL = LT f x−1df (2.20) x Vì v y: ∂L = (LT f x−1 ) = −λ T (2.21) x ∂f du =0 Do ñó -λ là ñ o hàm riêng c a L v i bi n ñi u khi n u là h ng s . ði u này nói lên tác d ng c a hàm ch tiêu ch t lư ng v i bi n ñi u khi n không ñ i khi ñi u ki n ràng bu c thay ñ i . Như là m t cách th ba ñ tìm ñư c (2.14), ta phát tri n thêm ñ s d ng cho các phân tích trong nh ng ph n sau. K t h p ñi u ki n ràng bu c và hàm ch tiêu ch t lư ng ñ thành l p hàm Hamilton . H ( x, u, λ ) = L( x, u ) + λT f ( x, u ) (2.22) V i λ ∈ R là th a s Lagrange chưa xác ñ nh . Mu n ch n x , u , λ ñ có n ñư c ñi m d ng , ta ti n hành các bư c sau . ð bi n thiên c a H theo các ñ bi n thiên c a x , u , λ ñư c vi t như sau : (2.23) T T T dH = H x dx + H u du + H λ dλ ∂H Lưu ý r ng : = f ( x, u ) (2.24) Hλ = ∂λ Gi s chúng ta ch n các giá tr c a u th a mãn: Hλ = 0 (2.25) Trang 132
Chương 2: ði u khi n t i ưu Sau ñó ta xác ñ nh x v i giá tr c a u ñã có b ng phương trình ñi u ki n ràng bu c f ( x, u ) = 0 . Trong trư ng h p này hàm Hamilton tương ñương v i hàm ch tiêu ch t lư ng: (2.26) H f =0 = L Nh c l i : n u f = 0 , ta s tìm ñư c dx theo du t (2.10) . Ta không nên xét m i quan h gi a du và dx ñ thu n ti n trong vi c ch n λ sao cho : Hx = 0 (2.27) ð o hàm (2.22) theo x: ∂H = L x + f xT λ = 0 (2.28) ∂x hay λT = − LT f x−1 . x N u gi nguyên (2.25) và (2.27) t (2.23): (2.29) T dL = dH = H u du Vì H = L, ñ có ñư c ñi m d ng ta ph i áp ñ t ñi u ki n: Hu = 0 (2.30) Tóm l i , ñi u ki n c n ñ có ñư c ñi m c c ti u c a L(x,u) th a mãn ñi u ki n ràng bu c f(x,u) = 0 g m có : ∂H = f =0 (2.31a) ∂λ ∂H = L x + f xT λ = 0 (2.31b) ∂x ∂H = Lu + f uT λ = 0 (2.31c) ∂u V i H ( x, u, λ ) xác ñ nh b i (2.22) . Cách thư ng dùng là t 3 phương trình ñã cho xác ñ nh x , λ , và u theo th t tương ng . So sánh 2 phương trình (2.31b) và (2.31c) ta th y chúng tương ng v i 2 phương trình (2.17) và (2.18) . Trong nhi u ng d ng , chúng ta không quan tâm ñ n giá tr c a λ , tuy nhiên ta v n ph i ñi tìm giá tr c a nó vì ñó là m t bi n trung gian cho phép chúng ta xác ñ nh các ñ i lư ng c n tìm là u , x và giá tr nh nh t c a L . Ưu ñi m c a th a s Lagrange có th tóm t t như sau : trên th c t , hai ñ i lư ng dx và du không ph i là hai ñ i lư ng bi n thiên ñ c l p v i nhau , theo (2.10) . B ng cách ñưa ra m t th a s b t ñ nh λ , chúng ta ch n λ sao cho dx và du có th ñư c xem là hai ñ i lư ng bi n thiên ñ c l p v i nhau . L y ñ o hàm riêng c a H l n lư t theo các bi n như trong (2.31) , như th ta s có ñư c ñi m d ng . Trang 133
Chương 2: ði u khi n t i ưu Khi ñưa ra th a s Lagrange , chúng ta có th thay th bài toán tìm giá tr nh nh t c a L(x,u) v i ñi u ki n ràng bu c f(x,u) = 0 , thành bài toán tìm giá tr nh nh t c a hàm Hamilton H(x,u,λ) không có ñi u ki n ràng bu c . ði u ki n ñã (2.31) xác ñ nh m t ñi m d ng . Ta s ti p t c ch ng minh ñây là ñi m c c ti u như ñã th c hi n trong ph n trư c . Vi t chu i Taylor m r ng cho ñ bi n thiên c a L và f như sau : L L xu   dx   dx  1 [ ] [ ] dL = LT LT   + dx T du T  xx    + O(3) (2.32) x u du  2  Lux Luu  du  f f xu   dx   dx  1 [ ] df = [ f x f u ]  + dx T du T  xx + O(3) (2.33) f uu  du  du  2  f ux   V i: ∂2 f ∆ f xu = ∂u∂x ð ñưa ra hàm Hamilton , ta s d ng các phương trình sau : H H xu   dx  dL  T  dx  1 [1 λT ] [ ] [ ] du T  xx + O(3) T H u   + dx T  = Hx H uu  du  du  2  H ux  df    (2.34) Bây gi , ñ có ñư c ñi m d ng ta c n có f = 0 , và ñ ng th i thành ph n th nh t c a dL b ng 0 v i m i s bi n thiên c a dx và du . Vì f = 0 nên df = 0 , và ñi u này ñòi h i H x = 0 và H u = 0 như trong (2.31) . ð tìm ñi u ki n ñ cho ñi m c c ti u , chúng ta xét ñ n thành ph n th hai . ð u tiên , ta c n xem m i quan h gi a dx và du trong (2.34) . Gi s r ng chúng ta ñang ñi m c c tr nên H x = 0 , H u = 0 và df = 0 . T (2.10) suy ra: dx = − f x−1 f u du + O(2) (2.35) Thay vào (2.34) ta ñư c : H xu  − f x−1 f u  H 1 [ ] dL = du T − f uT f x−T I  xx  du + O(3) (2.36)  2  H ux H uu   I  ð ñ m b o ñây là ñi m c c ti u , dL trong (2.36) ph i dương v i m i s bi n thiên c a du . ði u này ñư c ñ m b o n u như ma tr n u n v i f luôn b ng 0 là xác ñ nh dương . H xu  − f x−1 f u  H ∆ [ ] Luu = Luu f = − f uT f x−T I  xx f   (2.37)  H ux H uu   I  = H uu − f uT f x−T H xu − H ux f x−1 f u + f uT f x−T H xx f x−1 f u Trang 134
Chương 2: ði u khi n t i ưu Lưu ý r ng n u ñi u ki n ràng bu c f ( x, u ) = 0 v i m i x và u thì (2.37) ñư c rút l i thành Luu phương trình (2.3) . N u (2.37) là xác ñ nh âm ( ho c không xác ñ nh ) thì ñi m d ng s là ñi m c c ñ i ( ho c ñi m yên ng a ) . 2.1.3 Ví d T i ưu hóa không có ñi u ki n ràng bu c Ví d 2.1 : Không gian toàn phương . Cho u ∈ R 2 và : 1 q q12  u + [s1 s 2 ]u L(u ) = u T  11 (1) 2 q12 q 22   ∆1 (2) = u T Qu + S T u 2 ði m c c tr ñư c xác ñ nh b i : Lu = Qu + S = 0 (3) u ∗ = −Q −1 S (4) v i u* dùng ñ ch bi n ñi u khi n t i ưu. Lo i c a ñi m c c tr ñư c xác ñ nh b ng cách xét ma tr n hessian (5) Luu = Q 2 ði m u* là c c ti u n u Luu > 0 ( q11 > 0 và q11 q 22 − q12 > 0 ) . Là ñi m c c 2 ñ i n u Luu < 0 ( q11 < 0 và q11 q 22 − q12 > 0 ) . N u Q < 0 , thì u* là ñi m yên ng a . N u Q = 0 , thì u* là ñi m kỳ d , chúng ta không th xác ñ nh ñư c ñó là c c ti u hay c c ñ i t Luu . B ng cách thay (4) vào (2) ta s tìm ñư c giá tr c a hàm ch tiêu ch t lư ng như sau : 1 1 ∆ L* = L(u * ) = S T Q −1QQ −1 S − S T Q −1 S = − S T Q −1 S (6) 2 2 Gi s cho L như sau: 1 T 1 1  u + [0 1]u (7) L= u 2 1 2   Khi ñó giá tr u t i ưu s là: 2 − 1 0  1  u * = − (8)   =   1 1  1 − 1 Trang 135
Chương 2: ði u khi n t i ưu là m t c c ti u , vì Luu > 0 . T (6) giá tr nh nh t c a L là L* = -1/2 . Các ñư ng ñ ng m c c a L(u) trong (7) ñư c v trong Hình 2.4 , v i u = [u1 u2]T . Các mũi tên là gradient .  u + u2  Lu = Qu + S =  1 (9)  u1 + 2u 2 + 1 Lưu ý r ng gradient luôn luôn vuông góc v i các ñư ng ñ ng m c và có hư ng là hư ng tăng L(u) . Chúng ta dùng d u “*” ñ ch giá tr t i ưu c a u và L c n tìm . Tuy nhiên ta thư ng b qua d u “*” . Hình 2.4 : Các ñư ng ñ ng m c và vector gradient . Ví d 2.2 : T i ưu hóa b ng tính toán vô hư ng . Ph n trên chúng ta ñã ñ c p phương pháp gi i bài toán t i ưu b ng cách s d ng các vector và gradient . Sau ñây ta s ti p c n bài toán v i m t cách nhìn khác , xem chúng như là nh ng ñ i lư ng vô hư ng . ð ch ng minh , ta xét (7) d ng: 1 L(u1 , u 2 ) = u12 + u1u 2 + u 2 + u 2 2 (1) 2 V i u1 , u 2 là các ñ i lư ng vô hư ng . ði m c c tr xu t hi n khi ñ o hàm riêng c a L theo t t c các ñ i s ph i b ng 0 : ∂L = u1 + u 2 = 0 (2a) ∂u1 Trang 136
Chương 2: ði u khi n t i ưu ∂L = u1 + 2u 2 + 1 = 0 (2b) ∂u 2 Gi i h phương trình trên ta ñư c : u1 = 1, u 2 = −1 (3) V y , ñi m c c tr là (1 ,-1) . Bi u th c (1) là m t d ng m r ng c a bi u th c (7) trong ví d 2.1 , như v y chúng ta v a tìm ñư c m t k t qu tương t b ng m t cách khác . T i ưu hóa có ñi u ki n ràng bu c Ví d 2.3 : Không gian toàn phương v i ñi u ki n ràng bu c tuy n tính . Gi s hàm ch tiêu ch t lư ng ñư c cho b i ví d 2.1 v i các ñ i lư ng vô hư ng u1 , u 2 ñư c thay th b ng x, u : 1 1  x   x 1 L( x, u ) = [x u ]  u  + [0 1]u  (1) 1 2   2  V i ñi u ki n ràng bu c : f ( x, u ) = x − 3 = 0 (2) Hàm Hamilton s là : 1 H = L + λT f = x 2 + xu + u 2 + u + λ ( x − 3) (3) 2 v i λ là m t ñ i lư ng vô hư ng . ði u ki n ñ có ñi m d ng theo (2.31) là : Hλ = x − 3 = 0 (4) Hx = x +u + λ = 0 (5) H u = x + 2u + 1 = 0 (6) Gi i (4) , (5) , (6) ta ñư c : x = 3 , u = -2 , λ = -1 . ði m d ng là : (x, u )∗ = (3,−2) (7) ð xác ñ nh (7) là ñi m c c ti u , tìm ma tr n u n theo (2.3) : Luu = 2 (8) f Luu >= 0 , vì th ( x, u ) = (3,−2 ) là ñi m c c ti u . ∗ f Các ñư ng ñ ng m c c a L(x,u) và ñi u ki n ràng bu c (2) ñư c v trong Hình 2.5 . Grad c a f(x,u) trong h t a ñ (x,u) ñư c vi t như sau:  f x  1 (9)  f  = 0   u   ñư c v trong Hình 2.5 . Và grad c a L(x,u) : Trang 137
Chương 2: ði u khi n t i ưu  Lx   x + u  (10)  L  =  x + 2u + 1  u   T i ñi m c c ti u (3,-2) , grad L(x,u) s có giá tr :  L x  1 (11)  L  = 0   u   C n lưu ý r ng gradf và gradL tương ñương v i nhau t i ñi m d ng . Có nghĩa là ñi m c c ti u xu t hi n khi ñi u ki n ràng bu c (2) là ñư ng ti p tuy n c a các ñư ng ñ ng m c c a L. Di chuy n hư ng d c theo ñư ng th ng f = 0 s làm tăng giá tr c a L . Ta tìm ñư c giá tr c a L t i ñi m c c ti u b ng cách thay x = 3, u = -2 vào (1) , ta ñư c L*=0,5 . Vì λ = -1 , gi nguyên giá tr u = -2 , thay ñ i ñi u ki n ràng bu c df ( d ch chuy n ñư ng th ng trong Hình 2.5 v phía ph i ) s làm tăng L(x,u) v i dL = -λdf = df . Ví d 2.4 : Hàm ch tiêu ch t lư ng d ng toàn phương v i ñi u ki n ràng bu c tuy n tính - Trư ng h p vô hư ng . Xét hàm ch tiêu ch t lư ng d ng toàn phương : 1  x2 y 2  L ( x, u ) =  2 + 2  (1) 2a b   V i ñi u ki n ràng bu c tuy n tính : f ( x, u ) = x + mu − c (2) 2 Các ñư ng ñ ng m c c a L(x,u) là nh ng ellip; n u L(x,u) = l /2, thì bán kính tr c chính và bán kính tr c ph là al và bl . ði u ki n ràng bu c f(x,u) là m t h các ñư ng th ng ch a thông s c . Xem Hình 2.6 ( lưu ý r ng u là bi n ñ c l p , v i x ñư c xác ñ nh b i f(x,u) = 0 ) . Hàm Hamilton là : 1  x2 u2  H =  2 + 2  + λ ( x + mu − c) (3) 2a b   Và ñi u ki n ñ có ñi m d ng : H λ = x + mu − c = 0 (4) x Hx = 2 + λ = 0 (5) a u H u = 2 + λm = 0 (6) b Trang 138
Chương 2: ði u khi n t i ưu Hình2.5 : Các ñư ng ñ ng m c c a L(x,u) và ñi u ki n ràng bu c f(x,u) . Hình 2.6 : Các ñư ng ñ ng m c c a L(x,u) và ñi u ki n ràng bu c f(x,u). ð gi i h phương trình này , trư c h t ta s d ng phương trình (6) ñ ñưa ra bi n ñi u khi n t i ưu theo th a s Lagrange . u = −b 2 mλ (7) Bây gi thay phương trình (7) vào (4) ñ kh u , k t h p v i (5) và ñư c vi t l i :  1 − b 2 m 2   x  c  1  (8) = 1  λ  0  2   a  Gi i ra ta ñư c giá tr c a ñi m d ng : a 2c (9) x= 2 a + b2m2 Trang 139
Chương 2: ði u khi n t i ưu c (10) λ=− a + b2m2 2 Thay (9) , (10) vào (7) , ta có ñư c giá tr u t i ưu : b 2 mc (11) u= 2 a + b2m2 ð xác ñ nh ñi m d ng là c c ti u , dùng (2.37) ñ tìm ra ma tr n u n : 1 m2 (12) f Luu = 2 + 2 b a Luu > 0 vì v y ta tìm ñư c m t ñi m c c ti u . f Thay (9) và (11) vào (1) ta ñư c giá tr t i ưu c a hàm ch tiêu ch t lư ng : c2 1 L* = (13) 2 a2 + b2m2 ð ki m ch ng (2.21) , lưu ý r ng: ∂L* ∂L* (14) = −λ = ∂f du =0 ∂c Gradf trong mi n (u,x) là :  f u  m (15)  f  = 1  x   ñư c bi u di n trong Hình 2.6 . GradL là : u  Lu   b 2  (16) L  =  x   x   a2  và t i ñi m d ng (11) , (9) s có giá tr : *  Lu  m  c (17) L  =  1  2 2 2   a +b m  x ði u này tương ng v i (15) , vì v y ñi m d ng xu t hi n khi f(x,u) = 0 là ñư ng ti p tuy n v i m t ñư ng ñ ng m c c a L(x,u) . Ví d 2.5 : Hàm ch tiêu ch t lư ng d ng toàn phương v i ñi u ki n ràng bu c tuy n tính . Bây gi ta t ng quát hóa ví d 2.4 v i vector x ∈ R n , u ∈ R m , f ∈ R n , λ ∈ Rn . Xét hàm ch tiêu ch t lư ng d ng toàn phương: Trang 140
Chương 2: ði u khi n t i ưu 1T 1 (1) x Qx + u T Ru L= 2 2 v i ñi u ki n ràng bu c tuy n tính : (2) f = x + Bu + c = 0 v i Q , R và B là các ma tr n , c là vector n hàng . Gi s Q ≥ 0 và R > 0 ( v i Q , R là ma tr n ñ i x ng ) . Các ñư ng ñ ng m c c a L(x,u) là các ñư ng ellip trong không gian , và f(x,u)=0 là m t ph ng c t ngang qua chúng . ði m d ng xu t hi n khi gradf và gradL song song v i nhau . Hàm Hamilton là : 1 1 H = x T Qx + u T Ru + λT ( x + Bu + c) (3) 2 2 và các ñi u ki n ñ có ñi m d ng là : H λ = x + Bu + c = 0 (4) H x = Qx + λ = 0 (5) H u = Ru + B λ = 0 (6) T ð gi i các phương trình trên , ñ u tiên ta dùng ñi u ki n (6) ñ tìm u theo λ: u = − R −1 B T λ (7) T (5) ta có : (8) λ = −Qx K t h p v i (4) ta ñư c : (9) λ = QBu + Qc dùng k t qu này thay vào (7) cho ta : u = − R −1 B T (QBu + Qc) (10) ( ) hay : −1 T −1 T I + R B QB u = − R B Qc (R + B QB )u = − B (11) T T Qc T T Vì R > 0 và B QB ≥ 0 , chúng ta có th tìm ngh ch ñ o c a (R + B QB) và vì th giá tr u t i ưu là : u = −( R + B T QB) −1 B T Qc (12) So sánh k t qu này v i (11) trong ví d 2.4 . Thay (12) vào (4) và (9) cho ta giá tr tr ng thái t i ưu và th a s Lagrange t i ưu : ) ( −1 ( ) (13) x = − I − B R + BT QB BT Q c B Q) c ( −1 λ = Q − QB ( R + BT QB ) (14) T B ng b ñ c a ngh ch ñ o ma tr n : Trang 141
Chương 2: ði u khi n t i ưu λ = (Q −1 + BR −1 B T ) c −1 (15) n u Q ≠ 0 . Các k t qu trên s rút l i thành k t qu c a ví d 2.4 trong trư ng h p vô hư ng . ð xác ñ nh bi n ñi u khi n (12) là m t c c ti u , ta s d ng (2.37) ñ xác ñ nh ma tr n u n là xác ñ nh dương v i giá tr c a R và Q ñư c gi i h n . (16) Luu = R + B T QB f S d ng (12) và (13) th vào (1) ta có ñư c giá tr t i ưu : [ ] 1 ( ) −1 L* = c T Q − QB R + B T QB B T Q c (17) 2 1 L* = c T λ (18) 2 ∂L * (19) =λ Vì th : ∂c Ví d 2.6 : Bài toán v i nhi u ñi u ki n ràng bu c . Tìm kho ng cách nh nh t gi a parabol : y = ax 2 + bx + d (1) v i ñư ng th ng : (2) y = x+c Xem Hình 2.7 . Trong bài toán này s có hai ñi u ki n ràng bu c : f1 ( x1 , y1 ) = y1 − ax12 − bx1 − d = 0 (3) f 2 ( x2 , y 2 ) = y 2 − x2 − c = 0 Và : (4) v i ( x1 , y1 ) là 1 ñi m trên parabol và ( x 2 , y 2 ) là 1 ñi m trên ñư ng th ng . Chúng ta ch n hàm ch tiêu ch t lư ng là m t n a c a bình phương kho ng cách gi a 2 ñi m này . 1 1 L( x1 , x 2 , y1 , y 2 ) = ( x1 − x 2 ) 2 + ( y1 − y 2 ) 2 (5) 2 2 ð gi i bài toán này , ta x lý b ng cách ñ t : ∆ f  ∆ x  ∆y  f = 1  , x = 1  , u = 1  (6)  f2   x2   y2  và s d ng cách ti p c n vector ; tuy nhiên , s k t h p gi a m t ñi u ki n ràng bu c tuy n tính và m t ñi u ki n phi tuy n s làm ph c t p thêm bài toán . Thay vào ñó ta s s d ng các ñ i lư ng vô hư ng . Trang 142
Chương 2: ði u khi n t i ưu Hình 2.7 : Bài toán v i nhi u ñi u ki n ràng bu c . ðưa ra m t th a s Lagrange cho m i ñi u ki n ràng bu c , hàm Hamilton là : 1 1 H = ( x1 − x 2 ) 2 + ( y1 − y 2 ) 2 + λ1 ( y1 − ax12 − bx1 − d ) + λ 2 ( y 2 − x 2 − c) (7) 2 2 Khi ñó , ñ có ñi m d ng ta c n có : H x1 = x1 − x 2 − 2aλ1 x1 − bλ1 = 0 (8) H x2 = − x1 + x2 − λ2 = 0 (9) H y1 = y1 − y 2 + λ1 = 0 (10) H y2 = − y1 + y 2 + λ2 = 0 (11) H λ1 = y1 − ax12 − bx1 − d = 0 (12) H λ2 = y 2 − x2 − c = 0 (13) Gi i (12) ñ có ñư c y1 như sau : y1 = ax12 + bx1 + d (14) T (9) và (11) , ta có : (15) λ2 = x 2 − x1 = y1 − y 2 và s d ng (14) v i y 2 = x 2 + c t (13) có ñư c k t qu sau : Trang 143
Chương 2: ði u khi n t i ưu x 2 − x1 = ax12 + bx1 + d − x 2 − c (16) Khi ñó : 1 ( ) ax12 + (b + 1) x1 + d − c (17) x2 = 2 Theo (10) và (11) , λ1 = -λ2 , v y t (15) và (17) ta có : λ1 = x1 − x 2 1 ( ) λ1 = − ax12 + (b − 1) x1 + d − c (18) 2 Cu i cùng , chú ý r ng (8) là : (2ax1 + (b − 1))λ1 = 0 (19) hay : ( ) (2ax1 + (b − 1)) ax12 + (b − 1) x1 + d − c = 0 (20) * Phương trình b c 3 (20) ñư c gi i ñ có giá tr t i ưu x1 t giá tr a, b, c, d cho trư c . N u ñư ng th ng c t ngang qua parabol thì giao ñi m s là k t qu t i ưu ( khi ñó λ1=λ2=0 ) ; ngư c l i , s có ch m t c p g n nhau nh t (x1,x2) , (y1,y2) . M t khi tìm ñư c x1 thì ta s tìm ñư c x2 , y1 và y2 l n lư t theo các phương trình (17) , (14) và (15) . Thay các giá tr t i ưu này vào (5) s cho chúng ta kho ng cách ng n nh t là 2 L * . Ví d 2.7 : Moät con taøu ñang di chuyeån vôùi vaän toác 10 haûi lyù moät giôø theo phöông hôïp vôùi phöông baéc moät goùc 30o (laáy theo chieàu kim ñoàng hoà) ñeå ñeán moät hoøn ñaûo. Giaû söû raèng taïi thôøi ñieåm t = 0, con taøu ñang ôû vò trí caùch hoøn ñaûo 30 haûi lyù veà phía Baéc vaø 20 haûi lyù veà phía ñoâng (töùc laø con taøu leäch khoûi phöông di chuyeån ban ñaàu). Baøi toaùn ñaët ra laø xaùc ñònh ñieåm gaàn nhaát treân phöông di chuyeån ban ñaàu maø con taøu caàn trôû veà ñeå ñi tôùi hoøn ñaûo. Tìm khoaûng caùch töø hoøn ñaûo ñeán ñieåm gaàn nhaát vaø thôøi gian ñeán ñöôïc ñieåm ñoù. G i i: N Caùch 1: giaûi baèng phöông phaùp hình hoïc: My M 30 Trong ∆ MMxO: MM x 30 = 1.5 ⇒ ∠MOM x = 56.30 300 tg∠MOM x = = OM x 20 OM = MM x + OM x = 20 2 + 30 2 = 36.0555 2 2 Mx O 20 Trang 144