A Simplex-Like Algorithm for Interval Linear Systems

1 downloads 0 Views 711KB Size Report
https://hal.inria.fr/inria-00073536. Submitted on 24 May 2006. HAL is a ..... an approximation of the inverse of A. Moreover, let A = (v vv. \. Ly. L2 ... L{. \}}}.
A Simplex-Like Algorithm for Interval Linear Systems Olivier Beaumont

To cite this version: Olivier Beaumont. A Simplex-Like Algorithm for Interval Linear Systems. [Research Report] RR-3153, 1997.

HAL Id: inria-00073536 https://hal.inria.fr/inria-00073536 Submitted on 24 May 2006

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L’archive ouverte pluridisciplinaire HAL, est destin´ee au d´epˆot et `a la diffusion de documents scientifiques de niveau recherche, publi´es ou non, ´emanant des ´etablissements d’enseignement et de recherche fran¸cais ou ´etrangers, des laboratoires publics ou priv´es.

INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE

A Simplex-Like Algorithm for Interval Linear Systems. O. Beaumont

N ˚ 3153 Avril 1997 ` THEME 4

ISSN 0249-6399

apport de recherche

 

   

!#"%$&' ()+*,-$. (*+/021% 341 5

6 79;:=?A@CBEDGF 8 H I%KMJ L;KON;PRQ=STLUWVYXGZ[S]\_^.K Z(\a`WZbSTL;SYc[XGZbST\a^ d K,cbe=cfZ KJ L&KMc g \_L&`WV]KhiKMc jlk \amnK Zpo(VYX d ST^ q !X `r`\ k Z d K k KMg I K k g I Ks^taurvxw!uyPzo({ k STV|vM}_}~€PvM‚`rXaƒaKMc

„2…'†M‡ˆG‰WŠa‡a‹yŒ ^yZ I SYc+`rX!`K kxŽ Kc I \ sI \  S‘Z,SYc+`\_c[cfST’WVTK&Zb\ UrcbK“g \_^{_K h”`\aVTe I K d k X•–\ k cb\aVT{=S]^rƒ VTST^WKMX k ST^ZbK k {GX!VcfeicfZbKML&c  S]Z I \_UiZOUcfST^Wƒ2` k Kxg \a^ d S‘Z[S]\_^W^WST^Wƒr—O˜™Kš k cfZOc I \ ›I \  Zb\ d K k ST{aK Ž • k \aL Xa^AK^g VT\_cbU k K2\!•œpOžfŸ  p¡£¢MŸ ¤[¡¦¥ Ž X§`\aVTe I K d k \_^ pI STg I g \_^_ZX!ST^rc+Z I K.g \a^={_K h I UrV]V¨\!•(Z I Kcb\aVTUiZ[S]\_^ cbK Zx— H I KM^ Ž X©cfSTL&`WVTK h=ª«VTST¬aK;L;KZ I \ d K^rXa’WVTKMcOUrcOZb\šr^ d X^WK  \aUiZ[K k ST^rg VTUrcbS]\_^%—­ \ k KM\G{aK kMŽ Z I K g\a^rcfZ k X!ST^Z[c¨\a’WZ[X!ST^WK d L“X®e’K+UrcbK d Zb\“g \_L&`WUiZbK,Xa^ST^W^WK k ST^rg VTUrcbS]\_^\a•¨œpOžfŸ  p¡£¢MŸ ¤[¡¦¥ — ¯ °±C²[³´ ˆGµ'†a‹ ST^ZbK k {GX!VV]ST^WKMX k cfeicfZbKL“c Ž cbS]L&`WVTK hXaV]ƒ_\ k S‘Z I L

¶Y·;¹b¸ º«»x¼ ¹¾¸ '½ ¿Tº£À«Á  ÃĦŮÆÈÇÈÉ Ç«ÊMÇ'Ë®ÌlÍ%ÌnΣϮÌnЫΣÏxÌlÌnžĦÅxÑTÒÓÐÈÔ(Õ[Ç«É×ÖÊ®ÌÌØÇÙÚMÆÈÇfÌfÛ ÔpÌnÆ%Ü'Ý£ÌfÞ Õ[ÇÈÒÓÉ×ЫÌnÆ%ߍÌnÔ(ÕÉ×ÝYà ÒÓá!ÌfÕÊ®ÔpÒÓÅ!â'É×ЫÉ×Æ«Õxã ÑTÐ

Unit´e de recherche INRIA Rennes IRISA, Campus universitaire de Beaulieu, 35042 RENNES Cedex (France) T´el´ephone : (33) 02 99 84 71 00 – T´el´ecopie : (33) 02 99 84 71 71

$ *, Ó 4  C0™  

© ™.  ( 1' |  (1 1 021  3  (1#M$

( *, | (1#*() M$;' ()+*,M (1 5 ° † ° ‹  X!^rc“gK Z2X k Z[STgV]K Ž ^W\_Urc“L&\a^Z k \a^c;g\aL&L&K^Z“S]V3KMcfZ&`\_c[cbS]’WVTK d UiZbSTVTSTcbK k d Kxc&`\_V]eª I%KJ d k KMcg \_^{_K hiKMc`\_U k;k Kx cf\_U d k K“V]Kxc&cfeicfZ KMJ L;Kxc,VTS]^ KM XaS k KMc d  ST^ZbK k {GX!VTVTKMc;cbXa^rcUWZbSTV]SYcfK k d K` kKx g \a^Wª d S‘Z[S]\_^W^WKL&KM^_Zx—3\aUrcL&\a^Z k \a^rcOZ[\aUiZ d  Xa’\ k d g\aL&L&K^Z,Z k \aUr{aK kxŽ¨X.J `X k ZbS k d  UW^AcfU k ª£K^rcbKL’WVTK d K œ OžfŸ  p¡£¢MŸ ¤[¡–¥ Ž UW^8`\_V]e I%KJ d k K UWS g \a^Z[S]KM^_Z+V  KM^rcbKL’rV]Kcb\aVTUiZbST\a^— (\aUc+UiZbSTVTSTcb\a^rc+XaV]\ k cOUW^AXaV‘ª ƒ_\ k S‘Z I L&K d K2Zne=`K.cfSTL&`WVTK hiK`\aU k d K  ZbK k L&ST^WK k UW^-^W\aUW{_KMXaUAcfU k ª£K^rcbKL’WV]K d K2VYX™cf\_V]UiZ[S]\_^%—  K `rV]Urc Ž V]Kxcg \_^Z k X!ST^_Z[KMc€\_’iZbKM^UrKMc,`KUr{aK^Z K  Z k K&UiZbSTVTSTc K Kxc+`\aU k d K Z[K k L;ST^WK k UW^Acb\aUrcfª«KM^rcbKL’rV]K d K œpOžfŸ  p¡£¢MŸ ¤[¡–¥Ó— !s´ ‡x† ² #Š " ° ‹ cbeicnZ KJ L&Kxc VTST^ KM XaS k KMc d  ST^_Z[K k {GX!VTV]Kxc Ž XaV]ƒ_\ k S‘Z I L&K d U cbS]L&`WVTK hiK

      !"#$ %&(' #) *'+,(-. !/-%0 u

"%$;'|  ¨ ;$ o VT\!Z€\!•  \ k ¬ I XacO’KK^ d \a^rKST^y\ k d K k Zb\.cf\_V]{_KST^_Z[K k {GX!V'V]ST^WKMX k cfeicfZbKL“c;Ÿ  (¡3254 Ÿ ¤[¡(Ÿ 6®¡£— q \ I ^ Xa^ d87 k KST^W\G{=STg I Ÿ]v:9®¡ Ž Ÿ Na¡ I X_c` k \G{aK d Z I X!ZlZ I KOgMX!VYg UWVYXGZ[S]\_^“\!•œpOžfŸ  p¡£¢MŸ ¤[¡¦¥ Ž Z I KOcfL“X!VTVTKMcfZ ’\®h;Z I XGZ g\a^Z[XaS]^c3X!VTV'Z I K&cb\aVTUiZbST\a^rc¾\!•3Ÿ  p¡;2=¾^§Z I K;\!Z I K kOI X!^ d Ž cfKM{aK k XaV XaV]ƒ_\ k S‘Z I L“c¨\a’iZX!ST^ƒ_\\ d k KxcfUrV‘Zc Ž Kxcf`KMgSTXaV]VTe pI KM^Z I K d SYX!L&K Z[K k \a• Ÿ  p¡SYcpcbL&XaV]V Ÿ ‚®¡ Ž Ÿ 6!¡«— Œ ^Z I SYcl`X!`K kxŽ_ K3` k \a`\_cbK(X,^rK  X!VTƒa\ k S]Z I L pI SYg I SYc’rX_cfK d \a^“VTST^WKMX k ` k \aƒ k XaL&L;ST^Wƒ— Œ Zlg\a^rcbSTcfZ[c \a•lX!^ S]ZbK k X!ZbST{aKc[g I KL&K pI STg I g \a^cfS d K k c(X_c3ST^W`WUiZ€X!^§KM^rg VT\_cbU k K,\!•|Z I Kcb\aVTUiZ[S]\_^ \!•(œ OžbŸ  (¡«¢xŸ ¤¡¦¥ Xa^ d k KZbU k ^rcpX.ž–UrcbUrXaV]VTe2’K ZbZbK k ¥¨K^g VT\_cbU k K_— H I STc;X!VTƒa\ k S]Z I L g \_^={aK k ƒaKxc€Zb\  X k d X”cbUW`K k cfKZ\a•pZ I K.g \_^{_K h I UWV]V \a•pZ I K©UW^WS]ZbK d cb\aVTUiZbST\a^cbK Z OžbŸ  (¡«¢xŸ ¤¡¦¥ — Œ ^©Z I STc `X!`K kxŽW K€UrcbKOZ I K€•–\_V]VT\  ST^Wƒ&^W\!ZXGZbST\a^%c ? OžbŸ  (¡«¢xŸ ¤¡¦@¥ 4BA:2 D¢ CiB   E8Ÿ  (¡«¢ CW¤ E8Ÿ ¤¡«¢(  2F4 H¤ G I žfŸ  p¡£¢MŸ ¤[¡–J¥ 4LKMiž£OžfŸ  p¡£¢MŸ ¤[¡–¥f¥Ó¢ pI K k >K KM d KM^W\!Z[KMc Z I K,g \_^={aK h I UWVTV ˜™K,X!VYcb\ d KM^W\!Z[K+’=e§œ OžbŸ  (¡«¢xŸ ¤[¡–¥ Z I K+S]^ZbK k {GXaV I UrV]V%\a• Z I K,cb\aVTUiZbST\a^%— H I K3ST^_Z[K k {GX!VV]ST^WKMX k cbeicnZ[KL“c  K¾g\a^rcbS d K k S]^“Z I SYc `X!`K k X k K d K š^WK d ’=e;Z I K3•–\_V]VT\  S]^rƒ,^W\!ZXGZ[S]\_^rc%? Ÿ  pN¡ 4›Ÿ P   O5QE ,¢[S   RTQ% p¡ Ÿ ¤[¡)4›Ÿ ¤ ¢U¤b¡V4›Ÿ ¤WOXQ ¤®¢[¤YRSQl¤[¡£¢Ql¤Z4 U@¤ O ¤

1

[

 `_”$ * ™;Ó0.  |;$ |.*O p$;*,$.1 |. p;$&)

  |.  $   1' | ;$ 1 p Œ Z SYcl¬=^W\  ^“Z I XGZ¨Z I KOcb\aVTUiZbST\a^2\!•EZ I K¾` k \a’rV]KML Ÿ  (¡32F4sŸ ¤¡SYc¨S]^ƒaK^rK k X!Vr^r\!Zpg \a^={_K h%Ÿ 6x¡£—|o3c • X k Xac  K€X k K3\a^rV]e“ST^ZbK k KMcfZbK d ST^”œpOžfŸ  p¡£¢MŸ ¤[¡–¥ ŽW K¾L“X®e“UrcbK I žfŸ  p¡£¢MŸ ¤[¡–¥¨Xac¨S]^Z[K k L;K d SYXGZbK€cbK ZM@— =¾KZfZbVTS%X!^ d jlk XaƒaK k Ÿ×~M¡` k \G{aK d Z I K€•–\aVTV]\  ST^Wƒ;Z I K\ k KLa? bdc °=´ ˆ ° fe gž CWh  E8Ÿ  p¡£i¢ Cr¤ E8Ÿ ¤[¡«¢¾D  2j4 ¤ ¥ kml D  2nO ¤ lpoPQE/  l 2qlHRTQ ¤ H I STclK hi` k KMc[cfST\a^ d \=KMcl^W\!Z¨V]KxX d“d S k KMgÓZ[V]eZb\X,`\aVTe I K d k \a^ Ž ’KxgXaUrcfK3\!•CZ I KOX!’rcb\aVTUiZ[K({GX!VTUWK ŽpI SYg I Ur^ d K k VTS]^WKxc¨Z I K€• X_gÓZpZ I XGZpZ I K,cf\_V]UWZbST\a^.cfKZ SYcpUrcbUrXaV]VTe^W\aZ3X;g\a^={aKh2cfKZŸ 6G¡«— r K ZUc%Xac[cfUrL;KZ I XGZ  K¬=^W\  X!^€KM^rg VT\_cbU k K|\!•Eœ OžbŸ  (¡«¢xŸ ¤[¡–¥ Ž®pI SYg I L“X®e3’Kl\a’WZ[X!ST^WK d Ž •–\ k ST^rcnZX!^rgK Ž ’=e“Z I K,X!VTƒa\ k S]Z I L ` k \a`\_cbK d ’e q UWL&` Ÿ ‚!¡«—|˜™K,ƒaST{aK€X;š k cnZ k KMcbUWV]Z pI STg I ` k \G{=S d KxcpX  X®e“Z[\&ƒaK Z k S d \a• Z I K,X!’rcb\aVTUiZ[KO{GXaV]UWKxc— s °  AX‰ e "t œ OžbŸ  (¡«¢xŸ ¤[¡–v¥ u Ÿ 2 ¢ 2rU ¡

\

]  ^

Í Í ÅYw:xHy{zx

0 ' #  N y

x x __i



__ x

S]ƒ_U k K;v?Q=S‘Z[UrXGZ[S]\_^ p I K ^.’\!Z I 2 X!^ d

'# J 

i

2U 

X k K¾`\_cbS‘Z[S]{_K

4 l 2 U 2U DlOO 2 l 2 l '+# 4 2 U l 2 2U Zl+OO 22 l 2 U l  O *+ # %#  { 2  i '+&+    ! 2  %# +. !->   2 ™ E Ÿ 2 ¢ 2 U ¡£¢Hl 2 lo  2 R  0



ˆ ´%´ ‹

2U o 9 Ž Z I K^ 4 O+va¢  D4L9 Ž Xa^ d %2 R 4 OZ2F4 l 2ql‘—  Q=KMg\a^ d gMXacbK ?2  9 Ž Z I KM^ 4›v_¢  D4L9 Ž X!^ d %2 R  4 2F4 l 2ql‘—  H I S k d gX_cfK? 2 U %2 o 9 Ž Z I K^ 4 !!# "'"%$& !! "" ¢  4 &! ("*!& ")! ! " " —    & ( ! o , o 2 U Ž Z I K^  ,R  O l , l4 ! " & ! " "/. ž 2 U O0,¥  9W— + S]•Y9 Z I Urc Ž l ,)lpo  ,R  — o , o 9 Ž Z I K^ ,R  DO l , l4 ! "*( &! " ! " . ž1, O@# * A@%, B 4h=A žt2 ¢%:W¥Ó¢2 8E Ÿ 2 ¢ 42 U ¡£A ¢ 9žt2 Co ¢):W:F¥ 4 .*+ ?,8 .@A>-6B,8CBC,DA2EJ#) * @* %    $ !#(- +# ,