Todos estamos familiarizados con los muchos métodos computacionales para resolver el sistema lineal estándar.
A x = b .
Sin embargo, tengo curiosidad si hay algún método computacional "estándar" para resolver un sistema lineal más general (dimensión finita) de la forma
L A = B ,
donde, por ejemplo, es una matriz , es una matriz , y es un operador lineal que toma matrices matrices , lo
que no implica vectorización las matrices , es decir, convertir todo a la forma estándar .
m 1 × n 1 B m 2 × n 2 L m 1 × n 1 m 2 × n 2UNAmetro1× n1simetro2× n2Lmetro1× n1metro2× n2A x = b
La razón por la que pregunto es que necesito resolver la siguiente ecuación para :tu
( R∗R + λ I) u = f
donde es la transformación de radón 2d, es contigua, y tanto como son matrices (imágenes) 2d. Es posible vectorizar esta ecuación, pero es un dolor, especialmente si vamos a 3D.
R ∗ u fRR∗tuF
Más en general, ¿qué pasa con las matrices ? Por ejemplo, resolver donde y son matrices 3D (necesitaré hacer esto con Radon transform en algún punto también).L A = B A Bn DL A = BUNAsi
Gracias de antemano, y siéntase libre de enviarme a otro StackExchange si lo necesita.