HELP with creating matrix for Jacobi and Gauss-Seidel method problem!!

Question

equinox 2016-11-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/311913-help-with-creating-matrix-for-jacobi-and-gauss-seidel-method-problem

评论： Dang Khoa 2018-8-13

How can i create this matrix in matlab?

aij= [2i when j=i and i=1,2,...,80

0.5i when {j=i+2 and i=1,2,...,78 AND j=i-2 and i=3,4,...,80

0.25i when {j=i+4 and i=1,2,...,76 AND j=i-4 and i=5,6,...,80

0 otherwise]

and those of b are bi = π, for each i = 1,2,...,80.

I will need to use Jacobi and Gauss-seidel method to solve the linear system Ax = b to within 10−5 in the l∞ norm with the given matrix entries above. I am just stuck with how to create the matrix. I have spent so long reading the book and other sites but still have no idea where to even begin with the entries of this matrix so any help is appreciated!

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Dang Khoa 2018-8-13

After that, how can you create the matrix ? can you give me a code

请先登录，再进行评论。

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

Answer 1

Daniel kiracofe 2016-11-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/311913-help-with-creating-matrix-for-jacobi-and-gauss-seidel-method-problem#answer_243116

在 MATLAB Online 中打开

the brute force approach would be something like this.

for i = 1:80
  for j = 1:80
    if (i==j)
       a(i,j) = 2*i;
    elseif ( other conditions)
       a(i,j) = whatever
    etc.
  end 
end

There are quicker ways to do it, but this will be easiest to understand. for 80x80 matrix, speed will not be an issue.

2 个评论
显示无隐藏无

equinox 2016-11-15

Thank you for the help! i was able to construct the matrix using this approach. However now i am a little confused on using the Jacobi and Gauss-seidel methods to solve the linear system Ax=b. Any tips or suggestions?