趙翠旗，孔會蘭

（晉城供電分公司，山西 晉城 048000）

現(xiàn)代社會中，電能是一種最為廣泛使用的能源，其應(yīng)用程度是一個國家發(fā)展水平的主要標(biāo)志之一。隨著科學(xué)技術(shù)的發(fā)展、微電子器件與電力電子技術(shù)的廣泛應(yīng)用，人們對電能質(zhì)量的要求越來越高；同時由于擾動性負(fù)荷（如非線性、沖擊性或不對稱負(fù)荷）接入電力系統(tǒng)及其他擾動源（系統(tǒng)短路故障）的存在，造成了大量的電能質(zhì)量問題，諧波污染、三相不平衡度、跌落和閃變也越來越嚴(yán)重。電能質(zhì)量問題不僅對電網(wǎng)的安全運(yùn)行不利，嚴(yán)重干擾電網(wǎng)的穩(wěn)定經(jīng)濟(jì)運(yùn)行，還能造成某些對電能質(zhì)量要求較高的電力用戶（如紡織行業(yè)、機(jī)械電子制造業(yè)等）嚴(yán)重的經(jīng)濟(jì)損失。

電能質(zhì)量問題的產(chǎn)生可能來源于供電方的輸配電系統(tǒng)、用戶端的不合理用電及雷電等自然現(xiàn)象。只有對電能質(zhì)量進(jìn)行有效地監(jiān)測，才會對問題的產(chǎn)生和影響有清楚的認(rèn)識，為電能質(zhì)量的改善、供用電雙方的協(xié)調(diào)和供用電市場的規(guī)范提供真實(shí)依據(jù)，以便采取有效的解決措施。為保護(hù)電力系統(tǒng)安全和用戶的用電安全，迫切需要對電能質(zhì)量進(jìn)行監(jiān)測和分析，以提供整改方案、加強(qiáng)防范措施、限制強(qiáng)干擾源（如諧波源），從而保護(hù)電力系統(tǒng)的安全、可靠、經(jīng)濟(jì)運(yùn)行，并保護(hù)電力用戶的合法權(quán)益。

法國學(xué)者M(jìn). Grimaud在地形學(xué)的基礎(chǔ)上提出了一種極值點(diǎn)評價測度——Dyn測度（Dynamics），以解決噪聲環(huán)境下的極值點(diǎn)提取問題。Dyn 測度能夠反映信號極值點(diǎn)的結(jié)構(gòu)特征，并可識別信號的畸變極值點(diǎn)。這一特點(diǎn)使其非常適合電能質(zhì)量擾動檢測。本文在此基礎(chǔ)上，利用MATLAB編寫程序，對電能質(zhì)量一般性周期性擾動進(jìn)行分析，對擾動信號進(jìn)行實(shí)時定位，諧波擾動也可以看做是周期性擾動的特例，對于諧波擾動，更關(guān)心的是它的測量指標(biāo)、諧波畸變率和各次諧波含有率，所以此處不再選擇諧波擾動信號進(jìn)行分析。

本文采用的是MATLAB軟件，其名稱源自Matrix Laboratory，它是由美國MathWorks公司推出的一種科學(xué)計算軟件，專門以矩陣的形式處理數(shù)據(jù)。MATLAB將高性能的數(shù)值計算和可視化集成在一起，并提供了大量的內(nèi)置函數(shù)，從而使其廣泛地應(yīng)用于科學(xué)計算控制系統(tǒng)、信息處理等領(lǐng)域的分析、仿真和設(shè)計工作中。而且利用MATLAB產(chǎn)品的開放式結(jié)構(gòu)，用戶可以容易地對MATLAB的功能進(jìn)行擴(kuò)充，從而在不斷深化對問題認(rèn)識的同時，逐步完善MATLAB產(chǎn)品，以提高產(chǎn)品自身的競爭力。

1 MATLAB編程的優(yōu)點(diǎn)

用MATLAB來解算問題要比用C、FORTRAN等語言簡捷得多。MATLAB是國際公認(rèn)的優(yōu)秀數(shù)學(xué)應(yīng)用軟件之一，其優(yōu)點(diǎn)如下：

（1）以矩陣和數(shù)組為基礎(chǔ)的運(yùn)算，不需要預(yù)先定義變量和矩陣（包括數(shù)組）的維數(shù)，可以方便地進(jìn)行矩陣的算數(shù)運(yùn)算、關(guān)系運(yùn)算和邏輯運(yùn)算等。而且MATLAB有特殊矩陣專門的庫函數(shù)，可以高效地求解諸如信號處理、圖像處理、控制等問題。

（2）語言簡潔，使用方便。MATLAB程序書寫形式自由，函數(shù)名和表達(dá)式更加接近書寫計算公式的表達(dá)式。

（3）強(qiáng)大的科學(xué)計算機(jī)數(shù)據(jù)處理能力。比起 C語言，在計算要求相同的情況下，使用MATLAB編程，工作量會大大減少。MATLAB的這些函數(shù)集包括從最簡單最基本的函數(shù)到諸如矩陣、特征向量、快速傅里葉變換。

（4）強(qiáng)大的圖形處理功能。MATLAB具有強(qiáng)大的以圖形化顯示矩陣和數(shù)組的能力，同時它能給這些圖形增加注釋，對圖形進(jìn)行標(biāo)注和打印。

（5）應(yīng)用廣泛的模塊集合—工具箱。MATLAB的一個重要特色就是具有一套程序擴(kuò)展系統(tǒng)和一組稱之為工具箱的特殊應(yīng)用子程序，每一個工具箱都是為某一類學(xué)科專業(yè)和應(yīng)用而定制的。目前，MATLAB已經(jīng)把工具箱延伸到了科學(xué)研究和工程應(yīng)用的諸多領(lǐng)域，如數(shù)據(jù)采集、數(shù)據(jù)庫接口、概率統(tǒng)計、樣條擬合、優(yōu)化算法、神經(jīng)網(wǎng)絡(luò)、小波分析、電力系統(tǒng)仿真等都在工具箱家族有了一席之地。

（6）可擴(kuò)充性強(qiáng)，具有方便實(shí)用的應(yīng)用程序接口和發(fā)布平臺。MATLAB不僅有著豐富的庫函數(shù)，在進(jìn)行復(fù)雜的數(shù)學(xué)運(yùn)算時可以直接調(diào)用。而且用戶還可以根據(jù)需要，方便地在MATLAB環(huán)境中調(diào)用其他用C語言編寫的代碼，也可以在C語言程序中調(diào)用 MATLAB計算引擎來執(zhí)行 MATLAB代碼。新版本的MATLAB可以利用MATLAB編譯器和C/C＋＋數(shù)據(jù)庫和圖形庫，將自己的MATLAB程序自動轉(zhuǎn)換為獨(dú)立于MATLAB運(yùn)行的 C/C＋＋代碼。允許用戶編寫可以和 MATLAB進(jìn)行交互的C/C＋＋語言代碼。另外，MATLAB網(wǎng)頁服務(wù)程序還容許在Web應(yīng)用中使用自己的數(shù)學(xué)和圖形程序。

（7）源程序的開放性。除內(nèi)部函數(shù)外，所有MATLAB的核心文件和工具箱文件都是可讀可改的源文件，用戶可通過對源文件的修改以及加入自己的文件構(gòu)成新的工具箱。

MATLAB產(chǎn)品族可以用來進(jìn)行以下工作：①數(shù)值分析；②數(shù)值和符號計算；③工程與科學(xué)繪圖；④控制系統(tǒng)的設(shè)計與仿真；⑤數(shù)字圖像處理；⑥數(shù)字信號處理；⑦通訊系統(tǒng)設(shè)計與仿真；⑧財務(wù)與金融工程。

編寫MATLAB程序是在*.m格式文件中進(jìn)行的，其基本操作是：執(zhí)行“File”菜單下的“New”命令，在彈出的對話框中，選擇“M－File”子選項，生成工作環(huán)境，進(jìn)行程序編寫，程序編寫完成，運(yùn)行時點(diǎn)擊工具欄中的“Save and run”命令。

2 簡述MATLAB函數(shù)與指令

在進(jìn)行程序編寫時，會用到不同的MATLAB函數(shù)與指令來實(shí)現(xiàn)程序的功能，這里僅對幾個比較重要的函數(shù)和指令進(jìn)行描述。

2.1 數(shù)據(jù)載入指令

將 txt文檔中的數(shù)據(jù)載入 MATLAB中可以用的命令有很多，本文程序中用到的是load命令和importdata命令。

2.2 畫圖指令

在MATLAB中進(jìn)行數(shù)據(jù)可視化，使用最頻繁的繪制函數(shù)就是plot函數(shù)，該函數(shù)能夠?qū)⑾蛄炕蛘呔仃囍械臄?shù)據(jù)繪制在圖形窗體中，并且可以指定不同的線型和色彩。同一個plot函數(shù)不僅能繪制一條曲線，還可以一次繪制多條曲線。

plot函數(shù)的基本使用格式為：

繪制一條曲線：plot(xdata,ydata,'color_linestyle_marker')；繪制多條曲線：plot(xdata1,ydata2,'clm1',xdata2,ydata2,'clm2',......)；若在繪制曲線時沒有指定曲線的色彩、線型和標(biāo)識符，則MATLAB使用默認(rèn)的設(shè)置。

2.3 選擇結(jié)構(gòu)

當(dāng)判斷某一條件是否滿足，根據(jù)判斷的結(jié)果來選擇不同的解決問題的方法時，就需要使用選擇結(jié)構(gòu)。和C語言類似，MATLAB的條件判斷可以使用if語句或者switch語句。本文采用的是if語句。

2.4 for循環(huán)結(jié)構(gòu)

在解決很多問題時需要使用循環(huán)結(jié)構(gòu)，MATLAB中包含兩種循環(huán)結(jié)構(gòu)：一種是循環(huán)次數(shù)不確定的 while循環(huán)；另一種是循環(huán)次數(shù)確定的for循環(huán)。本文中主要用到的是for循環(huán)結(jié)構(gòu)，使用for語句構(gòu)成循環(huán)結(jié)構(gòu)是最靈活、簡便的方法，在MATLAB中，for循環(huán)的基本結(jié)構(gòu)如下：

for index=start:increment:end

MATLAB語句end

其中，index的取值取決于start和end的值。MATLAB的循環(huán)結(jié)構(gòu)也可以進(jìn)行嵌套使用，使用嵌套的循環(huán)需要注意 for關(guān)鍵字和end關(guān)鍵字之間的配對使用。

3 Dyn測度

在了解Dyn測度定義之前，首先要了解兩點(diǎn)之間路徑的概念。設(shè)f（t）為定義在有限區(qū)間I上的信號，且滿足條件f（t）≥0，可以將各采樣點(diǎn)的數(shù)值看作是各點(diǎn)的高度。兩點(diǎn)之間的路徑：設(shè)m，n為f上不同的兩點(diǎn)，則f上這兩點(diǎn)之間的部分稱為路徑P（m，n），其中，P1為m，PN為n。

（1）路徑的Dyn測度：路徑P（m，n）的Dyn測度定義為路徑上最高點(diǎn)和最低點(diǎn)的高度差。

（2）極小點(diǎn)的Dyn測度：設(shè)M為f的一個極小點(diǎn)，如果存在比其更低的極小點(diǎn)，則極小點(diǎn)M的Dyn測度等于由點(diǎn)M通向同高度點(diǎn)的所有路徑中最小的路徑Dyn測度。極小點(diǎn)M兩側(cè)各有一個或多個比點(diǎn)M更低的極小點(diǎn)時，點(diǎn)M兩側(cè)也一定存在兩個與點(diǎn)M等高度的點(diǎn)N1和N2。由極小點(diǎn)M通向點(diǎn)N1的路徑記為P1，由極小點(diǎn)M通向點(diǎn)N2的路徑記為P2。極小點(diǎn)M的Dyn測度等于路徑P1的Dyn測度和路徑P2的Dyn測度中較小者。路徑P1的Dyn測度小于路徑P2的Dyn測度，所以極小點(diǎn)M的Dyn測度應(yīng)該等于路徑P1的Dyn測度。即：

值得注意的是，極小點(diǎn)M的Dyn測度與路徑或路徑P2的長度無關(guān)。極小點(diǎn)M的一側(cè)存在比點(diǎn)M更低的極小點(diǎn)時，意味著只存在路徑P1或路徑P2，則極小點(diǎn)M的Dyn測度等于路徑P1或路徑P2的Dyn測度。極小點(diǎn)M為信號的最小點(diǎn)時，設(shè)置其Dyn測度等于信號最高點(diǎn)和最低點(diǎn)的高度差。這樣，最低點(diǎn)的 Dyn測度比其它極小點(diǎn)的 Dyn測度都大。極大點(diǎn)的 Dyn測度同極小點(diǎn)Dyn測度的定義相仿。

Dyn測度算法能夠快速檢測信號的所有極值點(diǎn)，并計算各極值點(diǎn)的Dyn測度，其基于溢流仿真的過程描述如下：

設(shè)信號f（t）在定義區(qū)間I＝[a，b]上含有多個極小值和極大值，hmin和 hmax分別為信號的最小值和最大值。某時刻的溢流水平為 h，對應(yīng)的溢流區(qū)域 C（h）定義為滿足 f（t）≤h的所有點(diǎn)t的集合，即：

隨著溢流水平的遞增，溢流區(qū)域不斷擴(kuò)大，直到溢流水平達(dá)到h＝hmax時，溢流過程結(jié)束，此時的溢流區(qū)域為信號的定義區(qū)間，即：

溢流過程是一個迭代過程。迭代步長為dh，假設(shè)第k步的溢流水平為hk，對應(yīng)的溢流區(qū)域為C（hk）；第k＋l步的溢流水平為 hk＋1＝hk＋dh，對應(yīng)的溢流區(qū)域為 C（hk＋1），則：

如果D為C（hk＋1）的一個連通子集，D與C（hk）可能有下面3種關(guān)系：

（1）D∩C（hk）為空集：這種情況下，D是f（t）的一個新的極小點(diǎn)，該新的極小點(diǎn)記為Mj，其高度為hk＋1。

（2）D∩C（hk）不為空集，并且為一個以上的連通成分：這是沒有出現(xiàn)新的極小點(diǎn)，也沒有出現(xiàn)新的極大點(diǎn)。

（3）D∩C（hk）不為空集，并且為一個以上的連通成分：這種情況下，出現(xiàn)了新的極大點(diǎn)，且該極大點(diǎn)的高度為 hk＋1。該極大點(diǎn)兩側(cè)的兩個極小點(diǎn)中，較高極小點(diǎn)Mj的Dyn測度為：

其中，alt（Mj）表示Mj點(diǎn)的高度。

迭代過程結(jié)束時，除最低極小點(diǎn)M外，其他極小點(diǎn)的Dyn測度全部得到。而最低極小點(diǎn)M的測度為：

如圖1所示，信號發(fā)生了周期性畸變，信號的Dyn測度檢測結(jié)果如圖1（b）所示，不同位置出現(xiàn)了幅值不等的Dyn測度，但其分布有一定的規(guī)律性：一部分Dyn測度的幅值較大，保持在2.5左右，且均勻分布在每個周期；另一部分Dyn測度的幅值較小，也均勻分布在每個周期。比較兩圖可知，前者幅值較大的Dyn測度對應(yīng)信號的峰（或谷）點(diǎn)；而幅值較小的Dyn測度對應(yīng)信號的畸變點(diǎn)。

圖1 周期性擾動信號檢測結(jié)果

圖2 調(diào)節(jié)閾值后的周期性擾動的檢測結(jié)果

取閾值為信號峰/谷點(diǎn)的Dyn測度的15%，得到如圖2所示的檢測結(jié)果。比較圖1和圖2可知，通過選擇合適的閾度，可以消除峰谷點(diǎn)，保留畸變點(diǎn)，達(dá)到實(shí)時檢測周期性擾動信號畸變點(diǎn)的目的。

通過測試可以知道，基于Dyn測度的檢測方法能有效地檢測多種電能質(zhì)量擾動（暫態(tài)擾動、電壓暫降和周期性擾動）；而且算法簡單、計算量小，其應(yīng)用和推廣對電能質(zhì)量監(jiān)測有重要意義。

1 李瑞峰.10 kV電網(wǎng)電能質(zhì)量在線監(jiān)測系統(tǒng)的研究與應(yīng)用[D].華南理工大學(xué)，2006

2 袁之泉.電能質(zhì)量分析系統(tǒng)研究[D].東北大學(xué)，2005

3 宋曉芳.電能質(zhì)量分析技術(shù)研究與實(shí)現(xiàn)[D].南京理工大學(xué)，2005

4 劉應(yīng)梅.電能質(zhì)量擾動檢測和分析的研究[D].中國電力科學(xué)研究院，2003

5 曹暉.基于DSP的電能質(zhì)量分析系統(tǒng)的開發(fā)[D].華中科技大學(xué)，2004

6 許康平.電能質(zhì)量監(jiān)測裝置的研究[D].東南大學(xué)，2007