Matlab 变压器Tt、T2、T3：SFT-40000/110，额定变比10.5/121kV,APg=215kW，V=10.5%，AP。=35kW，10=0.6%，线路DB段：l=75km， n=0.21.0/km，x=0.420/km；线路CB段：l=20km， r=0.152/km，X1=0.4.2/km;线路DC段：l=40km，i=0.220/km,x1=0.410/km；各段线路的对地导纳均可略去不计。节点E为平衡节点，电压取10.5kV，其余为PQ节点。负荷功率如表1.1所示。(Sp=100MVA，Va=Vao) (1) 列写节点导纳矩阵；

根据节点导纳矩阵的定义，我们需要先列出节点电流和电压之间的关系式，即：

I = YV

其中，I 和 V 分别表示节点电流和电压，Y 表示节点导纳矩阵。根据电力系统基本知识，我们可以将节点电流表示为：

I = Y (V - Vp)

其中，Vp 表示节点电压相角的余弦值。

为了方便计算，我们可以将节点按照顺序编号，从节点 1 开始。根据题目的条件，我们可以列出节点导纳矩阵的表达式：

Y11 = (1/RDB + 1i/XDB) + (1/RDC + 1i/XDC) + (1/R1 + 1i/X1)
Y12 = -1i/XDB
Y13 = -1i/XDC
Y14 = 0
Y15 = 0
Y16 = -1/RT2
Y17 = 0
Y18 = 0
Y19 = -1/RT3

Y21 = -1i/XDB
Y22 = (1/RDB + 1i/XDB) + (1/R1 + 1i/X1)
Y23 = 0
Y24 = -1/RT1
Y25 = 0
Y26 = 0
Y27 = 0
Y28 = 0
Y29 = 0

Y31 = -1i/XDC
Y32 = 0
Y33 = (1/RDC + 1i/XDC) + (1/R1 + 1i/X1)
Y34 = 0
Y35 = -1/RT4
Y36 = 0
Y37 = 0
Y38 = 0
Y39 = 0

Y41 = 0
Y42 = -1/RT1
Y43 = 0
Y44 = 1/RT1 + 1/RT2 + 1/RT4
Y45 = -1/RT4
Y46 = 0
Y47 = 0
Y48 = 0
Y49 = 0

Y51 = 0
Y52 = 0
Y53 = -1/RT4
Y54 = -1/RT4
Y55 = 1/RT4 + 1/RT5
Y56 = 0
Y57 = 0
Y58 = 0
Y59 = 0

Y61 = -1/RT2
Y62 = 0
Y63 = 0
Y64 = 0
Y65 = 0
Y66 = 1/RT2 + 1/RT3
Y67 = -1/RT3
Y68 = 0
Y69 = 0

Y71 = 0
Y72 = 0
Y73 = 0
Y74 = 0
Y75 = 0
Y76 = -1/RT3
Y77 = 1/RT3 + 1/RT6
Y78 = 0
Y79 = 0

Y81 = 0
Y82 = 0
Y83 = 0
Y84 = 0
Y85 = 0
Y86 = 0
Y87 = 0
Y88 = 1/RT5 + 1/RT6
Y89 = -1/RT5

Y91 = -1/RT3
Y92 = 0
Y93 = 0
Y94 = 0
Y95 = 0
Y96 = 0
Y97 = 0
Y98 = -1/RT5
Y99 = 1/RT3 + 1/RT5

其中，RDB、XDB、R1、X1、RDC、XDC、RT1、RT2、RT3、RT4、RT5、RT6 分别表示题目中的各个参数，需要根据实际情况进行替换。

需要注意的是，节点 1、2、3 分别对应于 Tt、T2、T3，因此它们的导纳矩阵中需要考虑变压器的影响。节点 4、5、6、7、8、9 分别对应于 DB、CB、DC、E、PQ1、PQ2。