Résolution d'un système linéaire avec des arguments matriciels

Nous connaissons tous les nombreuses méthodes de calcul pour résoudre le système linéaire standard

$$Ax=b.$$ Cependant, je suis curieux de savoir s'il existe des méthodes de calcul "standard" pour résoudre un système linéaire plus général (de dimension finie) de la forme

$$LA=B,$$ où, disons,$A$ est unematrice$m_1\times n_1$ ,$B$ est unematrice$m_2\times n_2$ et$L$ est un opérateur linéaire prenant$m_1\times n_1$ matrices en$m_2\times n_2$ matrices, ce quin'implique pas la vectorisation des matrices, c'est-à-dire tout convertir en la forme standard$Ax=b$ .

La raison pour laquelle je demande est que je dois résoudre l'équation suivante pour $u$ :

$$(R^*R+\lambda I)u=f$$ où$R$ est la transformée de Radon 2d,$R^*$ son adjoint, et$u$ et$f$ sonttous les deux destableaux (images). Il est possible de vectoriser cette équation, mais c'est pénible, surtout si on passe en 3D.

Plus généralement, qu'en est-il des tableaux $nD$ ? Par exemple, résoudre $LA=B$ où $A$ et $B$ sont des tableaux 3D (je devrai également le faire avec la transformation de Radon à un moment donné).

Merci d'avance et n'hésitez pas à m'envoyer vers un autre StackExchange si vous en ressentez le besoin.

Oui, vous avez raison et cela fonctionnera bien lorsque vous passerez à la 3D. La partie la plus simple, vraiment, est le produit intérieur --- il suffit de faire un produit scalaire standard sur l'équivalent, R n déroulé$\mathbb{R}^n$$\langle y, x\rangle\triangleq \mathop{\textrm{Re}}(y^Hx)$

Une chose à laquelle vous devez faire attention lorsque vous implémentez CG (ou des approches itératives similaires) avec des opérateurs linéaires généraux est d'implémenter correctement l' adjoint de votre opérateur linéaire. Autrement dit, les gens obtiennent souvent correctement, mais commettent une erreur en mettant en œuvre .$y=F(x)$$z=F^*(y)$

Je recommande d'implémenter un test simple qui tire parti de l'identité suivante: pour tout et conforme , Donc, ce que vous faites est de générer des valeurs aléatoires de et , de les exécuter à travers vos opérations directes et adjointes, respectivement, et de calculer les deux produits internes ci-dessus. Assurez-vous qu'ils correspondent à une précision raisonnable et répétez plusieurs fois.$x$$y$

$$\langle y, F(x) \rangle = \langle F^*(y), x \rangle.$$ $x$$y$

EDIT: que faites-vous si votre opérateur linéaire est supposé être symétrique? Eh bien, vous devez également vérifier cette symétrie. Utilisez donc le même test, notant simplement que --- appliquez la même opération à et . Bien sûr, l'OP a à la fois un opérateur asymétrique et un opérateur symétrique à gérer ...$F=F^*$$x$$y$

$R^*R+\lambda I$$r_k^Tr_k$$p_k^TAp_k$

$$\langle A,B\rangle=\sum_{i,j}A_{ij}B_{ij}$$

Je soupçonne que cela se passera très bien lorsque je passerai aux tableaux 3D, bien que je n'aie pas encore vu le produit interne Frobenius défini sur les tableaux 3D (je travaillerai en supposant que je peux à nouveau simplement résumer le produit ponctuel).

Je serais toujours intéressé par des méthodes plus générales si quelqu'un en connaissait!