de-CH
utf-8
math math-format
Matrix-Vektor-Iteration
mvit-01-01
multiple
163840
randRange(1,4) randRange(1,5) randRange(0,1) pow(-1,SN)*N/(N+K) randRangeNonZero(-8,-8) BZ/(N+K) fractionReduce(BZ,N+K) L1-2*BN fractionReduce(pow(-1,SN)*N-2*BZ,N+K) randRangeNonZero(-8,8) 2*AN+4*BN-2*D fractionReduce(2*(pow(-1,SN)*N-2*BZ)+4*(BZ)-2*D*(N+K),N+K) randRangeNonZero(-8,8)

Gegeben sei eine Folge von Vektoren \left(v_n \right)_n mit v_{n+1} =A \cdot v_n und A = \begin{pmatrix} A & B\\ C & D \end{pmatrix}.

Bestimmen Sie den Eintrag {\color{red}y}, sodass die die Folgen mit dem Startvektor {\color{orange}v_0} = \begin{pmatrix} X \\ {\color{red}y} \end{pmatrix} gegen den Nullvektor konvergiert.

a \color{red} y = 2 * X

Die Matrix hat einen Eigenwert {\color{blue}\lambda} mit |{\color{blue}\lambda}| < 1 .

Dies ist {\color{blue}\lambda}= fractionReduce(pow(-1,SN)*N,N+K).

Wir suchen {\color{red}y} mit v_1 = A \cdot {\color{orange}v_0} = {\color{blue}\lambda} {\color{orange}v_0},

mit anderen Worten der Startvektor {\color{orange}v_0} ist ein Eigenvektor zum Eigenwert {\color{blue}\lambda} = fractionReduce(pow(-1,SN)*N,N+K).

Dann ist weiter v_n = A^n \cdot {\color{orange}v_0} = {\color{blue}\lambda}^n {\color{orange}v_0} und mit |{\color{blue}\lambda}| < 1 konvergiert die Zahlenfolge \left({\color{blue}\lambda}^n \right)_n gegen Null.

Wir rechnen \begin{pmatrix} A & B\\ C & D \end{pmatrix} \cdot \begin{pmatrix} X \\ {\color{red}y} \end{pmatrix} =\begin{pmatrix} negParens(A) \cdot negParens(X) + negParens(B) \cdot {\color{red}y} \\ negParens(C) \cdot negParens(X) + negParens(D) \cdot {\color{red}y} \end{pmatrix} = \begin{pmatrix} negParens(fractionReduce(X*(pow(-1,SN)*N-2*BZ),N+K)) + negParens(B) \cdot {\color{red}y} \\ negParens(fractionReduce(X*(2*(pow(-1,SN)*N-2*BZ)+4*(BZ)-2*D*(N+K)),N+K)) + negParens(D) \cdot {\color{red}y} \end{pmatrix}

und suchen damit {\color{red}y} mit \begin{pmatrix} negParens(fractionReduce(X*(pow(-1,SN)*N-2*BZ),N+K)) + negParens(B) \cdot {\color{red}y} \\ negParens(fractionReduce(X*(2*(pow(-1,SN)*N-2*BZ)+4*(BZ)-2*D*(N+K)),N+K)) + negParens(D) \cdot {\color{red}y} \end{pmatrix} = {\color{blue}negParens(fractionReduce(pow(-1,SN)*N,N+K))} \cdot \begin{pmatrix} X \\ {\color{red}y} \end{pmatrix}.

Mit der ersten Zeile negParens(fractionReduce(X*(pow(-1,SN)*N-2*BZ),N+K)) + negParens(B) \cdot {\color{red}y} = fractionReduce(X*pow(-1,SN)*N,N+K) können wir dann {\color{red}y} = 2*X bestimmen.

Damit gilt für n \to \infty für die Koordinaten der Vektoren v_n = \begin{pmatrix} \left(fractionReduce(pow(-1,SN)*N,N+K)\right)^n \cdot negParens(X) \\ \left(fractionReduce(pow(-1,SN)*N,N+K)\right)^n \cdot negParens(2*X)\end{pmatrix} \to \begin{pmatrix} 0 \\ 0 \end{pmatrix}.