现在时间是:
当前位置:首 页 >> 矩阵计算>> 教学区>> 文章列表

130213011 雷航 实验二

作者:雷航   发布时间:2015-06-23 16:23:38   浏览次数:143

function x=Gauss(A,b)

n=size(A,1)

    tol=1e-16;

    x=zeros(n,1);

    for k=1:n-1

        if A(k,k)>tol

            A(k+1:n,k)=A(k+1:n,k)/A(k,k);

            A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);

            b(k+1:n)=b(k+1:n)-A(k+1:n,k)*b(k);

        else

            stop

        end

    end

    for j=n:-1:2

        b(j)=b(j)/A(j,j);

        b(1:j-1)=b(1:j-1)-b(j)*A(1:j-1,j);

    end

    x(k)=b(k)/A(k,k);

 

a=ones(1,83);

b=ones(1,84);

c=ones(1,83);

D=diag(a,1)+6*diag(b,0)+8*diag(c,-1);

B=[7;15*ones(82,1);14];

E=inv(D);

y=E*B;

x=gauss(D,B)

err=y-x

运行结果:

x =

  1.0e+008 *

         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
         0
   -1.0066
         0







上一篇:没有了    下一篇:没有了

Copyright ©2019    计算数学达人 All Right Reserved.

技术支持:自助建站 | 领地网站建设 |短信接口 |燕窝 版权所有 © 2005-2019 lingw.net.粤ICP备16125321号 -5