电力系统节点优化

是包含了半动态优化和静态优化，并且有着详细的中文注释，代码刚通过课设，保证正确。学生党，希望和各位一起交流

function [lr] = LR_my（ y ）
%对节点导纳矩阵进行三角分解
%   节点导纳矩阵为y
%   下三角矩阵为L
%   上三角单位矩阵为R
[numrow]=size（y）;
l=zeros（numrow）;
r=zeros（numrow）;
for i=1:num
    for j=1:i
        a=0;
        for k=1:j-1;
            a=a+l（ik）*r（kj）;
        end
            l（ij）=y（ij）-a;
    end
    for j=i+1:num
        a=0;
        for k=1:i-1;
            a=a+l（ik）*r（kj）;
        end
        r（ij）=（y（ij）-a）/l（ii）;
    end
    for m=1:num
        r（mm）=1;
    end  
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      10097  2018-01-16 21:59  node\Brand.xlsx

     文件        525  2018-01-07 12:15  node\LR_my.m

     文件       2785  2018-01-17 10:53  node\node_admit.m

     文件        491  2018-01-07 12:25  node\patent_U.m

     文件       3303  2018-01-17 23:39  node\semi.m

     文件       1977  2018-01-17 23:41  node\static.m

     目录          0  2018-01-26 16:59  node

----------- ---------  ---------- -----  ----

                19178                    7