3节电力系统状态估计(WLS算法).docxVIP

3节电系统状态估计报告【任务说明】3节点系统主接线图任务：采用最小二乘状态估计算法,所有量测的权重都取1.0,编写状态估计程序(C/Matlab)。按量测类型，列出量测方程（每一类写出一个方程）画出程序流程提交源程序，程序中每个函数的作用提交计算的输出结果(屏幕拷贝)系统参数：功率基值：100MW电压基值：230 kV线路阻抗参数（标么值）：线路首节点末节点R(标么值)X(标么值)半充电电纳Line1-2120.017 0.0920Line1-3130.00850.072 0Line2-3230. 01190.10080线路量测(流出母线为正)：线路首端量测首端真值末端量测末端真值Line1-2-61.3+j1.2-63.32+j1.2460+j2.463.99+j2.35Line1-3-46.7-j14.8-47.69-j14.7945.9+j16.547.9+j16.54Line2-324-j6.624.93-j6.6-24+j7.2-24.86+j7.24母线电压量测：母线电压量测真值(幅值/角度)1232233/02234.56235.56/3.233236.46236.46/1.82负荷量测(流出母线为正)：负荷量测真值Load1111+j13.50111.02+j13.55发电量测(流入母线为正)：发电机量测真值unit288-j4.2488.92-j4.24unit323+j2423.04+j23.78注：量测存在误差【数据预处理】首先根据基值将已知的量测值均转换为标幺值，并将功率值转换为流入量，得到如下数据：线路导纳参数（标么值）：线路首节点末节点G(标么值)B(标么值)半充电电纳Line1-2121.9422-10.51070Line1-3131.6171-13.698 0Line2-3231.1551-9.78430线路注入功率量测(标幺值)：线路首端量测首端真值末端量测末端真值Line1-20.613-j0.0120.6332-j0.0124-0.60-j0.024-0.6399-j0.0235Line1-30.467+j0.1480.4769+j0.1479-0.459-j0.165-0.479-j0.1654Line2-3-0.24+j0.066-0.2493+j0.0660.24-j0.0720.2486-j0.0724负荷点注入功率量测(标幺值)：负荷量测真值Load1-1.11-j0.1350-1.1102-j0.1355发电机节点注入量测(流入母线为正)：发电机量测真值unit20.88-j0.04240.8892-j0.0424unit30.23+j0.240.2304+j0.2378母线电压量测（标幺值）：母线电压量测真值(幅值/角度)11.00871.0130/021.01981.0242/3.2331.02811.0281/1.82【量测方程】选择节点1的电压相角为参考，为0度，以vi表示误差值。节点1电压量测方程：Vi=Vi+v1即1.0087=V1+v12）1-3支路1号节点处注入有功功率功率：Pij=Vi2gij-ViVj(gijcos+bijsin)+v20.613=V12g13-V1V3(g13cos+b13sin)+v2即0.613=-1.6171V12-V1V3(-1.6171cos+13.698sin)+v23）1号节点注入功率：Pi=Vi2Gii+Gijcos+Bijsin+v3P1=V12G11+G1jcos+B1jsin+v3即-1.11=3.5613V12+V1V2(-1.9442cos-10.5107sin)+V1V3(-1.6171 cos-13.698 sin)+v3【流程图】遥测数据给定V0，，k=0求解AX=b,得X计算H(V(k),)和h(V(k),)A=HTR-1H, b=HTR-1(Z-h)k=k+1X(k+1)=X(k)+Xmax|X|NY结束【计算结果】其中iterations为迭代次数，可见本例的迭代次数为4，收敛较快，状态估计得到的节点1、2、3电压分别为：234.01【程序说明】1、计算h矩阵的函数cal_hfunction h=cal_h(V,th0,B,G) %其中，V为节点电压估计值，th0为节点电压相角估计%值，B为节点电导矩阵，G为节点电纳矩阵b=-B; %线路电导矩阵g=-G; %线路电纳矩阵P=zeros(3,1); %初始化，节点注入功率Q=zeros(3,1);PP=zeros(3,3); %线路注入功率QQ=PP;th=[0;th0]; %节点1的电压相角为0for i=1:3 P_P=0; Q_Q=0;for j=1:3if(j~=i)P_P=P_P+V(i)*V(j)*(G(i,j)*cos(th(i