matlab导纳矩阵形成函数

function [ Y ] = lfybus（ zdata ）
%UNtitleD2 此处显示有关此函数的摘要
%   此处显示详细说明
nl=zdata（:1）;nr=zdata（:2）;R=zdata（:3）;X=zdata（:4）;B=zdata（:5）;A=zdata（:6）;
nbr=length（zdata（:1））;nbus=max（max（nl）max（nr））;
Z=R+j*X;
y=ones（nbr1）./Z;
Y=zeros（nbusnbus）;
for k=1:nbr;