Articulo de referencia

Producto estrella de Redheffer

En matemáticas , el producto estrella de Redheffer es una operación binaria sobre operadores lineales que surge al resolver sistemas acoplados de ecuaciones lineales . Fue intro...

En matemáticas , el producto estrella de Redheffer es una operación binaria sobre operadores lineales que surge al resolver sistemas acoplados de ecuaciones lineales . Fue introducido por Raymond Redheffer en 1959 [ 1 ] y posteriormente se ha adoptado ampliamente en métodos computacionales para matrices de dispersión . Dadas dos matrices de dispersión de diferentes dispersores lineales, el producto estrella de Redheffer produce la matriz de dispersión combinada que se genera cuando algunos o todos los canales de salida de un dispersor se conectan a las entradas de otro dispersor.

Definición

SuponerA,B{\displaystyle A,B}son las matrices de bloquesA=(A11A12A21A22){\displaystyle A={\begin{pmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{pmatrix}}} y B=(B11B12B21B22){\displaystyle B={\begin{pmatrix}B_{11}&B_{12}\\B_{21}&B_{22}\end{pmatrix}}}, cuyos bloquesAij,Bkl{\displaystyle A_{ij},B_{kl}}tienen la misma forma cuando ij=kl{\displaystyle ij=kl}. El producto estrella de Redheffer se define entonces por: [ 1 ]

AB=(B11(IA12B21)1A11B12+B11(IA12B21)1A12B22A21+A22(IB21A12)1B21A11A22(IB21A12)1B22){\displaystyle A\star B={\begin{pmatrix}B_{11}(I-A_{12}B_{21})^{-1}A_{11}&B_{12}+B_{11}(I-A_{12}B_{21})^{-1}A_{12}B_{22}\\A_{21}+A_{22}(I-B_{21}A_{12})^{-1}B_{21}A_{11}&A_{22}(I-B_{21}A_{12})^{-1}B_{22}\end{pmatrix}}} ,

suponiendo que(IA12B21),(IB21A12){\displaystyle (I-A_{12}B_{21}),(I-B_{21}A_{12})}son invertibles , dondeI{\displaystyle I}es una matriz identidad conforme aA12B21{\displaystyle A_{12}B_{21}}oB21A12{\displaystyle B_{21}A_{12}}, respectivamente. Esto se puede reescribir de varias maneras utilizando la llamada identidad de empuje.(IAB)A=A(IBA)A(IBA)1=(IAB)1A{\displaystyle (I-AB)A=A(I-BA)\iff A(I-BA)^{-1}=(I-AB)^{-1}A}.

La definición de Redheffer se extiende más allá de las matrices a los operadores lineales en un espacio de Hilbert.H{\displaystyle {\mathcal {H}}}. [ 2 ] . Por definición,Aij,Bkl{\displaystyle A_{ij},B_{kl}}son endomorfismos lineales deH{\displaystyle {\mathcal {H}}}, haciendoA,B{\displaystyle A,B}endomorfismos lineales deHH{\displaystyle {\mathcal {H}}\oplus {\mathcal {H}}}, dónde{\displaystyle \oplus }es la suma directa . Sin embargo, el producto estrella sigue teniendo sentido siempre que las transformaciones sean compatibles, lo cual es posible cuandoAL(HγHα,HαHγ){\displaystyle A\in {\mathcal {L(H_{\gamma }\oplus H_{\alpha },H_{\alpha }\oplus H_{\gamma })}}} yBL(HαHβ,HβHα){\displaystyle B\in {\mathcal {L(H_{\alpha }\oplus H_{\beta },H_{\beta }\oplus H_{\alpha })}}} de modo queABL(HγHβ,HβHγ){\displaystyle A\star B\in {\mathcal {L(H_{\gamma }\oplus H_{\beta },H_{\beta }\oplus H_{\gamma })}}}.

Propiedades

Existencia

(IA12B21)1{\displaystyle (I-A_{12}B_{21})^{-1}}existe si y solo si (IB21A12)1{\displaystyle (I-B_{21}A_{12})^{-1}}existe. [ 3 ] Por lo tanto, cuando existe cualquiera de ellos, también existe el producto estrella de Redheffer.

Identidad

La identidad de la estrella es la identidad enHH{\displaystyle {\mathcal {H}}\oplus {\mathcal {H}}}, o(I00I){\displaystyle {\begin{pmatrix}I&0\\0&I\end{pmatrix}}}. [ 2 ]

Asociatividad

El producto estrella es asociativo, siempre que se definan todas las matrices relevantes. [ 3 ] Por lo tantoABdo=(AB)do=A(Bdo){\displaystyle A\star B\star C=(A\star B)\star C=A\star (B\star C)}.

Adjunto

Siempre que exista cualquiera de los lados, el adjunto de un producto estrella de Redheffer es(AB)=BA{\displaystyle (A\star B)^{*}=B^{*}\star A^{*}}. [ 2 ]

Inverso

SiB{\displaystyle B}es la matriz inversa izquierda deA{\displaystyle A}de tal manera que BA=I{\displaystyle BA=I},A22{\displaystyle A_{22}}tiene un inverso derecho, y AB{\displaystyle A\star B}existe, entoncesAB=I{\displaystyle A\star B=I}. [ 2 ] De manera similar, siB{\displaystyle B}es la matriz inversa izquierda deA{\displaystyle A}de tal manera queBA=I{\displaystyle BA=I},A11{\displaystyle A_{11}}tiene un inverso derecho, y BA{\displaystyle B\star A}existe, entoncesBA=I{\displaystyle B\star A=I}.

Además, siAB=I{\displaystyle A\star B=I}yA22{\displaystyle A_{22}}tiene una inversa izquierda entoncesBA=I{\displaystyle BA=I}.

La inversa de la estrella es igual a la inversa de la matriz y ambas se pueden calcular con inversión de bloques como [ 2 ].

(A11A12A21A22)1=((A11A12A221A21)1(A21A22A121A11)1(A12A11A211A22)1(A22A21A111A12)1){\displaystyle {\begin{pmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{pmatrix}}^{-1}={\begin{pmatrix}(A_{11}-A_{12}A_{22}^{-1}A_{21})^{-1}&(A_{21}-A_{22}A_{12}^{-1}A_{11})^{-1}\\(A_{12}-A_{11}A_{21}^{-1}A_{22})^{-1}&(A_{22}-A_{21}A_{11}^{-1}A_{12})^{-1}\end{pmatrix}}}.

Derivación a partir de un sistema lineal

El sistema acoplado de ecuaciones, con flechas que etiquetan las entradas y salidas de cada matriz.

El producto estrella surge de la resolución de múltiples sistemas de ecuaciones lineales que comparten variables comunes. A menudo, cada sistema lineal modela el comportamiento de un subsistema en un proceso físico y, al conectar los múltiples subsistemas en un todo, se pueden eliminar las variables compartidas entre ellos para obtener el sistema lineal general. Por ejemplo, sea{incógnitai}i=16{\displaystyle \{x_{i}\}_{i=1}^{6}}ser elementos de un espacio de HilbertH{\displaystyle {\mathcal {H}}}de tal manera que [ 4 ]

(incógnita3incógnita6)=(A11A12A21A22)(incógnita5incógnita4){\displaystyle {\begin{pmatrix}x_{3}\\x_{6}\end{pmatrix}}={\begin{pmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{pmatrix}}{\begin{pmatrix}x_{5}\\x_{4}\end{pmatrix}}}

y

(incógnita1incógnita4)=(B11B12B21B22)(incógnita3incógnita2){\displaystyle {\begin{pmatrix}x_{1}\\x_{4}\end{pmatrix}}={\begin{pmatrix}B_{11}&B_{12}\\B_{21}&B_{22}\end{pmatrix}}{\begin{pmatrix}x_{3}\\x_{2}\end{pmatrix}}}

La "infraestructura" de uno de los sistemas de ecuaciones de Redheffer.

dando lo siguiente4{\displaystyle 4}ecuaciones en6{\displaystyle 6}variables:

incógnita3=A11incógnita5+A12incógnita4incógnita6=A21incógnita5+A22incógnita4incógnita1=B11incógnita3+B12incógnita2incógnita4=B21incógnita3+B22incógnita2{\displaystyle {\begin{aligned}x_{3}&=A_{11}x_{5}+A_{12}x_{4}\\x_{6}&=A_{21}x_{5}+A_{22}x_{4}\\x_{1}&=B_{11}x_{3}+B_{12}x_{2}\\x_{4}&=B_{21}x_{3}+B_{22}x_{2}\end{aligned}}}.

Sustituyendo la primera ecuación en la última, obtenemos:

incógnita4=(IB21A12)1(B21A11incógnita5+B22incógnita2){\displaystyle x_{4}=(I-B_{21}A_{12})^{-1}(B_{21}A_{11}x_{5}+B_{22}x_{2})}.

Sustituyendo la última ecuación en la primera, obtenemos:

incógnita3=(IA12B21)1(A11incógnita5+A12B22incógnita2){\displaystyle x_{3}=(I-A_{12}B_{21})^{-1}(A_{11}x_{5}+A_{12}B_{22}x_{2})}.

Eliminando incógnita3,incógnita4{\displaystyle x_{3},x_{4}} sustituyendo las dos ecuaciones anteriores en las deincógnita1,incógnita6{\displaystyle x_{1},x_{6}}da como resultado que el producto estrella de Redheffer sea la matriz tal que: [ 1 ]

El producto estrella elimina las variables compartidas en este sistema de ecuaciones acopladas.

(incógnita1incógnita6)=(AB)(incógnita5incógnita2){\displaystyle {\begin{pmatrix}x_{1}\\x_{6}\end{pmatrix}}=(A\star B){\begin{pmatrix}x_{5}\\x_{2}\end{pmatrix}}}.

Conexión con matrices de dispersión

La estructura de la matriz de dispersión sigue una convención diferente a la de Redheffer, que consiste en intercambiar y renombrar varias cantidades. La ventaja es que ahora los subíndices de la matriz S identifican los puertos de entrada y salida, así como los índices de los bloques.

Muchos procesos de dispersión adoptan una forma que motiva una convención diferente para la estructura de bloques del sistema lineal de una matriz de dispersión. Típicamente, un dispositivo físico que realiza una transformación lineal en las entradas, como los medios dieléctricos lineales en ondas electromagnéticas o en la dispersión cuántica, puede encapsularse como un sistema que interactúa con el entorno a través de varios puertos, cada uno de los cuales acepta entradas y devuelve salidas. Es convencional utilizar una notación diferente para el espacio de Hilbert,Hi{\displaystyle {\mathcal {H}}_{i}}, cuyo subíndice etiqueta un puerto en el dispositivo. Además, cualquier elemento,doi±Hi{\displaystyle c_{i}^{\pm }\in {\mathcal {H}}_{i}}, tiene un superíndice adicional que indica la dirección del viaje (donde + indica el movimiento del puerto i al i+1 y - indica el movimiento inverso).

La notación equivalente para una transformación de Redheffer, RL(H1H2,H2H1){\displaystyle R\in {\mathcal {L(H_{1}\oplus H_{2},H_{2}\oplus H_{1})}}}, utilizado en la sección anterior es

(do2+do1)=(R11R12R21R22)(do1+do2){\displaystyle {\begin{pmatrix}c_{2}^{+}\\c_{1}^{-}\end{pmatrix}}={\begin{pmatrix}R_{11}&R_{12}\\R_{21}&R_{22}\end{pmatrix}}{\begin{pmatrix}c_{1}^{+}\\c_{2}^{-}\end{pmatrix}}} .

La acción de la matriz S , SL(H1H2,H1H2){\displaystyle S\in {\mathcal {L(H_{1}\oplus H_{2},H_{1}\oplus H_{2})}}}, se define con un cambio adicional en comparación con la definición de Redheffer: [ 5 ]

(do1do2+)=(S11S12S21S22)(do1+do2){\displaystyle {\begin{pmatrix}c_{1}^{-}\\c_{2}^{+}\end{pmatrix}}={\begin{pmatrix}S_{11}&S_{12}\\S_{21}&S_{22}\end{pmatrix}}{\begin{pmatrix}c_{1}^{+}\\c_{2}^{-}\end{pmatrix}}} ,

entonces S=(0II0)R{\displaystyle S={\begin{pmatrix}0&I\\I&0\end{pmatrix}}R} . Tenga en cuenta que para que se definan las matrices identidad fuera de la diagonal, necesitamosH1,H2{\displaystyle {\mathcal {H_{1},H_{2}}}}sea ​​el mismo espacio de Hilbert subyacente. (El subíndice no implica ninguna diferencia, sino que es solo una etiqueta para fines contables).

El producto estrella,S{\displaystyle \star _{S}}, para dos matrices S,A,B{\displaystyle A,B}, viene dado por [ 5 ]

La "infraestructura" de un par acoplado de matrices de dispersión en un producto estrella.

ASB=(A11+A12(IB11A22)1B11A21A12(IB11A22)1B12B21(IA22B11)1A21B22+B21(IA22B11)1A22B12){\displaystyle A\star _{S}B={\begin{pmatrix}A_{11}+A_{12}(I-B_{11}A_{22})^{-1}B_{11}A_{21}&A_{12}(I-B_{11}A_{22})^{-1}B_{12}\\B_{21}(I-A_{22}B_{11})^{-1}A_{21}&B_{22}+B_{21}(I-A_{22}B_{11})^{-1}A_{22}B_{12}\end{pmatrix}}} ,

dóndeAL(H1H2,H1H2){\displaystyle A\in {\mathcal {L(H_{1}\oplus H_{2},H_{1}\oplus H_{2})}}} yBL(H2H3,H2H3){\displaystyle B\in {\mathcal {L(H_{2}\oplus H_{3},H_{2}\oplus H_{3})}}}, entoncesASBL(H1H3,H1H3){\displaystyle A\star _{S}B\in {\mathcal {L(H_{1}\oplus H_{3},H_{1}\oplus H_{3})}}}.

Propiedades

Estos son análogos de las propiedades de{\displaystyle \star }paraS{\displaystyle \star _{S}} La mayoría de ellas se derivan de la correspondencia. J(AB)=(JA)S(JB){\displaystyle J(A\star B)=(JA)\star _{S}(JB)}. J{\displaystyle J}, el operador de intercambio , es también la identidad estrella de la matriz S definida a continuación. Para el resto de esta sección,A,B,do{\displaystyle A,B,C}son matrices S.

Existencia

ASB{\displaystyle A\star _{S}B}existe cuando (IA22B11)1{\displaystyle (I-A_{22}B_{11})^{-1}} o (IB11A22)1{\displaystyle (I-B_{11}A_{22})^{-1}} existir.

Identidad

La identidad estelar de la matriz S,J{\displaystyle J}, es J=(0II0){\displaystyle J={\begin{pmatrix}0&I\\I&0\end{pmatrix}}}. Esto significaJSS=SSJ=S{\displaystyle J\star _{S}S=S\star _{S}J=S}

Asociatividad

Asociatividad deS{\displaystyle \star _{S}}Se deduce de la asociatividad de{\displaystyle \star }y de multiplicación de matrices .

Adjunto

De la correspondencia entre{\displaystyle \star }yS{\displaystyle \star _{S}}y el adjunto de{\displaystyle \star }, tenemos eso (ASB)=J(BSA)J{\displaystyle (A\star _{S}B)^{*}=J(B^{*}\star _{S}A^{*})J}

Inverso

La matrizΣ{\displaystyle \Sigma }que es el producto estrella de la matriz S inverso de S{\displaystyle S}en el sentido de queΣSS=SSΣ=J{\displaystyle \Sigma \star _{S}S=S\star _{S}\Sigma =J} esJS1J{\displaystyle JS^{-1}J}dóndeS1{\displaystyle S^{-1}}es la inversa de matriz ordinaria yJ{\displaystyle J}es como se define anteriormente.

Conexión a matrices de transferencia

Las matrices de transferencia tienen una "estructura de tuberías" diferente a la de las matrices de dispersión. Conectan un puerto con otro, en lugar de conectar las entradas de todos los puertos con las salidas de todos los puertos.

Obsérvese que una matriz de dispersión puede reescribirse como una matriz de transferencia ,T{\displaystyle T}, con acción (do2+do2)=T(do1+do1){\displaystyle {\begin{pmatrix}c_{2}^{+}\\c_{2}^{-}\end{pmatrix}}=T{\begin{pmatrix}c_{1}^{+}\\c_{1}^{-}\end{pmatrix}}}, donde [ 6 ]

T=(T++T+T+T)=(S21S22S121S11S22S121S121S11S121){\displaystyle T={\begin{pmatrix}T_{\scriptscriptstyle ++}&T_{\scriptscriptstyle +-}\\T_{\scriptscriptstyle -+}&T_{\scriptscriptstyle --}\end{pmatrix}}={\begin{pmatrix}S_{21}-S_{22}S_{12}^{-1}S_{11}&S_{22}S_{12}^{-1}\\-S_{12}^{-1}S_{11}&S_{12}^{-1}\end{pmatrix}}} .

Aquí los subíndices relacionan las diferentes direcciones de propagación en cada puerto. Como resultado, el producto estrella de matrices de dispersión

(do3+do1)=(SASB)(do1+do3){\displaystyle {\begin{pmatrix}c_{3}^{+}\\c_{1}^{-}\end{pmatrix}}=(S^{A}\star S^{B}){\begin{pmatrix}c_{1}^{+}\\c_{3}^{-}\end{pmatrix}}} ,

es análogo a la siguiente multiplicación matricial de matrices de transferencia [ 7 ]

(do3+do3)=(TATB)(do1+do1){\displaystyle {\begin{pmatrix}c_{3}^{+}\\c_{3}^{-}\end{pmatrix}}=(T^{A}T^{B}){\begin{pmatrix}c_{1}^{+}\\c_{1}^{-}\end{pmatrix}}} ,

dóndeTAL(H1H1,H2H2){\displaystyle T^{A}\in {\mathcal {L(H_{1}\oplus H_{1},H_{2}\oplus H_{2})}}} yTBL(H2H2,H3H3){\displaystyle T^{B}\in {\mathcal {L(H_{2}\oplus H_{2},H_{3}\oplus H_{3})}}}, entoncesTATBL(H1H1,H3H3){\displaystyle T^{A}T^{B}\in {\mathcal {L(H_{1}\oplus H_{1},H_{3}\oplus H_{3})}}}.

Generalizaciones

Redheffer generalizó el concepto de producto estrella de varias maneras:

Biyecciones arbitrarias

Si existe una biyecciónMETROL{\displaystyle M\leftrightarrow L}dado por L=F(METRO){\displaystyle L=f(M)}Entonces, un producto estrella asociativo puede definirse mediante: [ 7 ]

AB=F1(F(A)F(B)){\displaystyle A\star B=f^{-1}(f(A)f(B))}.

El producto estrella particular definido por Redheffer anteriormente se obtiene de:

F(A)=((IA)+(I+A)J)1((AI)+(A+I)J){\displaystyle f(A)=((I-A)+(I+A)J)^{-1}((A-I)+(A+I)J)}

dóndeJ(incógnita,y)=(incógnita,y){\displaystyle J(x,y)=(-x,y)}.

Producto de 3 estrellas

También se puede definir un producto estrella para matrices de 3x3. [ 8 ]

Aplicaciones a matrices de dispersión

En física , el producto estrella de Redheffer aparece al construir una matriz de dispersión total a partir de dos o más subsistemas. Si el sistemaA{\displaystyle A}tiene una matriz de dispersiónSA{\displaystyle S^{A}}y sistema B{\displaystyle B}tiene matriz de dispersiónSB{\displaystyle S^{B}}, entonces el sistema combinado AB{\displaystyle AB}tiene matriz de dispersiónSAB=SASB{\displaystyle S^{AB}=S^{A}\star S^{B}}. [ 5 ]

teoría de la línea de transmisión

Muchos procesos físicos, incluyendo la transferencia radiativa , la difusión de neutrones, la teoría de circuitos y otros, se describen mediante procesos de dispersión cuya formulación depende de la dimensión del proceso y la representación de los operadores. [ 6 ] Para problemas probabilísticos, la ecuación de dispersión puede aparecer como una ecuación de tipo Kolmogorov .

Electromagnetismo

El producto estrella de Redheffer se puede utilizar para resolver la propagación de campos electromagnéticos en medios estratificados y multicapa. [ 9 ] Cada capa en la estructura tiene su propia matriz de dispersión y la matriz de dispersión total de la estructura se puede describir como el producto estrella entre todas las capas. [ 10 ] Un programa de software libre que simula el electromagnetismo en medios estratificados es el Stanford Stratified Structure Solver .

Interfaces de semiconductores

Los modelos cinéticos de interfaces semiconductoras consecutivas pueden utilizar una formulación de matriz de dispersión para modelar el movimiento de electrones entre los semiconductores. [ 11 ]

Factorización en grafos

En el análisis de operadores de Schrödinger en grafos, la matriz de dispersión de un grafo se puede obtener como un producto estrella generalizado de las matrices de dispersión correspondientes a sus subgrafos. [ 12 ]

Referencias

  1. 1 2 3 Redheffer, Raymond (1959). "Desigualdades para una ecuación de Riccati matricial" . Journal of Mathematics and Mechanics . 8 (3): 349– 367. ISSN 0095-9057 . JSTOR 24900576 .  
  2. 1 2 3 4 5 Redheffer, RM (1960). "Sobre una cierta transformación fraccionaria lineal" . Journal of Mathematics and Physics . 39 ( 1–4 ): 269–286 . doi : 10.1002/sapm1960391269 . ISSN 1467-9590 . 
  3. 1 2 Mistiri, F. (1986-01-01). "El producto estrella y sus propiedades algebraicas" . Journal of the Franklin Institute . 321 (1): 21– 38. Bibcode : 1986FrInJ.321...21M . doi : 10.1016/0016-0032(86)90053-0 . ISSN 0016-0032 . 
  4. Liu, Victor. "Sobre matrices de dispersión y el producto estrella de Redheffer" (PDF) . Consultado el 26 de junio de 2021 .
  5. 1 2 3 Rumpf, Raymond C. (2011). "Formulación mejorada de matrices de dispersión para métodos semianalíticos que es consistente con la convención" . Progress in Electromagnetics Research B. 35 : 241–261 . doi : 10.2528 /PIERB11083107 . ISSN 1937-6472 . 
  6. 1 2 Redheffer, Raymond (1962). "Sobre la relación de la teoría de líneas de transmisión con la dispersión y la transferencia" . Journal of Mathematics and Physics . 41 ( 1– 4): 1– 41. doi : 10.1002/sapm19624111 . ISSN 1467-9590 . 
  7. 1 2 Redheffer, Raymond (1960). "Nota complementaria sobre ecuaciones de Riccati matriciales" . Journal of Mathematics and Mechanics . 9 (5): 745–7f48. ISSN 0095-9057 . JSTOR 24900784 .  
  8. Redheffer, Raymond M. (1960). "El problema de difusión de Mycielski-Paszkowski" . Journal of Mathematics and Mechanics . 9 (4): 607– 621. ISSN 0095-9057 . JSTOR 24900958 .  
  9. Ko, ​​DYK; Sambles, JR (1988-11-01). "Método de matriz de dispersión para la propagación de radiación en medios estratificados: estudios de reflexión total atenuada de cristales líquidos" . JOSA A. 5 ( 11): 1863– 1866. Bibcode : 1988JOSAA...5.1863K . doi : 10.1364/JOSAA.5.001863 . ISSN 1520-8532 . 
  10. Whittaker, DM; Culshaw, IS (1999-07-15). "Tratamiento de matriz de dispersión de estructuras fotónicas multicapa con patrones" . Physical Review B. 60 ( 4): 2610– 2618. Bibcode : 1999PhRvB..60.2610W . doi : 10.1103/PhysRevB.60.2610 .
  11. Gosse, Laurent (2014-01-01). "Productos de Redheffer y aproximación numérica de corrientes en modelos cinéticos de semiconductores unidimensionales" . Multiscale Modeling & Simulation . 12 (4): 1533– 1560. doi : 10.1137/130939584 . ISSN 1540-3459 . 
  12. Kostrykin, V.; Schrader, R. (2001-03-22). "El producto estrella generalizado y la factorización de matrices de dispersión en grafos" . Journal of Mathematical Physics . 42 (4): 1563– 1598. arXiv : math-ph/0008022 . Bibcode : 2001JMP....42.1563K . doi : 10.1063/1.1354641 . ISSN 0022-2488 . S2CID 6791638 .