你能帮助这个错误吗?
指数在位置2超出数组边界。
错误91行。
% |汽车| |型垂直地震剖面|θ| PGi | QGi | PLi | QLi |
1.06 busdata1 = [1 1 0 0 0 0 0
2 2 0 0 0 0 0.8 0.25);
% IEEE-30公交系统线数据. .
X R %从| | | | |
% |汽车|汽车|聚氨酯|聚氨酯|
linedata1 = (1 2 0.0 - 0.15);
% b总线的数据行(并联导纳)的形成。
% b总线返回. .
fb = linedata1 (: 1);
结核病= linedata1 (:, 2);
nbus = max (max(神奇动物),max (tb));
%。公共汽车……
nbranch =长度(神奇动物);
%。树枝……
bbus = 0 (nbus nbus);
为
k = 1: nbranch
bbus (tb (k), facebook (k)) = bbus (fb (k),结核(k));
结束
%程序形成导纳和阻抗总线形成....
%与变压器抽头设置. .
%返回ybus
% linedata3打电话”。m”行数据……
fb = linedata1 (: 1);
%从总线数字…
结核病= linedata1 (:, 2);
%总线数字…
r = linedata1 (:, 3);
%的阻力,R……
x = linedata1 (:, 4);
%电抗X……
z = r + * x;
% Z矩阵…
y = 1. / z;
%获得每个元素的倒数……
nbus = max (max(神奇动物),max (tb));
%。公共汽车……
nbranch =长度(神奇动物);
%。树枝……
ybus = 0 (nbus nbus);
%初始化YBus……
%的形成了对角元素……
为
k = 1: nbranch
ybus (fb (k),结核(k)) = ybus (fb (k),结核(k)) - y (k);
ybus (tb (k), facebook (k)) = ybus (fb (k),结核(k));
结束
%的对角元素....
为
m = 1: nbus
为
n = 1: nbranch
如果
fb (n) = = m
ybus (m m) = ybus (m m) + y (n);
elseif
结核病(n) = = m
ybus (m m) = ybus (m m) + y (n);
结束
结束
结束
%程序牛顿负载流分析
% 1号公共汽车被认为是希隆斯克汽车. .
%得到Y-Bus . .
busdata = busdata1 ();
%获得总线数据. .
总线= busdata1 (: 1);
%的巴士数量. .
类型= busdata1 (:, 2);
%的总线1-Slack、2-PV 3-PQ . .
V = busdata1 (: 3);
指定电压% . .
德尔= busdata1 (: 4);
%电压角度. .
Pg = busdata1 (: 5);
% PGi . .
路上= busdata1 (: 6);
% QGi . .
Pl = busdata1 (: 7);
% PLi . .
Ql = busdata1 (: 8);
% QLi . .
nbus = max(总线);
%没有。的巴士. .
P = Pg - Pl;
%π= PGi - PLi . .
Q =路上- Ql;
%气= QGi - QLi . .
托尔= 10;
%公差保持在高价值。
Iter = 1;
%迭代开始
Psp = P;
Qsp =问;
G =实际(y);
%电导。
B =图像放大(y);
%电纳. .
pq =找到(type = = 1);
%指数PQ巴士. .
npq =长度(pq);
%的PQ巴士. .
而
(Tol > 1 e-5)
%迭代开始. .
P = 0 (nbus, 1);
Q = 0 (nbus, 1);
%计算P和Q
为
我= 1:nbus
为
k = 1: nbus
P (i) = P (i) + V (i) * V (k) * (G (i (k) * cos (del (i)德尔(k)) + B (i (k) * sin (del (i)德尔(k)));
问(我)=问(我)+ V(我)* (k) * (G (i (k) * sin (del (i)德尔(k) - B (i (k) * cos (del (i)德尔(k)));
结束
结束
