解上三角形方程组

本文介绍了一个MATLAB函数,用于求解上三角形线性方程组。该函数首先判断系数矩阵A和增广矩阵B的秩是否相等,并据此决定方程组是否有解。若秩相等且等于方程数,则方程组有唯一解;若不相等,则方程组无解;若秩相等但小于方程数,则方程组有无穷多个解。通过回代法逐步计算出解。
function [RA,RB,n,X]=shangsan(A,b)
%解上三角形方程组
B=[A b];n=length(b);RA=rank(A);RB=rank(B);zhica=RB-RA;
if(zhica>0)
    disp('RA~=RB,无解');
end


if(RA==RB)
    if(RA==n)
        disp('有唯一解');
        X=zeros(n,1);X(n)=b(n)/A(n,n);
        for k=n-1:-1:1
            X(k)=(b(k)-sum(A(k,k+1:n)*X(k+1:n)))/A(k,k);
        end
    else
        disp('无穷解');
    end
end
end
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值