Need help with Jacobi iteration

4 次查看（过去 30 天）
Osame Amayo
Osame Amayo 2017-11-29
评论： Osame Amayo 2017-12-6
Currently attempting to write a jacobi iteration that will give me the same results as A\b. I am fairly new to matlab I am getting close but a few numbers are off. Is there anything that I am doing wrong? Also my professor said I don't need D but without D I can not get D_inv what am I suppose to do?
d=diag(A);
D=diag(d);
D_inv=inv(D);
E=A-D;
x=rand(n,1);
T=-D_inv*E;
C=D_inv*b;
j-1:N
x=T*x+C;
disp(x);
  4 个评论
显示 2更早的评论
Elizabeth Reese
Elizabeth Reese 2017-12-6
You do not need to form the full D or D_inv. They are both just diagonal matrices. Think about how multiplying by a diagonal matrix changes the original matrix or vector. If you leave everything in terms of d, then you can get matrix-vector multiplication instead of matrix-matrix multiplication.
Osame Amayo
Osame Amayo 2017-12-6
I'm sorry I am a little confused leaving everything in terms of d. If I leave everything in terms of d, how would I complete the function?

请先登录，再进行评论。

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

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by