MHB What are the values of α for any initial value of x(0) in a linear system?

AI Thread Summary
The discussion focuses on determining the values of α for which a linear system converges, specifically in the context of numerical stability. It is established that for convergence, the condition ||αA|| < 1 must be satisfied, where A is a matrix and ||·|| denotes the matrix norm. This implies that |α| must be less than 1/||A||. The conversation also references the solution of a first-order difference equation and how it relates to the convergence of the iteration. Overall, the key point is that the choice of α directly impacts the stability and convergence of the system.
natalia
Messages
6
Reaction score
0
Considering the following

View attachment 2409
find the values for α for any initial value of x(0).

Any help will be useful, thank you!
 

Attachments

  • cap2_CalculNumeric.png
    cap2_CalculNumeric.png
    797 bytes · Views: 83
Mathematics news on Phys.org
natalia said:
Considering the following

View attachment 2409
find the values for α for any initial value of x(0).

Any help will be useful, thank you!

Welcome to MHB natalia! :)

Your problem needs a little more context.

I'm going to venture out and guess that this is about numerical stability.
The question might then be which values of $\alpha$ will reduce a worst case error in $x^{(0)}$.
In that case we can say that $\alpha$ times the norm of the matrix must have an absolute value that is smaller than 1.

Which context can you provide for this problem?
 
I need to find the values of $\alpha$ for which the iteration converges.
 
Last edited:
natalia said:
Considering the following

View attachment 2409
find the values for α for any initial value of x(0).

Any help will be useful, thank you!

Wellcome on MHB natalia!...

In...

http://mathhelpboards.com/discrete-mathematics-set-theory-logic-15/difference-equation-tutorial-draft-part-i-426-post2494.html

... it has been explained that the first order difference equation...$\displaystyle x_{n+1} = a\ x_{n} + b\ (1)$ ... has solution...

$\displaystyle x_{n} = x_{0}\ a^{n} + b\ \frac{1 - a^{n}}{1 - a}\ (2)$

In Your case I suppose that $\overrightarrow {x}_{n} $ and $\displaystyle \overrightarrow {b}$ are vectors of dimesion 2, $\alpha$ is a scalar and $\displaystyle A = \left | \begin{matrix} 2 & 1 \\ 1 & 2 \end{matrix} \right |$ a 2 x 2 Matrix, so that is...

$\displaystyle \overrightarrow {x}_{n+1} = \overrightarrow {b} + \alpha\ A\ \overrightarrow {x}_{n}\ (3)$

... and the solution of (3) is...

$\displaystyle \overrightarrow {x}_{n} = \overrightarrow {x}_{0}\ \alpha^{n}\ A^{n} + \overrightarrow {b}\ (I- \alpha^{n}\ A^{n})\ (I- \alpha\ A)^{-1}\ (4)$

... where I is the 2 x 2 identity matrix...

Kind regards

$\chi$ $\sigma$
 
Last edited:
natalia said:
I need to find the values of $\alpha$ for which the iteration converges.

Aha!

Actually, that turns out the same: you need $\alpha$ such that
$$||\alpha A || < 1$$
where $|| \cdot ||$ represents the matrix norm.
This is the same as:
$$|\alpha| \cdot || A || < 1$$
$$|\alpha| < \frac 1 {|| A ||}$$To explain, suppose the sequence converges (or should converge) to $y$.

Let $\Delta y^{(k)} = x^{(k)} - y$, which represents how far you are from the limit.
And let's call your matrix $A$ for ease of notation.

Then:
\begin{aligned}
x^{(k+1)} &= b + \alpha A x^{(k)} \\

&= b + \alpha A (y + \Delta y^{(k)}) \\

&= (b + \alpha A y) + \alpha A \Delta y^{(k)} \\

&= y + \alpha A \Delta y^{(k)}
\end{aligned}

That means that:
$$\Delta y^{(k+1)} = \alpha A \Delta y^{(k)}$$

In a worst case scenario this will only converge if
$$||\alpha A || < 1$$
where $|| \cdot ||$ represents the matrix norm.
 
Thank you very much, I like Serena :), your explanations are really clear.
 
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...
Fermat's Last Theorem has long been one of the most famous mathematical problems, and is now one of the most famous theorems. It simply states that the equation $$ a^n+b^n=c^n $$ has no solutions with positive integers if ##n>2.## It was named after Pierre de Fermat (1607-1665). The problem itself stems from the book Arithmetica by Diophantus of Alexandria. It gained popularity because Fermat noted in his copy "Cubum autem in duos cubos, aut quadratoquadratum in duos quadratoquadratos, et...
Thread 'Imaginary Pythagorus'
I posted this in the Lame Math thread, but it's got me thinking. Is there any validity to this? Or is it really just a mathematical trick? Naively, I see that i2 + plus 12 does equal zero2. But does this have a meaning? I know one can treat the imaginary number line as just another axis like the reals, but does that mean this does represent a triangle in the complex plane with a hypotenuse of length zero? Ibix offered a rendering of the diagram using what I assume is matrix* notation...
Back
Top