Conosciamo tutti i numerosi metodi computazionali per risolvere il sistema lineare standard
A x = b .
Tuttavia, sono curioso di sapere se esistono metodi computazionali "standard" per risolvere un sistema lineare più generale (dimensione finita) della forma
A m 1 × n 1 B m 2 × n 2 L m 1 × n 1 m 2 × n 2 A x = b
L A = B ,
dove, diciamo, è una matrice , è una e è un operatore lineare che porta matrici alle matrici , che
non comporta la vettorializzazione le matrici , ovvero la conversione di tutto nella forma standard .
UNm1× n1Bm2× n2Lm1× n1m2× n2A x = b
Il motivo che mi chiedo è che ho bisogno di risolvere la seguente equazione per :u
R R ∗ u f
( R*R + λ I) u = f
dove è la trasformata Radon 2d, suo adiacente, e sia che sono matrici 2d (immagini). È possibile vettorializzare questa equazione, ma è una seccatura, specialmente se si passa al 3D.
RR*uf
Più in generale, che dire delle matrici ? Ad esempio, risolvendo dove e sono array 3D (dovrò farlo anche con la trasformazione del Radon ad un certo punto).L A = B A Bn DL A = BUNB
Grazie in anticipo e sentiti libero di spedirmi in un altro StackExchange se ne senti il bisogno.