matlab 改进的平方根法 LDL 解线性方程组 数值分析

%平方根法LDL
%author LijiaYi（foddcus） FaFu university 2022.6
%email：foddcus@163.com


%input输入矩阵，包括系数和参数

%output 输出答案，对应X的各个值
%xnum,ynum 输入矩阵的尺寸
%Lneed 消元中，主式需要乘的变量
%keepnum 保留小数位数
%errorF（各个多项式的误差）

clear all%启用模块化后需要注释该句
input=[2,-1,1,4;
       -1,2,3,5;
       1,3,1,6];%输入矩阵
oringinData=input;
keepnum=-10;%设置保留n位小数,记得加负号表示小数位
[ynum,xnum]=size(input);%获取输入矩阵大小
d(1,1)=input(1,1);
for j=2:ynum
    u(j,1)=input(j,1);
    l(j,1)=u(j,1)/d(1,1);
end
for k=2:ynum
    for i=k+1:ynum
    d(k,1)=input(k,k)-sum(u(k,:).*l(k,:));
    u(i,k)=input(j,k)-sum(u(i,:).*l(k,:));
    l(i,k)=u(i,k)/d(k,1);
    end
end
d(ynum,1)=input(ynum,ynum)-sum(u(ynum,:).*l(ynum,:));

b(1,1)=input(1,xnum);
for i=2:ynum
    b(i,1)=input(i,xnum)-sum(l(i,:).*b');
end
output=zeros(ynum,1);
output(ynum,1)=b(ynum,1)/d(ynum,1)
for i=2:ynum
    output(ynum+1-i,1)=b(ynum+1-i,1)/d(ynum+1-i,1)-sum(l(:,ynum+1-i).*output)
end


%验算
for i=1:ynum
errorF(i,1)=sum(output.*input(i,1:end-1)')-input(i,xnum);
end