Web Analytics
Privacy Policy Cookie Policy Terms and Conditions Identità di Woodbury - Wikipedia

Identità di Woodbury

Da Wikipedia, l'enciclopedia libera.

In matematica, particolarmente in algebra lineare, la identità matriciale di Woodbury per matrici di aspetto n × n è data dalla seguente formula:

\left( \textrm{A}+\textrm{UC}\textrm{V}\right) ^{-1}= \textrm{A}^{-1}-\textrm{A}^{-1}\textrm{U}\left( \textrm{C}^{-1}+\textrm{VA}^{-1}\textrm{U}\right) ^{-1}\textrm{VA}^{-1}.

In essa A e UCV sono matrici di aspetto n × n, mentre C è una matrice quadrata che può avere aspetto diverso r × r; conseguentemente U ha aspetto n × r e V aspetto r × n.

Indice

[modifica] Applicazioni

Questa identità risulta utile in elaborazioni numeriche nelle quali \,A^{-1}\, è già costruita e risulta necessario costruire \,(A+UCV)^{-1}\,. Essa risulta vantaggiosa quando A è numericamente preponderante rispetto a UCV che può considerarsi una sua piccola perturbazione. Se inoltre l'estensione di C è molto più piccola di quella di A (r molto minore di n), possiamo ottenere il secondo membro invertendo solo due matrici di estensione ridotta.

Questo procedimento si applica, ad es., nel filtro di Kalman e in altri metodi di stima dei minimi quadrati, per rimpiazzare la soluzione parametrica, che richiede l'inversione di una matrice quadrata di un ordine dato dalle dimensioni di un vettore di stato, con una soluzione basata su equazioni di condizione. Nel caso del filtro di Kalman questa matrice ha le dimensioni del vettore delle osservazioni, fino a ridursi a 1 nel caso si elabori una sola nuova osservazione per turno. L'utilizzo della identità di Woodbury spesso accelera significativamente calcoli del filtro che si devono effettuare in tempo reale.

[modifica] Dimostrazione

Consideriamo la seguente equazione:

\begin{bmatrix} \textrm{A} & \textrm{U}\\ \textrm{V} & -\textrm{C}^{-1} \end{bmatrix} \begin{bmatrix} \textrm{D}_{11} & \textrm{D}_{12}\\ \textrm{D}_{21} & \textrm{D}_{22} \end{bmatrix} = \begin{bmatrix} \textrm{I} & 0\\ 0 & \textrm{I} \end{bmatrix}

Trascriviamo questa mediante quattro equazioni:

\begin{matrix} \textrm{AD}_{11}+\textrm{UD}_{21} & = & \textrm{I} \\ \textrm{AD}_{12}+\textrm{UD}_{22} & = & 0 \\ \textrm{VD}_{11}-\textrm{C}^{-1}\textrm{D}_{21} & = & 0 \\ \textrm{VD}_{12}-\textrm{C}^{-1}\textrm{D}_{22} & = & \textrm{I} \end{matrix}

Di queste, solo la prima e la terza ci occorrono.

Addizioniamo la terza alla prima dopo aver moltiplicato per UC:

\left( \textrm{A}+\textrm{UC}\textrm{V}\right) \textrm{D}_{11}=\textrm{I}\, \Rightarrow \, \textrm{D}_{11}=\left( \textrm{A}+\textrm{UC}\textrm{V}\right) ^{-1}

A questo punto si sottrae la prima dalla terza dopo aver moltiplicato per VA - 1:

\left( -\textrm{C}^{-1}-\textrm{VA}^{-1}\textrm{U}\right) \textrm{D}_{21}=-\textrm{VA}^{-1}\, \Rightarrow \, \textrm{D}_{21}=\left( \textrm{C}^{-1}+\textrm{VA}^{-1}\textrm{U}\right) ^{-1}\textrm{VA}^{-1}

Torniamo a sostituire nella prima equazione:

\textrm{AD}_{11}+\textrm{U}\left( \textrm{C}^{-1}+\textrm{VA}^{-1}\textrm{U}\right) ^{-1}\textrm{VA}^{-1}=\textrm{I}\, \Rightarrow \textrm{D}_{11}=\textrm{A}^{-1}-\textrm{A}^{-1}\textrm{U}\left( \textrm{C}^{-1}+\textrm{VA}^{-1}\textrm{U}\right) ^{-1}\textrm{VA}^{-1}

Ora abbiamo due differenti espressioni per la sottomatrice D11 che dovrebbero essere identiche. Otteniamo così:

\left( \textrm{A}+\textrm{UC}\textrm{V}\right) ^{-1}=\textrm{A}^{-1}-\textrm{A}^{-1}\textrm{U}\left( \textrm{C}^{-1}+\textrm{VA}^{-1}\textrm{U}\right) ^{-1}\textrm{VA}^{-1}.

Con questa si completa la dimostrazione.

[modifica] Vedi anche:

[modifica] Rinvio esterno

Altre lingue
THIS WEB:

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - be - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - closed_zh_tw - co - cr - cs - csb - cu - cv - cy - da - de - diq - dv - dz - ee - el - eml - en - eo - es - et - eu - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gd - gl - glk - gn - got - gu - gv - ha - haw - he - hi - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mg - mh - mi - mk - ml - mn - mo - mr - ms - mt - mus - my - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - rm - rmy - rn - ro - roa_rup - roa_tara - ru - ru_sib - rw - sa - sc - scn - sco - sd - se - searchcom - sg - sh - si - simple - sk - sl - sm - sn - so - sq - sr - ss - st - su - sv - sw - ta - te - test - tet - tg - th - ti - tk - tl - tlh - tn - to - tokipona - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu

Static Wikipedia 2008 (no images)

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - bcl - be - be_x_old - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - co - cr - crh - cs - csb - cu - cv - cy - da - de - diq - dsb - dv - dz - ee - el - eml - en - eo - es - et - eu - ext - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gan - gd - gl - glk - gn - got - gu - gv - ha - hak - haw - he - hi - hif - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kaa - kab - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mdf - mg - mh - mi - mk - ml - mn - mo - mr - mt - mus - my - myv - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - quality - rm - rmy - rn - ro - roa_rup - roa_tara - ru - rw - sa - sah - sc - scn - sco - sd - se - sg - sh - si - simple - sk - sl - sm - sn - so - sr - srn - ss - st - stq - su - sv - sw - szl - ta - te - tet - tg - th - ti - tk - tl - tlh - tn - to - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu -

Static Wikipedia 2007:

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - be - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - closed_zh_tw - co - cr - cs - csb - cu - cv - cy - da - de - diq - dv - dz - ee - el - eml - en - eo - es - et - eu - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gd - gl - glk - gn - got - gu - gv - ha - haw - he - hi - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mg - mh - mi - mk - ml - mn - mo - mr - ms - mt - mus - my - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - rm - rmy - rn - ro - roa_rup - roa_tara - ru - ru_sib - rw - sa - sc - scn - sco - sd - se - searchcom - sg - sh - si - simple - sk - sl - sm - sn - so - sq - sr - ss - st - su - sv - sw - ta - te - test - tet - tg - th - ti - tk - tl - tlh - tn - to - tokipona - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu

Static Wikipedia 2006:

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - be - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - closed_zh_tw - co - cr - cs - csb - cu - cv - cy - da - de - diq - dv - dz - ee - el - eml - en - eo - es - et - eu - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gd - gl - glk - gn - got - gu - gv - ha - haw - he - hi - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mg - mh - mi - mk - ml - mn - mo - mr - ms - mt - mus - my - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - rm - rmy - rn - ro - roa_rup - roa_tara - ru - ru_sib - rw - sa - sc - scn - sco - sd - se - searchcom - sg - sh - si - simple - sk - sl - sm - sn - so - sq - sr - ss - st - su - sv - sw - ta - te - test - tet - tg - th - ti - tk - tl - tlh - tn - to - tokipona - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu