%%?Power?Grid?State?Estimation?（Test?1）
clc;????%?Clear?windows
format

%%?使用變量說明
%?GATE??????對不良數據檢測與辨識的調用開關
%?JM????????min?（J）?殘差平方和
%?uij，k??用于循環的變量
%?RE????????residual?error（殘差矢量）
%?VA????????variance（方差）
%?FN????????標準化殘差
%?BDT???????Bad?Data‘s?Type?（壞數據類型）
%?BDP???????Bad?Data‘s?Position?（壞數據位置）
%?F?????????狀態估計方程計算出的值

%%?Read?power?flow?data?from?document?“iSE30Bus.txt“
fp=fopen（‘iSE30Bus1.txt‘‘r‘）;???%?1號為有量測誤差數據，2號為無誤差數據
fscanf（fp‘%s‘1）;???????????????%?read?one?sting
node=fscanf（fp‘%d‘1）;??????????%?read?number?of?nodes
nzl=fscanf（fp‘%d‘1）;???????????%?read?number?of?branches
kPQ=fscanf（fp‘%d‘1）;???????????%?read?number?of?known?PQ
kBus=fscanf（fp‘%d‘1）;??????????%?read?number?of?known?Bus
kV=fscanf（fp‘%d‘1）;????????????%?read?number?of?known?Vamp

fscanf（fp‘%s‘1）;???????????????%?read?a?string
ReV=fscanf（fp‘%f‘1）;???????????%?read?number?of?Real?Vamp
ImV=fscanf（fp‘%f‘1）;???????????%?read?number?of?Imag?Vamp

fscanf（fp‘%s‘1）;???????????????%?read?a?string
zl=zeros（nzl4）;?????????????????%?確定節點導納矩陣的規格，并保存i??j??R??X的數據
for?i=1:nzl
????for?j=1:4
??????zl（ij）=fscanf（fp‘%f‘1）;
????end??????????????????????????%?read?branch?parameters
end

fscanf（fp‘%s‘1）;???????????????%?read?one?string
N=zeros（kPQ3）;??????????????????%?save:確定所有節點的功率值，保存i?P?Q
for?i=1:kPQ???
????for?j=1:3
????????N（ij）=fscanf（fp‘%f‘1）;???%?read?generator?parameters
????end
end

PQ=zeros（2*kPQ1）;???????????????%?將節點功率轉換為列向量的形式
for?i=1:kPQ
????PQ（2*i-1）=N（i2）;
????PQ（2*i）=N（i3）;
end

fscanf（fp‘%s‘1）;???????????????%?read?one?string
Bus=zeros（kBus4）;???????????????%?save?known?Bus:?i??j??Pij??Qij
for?i=1:kBus
????for?j=1:4
??????Bus（ij）=fscanf（fp‘%f‘1）;
????end?????????????????????????
end

fscanf（fp‘%s‘1）;???????????????%?read?one?string
Vam=zeros（kV2）;?????????????????%?save?known?Vamp:?i?Vamp
for?i=1:kV
????for?j=1:2
??????Vam（ij）=fscanf（fp‘%f‘1）;
????end?????????????????????????
end
fclose（fp）;

%%?形成節點導納矩陣?G+jB
G=zeros（nodenode）;
B=zeros（nodenode）;
for?k=1:nzl
????i=zl（k1）;
????j=zl（k2）;
????h=zl（k3）*zl（k3）+zl（k4）*zl（k4）;
????f=zl（k3）/h;???????????????%?branch?conductance
????g=-zl（k4）/h;??????????????%?branch?susceptance
????G（ii）=G（ii）+f;
????G（jj）=G（jj）+f;
????B（ii）=B（ii）+g;
????B（jj）=B（jj）+g;
????G（ij）=G（ij）-f;
????G（ji）=G（ji）-f;
????B（ij）=B（ij）-g;
????B（ji）=B（ji）-g;
end
%?G???????????????????????????????%?show
%?B???????????????????????????????%?show


%%?State?Estimation
GATE=0;
zcs=0;
?for?u=1:2*kPQ+2*kBus+kV-2*node+2
???
?????if?GATE==1
?????????break;
?????????JM=RE‘*RE;?????????????????????????%?不良數據檢測與辨識
?????????JM
?????????if?JM<0.00003??????????????????????%?不良數據檢測
?????????????break;
?????????else
?????????????VA=JM/（2*kPQ+2*kBus+kV-1）;?????%?VA：殘差的方差
?????????????VA1=1/sqrt（VA）;
?????????????FN=RE*VA1;?????????????????????%?FN：經過標準化處理的殘差向量
?????????????FN1=zeros（2*kPQ+2*kBus+kV2）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????15428??2019-02-02?18:39??SEF.m

?????文件???????3121??2018-05-05?23:19??iSE30Bus1.txt

?????文件???????3101??2018-05-01?17:54??iSE30Bus2.txt

-----------?---------??----------?-----??----

????????????????21650????????????????????3