????????歡迎來(lái)到本博客????????

??博主優(yōu)勢(shì)：??????博客內(nèi)容盡量做到思維縝密，邏輯清晰，為了方便讀者。

??座右銘：行百里者，半于九十。

??????本文目錄如下：??????

目錄

??1 概述

??2 運(yùn)行結(jié)果

??3?參考文獻(xiàn)

??4 Matlab代碼實(shí)現(xiàn)

---

??1 概述

文獻(xiàn)來(lái)源：

文獻(xiàn)分享閱讀：

我用夸克網(wǎng)盤(pán)分享了「【狀態(tài)估計(jì)】電力系統(tǒng)狀態(tài)估計(jì)的虛假數(shù)據(jù)注入攻擊建模與對(duì)策.pdf」，點(diǎn)擊鏈接即可保存。打開(kāi)「夸克APP」在線(xiàn)查看，支持多種文檔格式轉(zhuǎn)換。
鏈接：夸克網(wǎng)盤(pán)分享
提取碼：4vRm?

虛假數(shù)據(jù)注入攻擊（FDIA）已被證明是對(duì)抗電力系統(tǒng)的嚴(yán)重安全挑戰(zhàn)之一。這正成為電力公司日益關(guān)注的問(wèn)題，最近引起了電力系統(tǒng)研究人員和工程師的注意。因此，電力系統(tǒng)運(yùn)行和規(guī)劃中的狀態(tài)估計(jì)是監(jiān)測(cè)和控制系統(tǒng)通過(guò)電表測(cè)量和電力系統(tǒng)拓?fù)涔烙?jì)電力系統(tǒng)最佳狀態(tài)的重要和必不可少的工具。因此，本文提出了避免系統(tǒng)內(nèi)不必要的完全停電的建模和對(duì)策。本文還考慮了輸電和配電系統(tǒng)對(duì)FDIA的脆弱性。本文的結(jié)果可為制定針對(duì)電力系統(tǒng)網(wǎng)絡(luò)漏洞的必要保護(hù)對(duì)策奠定基礎(chǔ)。

原文摘要：

Abstract:

False Data Injection Attacks (FDIA) has been shown to be one of the serious security challenges combating power systems. This is becoming a grown concern to power utilities and has drawn the attention of power system researchers and Engineers in recent times. State estimation in power system operation and planning is therefore an important and an essential tool for monitoring and controlling the system to estimate the best state of the power system through meter measurements and power system topologies. This paper therefore presents the modeling and countermeasures for avoiding unnecessary total blackout within the system. The vulnerability of both the Transmission and the distribution power system to FDIA is also considered in this paper. The outcomes of this paper could serve as a basis for the development of necessary protective countermeasures against vulnerabilities within power system networks.

確保任何系統(tǒng)的數(shù)據(jù)完整性和安全性對(duì)于安全運(yùn)行和保護(hù)系統(tǒng)免受操縱至關(guān)重要。隨著傳統(tǒng)電網(wǎng)被智能電網(wǎng)取代，這就構(gòu)成了電網(wǎng)容易受到網(wǎng)絡(luò)攻擊的情況[1]。智能電網(wǎng)作為國(guó)家關(guān)鍵基礎(chǔ)設(shè)施，是一個(gè)將網(wǎng)絡(luò)和物理系統(tǒng)集成在一起的系統(tǒng)，該系統(tǒng)的不同部分之間不斷通信。因此，確保數(shù)據(jù)完整性并保護(hù)系統(tǒng)免受操縱非常重要。物聯(lián)網(wǎng)（IoT）的應(yīng)用使通信鏈路和整個(gè)電力系統(tǒng)容易受到攻擊。電網(wǎng)不能被描述為傳統(tǒng)的網(wǎng)絡(luò)系統(tǒng)，而是一個(gè)信息物理系統(tǒng)，其特征是網(wǎng)絡(luò)和物理系統(tǒng)之間的通信或集成[2]。智能電網(wǎng)帶來(lái)的進(jìn)步在自動(dòng)化、電網(wǎng)監(jiān)控和管理、智能計(jì)量基礎(chǔ)設(shè)施和客戶(hù)端管理領(lǐng)域值得稱(chēng)贊。由此產(chǎn)生的電力系統(tǒng)的通信網(wǎng)絡(luò)更容易受到入侵者的攻擊和操縱[2]。

智能電網(wǎng)的一個(gè)重要元素是狀態(tài)估計(jì)。狀態(tài)估計(jì)的輸入僅限于負(fù)載總線(xiàn)上的P，Q注入和電壓控制總線(xiàn)上的P，V值。需要狀態(tài)估計(jì)有兩個(gè)原因。首先，來(lái)自遠(yuǎn)程終端單元（RTU）和傳感器等現(xiàn)場(chǎng)設(shè)備的測(cè)量被認(rèn)為是有噪聲的，并且包含一些誤差。其次，在潮流的情況下，測(cè)量值包含一些變量，例如P，Q線(xiàn)流，這些變量對(duì)于潮流計(jì)算不是必需的。因此，運(yùn)行狀態(tài)的最佳估計(jì)是通過(guò)狀態(tài)估計(jì)獲得的，并且該估計(jì)值在被接受之前需要進(jìn)行統(tǒng)計(jì)分析[3]。

電力系統(tǒng)由能源管理系統(tǒng)（EMS）或監(jiān)督控制和數(shù)據(jù)采集（SCADA）系統(tǒng)持續(xù)監(jiān)測(cè)和控制，這些系統(tǒng)需要狀態(tài)估計(jì)器的輸入來(lái)維持電力系統(tǒng)的運(yùn)行狀態(tài)。SCADA系統(tǒng)從現(xiàn)場(chǎng)設(shè)備（如網(wǎng)絡(luò)中的傳感器）收集測(cè)量值。傳感器測(cè)量值用作狀態(tài)估計(jì)器的輸入。狀態(tài)估計(jì)器根據(jù)這些測(cè)量值估計(jì)電力系統(tǒng)的真實(shí)狀態(tài)。狀態(tài)估計(jì)器的輸出是一些電壓幅度和相位角?？刂苿?dòng)作和能量管理系統(tǒng)依賴(lài)于狀態(tài)估計(jì)器的輸出，以實(shí)現(xiàn)電力系統(tǒng)的安全運(yùn)行。因此，狀態(tài)估計(jì)使用來(lái)自?xún)x表測(cè)量的數(shù)據(jù)，計(jì)算這些測(cè)量的最佳估計(jì)值，然后將結(jié)果用于控制網(wǎng)格。因此，狀態(tài)估計(jì)在電力系統(tǒng)監(jiān)測(cè)和控制中起著重要作用[3]?[4]。

狀態(tài)估計(jì)給出了電力系統(tǒng)真實(shí)狀態(tài)的估計(jì)，大多數(shù)控制動(dòng)作、研究和系統(tǒng)操作都依賴(lài)于狀態(tài)估計(jì)。對(duì)電力系統(tǒng)狀態(tài)估計(jì)的攻擊往往會(huì)危及電網(wǎng)。電力系統(tǒng)模型可以是交流或直流，其中交流模型同時(shí)考慮實(shí)際功率和無(wú)功功率，因此該模型的方程變?yōu)榉蔷€(xiàn)性。對(duì)于簡(jiǎn)化的直流模型，方程是線(xiàn)性的，不需要迭代。

電網(wǎng)對(duì)攻擊的敏感性是由于電網(wǎng)作為網(wǎng)絡(luò)物理系統(tǒng)的架構(gòu)，通信鏈路更容易受到網(wǎng)絡(luò)物理攻擊。配備智能電表的先進(jìn)計(jì)量系統(tǒng)基礎(chǔ)設(shè)施的引入在實(shí)現(xiàn)智能電網(wǎng)方面發(fā)揮了重要作用。然而，由于電網(wǎng)上有許多輸入或插座可用，系統(tǒng)就越脆弱[4]。

針對(duì)電網(wǎng)的攻擊之一是虛假數(shù)據(jù)注入攻擊（FDIA）。FDIA 是一種針對(duì)電力系統(tǒng)中狀態(tài)估計(jì)的數(shù)據(jù)完整性攻擊，攻擊者會(huì)破壞傳感器節(jié)點(diǎn)。大多數(shù)情況下，攻擊者操縱和誤導(dǎo)系統(tǒng)操作員發(fā)出不準(zhǔn)確的命令，如果不快速檢測(cè)和控制，會(huì)導(dǎo)致電力系統(tǒng)停電[6]。因此，F(xiàn)DIA的后果是重大設(shè)備故障和完全停電。通過(guò)破壞智能電表測(cè)量、傳感器或遠(yuǎn)程終端單元 （RTU） 或侵入網(wǎng)絡(luò)上的監(jiān)控通信?[7]，可能會(huì)注入虛假數(shù)據(jù)。

攻擊者可以插入通常不會(huì)被BDD檢測(cè)到的測(cè)量值[8]。了解電力系統(tǒng)的當(dāng)前配置可能會(huì)導(dǎo)致BDD的當(dāng)前實(shí)現(xiàn)無(wú)法檢測(cè)到攻擊，因?yàn)楣粽呦Ｍ⑷肱c估計(jì)狀態(tài)密切相關(guān)的數(shù)據(jù)。大多數(shù)BDD技術(shù)基于最小加權(quán)平方法。了解電源系統(tǒng)配置肯定會(huì)違反這種方法。許多研究表明，保護(hù)一些戰(zhàn)略性選擇的儀表測(cè)量值，以防止FDIA [4]，[9]，[10]。確定一組儀表測(cè)量值，以便攻擊者注入 FDA 所需的矢量受到限制是一種很好的做法。但是，這種方法不會(huì)完全使系統(tǒng)免受攻擊，并且必須始終提供受保護(hù)傳感器的測(cè)量值。對(duì)于實(shí)際系統(tǒng)，可能無(wú)法實(shí)現(xiàn)儀表免受攻擊的絕對(duì)保護(hù)。由于SCADA設(shè)備和網(wǎng)絡(luò)拓?fù)涓鞑幌嗤?，因此正確監(jiān)控系統(tǒng)以確定最合適的FDIA實(shí)施非常重要。因此，需要制定超越這些傳統(tǒng)方法的檢測(cè)方法，因?yàn)楣粽呖梢岳米钚《朔ㄔ试S的誤差范圍之間的測(cè)量而不會(huì)被檢測(cè)到[6]。

因此，考慮到電網(wǎng)作為國(guó)家關(guān)鍵基礎(chǔ)設(shè)施的重要性以及FDIA對(duì)電力系統(tǒng)的嚴(yán)重影響，確保數(shù)據(jù)完整性和電網(wǎng)保護(hù)非常重要[5]。網(wǎng)絡(luò)計(jì)算應(yīng)用帶來(lái)的自動(dòng)化導(dǎo)致電網(wǎng)容易受到網(wǎng)絡(luò)攻擊。就像傳統(tǒng)的網(wǎng)絡(luò)系統(tǒng)一樣，系統(tǒng)上可用的許多插座或連接鏈接往往會(huì)增加其漏洞[1]。針對(duì)狀態(tài)估計(jì)的虛假數(shù)據(jù)注入攻擊 （FDIA） 被證明是攻擊者通過(guò)向電力系統(tǒng)注入錯(cuò)誤或錯(cuò)誤儀表測(cè)量來(lái)破壞電力系統(tǒng)運(yùn)行的一種方式。

大多數(shù)系統(tǒng)參數(shù)（如系統(tǒng)穩(wěn)定性和控制）都取決于狀態(tài)估計(jì)。狀態(tài)估計(jì)參數(shù)是做出基于系統(tǒng)的物理系統(tǒng)穩(wěn)定性和控制決策所必需的。FDIA可能會(huì)誤導(dǎo)系統(tǒng)運(yùn)營(yíng)商根據(jù)注入的數(shù)據(jù)做出不準(zhǔn)確的決策。最近，已經(jīng)證明虛假數(shù)據(jù)注入（FDI）攻擊可以繞過(guò)當(dāng)今EMS / SCADA系統(tǒng)中的錯(cuò)誤數(shù)據(jù)檢測(cè)（BDD）而不會(huì)被發(fā)現(xiàn)[2] [8] [11]。

??2 運(yùn)行結(jié)果

Del = 180/pi*del;
E2 = [V Del]; % Bus Voltages and angles..
disp('-------- State Estimation ------------------');
disp('--------------------------');
disp('| Bus |    V   |  Angle  | ');
disp('| No  |   pu   |  Degree | ');
disp('--------------------------');
for m = 1:num
    fprintf('%4g', m); fprintf('  %8.4f', V(m)); fprintf('   %8.4f', Del(m)); fprintf('\n');
end
disp('---------------------------------------------');

?

??

主函數(shù)代碼：

% Power System State Estimation using Weighted Least Square Method..

num = 30; % IEEE - 14 or IEEE - 30 bus system..(for IEEE-14 bus system replace 30 by 14)...
ybus = ybusppg(num); % Get YBus..
zdata = zdatas(num); % Get Measurement data..
bpq = bbusppg(num); % Get B data..
nbus = max(max(zdata(:,4)),max(zdata(:,5))); % Get number of buses..
type = zdata(:,2); % Type of measurement, Vi - 1, Pi - 2, Qi - 3, Pij - 4, Qij - 5, Iij - 6..
z = zdata(:,3); % Measuement values..
fbus = zdata(:,4); % From bus..
tbus = zdata(:,5); % To bus..
Ri = diag(zdata(:,6)); % Measurement Error..
V = ones(nbus,1); % Initialize the bus voltages..
del = zeros(nbus,1); % Initialize the bus angles..
E = [del(2:end); V]; ? % State Vector..
G = real(ybus);
B = imag(ybus);

vi = find(type == 1); % Index of voltage magnitude measurements..
ppi = find(type == 2); % Index of real power injection measurements..
qi = find(type == 3); % Index of reactive power injection measurements..
pf = find(type == 4); % Index of real powerflow measurements..
qf = find(type == 5); % Index of reactive powerflow measurements..

nvi = length(vi); % Number of Voltage measurements..
npi = length(ppi); % Number of Real Power Injection measurements..
nqi = length(qi); % Number of Reactive Power Injection measurements..
npf = length(pf); % Number of Real Power Flow measurements..
nqf = length(qf); % Number of Reactive Power Flow measurements..

iter = 1;
tol = 5;

while(tol > 1e-4)
? ??
? ? %Measurement Function, h
? ? h1 = V(fbus(vi),1);
? ? h2 = zeros(npi,1);
? ? h3 = zeros(nqi,1);
? ? h4 = zeros(npf,1);
? ? h5 = zeros(nqf,1);
? ??
? ? for i = 1:npi
? ? ? ? m = fbus(ppi(i));
? ? ? ? for k = 1:nbus
? ? ? ? ? ? h2(i) = h2(i) + V(m)*V(k)*(G(m,k)*cos(del(m)-del(k)) + B(m,k)*sin(del(m)-del(k)));
? ? ? ? end
? ? end
? ??
? ? for i = 1:nqi
? ? ? ? m = fbus(qi(i));
? ? ? ? for k = 1:nbus
? ? ? ? ? ? h3(i) = h3(i) + V(m)*V(k)*(G(m,k)*sin(del(m)-del(k)) - B(m,k)*cos(del(m)-del(k)));
? ? ? ? end
? ? end
? ??
? ? for i = 1:npf
? ? ? ? m = fbus(pf(i));
? ? ? ? n = tbus(pf(i));
? ? ? ? h4(i) = -V(m)^2*G(m,n) - V(m)*V(n)*(-G(m,n)*cos(del(m)-del(n)) - B(m,n)*sin(del(m)-del(n)));
? ? end
? ??
? ? for i = 1:nqf
? ? ? ? m = fbus(qf(i));
? ? ? ? n = tbus(qf(i));
? ? ? ? h5(i) = -V(m)^2*(-B(m,n)+bpq(m,n)) - V(m)*V(n)*(-G(m,n)*sin(del(m)-del(n)) + B(m,n)*cos(del(m)-del(n)));
? ? end
? ??
? ? h = [h1; h2; h3; h4; h5];
? ??
? ? % Residue..
? ? r = z - h;
? ??
? ? % Jacobian..
? ? % H11 - Derivative of V with respect to angles.. All Zeros
? ? H11 = zeros(nvi,nbus-1);

? ? % H12 - Derivative of V with respect to V..?
? ? H12 = zeros(nvi,nbus);
? ? for k = 1:nvi
? ? ? ? for n = 1:nbus
? ? ? ? ? ? if n == k
? ? ? ? ? ? ? ? H12(k,n) = 1;
? ? ? ? ? ? end
? ? ? ? end
? ? end

? ? % H21 - Derivative of Real Power Injections with Angles..
? ? H21 = zeros(npi,nbus-1);
? ? for i = 1:npi
? ? ? ? m = fbus(ppi(i));
? ? ? ? for k = 1:(nbus-1)
? ? ? ? ? ? if k+1 == m
? ? ? ? ? ? ? ? for n = 1:nbus
? ? ? ? ? ? ? ? ? ? H21(i,k) = H21(i,k) + V(m)* V(n)*(-G(m,n)*sin(del(m)-del(n)) + B(m,n)*cos(del(m)-del(n)));
? ? ? ? ? ? ? ? end
? ? ? ? ? ? ? ? H21(i,k) = H21(i,k) - V(m)^2*B(m,m);
? ? ? ? ? ? else
? ? ? ? ? ? ? ? H21(i,k) = V(m)* V(k+1)*(G(m,k+1)*sin(del(m)-del(k+1)) - B(m,k+1)*cos(del(m)-del(k+1)));
? ? ? ? ? ? end
? ? ? ? end
? ? end
? ??
? ? % H22 - Derivative of Real Power Injections with V..
? ? H22 = zeros(npi,nbus);
? ? for i = 1:npi
? ? ? ? m = fbus(ppi(i));
? ? ? ? for k = 1:(nbus)
? ? ? ? ? ? if k == m
? ? ? ? ? ? ? ? for n = 1:nbus
? ? ? ? ? ? ? ? ? ? H22(i,k) = H22(i,k) + V(n)*(G(m,n)*cos(del(m)-del(n)) + B(m,n)*sin(del(m)-del(n)));
? ? ? ? ? ? ? ? end
? ? ? ? ? ? ? ? H22(i,k) = H22(i,k) + V(m)*G(m,m);
? ? ? ? ? ? else
? ? ? ? ? ? ? ? H22(i,k) = V(m)*(G(m,k)*cos(del(m)-del(k)) + B(m,k)*sin(del(m)-del(k)));
? ? ? ? ? ? end
? ? ? ? end
? ? end
? ??
? ? % H31 - Derivative of Reactive Power Injections with Angles..
? ? H31 = zeros(nqi,nbus-1);
? ? for i = 1:nqi
? ? ? ? m = fbus(qi(i));
? ? ? ? for k = 1:(nbus-1)
? ? ? ? ? ? if k+1 == m
? ? ? ? ? ? ? ? for n = 1:nbus
? ? ? ? ? ? ? ? ? ? H31(i,k) = H31(i,k) + V(m)* V(n)*(G(m,n)*cos(del(m)-del(n)) + B(m,n)*sin(del(m)-del(n)));
? ? ? ? ? ? ? ? end
? ? ? ? ? ? ? ? H31(i,k) = H31(i,k) - V(m)^2*G(m,m);
? ? ? ? ? ? else
? ? ? ? ? ? ? ? H31(i,k) = V(m)* V(k+1)*(-G(m,k+1)*cos(del(m)-del(k+1)) - B(m,k+1)*sin(del(m)-del(k+1)));
? ? ? ? ? ? end
? ? ? ? end
? ? end
? ??
? ? % H32 - Derivative of Reactive Power Injections with V..
? ? H32 = zeros(nqi,nbus);
? ? for i = 1:nqi
? ? ? ? m = fbus(qi(i));
? ? ? ? for k = 1:(nbus)
? ? ? ? ? ? if k == m
? ? ? ? ? ? ? ? for n = 1:nbus
? ? ? ? ? ? ? ? ? ? H32(i,k) = H32(i,k) + V(n)*(G(m,n)*sin(del(m)-del(n)) - B(m,n)*cos(del(m)-del(n)));
? ? ? ? ? ? ? ? end
? ? ? ? ? ? ? ? H32(i,k) = H32(i,k) - V(m)*B(m,m);
? ? ? ? ? ? else
? ? ? ? ? ? ? ? H32(i,k) = V(m)*(G(m,k)*sin(del(m)-del(k)) - B(m,k)*cos(del(m)-del(k)));
? ? ? ? ? ? end
? ? ? ? end
? ? end
? ??
? ? % H41 - Derivative of Real Power Flows with Angles..
? ? H41 = zeros(npf,nbus-1);
? ? for i = 1:npf
? ? ? ? m = fbus(pf(i));
? ? ? ? n = tbus(pf(i));
? ? ? ? for k = 1:(nbus-1)
? ? ? ? ? ? if k+1 == m
? ? ? ? ? ? ? ? H41(i,k) = V(m)* V(n)*(-G(m,n)*sin(del(m)-del(n)) + B(m,n)*cos(del(m)-del(n)));
? ? ? ? ? ? else if k+1 == n
? ? ? ? ? ? ? ? H41(i,k) = -V(m)* V(n)*(-G(m,n)*sin(del(m)-del(n)) + B(m,n)*cos(del(m)-del(n)));
? ? ? ? ? ? ? ? else
? ? ? ? ? ? ? ? ? ? H41(i,k) = 0;
? ? ? ? ? ? ? ? end
? ? ? ? ? ? end
? ? ? ? end
? ? end
? ??
? ? % H42 - Derivative of Real Power Flows with V..
? ? H42 = zeros(npf,nbus);
? ? for i = 1:npf
? ? ? ? m = fbus(pf(i));
? ? ? ? n = tbus(pf(i));
? ? ? ? for k = 1:nbus
? ? ? ? ? ? if k == m
? ? ? ? ? ? ? ? H42(i,k) = -V(n)*(-G(m,n)*cos(del(m)-del(n)) - B(m,n)*sin(del(m)-del(n))) - 2*G(m,n)*V(m);
? ? ? ? ? ? else if k == n
? ? ? ? ? ? ? ? H42(i,k) = -V(m)*(-G(m,n)*cos(del(m)-del(n)) - B(m,n)*sin(del(m)-del(n)));
? ? ? ? ? ? ? ? else
? ? ? ? ? ? ? ? ? ? H42(i,k) = 0;
? ? ? ? ? ? ? ? end
? ? ? ? ? ? end
? ? ? ? end
? ? end
? ??
? ? % H51 - Derivative of Reactive Power Flows with Angles..
? ? H51 = zeros(nqf,nbus-1);
? ? for i = 1:nqf
? ? ? ? m = fbus(qf(i));
? ? ? ? n = tbus(qf(i));
? ? ? ? for k = 1:(nbus-1)
? ? ? ? ? ? if k+1 == m
? ? ? ? ? ? ? ? H51(i,k) = -V(m)* V(n)*(-G(m,n)*cos(del(m)-del(n)) - B(m,n)*sin(del(m)-del(n)));
? ? ? ? ? ? else if k+1 == n
? ? ? ? ? ? ? ? H51(i,k) = V(m)* V(n)*(-G(m,n)*cos(del(m)-del(n)) - B(m,n)*sin(del(m)-del(n)));
? ? ? ? ? ? ? ? else
? ? ? ? ? ? ? ? ? ? H51(i,k) = 0;
? ? ? ? ? ? ? ? end
? ? ? ? ? ? end
? ? ? ? end
? ? end
? ??
? ? % H52 - Derivative of Reactive Power Flows with V..
? ? H52 = zeros(nqf,nbus);
? ? for i = 1:nqf
? ? ? ? m = fbus(qf(i));
? ? ? ? n = tbus(qf(i));
? ? ? ? for k = 1:nbus
? ? ? ? ? ? if k == m
? ? ? ? ? ? ? ? H52(i,k) = -V(n)*(-G(m,n)*sin(del(m)-del(n)) + B(m,n)*cos(del(m)-del(n))) - 2*V(m)*(-B(m,n)+ bpq(m,n));
? ? ? ? ? ? else if k == n
? ? ? ? ? ? ? ? H52(i,k) = -V(m)*(-G(m,n)*sin(del(m)-del(n)) + B(m,n)*cos(del(m)-del(n)));
? ? ? ? ? ? ? ? else
? ? ? ? ? ? ? ? ? ? H52(i,k) = 0;
? ? ? ? ? ? ? ? end
? ? ? ? ? ? end
? ? ? ? end
? ? end
? ??
? ? % Measurement Jacobian, H..
? ? H = [H11 H12; H21 H22; H31 H32; H41 H42; H51 H52];
? ??
? ? % Gain Matrix, Gm..
? ? Gm = H'*inv(Ri)*H;
? ??
? ? %Objective Function..
? ? J = sum(inv(Ri)*r.^2); ?
? ??
? ? % State Vector..
? ? dE = inv(Gm)*(H'*inv(Ri)*r);
? ? E = E + dE;
? ? del(2:end) = E(1:nbus-1);
? ? V = E(nbus:end);
? ? iter = iter + 1;
? ? tol = max(abs(dE));
end

CvE = diag(inv(H'*inv(Ri)*H)); % Covariance matrix..

Del = 180/pi*del;
E2 = [V Del]; % Bus Voltages and angles..
disp('-------- State Estimation ------------------');
disp('--------------------------');
disp('| Bus | ? ?V ? | ?Angle ?| ');
disp('| No ?| ? pu ? | ?Degree | ');
disp('--------------------------');
for m = 1:num
? ? fprintf('%4g', m); fprintf(' ?%8.4f', V(m)); fprintf(' ? %8.4f', Del(m)); fprintf('\n');
end
disp('---------------------------------------------');

??3?參考文獻(xiàn)

部分理論來(lái)源于網(wǎng)絡(luò)，如有侵權(quán)請(qǐng)聯(lián)系刪除。

[1]A. S. Alayande, N. Nwulu and A. E. Bakare, "Modelling and Countermeasures of False Data Injection Attacks Against State Estimation in Power Systems," 2018 International Conference on Computational Techniques, Electronics and Mechanical Systems (CTEMS), Belgaum, India, 2018, pp. 129-134, doi: 10.1109/CTEMS.2018.8769295.文章來(lái)源地址http://www.zghlxwxcb.cn/news/detail-420794.html

??4 Matlab代碼實(shí)現(xiàn)

到了這里，關(guān)于【狀態(tài)估計(jì)】電力系統(tǒng)狀態(tài)估計(jì)的虛假數(shù)據(jù)注入攻擊建模與對(duì)策（Matlab代碼實(shí)現(xiàn)）的文章就介紹完了。如果您還想了解更多內(nèi)容，請(qǐng)?jiān)谟疑辖撬阉鱐OY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章，希望大家以后多多支持TOY模板網(wǎng)！