plot command is not work in Gauss Seidel Method

Question

suketu vaidya 2020-12-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/678183-plot-command-is-not-work-in-gauss-seidel-method

回答： suketu vaidya 2020-12-4

采纳的回答： John D'Errico

A = [

1 0 -1 0 0 0 0

-1 1 1 0 0 0 0

0 -1 2 -1 0 0 0

0 -1 0 2 0 1 0

0 -1 -1 0 2 1 1

0 0 1 -1 0 2 -1

0 0 0 -1 1 -1 2

];

b = [0 0 0 0 0 20 0]';

x=[1 1 1 1 1 1 1];

n=size(x,1);

normVal=Inf;

tol=1e-5; itr=0;

%%Gauss Seidel Method

while normVal>tol

x_old=x;

for i=1:n

sigma=0;

for j=1:i-1

sigma=sigma+A(i,j)*x(j);

end

for j=i+1:n

sigma=sigma+A(i,j)*x_old(j);

end

x(i)=(1/A(i,i))*(b(i)-sigma);

end

itr=itr+1;

normVal=norm(x_old-x);

end

%%

fprintf('system solution : \n%f\n%f\n%f\n%f in %d iterations',x,itr);

hold on

plot(itr,x, '+')

xlim([0 10])

ylim([0 2])

hold off

end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

John D'Errico 2020-12-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/678183-plot-command-is-not-work-in-gauss-seidel-method#answer_565163

在 MATLAB Online 中打开

When you tell us that something does not work, you need ot explain why you THINK it does not work. In fact, plot worked perfectly well there. At least, it tiried to work as well as it could.

itr
itr =
     2
     
x
x =
     0     1     1     1     1     1     1
     
plot(itr,x, '+')

So itr is the number 2. x is a vector. What did you expect it to do?

2 个评论
显示无隐藏无

suketu vaidya 2020-12-4

hello ,

i is not shosws full vactor x . i means x1 = 0 x2 = 1 x3=1 and so on . i want a graph that shows x1,x2,x3 in x axist and on y axis value of x1,x2,x3 and so on

suketu vaidya 2020-12-4