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

130213023 王登高--列主元

作者:130213023王登高   发布时间:2015-06-23 16:50:02   浏览次数:207

function [L,U,q]=guassliexuan(A)

[m,n]=size(A);

p=[];

tol=1e-16;

if m==n

    for k=1:n-1

        [r,c]=max(abs(A(k:n,k)));

        q(k)=c+k-1;

        A([k,q(k)],:)=A([q(k),k],:);

        if abs(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);

        else

            'A(k,k)=0!'

            break;

        end
    end
else

    'm~=n';

end

L=eye(n)+tril(A,-1);

U=triu(A);
 

A=[1,2,3;4,5,6;7,8,9];
[L,U]=guassliexuan(A)

L =

    1.0000         0         0
    0.1429    1.0000         0
    0.5714    0.5000    1.0000


U =

    7.0000    8.0000    9.0000
         0    0.8571    1.7143
         0         0    0.0000

 







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

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

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