Lý thuyết điều khiển tối ưu

87 915 1
Lý thuyết điều khiển tối ưu

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

Chng 1 : iu khin ti u Hc kì 1 nm hc 2005-2006 Chng 1 IU KHIN TI U Vài nét lch s phát trin lý thuyt điu khin . - Phng pháp bin phân c đin Euler_Lagrange 1766 . - Tiêu chun n đnh Lyapunov 1892 . - Trí tu nhân to 1950 . - H thng điu khin máy bay siêu nh 1955 . - Nguyên lý cc tiu Pontryagin 1956 . - Phng pháp quy hoch đng Belman 1957 . - iu khin ti u tuyn tính dng toàn phng LQR ( LQR : Linear Quadratic Regulator ) . - iu khin kép Feldbaum 1960 . - Thut toán di truyn 1960 . - Nhn dng h thng 1965 . - Logic m 1965 . - Lut điu khin h thng thích nghi mô hình tham chiu MRAS và b t chnh đnh STR 1970 ( MRAS : Model-Reference Adaptive System , STR : Self-Tuning Regulator ) . - H t hc Tsypkin 1971 . - Sn phm công nghip 1982 . - Lý thuyt bn vng 1985 . - Công ngh tính toán mm và điu khin tích hp 1985 . PGS.TS Nguyn Th Phng Hà 2 http://www.khvt.com 1.1 CHT LNG TI U 1.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 1 . Hình 1.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 : x 0 : tín hiu đu vào u : tín hiu điu khin x : tín hiu đu ra ε = x 0 – x : tín hiu sai lch 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 s mong mun x 0 , lng quá điu khin ( tr s cc đi x max 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 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 ca cht lng ti u khi lng thông tin ban đu thay đi ( Hình 1.2 ) . Chng 1 : iu khin ti u Trang 3 Hình 1.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 [u 1 ,u 2 ] , ta có đc giá tr ti u cc đi 1 J ∗ ca ch tiêu cht lng J ng vi tín hiu điu khin 1 u ∗ . Khi tín hiu điu khin u không b ràng buc bi điu kin 12 uuu ≤ ≤ , ta có đc giá tr ti u 21 JJ ∗ ∗ > ng vi 2 u ∗ . Nh vy giá tr ti u thc s bây gi là 2 J ∗ . Tng quát hn , khi ta xét bài toán trong mt min [ ] , mn uu nào đó và tìm đc giá tr ti u i J ∗ 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à () i JextremumJ ∗∗ = vi i J ∗ 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 : 0= ∂ ∂ u J • Xét giá tr đo hàm bc hai ca J theo u ti đim cc tr : 0 2 2 > ∂ ∂ u J : đim cc tr là cc tiu 0 2 2 < ∂ ∂ u J : đim cc tr là cc đi PGS.TS Nguyn Th Phng Hà 4 http://www.khvt.com 2. iu kin thành lp bài toá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 : 0 [(),(),] T JLxtuttdt= ∫ 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 i u và tín hiu ra x là góc quay ϕ ca trc đng c . Hình 1.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 : Chng 1 : iu khin ti u Trang 5 Mu c q d ki M M dt ω −= (1) d dt ϕ ω = (2) trong đó MM k C const=Φ= ; M q 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 ( 0 c M = ) thì : 2 2 Mu q d ki M dt ϕ = (3) Nu xét theo thi gian tng đi bng cách đt : / M q tk M τ = thì (3) có dng : 2 2 u d i d ϕ τ = (4) T đó ta có : 2 2 dx u d τ = (5) 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 . • 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 1u ≤ đ đ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à : 0 [(),(),] T JLxtuttdtT = = ∫ Rõ ràng t phng trình trên ta phi có [(),(),] 1Lxtutt = . 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 TdtJ 0 1 PGS.TS Nguyn Th Phng Hà 6 http://www.khvt.com • 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 : 0 00 [ (), (),] () TT T JLxtuttdt tdt ϕϕ ϕ ==−= ∫∫  Do đó [ (), (),] () ()Lxtutt t xt ϕ = =   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 : () 0 T Jxtdt= ∫  • Bài toán nng lng ti thiu : Tn hao nng lng trong h thng : 0 T uu QUidt= ∫ Da vào phng trình cân bng đin áp : uuu e UiRk ω = + và phng trình cân bng moment : Mu c q d ki M M dt ω −= Ta tính đc : 2 0 00 () TT ec uu T uu M kM QUidt Ridt k ϕϕ == −+ ∫∫  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 : 2 00 [(),(),] TT u JLxtuttdtidt== ∫ ∫ Mà dòng đin phn ng i u  đâ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 : 2 0 () T Jutdt= ∫ Chng 1 : iu khin ti u Trang 7 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 . 1.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 () 0 = uL đ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 m Ru ∈ . 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 : )3( 2 1 OduLduduLdL uu TT u ++= (1.1) Vi O(3) có th coi là s hng th 3 . Grad ca L theo u là mt vector m ct : ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ ∂∂ ∂∂ ∂∂ = ∂ ∂ Δ m u uL uL uL u L L / / / 2 1  (1.2) 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 ) : ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ ∂∂ ∂ = ∂ ∂ Δ ji uu uu L u L L 2 2 2 (1.3) L uu đ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ì : 0 = u L (1.4) Gi s đang  ti đim cc tr , có L u = 0 nh (1.4) .  đim cc tr tr thành đim cc tiu , chúng ta cn có : PGS.TS Nguyn Th Phng Hà 8 http://www.khvt.com )3( 2 1 OduLdudL uu T += (1.5) 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 L uu là xác đnh dng : 0> uu L (1.6) Nu L uu là xác đnh âm thì đim cc tr chính là đim cc đi ; còn nu L uu là không xác đnh thì đim cc tr chính là đim yên nga . Nu L uu là bán xác đnh thì chúng ta s xét đn thành phn bc cao hn trong (1.1) đ xác đnh đc loi ca đim cc tr . Nhc li : L uu 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 0= uu L , 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 ( ) uxL , , vi vector điu khin m Ru ∈ và vector trng thái n Rx ∈ . Bài toán đa ra là chn u sao cho hàm 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 . ( ) 0, = uxf (1.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 (1.7) , vì th f là mt h gm n phng trình vô hng , n Rf ∈ .  tìm điu kin cn và đ ca giá tr cc tiu , đng thi tha mãn () 0, =uxf , 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 . 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ó: 0=+= dxLduLdL T x T u (1.8) và: 0 = + = dxfdufdf xu (1.9) Chng 1 : iu khin ti u Trang 9 T (1.7) ta xác đnh đc x t giá tr u đã có, đ bin thiên dx đc xác đnh bi (1.9) t giá tr bin thiên du đã có . Nh vy , ma trn Jacobi f x không k d và : duffdx ux 1− −= (1.10) Thay dx vào (1.8) ta đc : duffLLdL ux T x T u )( 1− −= (1.11) o hàm riêng ca L theo u cha hng s f đc cho bi phng trình : () x T x T uu T ux T x T u df LffLffLL u L −− = −=−= ∂ ∂ 1 0 (1.12) vi () T x T x ff 1−− = . Lu ý rng : u dx L u L = ∂ ∂ =0 (1.13)  thành phn th nht ca dL bng không vi giá tr du tùy ý khi 0=df , ta cn có : 0=− − x T x T uu LffL (1.14) â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 (1.14) . Vit (1.8) và (1.9) di dng: 0= ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ du dx ff LL df dL ux T u T x (1.15) 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 kt qu [] T TT dudx . iu này ch xy ra nu ma trn h s ( )( ) mnn + × + 1 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: [ ] 0.1 = ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ux T u T x T ff LL λ (1.16) Hay: 0=+ x TT x fL λ (1.17) 0=+ u TT u fL λ (1.18) PGS.TS Nguyn Th Phng Hà 10 http://www.khvt.com Gii (1.17) ta đc λ : 1− −= x T x T fL λ (1.19) và thay vào (1.18) đ có đc (1.14) . Vector n R∈ λ đ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 (1.8) và (1.9) ta kh dx đ đc : dffLdL x T x 1− = (1.20) Vì vy: () λ −== ∂ ∂ − = T x T x du fL f L 1 0 (1.21) 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 thay đi . Nh là mt cách th ba đ tìm đc (1.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 và hàm ch tiêu cht lng đ tìm ra hàm Hamilton . ( ) ( ) ( ) uxfuxLuxH T ,,,, λλ += (1.22) Vi n R∈ λ là tha s Lagrange cha xác đnh . Mun chn x , u , λ đ có đ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 : λ λ dHduHdxHdH TT u T x ++= (1.23) Lu ý rng : ),( uxf H H = ∂ ∂ = λ λ (1.24) Gi s chúng ta chn các giá tr ca u tha mãn : 0 = λ H (1.25) 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 () 0, =uxf . Trong trng hp này hàm Hamilton tng đng vi hàm ch tiêu cht lng: [...]... trên nguyên lý t i u s khai c a M t chi n l c t i u có tính ch t không ph thu c vào nh ng quy t nh tr c ó ( ví d nh nh ng lu t i u khi n ) song các quy t nh còn l i ph i c u thành nên chi n l c t i u có liên quan v i k t qu c a nh ng quy t nh tru c ó Trang 33 PGS.TS Nguy n Th Ph 34 ng Hà Nguyên lý t i u c a Belman : “ B t k m t o n cu i cùng nào c a qu t i u c ng là m t qu ot i u” o Nguyên lý này gi... u i t a n i khi không bi t nguyên lý t i u và i theo chi u thu n T i a ta s so sánh chi phí khi i n b ho c d , và chúng ta quy t nh i n d Ti p t c nh v y ta s i n g ó không còn l a Trang 35 PGS.TS Nguy n Th Ph 36 ng Hà ch n nào khác là i n i qua h Toàn b chi phí cho ph + 4 + 2 = 9 và không ph i là t i u ng án này là 1 + 2 Cu i cùng chúng ta ch ra r ng nguyên lý t i u c a Belman giúp gi m s l ng... Lk xk , uk Jk xk 1 (1.58) 1 v i xk là tr ng thái th i i m k , và xk 1 c cho b i (1.56) Theo nguyên lý Belman thì t n hao t i u t th i i m k s là : J k xk min Lk xk , uk uk Jk 1 xk 1 * và lu t i u khi n t i u u k t i th i i m k là u k làm cho t n hao ti u (1.59) tc c Ph ng trình (1.59) chính là nguyên lý t i u cho h r i r c Vai trò quan tr ng c a nó là có th cho phép chúng ta t i u hóa cùng lúc t i... ng bay c a máy bay M t máy bay bay theo h ng t trái sang ph i nh Hình 1.9 qua các i m a, b, c… t ng tr ng cho các thành ph , và m c nhiên li u c n thi t hoàn t t m i ch ng ng Chúng ta s dùng nguyên lý t i u c a Belman gi i bài toán c c ti u hóa nhiên li u tiêu hao Li t kê các tr ng thái k t 0 n 4 trong quá trình ra quy t nh nh Hình 1.9 ( u m i tên và con s trong khung b c u có th ch a c n quan tâm)... chúng ta có bài toán t i thi u hóa n ng l ng tiêu hao v i tr ng thái cu i c nh , lu t i u khi n và các giá tr tr ng thái tìm ra lu t i u khi n ng v i m c tiêu hao nhiên li u t i thi u , ta s d ng nguyên lý t i u c a Belman , c b t u k N 4 Không có quy t nh nào c yêu c u ây do ó ta gi m k 3 N u x3 c th t f n c th f thì lu t i u khi n t i u là u3 1 và chi phí là 4 i u này hi n b ng cách t (4) phía trên... u là u 2 1 v i chi phí là 4 N u x2 g thì ch có m t s ch n l a duy nh t là u 2 chuy n là 6 1 v i chi phí di B ng cách l n l t gi m k và ti p t c so sánh các ph ng án i u khi n t i u c cho b i nguyên lý t i u , chúng ta có th i n vào các l a ch n còn l i ( u m i tên ) và chi phí t i u c th hi n trong Hình 1.9 D dàng nh n ra r ng chu i i u khi n c l a ch n là chu i t i u Chú ý r ng khi k = 0 , lu... , y1 là 1 i m trên parabol và x 2 , y 2 là 1 i m trên Chúng ta ch n hàm ch tiêu ch t l cách gi a 2 i m này L( x1 , x 2 , y1 , y 2 ) ng là m t n a c a bình ph 1 ( x1 2 x2 ) 2 gi i bài toán này , ta x lý b ng cách f f1 f2 , x x1 x2 1 ( y1 2 y2 ) 2 ng th ng ng kho ng (5) t: , u y1 y2 (6) 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... a qu t i u c ng là m t qu ot i u” o Nguyên lý này gi i h n xem xét trên m t s các ch tiêu t i u Nó ch ra r ng ph ng án t i u ph i c xác nh t tr ng thái cu i i ng c v tr c ó i u ki n áp d ng : nguyên lý t i u là m t ph ng pháp s , ch áp d ng c khi h th ng có phân c p i u khi n và ta bi t tr c s m tl i c xây d ng b ng th c nghi m Ví d n gi n sau s ch ra nh ng v n Bài toán m u ch t c a ph ng pháp này

Ngày đăng: 17/06/2015, 10:26

Từ khóa liên quan

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan