matlab由导纳求阻抗,关于利用矩阵稀疏技术求解节点导纳矩阵的MATLAB编程

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

求大神帮忙看看这个程序，解释一下！

%形成节点导纳矩阵

%考虑参考节点

%考虑互感支路

function [d_g,d_b,u_g,u_b]=calyy(itoi,noels_hgln,noelsnum_hgln,bs_noels,bs_cp,gb_hgln,gb_noels,gb_cp,lv_bs,lv_hgln,lv_noels,lv_cp,type_bs,iy,jy,remv_noels)

%lv_bs

d_g=zeros(1,lv_bs);

d_b=zeros(1,lv_bs);

u_g=zeros(1,length(jy));

u_b=zeros(1,length(jy));

type_noels=zeros(1,lv_noels);

for i=1:lv_hgln

%%支路自阻抗

lnnum=noelsnum_hgln(i);

lnnum;

gb_noels;

for ii1=1:lnnum

iln=noels_hgln(i,ii1);

[d_g,d_b,u_g,u_b]=getlnselfgb(iln,bs_noels,type_bs,itoi,gb_noels,iy,jy,d_g,d_b,u_g,u_b,lv_bs);

type_noels(iln)=1;

end

for ii1=1:lnnum

iln=noels_hgln(i,ii1);

for ii2=ii1+1:lnnum

zln=noels_hgln(i,ii2);

b=gb_hgln(i,ii1,ii2);

[d_g,d_b,u_g,u_b]=getlnhugb(iln,zln,bs_noels,type_bs,itoi,0,b,iy,jy,d