the loop is not converging it is giving me an infinite value(it is increasing)

1 次查看（过去 30 天）
Mahmoud Chawki
Mahmoud Chawki 2023-5-20
移动：Matt J 2023-5-20
clear variables
close all
clc
A = [1 1 2; 1 2 4; 1 2 5];
x = [1; 2; 3];
Error = 0.00001;
Diff= 1;
while Diff>Error
x1=A*x;
Diff=norm(x1-x);
x=x1;
end

请先登录，再回答此问题。

回答（1 个）

Torsten
Torsten 2023-5-20
移动：Matt J 2023-5-20
Ran in:
Why do you think your fixed-point iteration should converge ?
If the absolute values of the eigenvalues of A were all < 1, it would work. But it is not the case:
A = [1 1 2; 1 2 4; 1 2 5];
eig(A)
ans = 3×1
7.1842 0.5728 0.2430
  5 个评论
显示 3更早的评论
Torsten
Torsten 2023-5-20
移动：Matt J 2023-5-20
What do you mean by "make it work" ? It works as you have programmed it and gives you the correct result for the matrix A and the initial vector x: Divergence.
Mahmoud Chawki
Mahmoud Chawki 2023-5-20
移动：Matt J 2023-5-20
i found the answer, i need to normalize the new vector before using it again. thank you

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by