你能详细说明你想确切地了解什么？
此外，使用<插入一行代码>按钮对他人更容易阅读您的代码。
我需要帮助了解Load Flow Newton Raphson代码
8个视图(过去30天)
显示较旧的评论
格式短G.
disp ('表9.2页面＃337行数据为9.2'）
％总线到总线r x g b
％总充电MVAR Y / 2
Linedata =[1 2 0.01008, 0.05040, 3.815629， -19.078144, 10.25, 0.05125;
1 3 0.00744、0.03720、5.169561、-25.847809、7.75、0.03875;
2 4 0.00744、0.03720、5.169561、-25.847809、7.75、0.03875;
3 4 0.01272,0.06360,3.023705，-15.118528,12.75,0.06375]
disp (' table 9.3 page # 338 bus data for example 9.2 '）
% Bus PGi QGi PLi QLi V Delta总线类型
Busdata =[1 0, 0, 50, 30.99, 1.00, 0 1;
2 0、0、170、105.35、1.00、0 2;
3 0、0、200、123.94、1.00、0 2;
4 318、0、80、49.58、1.02、0 3]
％最后一列显示总线类型：1.slack总线2.PQ总线3.PV总线
％=================================================.==程序开始此处===================================================％
ss = i * linedata（:,8）;% Y / 2
y = leinata（:,5）+ i * leinata（:,6）;
总值= max（max（linedata（：，1）），max（Linedata（:,2）））;％总线总线
rotsbranches =长度（Linedata（：，1））;%。的分支
YBUS = Zeros（总排水沟，总禁止）;
w = 0;
u = 0;
为了n = 1：总空间％总值不含PV总线
如果Busdata（N，2）> 0
w = w + 1;
结尾
结尾
为了n = 2：总值％总值没有PQ总线
如果busdata (n, 2) = = 0
u = u + 1;
结尾
结尾
％------初始化YBUS ------
为了B = 1：截止银行
1) ybus (linedata (b)、(linedata (b, 2))) = - y (b);
ybus (linedata (b, 2)), (linedata (b - 1))) = ybus ((linedata (b - 1)), (linedata (b, 2)));
结尾
为了C = 1：总值
为了d = 1: totalbranches
如果Linedata (d,1) == c |
YBUS（C，C）= YBUS（C，C）+ Y（d）+ ss（d）;
结尾
结尾
结尾
disp ('表9.3页面＃338总线入场矩阵为9.2'）
ybus
ybusR = 0 (totalbuses totalbuses);
YBUSA = Zeros（总空间，全排放量）;
z =零（总排水管，8）;
busnumber = busdata (: 1);
PG = busdata (:, 2);
路上= busdata (: 3);
PL = BUSDATA（:,4）;
ql = busdata（:,5）;
V = busdata (: 6);
VV = V;
Ang = Busdata（:,7）;
类型= busdata (: 8);
s = (2 * totalbuses) -w-2;
j =零（s，s）;％初始化空雅各比矩阵
mm =零（s，1）;
一个= Ang;
p =（pg-pl）./ 100;在公共汽车上每单位有源电力
q =（qg-ql）./ 100;每单位无功功率在母线
tol = 1;
iter = 0;
kk =输入（“进入迭代容忍度”)；
%------------------------------ 矩形到极 --------------------------------
为了B = 1：截止银行
ybus的实部
1) ybusR (linedata (b)、(linedata (b, 2))) = abs (- y (b));
YBUSR（（LoData（B，2）），（Linedata（B，1）））= YBUSR（（Linedata（B，1）），（Linedata（B，2）））;
ybusR (b, b) = abs (-ybus (b, b));
YBUS的％角度
1) ybusA (linedata (b)、(linedata (b, 2))) =角(- y (b));
YBUSA（（LoData（B，2）），（LoData（B，1）））= YBUSA（（Linedata（B，1）），（Linedata（B，2）））;
ybusA (b, b) =角(-ybus (b, b));
结尾
ybusR;
ybusA;
而托尔>乐
％------------------------------------------------％
为了我= 2:totalbuses
YVV = 0;
为了n = 1：总空间
如果我~ = n
YVV = YVV +（YBUSR（i，n）* v（n）* v（i））* cos（（ybusa（i，n））+ busdata（n，7）-busdata（i，7））;％乘法导纳和电压
结尾
YVV;
结尾
电脑(i) = ((abs (V (i)) ^ 2) *实际(ybus(我)))+ YVV;%计算计算出的P。
deltap（i）= p（i）-pc（i）;
结尾
个人电脑;
dp = deltaP ';
为了我= 2：总空间
张毫米(1)= dp(我);
结尾
%------------------------ δQ -----------------------------%
为了i = 2：总值-W
如果Busdata（I，8）== 2
YVV = 0;
为了n = 1：总空间
如果我~ = n
YVV = YVV + (ybusR(我,n) * V (n) *(我))*罪((ybusA(我,n)) + busdata (n, 7) -busdata(我7));％乘法导纳和电压
结尾
YVV;
结尾
qc（i）= - （（abs（abs（v（i））^ 2）* imag（ybus（i，i））） - YVV;%计算计算出Q。
deltaq（i）= q（i）-qc（i）;
结尾
结尾
QC;
dq = deltaq';
为了我= 2:totalbuses-w
如果Busdata（I，8）== 2
毫米(totalbuses-2 +我,1)= dq(我);
别的
我张=;
结尾
结尾
J11 = 0 (totalbuses-1 totalbuses-1);
纺纱,= 0 (totalbuses-1 u);
J21 =零（U，TotalBuses-1）;
J22 = Zeros（U，U）;
％----------------------------------------------------％
为了我= 2：总空间
为了j = 2：总值
如果我~ = j
J11（I-1，J-1）= -Abs（（（ybusr（i，j）* v（i）* v（j）））* sin（（ybusa（i，j））+ busdata（j，7）-busdata（i，7）））;
结尾
如果我= =我
为了n = 1：总空间
如果n ~ =我
张J11(张)= J11(张、张)+ (abs (((ybusR(我,n) * V(我)* (n))) *罪((ybusA(我,n)) + busdata (n, 7) -busdata(我7))));
结尾
结尾
结尾
结尾
结尾
J11;
%---------------------------- 纺纱, -------------------------------%
为了我= 2：总空间
为了j = 2: totalbuses-w
如果Busdata（J，8）== 2
如果我~ = j
纺纱,(张,j - 1) = (abs (V (j)) * abs ((ybusR (i, j) * V (i)))) * cos ((ybusA (i, j)) + busdata (j, 7) -busdata(我7));
结尾
YV12 = 0;
如果我= =我
为了n = 1：总空间
如果n ~ =我
yv12 = yv12 + abs ((ybusR(我,n) * V (n))) * cos ((ybusA(我,n)) + busdata (n, 7) -busdata(我7));
结尾
结尾
J12（I-1，J-1）= ABS（v（i））*（2 * ABS（v（i））*真实（YBUS（I，I））+ YV12）;
结尾
结尾
结尾
结尾
纺纱,;
%---------------------------- J21 -------------------------------%
为了我= 2:totalbuses-w
如果Busdata（I，8）== 2
为了j = 2：总值
如果我~ = j
J21(张,j - 1) = abs (((ybusR (i, j) * V(我)* (j))) * cos ((ybusA (i, j)) + busdata (j, 7) -busdata(我7)));
结尾
如果我= =我
为了n = 1：总空间
如果n ~ =我
J21(张,j - 1) = J21 (j - 1张)+ (abs (((ybusR(我,n) * V(我)* (n))) * cos ((ybusA(我,n)) + busdata (n, 7) -busdata(我7))));
结尾
结尾
结尾
结尾
结尾
结尾
J21;
％----------------------------------------------------％
为了我= 2:totalbuses-w
为了j = 2: totalbuses-w
如果BusData（I，8）== 2 && BusData（J，8）== 2
如果我~ = j
J22（I-1，J-1）= - （ABS（v（j））* abs（（ybusr（i，j）* v（i）））））* sin（（ybusa（i，j））+Busdata（J，7）-BusData（I，7））;
结尾
YV12 = 0;
如果我= =我
J22（I-1，J-1）= - （J11（I-1，I-1）） - （2 *（ABS（v（i））^ 2）* imag（ybus（i，i）））;
结尾
结尾
结尾
结尾
J22;
%---------------------------------------------------------------------------------%
为了k = 1: totalbuses-1
为了l = 1: totalbuses-1
J (k, l) = J11 (k, l);
结尾
结尾
为了k = 1：你
为了l = 1: totalbuses-1
J的(l, k + totalbuses-1) =纺纱,(l, k);
结尾
结尾
为了k = 1：你
为了l = 1: totalbuses-1
J (k + totalbuses-1, l) = J21 (k, l);
结尾
结尾
为了k = 1：你
为了l = 1：你
J (l + totalbuses-1, k + totalbuses-1) = J22的(l, k);
结尾
结尾
毫米;
J;
格式短G.
k = inv（j）* mm;
为了n = 1：总值-1
(n) = (ANG (n + 1) + k (n));
sa = radtodeg(年代);
%-------------
Busdata（n + 1,7）= s（n）;
%-------------
结尾
Ang = Busdata（:,7）;
sa);
为了m = 2：u + 1
v（m）= abs（Busdata（m，6）+ k（全排放量-2 + m））;
% abs (VV)
busdata (m, 6) = V (m);
结尾
v;
tol = max（abs（abs（v） - abs（vv）））&& max（abs（abs（Ang） - ABS（AN）））;
iter = iter + 1;
VV = V;
一个= Ang;
结尾
j
it
％----------------------------------------------------％
％真实（Vacc'）
Z（1：总值，1）= BusData（：，1）;
z (1: totalbuses, 2) = busdata (:, 8);
z (1: totalbuses, 3) = abs (busdata (:, 6));
z（1：总值，4）= Radtodeg（Busdata（:,7））;
z (1: totalbuses, 5) = PG;
z (1: totalbuses 7) = busdata (:, 4);
z (1: totalbuses 8) = busdata (:, 5);
disp ('|公共汽车号||巴士类型||电压||角度|| MW（g）|| MVAR（g）|| MW（L）|| MVAR（L）|')；
格式短G.
Z.
答案（0）
翻译的
你也可以从以下列表中选择一个网站: