ويرد ادناه ج / ج + + وبرنامج شفره المصدر ان يحسب ybus فان وzbus معين من الشبكه الكهرباءيه من أي حجم. مدخلات النظام هو مكتوب على "rx.txt" وفيما يلى شكل العمود ؛ 'من الحافله' ، 'لحافله '،' ص '، و،' x '. 'من الحافله' و 'لحافله' يجب ان يكون موجبا على التوالي ابتداء من الصفر. حافله تحمل الصفر هو ان يكون البديل او تراخ في الحافله. معقدة ملف الرأس واستخدم لمساعدة ماتريكس العمليات. وحسب zbus كان مجرد عملية قلب ybus.

//-----------------------------
# pragma hdrstop
//----------------------------
# pragma argsused
# تشمل <stdio.h>
# تشمل <conio.h>
# تشمل "matrix.h"
# تشمل <complex.h>
باستخدام الامراض التي تنتقل بالاتصال الجنسي : : معقدة ؛
الرياضيات باستخدام المجرد ؛
typedef معقدة معقدة ؛
typedef مصفوفه مصفوفه ؛
الرئيسية الفراغ (الفراغ)
(
/ / تعلن المتغيرات
المركز الوطني للاستشعار عن دولي ، الاول ، ي ، matsize = 1 ؛
tmp3 مزدوجة ، tmp4 ؛
Tmp6 معقدة ؛
مصفوفه y ، ybus ، zbus ؛
* في الملف ؛
ي. Setsize (matsize ، matsize) ؛ / / مجموعة مؤقتة حجم المصفوفه
/ / للحصول على بيانات فرع
فى = fopen ( "rx.txt" ، و "ص") ؛ / / فتح الملف
اذا كانت الاجابه ب (== في لاغيه)
(printf ( "\ nrx.txt غير موجود") ؛
getch () ؛
)
المركز الوطني للاستشعار عن = ي أنا = = 0 ؛
ل) ؛ ؛)
(اذا (fscanf (في "د ٪" ، & الأول) == eof)
(كسر ؛)
اذا كان (ط> matsize)
= (matsize الاول ؛)
fscanf (في "د ٪" ، & ي) ؛
اذا كانت الاجابه ب (ي> matsize)
(matsize = ي ؛)
ي. Setsize (matsize +1 ، matsize +1) ؛
fscanf (في "المجله ٪" ، & tmp3) ؛
fscanf (في "المجله ٪" ، & tmp4) ؛
tmp6 = معقدة (tmp3 ، tmp4) ؛
ذ (ط ، ي) = 1.0/tmp6 ؛
المركز الوطني للاستشعار عن + + ؛
)
fclose (في) ؛
/ / نهاية فرع للوصول الى البيانات
/ / خلق ybus
Ybus.setsize (matsize +1 ، matsize +1) ؛
/ / مداخل قطر
ل(ط = 0 ؛ الاول <= matsize ؛ الاول + +)
(ل(ي = 0 ؛ ي <= matsize ؛ ي + +)
(Ybus (أنا ، أنا) = ybus (أنا ، أنا) + y (ط ، ي) ؛)
)
ل(ي = 0 ؛ ي <= matsize ؛ ي + +)
(ل(ط = 0 ؛ الاول <= matsize ؛ الاول + +)
(Ybus (ي ، ي) = ybus (ي ، ي) + y (ط ، ي) ؛)
)
/ / نهاية
/ / من مداخل قطر
المركز الوطني للاستشعار عن = 0 ؛
ل(ي = المركز الوطني للاستشعار ؛ ي <= matsize ؛ ي + +)
(ل(ط = 0 ؛ الاول <= matsize ؛ الاول + +)
(اذا (ى! = ط)
(Ybus (ط ، ي) = ybus (ي ، ط) =- 1،0 * y (ط ، ي) ؛)
)
المركز الوطني للاستشعار عن + + ؛
)
/ / نهاية
/ / y ايجاد نهاية حافله
/ / 0 حافله لا تدرج في المصفوفه ybus
/ / لانها تشكل اشارة حافله
ل(ي = 1 ؛ ي <= matsize ؛ ي + +)
(ل(ط = 1 ؛ الاول <= matsize ؛ الاول + +)
(Ybus ط - 1 ، 1 - ي) = ybus (ط ، ي) ؛
)
)
Ybus.setsize (matsize ، matsize) ؛
Zbus.setsize (matsize ، matsize) ؛
printf ( "\ n \ nybus مصفوفه \ n") ؛
محكمة <<ybus <<endl ؛
Zbus =! Ybus ؛ / / للحصول على عكس من ybus
printf ( "\ n \ nzbus مصفوفه \ n") ؛
محكمة <<zbus <<endl ؛
getch () ؛
)