下面是C / C + +程序和源代码，计算ybus和zbus某一特定电网络的任何大小。输入系统是写上“ rx.txt ”具有下列栏格式; '巴士' ， '巴士' ， ' ， ，的' X ' 。 '从巴士'和'巴士'必须有一个连续正整数开始，从零。 巴士零假设是要疲弱或摆动的巴士。 一个复杂的头文件被用来协助矩阵运算。 zbus计算简单的逆变该ybus 。

//-----------------------------
的＃ pragma hdrstop
//----------------------------
的＃ pragma argsused
＃包括 < stdio.h >
＃包括< conio.h >
＃包括“ matrix.h ”
＃包括< complex.h >
使用性病： ：复杂;
使用名字空间的数学;
typedef复杂，复杂的;
typedef矩阵矩阵;
无效的主体（无效）
（
/ /申报变数
CNT的诠释，我律政司， matsize = 1 ;
双tmp3 ， tmp4 ;
复杂的tmp6 ;
矩阵 y ， ybus ， zbus ;
*在文件;
耀setsize （ matsize ， matsize ） ; / /设置一个临时矩阵大小
/ /得到分行的数据
在= fopen （ “ rx.txt ” ， “住宅” ） ; / /开启档案
如果（在==空）
（ printf （ “ \ nrx.txt没有发现” ） ;
getch （ ） ;
）
CNT的= j = = 0 ;
（ ; ; ）
（如果（ fscanf （在， “ ％ D ”类， ＆ I的） ==的EOF ）
（打破; ）
如果（我> matsize ）
（ matsize = ; ）
fscanf （在“第％ d ” ， ＆ j ）条;
如果（十> matsize ）
（ matsize = j ; ）
耀setsize （ matsize一， matsize 1 ） ;
fscanf （在， “ ％的LF ” ， ＆ tmp3 ） ;
fscanf （在， “ ％的LF ” ， ＆ tmp4 ） ;
tmp6 =复杂（ tmp3 ， tmp4 ） ;
Y （下，我j ）条= 1.0/tmp6 ;
CNT的+ + ;
）
fclose （在） ;
/ /年底越来越分行的数据
/ /创建ybus
ybus.setsize （ matsize一， matsize 1 ） ;
/ /对角线作品
（一= 0 ;我< = matsize ，我+ + ）
（为（十= 0 ; j < = matsize ; J + +的）
（ ybus （我，我） = ybus （我，我） + Y （下，我j ）条; ）
）
为（十= 0 ; j < = matsize ; J + +的）
（ （一= 0 ;我< = matsize ，我+ + ）
（ ybus （十，十） = ybus （十，十） + Y （下，我j ）条; ）
）
/ /完
/ /小康对角线作品
CNT的= 0 ;
为（十= CNT的; j < = matsize ; J + +的）
（ （一= 0 ;我< = matsize ，我+ + ）
（如果（十！ =一）
（ ybus （一，十） = ybus （十，我） =- 1.0 * y （上，我j ）条; ）
）
CNT的+ + ;
）
/ /完
/ /年底创造y巴士
/ /巴士0是不包括在ybus矩阵
/ /因为它是参考巴士
为（十= 1 ; j < = matsize ; J + +的）
（ （一= 1 ，我< = matsize ，我+ + ）
（ ybus （我- 1 ，的J - 1 ） = ybus （一， j ）条;
）
）
ybus.setsize （ matsize ， matsize ） ;
zbus.setsize （ matsize ， matsize ） ;
printf （ “ \ n \ nybus矩阵\ n ” ） ;
法院< < ybus < < endl ;
zbus = ！ ybus ; / /得到逆的ybus
printf （ “ \ n \ nzbus矩阵\ n ” ） ;
法院< < zbus < < endl ;
getch （ ） ;
）