خب این کاری که انجام شده از به اصطلاح چندین برنامه براش نوشته شده، یکی runme هست که برای clear all و clc میاد و تمام این به اصطلاح variable هایی که در workspace هست رو پاک می کنه و close all force میاد و برای شما بلافاصله تمام پنجره ها رو می بنده،
یک basemva صد mva در نظر گرفتیم و یک accuracy یک صدمی رو براش در نظر گرفتیم که تا حد یک صدم mva می بایست که itretion ها را ادامه بده. busdata را برحسب دادها های مقاله در نظر گرفتیم
و linedata و cost و mvlimits رو براش در نظر گرفتیم. اینجا میاد و برنامه به اصطلاح lfybus رو فراخوانی می کنه
که اومدیم یک itretion اولیه براش 10 تا در نظر گرفتیم. (1-)j=sqrt(-1); i=sqrt هست، اینها رو براشون در نظر گرفتیم که بیاییم و قسمت های imagation یا به اصطلاح قسمت های موهومی را باهاش بسازیم. linedata رو میخونه و قسمت اول linedata رو در nl1 ، nr میشه دومین ستون R ، X ، Bc و a هر کدوم میان و linedata رو برابر خودشون قرار میدن.
nbr برابر طول linedata1 میشه و تعداد bus ها برابر ماکسیمم nl , nr قرار میگیره که بیشترین مقدارش که هر کدوم باشه برابر اون قرار خواهد گرفت و z هم که مشخصه،
z= R+j*x و y رو برای اولین گام میاییم برابر یک قرار میدیم و از اونجا میاییم و fbus رو شروع می کنیم تعداد رو برحسب 1 تا nbr میذاریم
و چنانچه a<0 شد a=1 قرار میدیم و Ybus رو برابر zero همه رو صفر در نظر می گیریم برای اولین بار تا بیاییم Ybus ها رو با استفاده از فرمول های فصل شش کتاب به اصطلاح هادی سعادت بیاییم و محاسبه کنیم
که می بینید همین گونه رفتار کرده و حالا دوباره اینجا میاییم و Ybus ها رو برحسب تعداد nbus میاییم و محاسبه می کنیم با استفاده از فرمول و بعد از اون میاییم و شکل هاش رو ترسیم می کنیم و plot می کنیم، بعد از این قسمت می رسیم به loadflow که برحسب Newton-Raphson هست و اینجا نوشته شده این loadflow برحسب Newton-Raphson .
ns , ng و یک سری پارامترهایی را اول برابر صفر قرار می دهیم و nbus رو برابر length busdata قرار میدیم و حالا میاییم طبق اون loadflow گفته شده توی درس به اصطلاح درس بررسی سیستم های قدرت هادی سعادت ، اینها نوشته شده و در آخرش plot میشه و خلاصه اینکه بعد از اینکه ما بیاییم و سیستم رو run کنیم.
تمام این برنامه میاد برای ما اجرا میشه و در نهایت می بینیم که نمودارهایی که به اصطلاح داده شده در این مقاله برای ما بدست میاد که حالا در اینجا هم به وضوح دیده می شوند. 
توضیح کامل برنامه run me
همه را پاک کن؛
CLC؛
basemva = 100؛ دقت = 0.01؛ maxiter = 10؛
٪ باس ولتاژ زاویه — وزنه —- ——- ژنراتور —– استاتیک مگاوار
٪ خیر مجله کد. درجه MW مگاوار MW مگاوار Qmin Qmax به + QC / -ql
busdata = [1 1 1.06 0.0 0 0 0 0 10 50 0
2 2 1.045 0.0 20 10 40 30 10 50 0
3 2 1.03 0.0 20 15 30 10 10 40 0
4 0 1.00 0.0 50 30 0 0 0 0 0
5 0 1.00 0.0 60 40 0 0 0 0 0
6 0 1.08 0.0 0 0 0 0 0 0 0]؛
کد٪ خط
٪ باس باس R X 1/2 B = 1 برای خطوط
٪ NL NR p.u. p.u. p.u. > 1 یا <1 TR. شیر در باس NL
linedata = [1 2 0.02 0.06 0.030 1
1 3 .08 .24 .025 1
2 3 0.06 0.18 0.020 1
2 4 0.06 0.18 0.020 1
2 5 0.04 0.12 0.015 1
3 4 0.01 0.03 0.010 1
4 5 0.08 0.24 0.025 1
1 6 0.08 0.24 0.025 1
3 6 0.08 0.20 0.035 1.1]؛
هزینه = [200 7.0 0.008
180 6.3 0.009
140 6.8 0.007]؛
mwlimits = [10 85
10 80
10 70].
lfybus٪ شکل ماتریس ادمیتانس باس
loadflownewton٪ راه حل جریان به روش نیوتن رافسون
٪ busout چاپ راه حل جریان برق بر روی صفحه نمایش
از دست دادن٪ به دست آوردن ضرایب فرمول از دست دادن
٪ محاسبه کل هزینه تولید $ / ساعت
دیسپاچینگ٪ به دست آوردن دیسپاچینگ بهینه از نسل
٪ dpslack تفاوت (مقدار مطلق) بین است
٪ نسل شل برنامه ریزی تعیین شده از
٪ معادله هماهنگی، و نسل شل،
٪ از راه حل جریان قدرت به دست آمده.
در حالی که dpslack> 0.01٪ تست را برای همگرایی
loadflownewton٪ محلول جریان برق جدید
ضرایب٪ از دست دادن از دست دادن به روز می شود
دیسپاچینگ٪ دیسپاچینگ بهینه از نسل. با جدید B-ضرایب
پایان
٪ از٪ چاپ راه حل جریان قدرت نهایی
DISP (تکنولوژی)
fprintf ( ‘حداکثر قدرت عدم تطابق =٪ G \ N’، maxerror)
fprintf ( ‘تعداد تکرارها =٪ G \ n \ n در، ITER)
سر = [ ‘باس ولتاژ زاویه —— —— بار — — نسل تزریق،
‘تعداد مجله. درجه MW مگاوار MW مگاوار مگاوار،
”]؛
DISP (سر)
برای n = 1: nbus
fprintf ( ‘٪ 5G، N)، fprintf (‘٪ 7.3f، VM (N))،
fprintf ( ‘٪ 8.3f، deltad (N))، fprintf (‘٪ 9.3f، PD (N))،
fprintf ( ‘٪ 9.3f، QD (N))، fprintf (‘٪ 9.3f، صفحه (N))،
fprintf ( ‘٪ 9.3f، g Q را (ها))، fprintf (‘٪ 8.3f \ N ‘، QSH (N))
پایان
fprintf ( ‘\ N’)، fprintf ( “مجموع”)
fprintf ( ‘٪ 9.3f، PDT)، fprintf (‘٪ 9.3f، Qdt)،
fprintf ( ‘٪ 9.3f، PGT)، fprintf (‘٪ 9.3f، Qgt)، fprintf ( ‘٪ 9.3f \ n \ n در، Qsht)
gencost هزینه٪ نسل با برنامه ریزی مطلوب GEN.
ترجمه برنامه پاور فلو
NS = 0؛ NG = 0؛ VM = 0؛ دلتا = 0؛ yload = 0؛ deltad = 0؛
nbus = طول (busdata (:، 1))؛
برای k = 1: nbus
N = busdata (K، 1)؛
KB (N) = busdata (K، 2)؛ VM (N) = busdata (K، 3)؛ دلتا (N) = busdata (K، 4)؛
PD (N) = busdata (K، 5)؛ QD (N) = busdata (K، 6). صفحه (N) = busdata (K، 7)؛ QG (N) = busdata (K، 8)؛
Qmin (N) = busdata (K، 9)؛ با Qmax (N) = busdata (K، 10)؛
QSH (N) = busdata (K، 11)؛
اگر VM (N) <= 0 VM (N) = 1.0؛ V (N) = 1 + J * 0؛
دلتا دیگری (N) = PI / 180 * دلتا (N)؛
V (N) = VM () COS (دلتا (N) + J * گناه (دلتا (N))) (N) *؛
P (N) = (صفحه (N) -Pd (N)) / basemva؛
Q (N) = (g Q را (N) -Qd (N) + QSH (N)) / basemva؛
S (N) = P (N) + J * Q (N)؛
پایان
پایان
برای k = 1: nbus
اگر KB (K) == 1، NS = NS + 1؛ دیگری، پایان
اگر KB (K) == 2 نانوگرم = NG + 1؛ دیگری، پایان
NGS (K) = NG؛
NSS (K) = NS؛
پایان
YM = ABS (Ybus)؛ T = زاویه (Ybus)؛
متر = 2 * nbus-NG-2 * NS؛
maxerror = 1؛ همگرا = 1؛
ITER = 0؛
٪ شروع از تکرار
پاک کردن یک DC J DX
در حالی که maxerror> = دقت و ITER <= maxiter٪ تست را برای حداکثر. عدم تطابق قدرت
برای i = 1: متر
برای k = 1: متر
(من، K) = 0؛ ٪ ماتریس ژاکوبین مقدار دهی اولیه
پایان، پایان
ITER = ITER + 1؛
برای n = 1: nbus
NN = N-NSS (N)؛
LM = nbus + N-NGS (N) -nss (N) -ns؛
J11 = 0؛ J22 = 0؛ J33 = 0؛ J44 = 0؛
برای i = 1: NBR
اگر NL (من) == N | NR (من) == N
اگر NL (من) == N، L = NR (من)؛ پایان
اگر NR (من) == N، L = NL (من)؛ پایان
J11 J11 = + VM (N) * VM (L) * YM (N، L) * گناه (T (N، L) -) دلتا (N) + دلتا (L)؛
J33 J33 = + VM (N) * VM (L) * YM (N، L) * COS (T (N، L) -) دلتا (N) + دلتا (L)؛
اگر KB (N) ~ = 1
J22 J22 = + VM (L) * YM (N، L) * COS (T (N، L) -) دلتا (N) + دلتا (L)؛
J44 J44 = + VM (L) * YM (N، L) * گناه (T (N، L) -) دلتا (N) + دلتا (L)؛
دیگری، پایان
اگر KB (N) ~ = 1 & KB (L) ~ = 1
LK = nbus + L-NGS (L) -nss (L) -ns؛
LL = L -nss (L)؛
٪ تخفیف diagonalelements از J1
A (NN، LL) -Vm = (N) * VM (L) * YM (N، L) * گناه (T (N، L) – دلتا (N) + دلتا (L))؛
اگر KB (L) == 0٪ عناصر خارج قطر J2
A (NN، LK) = VM (N) * YM (N، L) * COS (T (N، L) – دلتا (N) + دلتا (L))؛ پایان
اگر KB (N) == 0٪ عناصر خارج مورب از J3
A (LM، LL) -Vm = (N) * VM (L) * YM (N، L) * COS (T (N، L) -) دلتا (N) + دلتا (L)؛ پایان
اگر KB (N) == 0 و KB (L) == 0٪ عناصر خارج قطر J4
A (LM، LK) = -Vm (N) * YM (N، L) * گناه (T (N، L) – دلتا (N) + دلتا (L))؛ پایان
پایان دیگری
دیگری، پایان
پایان
PK = VM (N) ^ 2 * YM (N، N) * COS (T (N، N)) + J33؛
QK = -Vm (N) ^ 2 * YM (N، N) * گناه (T (N) N) – J11؛
اگر KB (N) == 1 P (N) = PK؛ Q (N) = QK؛ اتوبوس پایان٪ نوسان P
اگر KB (N) == 2 Q (N) = QK؛
اگر با Qmax (N) ~ = 0
Qgc = Q (N) * basemva + QD (N) – QSH (N)؛
اگر ITER <= 7٪ بین 2TH و 6 تکرار
اگر ITER> 2٪ مگاوار اتوبوس ژنراتور
اگر Qgc <Qmin (N)،٪ تست شده است. اگر نه در محدوده VM (N)
VM (N) = VM (N) + 0.01؛ ٪ در مراحل 0.01 PU به تغییر
ELSEIF Qgc> Qmax به (N)،٪ را از ژنراتور مگاوار در
VM (N) = VM (N) – 0.01؛ پایان٪ محدوده مشخص.
دیگری، پایان
دیگری، پایان
دیگری، پایان
پایان
اگر KB (N) ~ = 1
A (NN، NN) = J11؛ ٪ عناصر مورب از J1
DC (NN) = P (N) -PK؛
پایان
اگر KB (N) == 0
A (NN، LM) = 2 * VM (N) * YM (N، N) * COS (T (N، N)) + J22؛ ٪ عناصر مورب J2
A (LM، NN) = J33؛ ٪ عناصر مورب از J3
A (LM، LM) = -2 * VM (N) * YM (N، N) * گناه (T (N) N) – J44؛ ٪ مورب از عناصر J4
DC (LM) = Q (N) -Qk؛
پایان
پایان
DX = A \ DC ‘؛
برای n = 1: nbus
NN = N-NSS (N)؛
LM = nbus + N-NGS (N) -nss (N) -ns؛
اگر KB (N) ~ = 1
دلتا (N) = دلتا (N) + DX (NN)؛ پایان
اگر KB (N) == 0
VM (N) = VM (N) + DX (LM)؛ پایان
پایان
maxerror = حداکثر (ABS (DC))؛
اگر ITER == maxiter و maxerror> دقت
fprintf ( ‘\ nWARNING: راه حل تکرار شونده بعد از همگرا نیست’)
fprintf ( ‘٪ G’، ITER)، fprintf ( ‘تکرار است. \ n \ N’)
fprintf ( ‘مطبوعات را وارد کنید برای پایان دادن به تکرار و چاپ نتایج را \ n’)
همگرا = 0؛ مکث، دیگری، پایان
پایان
اگر همگرایی ~ = 1
تکنولوژی = ( ‘راه حل تکرار شونده آیا همگرا نیست’)؛ دیگری،
تکنولوژی = ( ‘جریان راه حل های قدرت نیوتن رافسون روش’)؛
پایان
. V = VM * COS (دلتا) + J * VM * گناه (دلتا).
deltad = 180 / PI * دلتا؛
من = SQRT (-1)؛
K = 0؛
برای n = 1: nbus
اگر KB (N) == 1
K = K + 1؛
S (N) = P (N) + J * Q (N)؛
صفحه (N) = P (N) * basemva + PD (N)؛
QG (N) = Q (N) * basemva + QD (N) – QSH (N)؛
PGG (K) = صفحه (N)؛
Qgg (K) = g Q را (N)؛ ٪ ژوئن 97
ELSEIF KB (N) == 2
K = K + 1؛
S (N) = P (N) + J * Q (N)؛
QG (N) = Q (N) * basemva + QD (N) – QSH (N)؛
PGG (K) = صفحه (N)؛
Qgg (K) = g Q را (N)؛ ٪ ژوئن 1997
پایان
yload (N) = (PD (N) – J * QD (N) + J * QSH (N)) / (basemva * VM (N) ^ 2)؛
پایان
busdata (:، 3) = VM ‘؛ busdata (:، 4) = deltad ‘؛
PGT = مجموع (PG)؛ Qgt = SUM (g Q را). PDT = مجموع (PD)؛ Qdt = مجموع (QD)؛ Qsht = SUM (QSH)؛
٪ روشن A DC DX J11 J22 J33 J44 QK دلتا LK LM LL
٪ روشن A DC DX J11 J22 J33 QK دلتا LK LM LL