matlab平面刚架矩阵分析法-1

函数  halfBandWidth  将对称矩阵化为半带宽矩阵

函数  halfBandChol  对半带宽矩阵进行乔列斯基分解

function h=halfBandWidth(a)
%将对称矩阵化为半带宽矩阵
%时间：2021/12/15/ 姓名：尹林博
n=length(a);
assert(issymmetric(a),...
    ['Arguement input error:  '...
    'Matrix must besymmetric.'])
V=1;
h=zeros(n,V);
for i=1:n
     for j=n:-1:i
         h_=0;
         if a(i,j)~=0
             h_=a(i,i:j);
             break
         end
     end
     if V<length(h_)
         V=length(h_);
     end
     h(i,V)=0;
     h(i,1:length(h_))=h_;
end

function s=halfBandChol(s)
%对半带宽矩阵进行乔列斯基分解
%时间：2021/12/16 姓名：尹林博
[n,v]=size(s);
for i=n:-1:1
    for j=v:-1:2
        if j+i-1<=n
           for k=j+1:v
               if k+i-1<=n
                  s(i,j)=s(i,j)-s(i,k)*s(k+i-1,1)*s(j+i-1,k-j+1);
               end
           end
           s(i,j)=