想请板上英雄帮实现程序

[半條龍]

:confused:小弟最近在写一段程式遇到个问题
在此跪求恳请板上英雄帮解
或一同帮想方法
恳请拜托 谢谢
-----------------------------------------------------------------------范例数据程序
clear all;
clc
table1 = [ 1 0 1 1
2 1 2 2
3 2 3 3
4 2 4 1 ];
% col1 支路号 col2 起始点 col3 终止点 col4 支路类型
col1 = table1(:,1);
col2 = table1(:,2);
col3 = table1(:,3);
col4 = table1(:,4);
table2 = [1001 1002 1003 ...
2001 2002 2003 ...
3001 3002 3003 ...
4001 4003]';
% col5 取编码的千位也就是支路号(用来知道V与 I 哪几个是属于哪个支路) col6 取编码的个位(可知道三相支路三点是否全通或是缺少哪点)
col5 = floor(table2 /1000);
col6 = mod(table2 ,10);
Vtr = [10;12;7;9;11;4;14;13;21;3;8];
-----------------------------------------------------------------------問
以Vtr与table2关系来看会看出
Vtr=[10 ----->1001
12 ----->1002
7 ----->1003
9 ----->2001
11 ----->2002
4 ----->2003
14 ----->3001
13 ----->3002
21 ----->3003
3 ----->4001
8] ----->4003
以table1表格支路号col1=3当范例
支路号col1=3 起始点col2=2 终止点col3=3
那起始点Vtr(VS)(看千位数=2)
找到 9 11 4 (VS1=个位数=1,VS2=个位数=2,VS3=个位数=3)
VS1 = 9 VS2 = 11 VS3 = 4
那终止点Vtr(VF))(看千位数=3)
找到 14 13 21 (VF1=个位数=1,VF2=个位数=2,VF3=个位数=3)
VF1 = 14 VF2 = 13 VF3 = 21

类型,col4编号 = 1 注入起始电 注入终止电 都是0

类型,col4编号 = 2
注入终止I
I = [ (VF1 - (VS1-VF1))
VF2 - (VS2-VF2))
VF3 - (VS2-VF2)) ]
注入起始I
I = [ (VS1 - (VS1-VF1))
VS2 - (VS2-VF2))
VS3 - (VS2-VF2)) ]

类型,col4编号 = 3
注入终止I
I = [ VF1 - (VS2-VF1) ---> 仔細看是 VS2
VF2 - (VS3-VF2) ---> 仔細看是 VS3
VF3 - (VS1-VF3) ] ---> 仔細看是 VS1
注入起始I
I = [ (VS1-VS3) - (VS1-VS2) + (VS1-VF3)
(VS1-VS2) - (VS2-VS3) + (VS2-VF1)
(VS2-VS3) - (VS1-VS3) + (VS3-VF2)]

------------------步骤分解:
(1)看table1支路号col1 = 1 支路类型编号col4 = 1找到 终止点col3 = 1 起始点col2 = 0
因为支路类型编号col4 = 1 (遇到1 I就是直接代0)
I=[0; 资料库ITR=[0;
0; 0;
0]; 0;]
因为起始点col2 = 0 所以不再动作
(2)看table1支路号col1 = 2 支路类型编号col4 = 2找到 终止点col3 = 2 起始点col2 = 1
因为支路类型编号col4 = 1
I = [ (VF1 - (VS1-VF1)) =[9 -(10-9)
VF2 - (VS2-VF2)) 11-(12-11)
VF3 - (VS2-VF2)) ] 4-(7-4) ]
ITR=[0;
0;
0;
8;
10;
1 ];
因为起始点col2 = 1 不是0所以刚刚终止点 = 1 的矩阵要多加注入起始I
I = [ (VS1 - (VS1-VF1)) =[10 -(10-9)
VS2 - (VS2-VF2)) 12-(12-11)
VS3 - (VS2-VF2)) ] 7-(7-4) ]
ITR=[0+9 =[9
0+11; 11
0+4; 4
8; 8
10; 10
1 ]; 1]
(3)看table1支路号col1 = 3 支路类型编号col4 = 3找到 终止点col3 = 3 起始点col2 = 2
因为支路类型编号col4 = 3
I = [ VF1-(VS2-VF1) ---> 仔細看 VS2 =[14-(11-14)
VF2-(VS3-VF2) ---> 仔細看 VS3 13-(4-13)
VF3-(VS1-VF3)] ---> 仔細看 VS1 21-(9-21)
ITR=[9
11
4
8
10
1
17
22
33]
因为起始点col2 = 2 不是0所以刚刚终止点 = 2 的矩阵要多加注入起始I
I = [ (VS1-VS3)-(VS1-VS2)+(VS1-VF3) =[(9-4)-(9-11))+(9-21)
(VS1-VS2)-(VS2-VS3)+(VS2-VF1) (9-11)-(11-4)+(11-14)
(VS2-VS3)-(VS1-VS3)+(VS3-VF2)] (11-4))-(9-4)+(4-13)]
ITR=[9 =[9
11 11
4 4
8+(-5) 3
10+(-12) -2
1+(-7) -6
17 17
22 22
33] 33]
(3)看table1支路号col1 = 4 支路类型编号col4 = 1找到 终止点col3 = 4 起始点col2 = 2
因为支路类型编号col4 = 1
I=[0; 因为col1 = 4而table2千位数=4只有4001,4003兩个
0]
ITR=[9
11
4
3
-2
-6
17
22
33
0
0]
因为起始点col2 = 2 不是0所以刚刚终止点 = 2 的矩阵要多加注入起始I
ITR=[9
11
4
3 +0 因为table2个位数只有1,3 所以相对只加入在 2001 2003
-2
-6 +0
17
22
33
0
0]
-----------------------------------------------------------------------答案
ITR
ITR=[9
11
4
3
-2
-6
17
22
33
0
0]