lup5run.m

分块五对角线性方程组的求解

%求解块五对角方程测试数据
function [x,err]=luP5Run(n)
a=[14,-4,0;-2,14,-4;0,-2,14];
b=[4,-2,0;2,4,-2;0,2,4];
c=[-5,1,0;-1,-5,1;0,-1,-5];
d=[-2,0.1,0;1,-2,0.1;0,1,-2];
e=[-2,4,0;5,-2,4;0,5,-2];
f=[3,4,5;6,7,8;9,10,11];

% a=[-3,1,0;1,-3,1;0,1,-3];
% b=zeros(3);
% c=zeros(3);
% d=zeros(3);
% e=zeros(3);
% f=[-3,0,0;0,-3,0;0,0,-3];


o=zeros(3);

A=cell(n,1);
B=cell(n,1);
C=cell(n,1);
D=cell(n,1);
E=cell(n,1);
F=cell(n,1);

for i=1:n
    A{i}=a;
    B{i}=b;
    C{i}=c;
    D{i}=d;
    E{i}=e;
    F{i}=f;
end

x=luP5(A,B,C,D,E,F,3,n);
if n==2
    mA=[a,c;b,a];
    mX=cell2mat(x);
    mF=cell2mat(F);
    err=norm(mA*mX-mF,inf);
end
if n==3
    mA=[a,c,e;b,a,c;d,b,a];
    mX=cell2mat(x);
    mF=cell2mat(F);
    mA*mX-mF
    err=norm(mA*mX-mF,inf);
end
if n==4
    mA=[a,c,e,o;b,a,c,e;d,b,a,c;o,d,b,a];
    mX=cell2mat(x);
    mF=cell2mat(F);
    mA*mX-mF
    err=norm(mA*mX-mF,inf);
end
if n==5
    mA=[a,c,e,o,o;b,a,c,e,o;d,b,a,c,e;o,d,b,a,c;o,o,d,b,a];
    mX=cell2mat(x);
    mF=cell2mat(F);
    mA*mX-mF
    err=norm(mA*mX-mF,inf);
end