ex6.m

数值分析, 同济大学教材<<现代数值数学和计算>>数值分析课程所有算法的matlab代码,所有例程均经过测试,私家珍藏

A=[ 4  -1 0  0  0 0  0  0  0  0; 
   -1  4 -1  0  0 0  0  0  0  0; 
    0  -1 4 -1  0 0  0  0  0  0; 
    0  0  -1 4 -1 0  0  0  0  0; 
    0  0  0  -1 4 -1 0  0  0  0; 
    0  0  0  0  -1 4 -1 0  0  0; 
    0  0  0  0  0  -1 4 -1 0  0; 
    0  0  0  0  0  0 -1 4  -1 0; 
    0  0  0  0  0  0  0 -1 4  -1; 
    0  0  0  0  0  0  0  0 -1 4;]; 
  
b = [7 5 -13 2 6 -12 14 -4 5 -5]'; 

[n,m] = size(A);

yupsilong = [1e-3 1e-4 1e-5 1e-6 1e-7 1e-8 1e-9 1e-10]; 

x_real = [2 1 -3 0 1 -2 3 0 1 -1];

for i=1:8 
    x0=ones(n,1)*10^(i-1); 
    x0'; 
    result_Jacobi=(Jacobi(A,b,x0,yupsilong(3))) 

    result_GS=(Gauss_Seidol(A,b,x0,yupsilong(3))) 

    result_SOR=(SOR(A,b,x0,yupsilong(3),1.1))

end 

x0=ones(n,1); 
for i=1:length(yupsilong) 
    (yupsilong(i))' 
  result_Jacobi=(Jacobi(A,b,x0,yupsilong(i)))'; 
  result_GS=(Gauss_Seidol(A,b,x0,yupsilong(i)))'; 
  result_SOR=(SOR(A,b,x0,yupsilong(i),1.1))'; 
end 

omiga=0.1:0.1:2 
for i=1:length(omiga) 
    omiga(i) 
  result_SOR=(SOR(A,b,x0,yupsilong(3),omiga(i)))'; 
end