的得到校正值為相應(yīng)的直流分量、各交流的 正弦分量、余弦分量和估計(jì)頻率的值，即
[0107] S106、依據(jù)校正處理后的直流分量UQ、正弦分量SUl、SU2、…、SUN、余弦分量CU1、 CU2、…、CUN、估計(jì)頻率f 1、f 2、…、f N的值，令j = 1，e = Ο，利用式⑴更新直流分量校正值增量 X0[l]，令下標(biāo)η分別取值為1，2，···，Ν，循環(huán)執(zhí)行式(2)，分別更新正弦分量校正值增量&[1]、 余弦分量校正值增量y n[ 1 ]、估計(jì)頻率校正值增量Ζη[ 1 ];
[0108] S107、利用式（7)獲得存檔用直流分量校正值增量f〇[m]，令下標(biāo)η分別取值為1， 2，···，Ν，循環(huán)執(zhí)行式(8)，分別獲得存檔用正弦分量校正值增量sufn[m]、余弦分量校正值增 量cufn[m]、估計(jì)頻率校正值增量ω fn[m];
[0111]其中，m從1開(kāi)始，在前8個(gè)采樣周期時(shí)間內(nèi)的采樣時(shí)刻2T，4T，6T，8T時(shí)分別獲得1組 存檔用校正值增量數(shù)據(jù)，m的值依次加1，分別為1、2、3、4,共獲得4組校正值增量數(shù)據(jù)和8Τ時(shí) 刻校正處理后獲得的直流分量UQ、正弦分量SUl、SU2、…、SUN、余弦分量CUl、CU2、…、CUN的值 作為算法第二階段的初始值；
[0? 12] S108、輸出數(shù)據(jù):依據(jù)校正處理后的直流分量UQ、正弦分量SU1、SU2、···、SUN、余弦分 量cm、cu2、…、CUN、估計(jì)頻率fhfy'fN的值，利用式(9)將直流分量輸出保存至u〇[m]、令 下標(biāo)η分別取值為1，2，···，Ν，循環(huán)執(zhí)行式(10)，
[0115] 將各交流成分的諧波或間諧波正弦分量、諧波或間諧波余弦分量、估計(jì)頻率、估計(jì) 幅值和相角依次輸出保存至sun[m]、cu n[m]、fn[m]、dn[m]和θη[ηι]。
[0116] 算法的第一階段完成后，接著執(zhí)行第二階段，此時(shí)，讀取離散時(shí)間序列中的第k(k >8)個(gè)數(shù)據(jù)，作為當(dāng)前處理的采樣數(shù)據(jù)us，依次執(zhí)行以下步驟：
[0117] S201、依據(jù)校正處理后獲得的直流分量uo,直流分量校正值增量fQ[l]、fQ[2]、f 0 [3]、f〇[4]，利用式（11 )獲得直流分量預(yù)估值.，.依據(jù)正弦分量SU1、SU2、…、SUN和正弦分量 校正值增量sufi[l]、sufi[2]、sufi[3]、sufi[4]，…，sufN[l ]、sufN[2]、sufN[3]、sufN[4]，余 弦分量CU1、CU2、…、CUN和余弦分量校正值增量cufi[l]、cufi[2]、cufi[3]、cufi[4]，…，cufN [1 ]、cufN[2]、cufN[3]、cufN[4]，估計(jì)頻率fi、f2、…、fN，估計(jì)頻率校正值增量 ω 負(fù)[1 ]、ω [2]、ωΜ3]、ωΜ4]，···，ω?·Ν[1]、ω?·Ν[2]、ω?·Ν[3]、ω?· Ν[4]，令n分別取值為 1，2，···，Ν，循 環(huán)執(zhí)行式(12)，獲得正弦分量預(yù)估值?丨、余弦分量預(yù)估值〇<、估計(jì)頻率預(yù)估值;

[0120] S202、依據(jù)直流分量預(yù)估值喊，利用式（13)獲得直流分量校正值增量#，運(yùn)用正 弦分量預(yù)估值<、《4、…、余弦分量預(yù)估值GMf、_f、·，·、βΜ:|，.估計(jì)頻率預(yù)估 值/Λ、…7,令η分別取值為1，2，"_，1循環(huán)執(zhí)行式(14)，
[0123] 獲得正弦分量校正值增量?jī)?cè)/Γ、余弦分量校正值增量⑶乂"、估計(jì)頻率校正值增量 (〇f：；
[0124] S203、依據(jù)校正處理后獲得的直流分量uo、直流分量校正值增量f〇[ 2 ]、f〇[ 3 ]、fo [4 ]、#，依據(jù)式（15)對(duì)直流分量UQ進(jìn)行校正處理，獲得直流分量校正值:;
[0125] 正弦分量sui、SU2、…、sun和正弦分量校正值增量sufi[2]、sufi[3]、sufi[4]，…， sufN[2]、sufN[3]、sufN[4]、《(//，余弦分量CUl、CU2、···、CUN和余弦分量校正值增量cufl[2]、 cufl[3]、cufi[4]，···，cufN[2]、cufN[3]、cufN[4]、《(//，估計(jì)頻率 fi、f2、…、fN，估計(jì)頻率校正 值增量 ω M2]、ω f\[3]、ω f\[4]，…，ω M2]、ω M3]、ω M4]、?紀(jì)、~，?#，令 下標(biāo)η分別取值為1，2, ···，Ν，循環(huán)執(zhí)行式（16)，先對(duì)正弦分量sun、余弦分量cun、估計(jì)頻率fn 進(jìn)行校正處理；分別獲得正弦分量校正值《<，余弦分量校正值和估計(jì)頻率校正值# 并對(duì)校正后的估計(jì)頻率./：進(jìn)行限幅處理以滿足條件《y'/丨 ;
[0128] S204、按照式（17)和（18)，令步驟S203中的得到校正值為相應(yīng)的直流分量、各交流 的正弦分量、余弦分量和估計(jì)頻率的值，即
[0131] S205、更新校正值增量，取i分別取值為2，3，4執(zhí)行式（19)更新直流分量校正值增 量&[1]、&[2]、&[3];當(dāng)1分別取值為2,3,4時(shí)，令下標(biāo)11分別取值為1，2，一，1循環(huán)執(zhí)行式 (20)更新正弦分量校正值增量sufi[l]、sufi[2]、sufi[3]，···，sufN[l]、sufN[2]、sufN[3]，以 及余弦分量校正值增量cufi[l]、cufi[2]、cufi[3]，···，CUfN[l]、CUfN[2]、CUfN[3]，以及估計(jì) 頻率校正值增量ωΜΙΚωΜ〗]、ωΜ3]，···，ω?· Ν[1]、ωΜ2]、cofN[3];
[0134]按式(21)更新直流分量校正值增量fo[4]，令下標(biāo)η分別取值為1，2，···，Ν，循環(huán)執(zhí) 行式（22)更新正弦分量校正值增量sufi[4]、suf2[4]，···，sufN[4]，余弦分量校正值增量 cufi[4]、cuf2[4]，…，cufN[4]和估計(jì)頻率校正值增量 ω ??]、ω f2[4]，…，ω fN[4];
[0137] S206、輸出數(shù)據(jù):依據(jù)校正處理后的直流分量UQ、正弦分量SU1、SU2、···、sun、余弦分 量cm、CU2、…、CUN、估計(jì)頻率f i、f2、…、fN的值，利用式(23)將直流分量輸出保存至uo[m]、令 下標(biāo)η分別取值為1，2，···，Ν，循環(huán)執(zhí)行式（24)，將各交流成分的電壓諧波或間諧波正弦分 量、電壓諧波或間諧波余弦分量、估計(jì)頻率、估計(jì)幅值和相角依次輸出保存至su n[m]、cun [m]、fn[m]、dn[m]和9n[m];
[0140]算法第一階段后，m的值為4,此處m從5開(kāi)始，每次輸出保存后加1。
[0141] 在執(zhí)行第一階段和第二階段的相關(guān)步驟時(shí)，依據(jù)預(yù)估一校正處理后的直流分量 U0、正弦分量SU1、SU2、…、SUN、余弦分量CU1、CU2、…、CUN，計(jì)算諧波或間諧波的估計(jì)頻率、估 計(jì)幅值和相角。為分析算法結(jié)果，將直流分量、正弦分量、余弦分量、諧波或間諧波的估計(jì)頻 率、估計(jì)幅值和相角保存到輸出系列的第m個(gè)位置，然后把位置號(hào)m增加1。
[0142] 在循環(huán)控制步驟中，先把位置號(hào)k增加1，再依據(jù)位置號(hào)k和離散時(shí)間序列中數(shù)據(jù)的 組數(shù)K的值判斷是否返回執(zhí)行數(shù)據(jù)預(yù)估一校正處理步驟。若k < K，返回執(zhí)行數(shù)據(jù)預(yù)估一校正 處理步驟;若k>K，表示離散時(shí)間序列中所有采樣數(shù)據(jù)都已處理完畢，應(yīng)該終止運(yùn)行，結(jié)束 離線分析過(guò)程。
[0143] 實(shí)施例2
[0144] 如圖2和圖3所示，為本發(fā)明實(shí)施例2的基于Adams算法的頻率自適應(yīng)諧波間諧波分 析方法的流程示意圖，實(shí)施例2的基于Adams算法的頻率自適應(yīng)諧波間諧波分析方法為在線 分析實(shí)施方式，特點(diǎn)是一邊采樣一邊分析計(jì)算，即每采樣一組數(shù)據(jù)，就進(jìn)行一次分析計(jì)算； 適用于需要對(duì)被測(cè)信號(hào)的每個(gè)采樣數(shù)據(jù)都進(jìn)行實(shí)時(shí)分析處理的情況。
[0145] 實(shí)施例2的基于Adams算法的頻率自適應(yīng)諧波間諧波分析方法包括主程序和定時(shí) 中斷服務(wù)程序。
[0146] 圖2所示的主程序包括初始化步驟和主循環(huán)步驟。在初始化步驟中，首先設(shè)定參數(shù) Τ、Ν、αι、α2、···、CX2N-1、Cl2N、V、r、d的值;設(shè)定直流分量UQ、正弦分量SUl、SU2、…、SUN、余弦分量 CU1、CU2、…、CUN、估計(jì)頻率fl、f2、…、ft的初值;接著設(shè)定直流分量校正值增量XQ [ 1 ]，正弦分 量校正值增量11[1]、12[1]、'"、1〃[1]、余弦分量校正值增量71[1]、72[1]、'"、5^[1]、估計(jì)頻 率校正值增量 Z1[1]、Z2[1]、…、ZN[1]的初值;設(shè)定當(dāng)前要采集的數(shù)據(jù)在離散時(shí)間序列中的 位置號(hào)k為1，設(shè)定當(dāng)前校正處理后的數(shù)據(jù)在輸出序列中的位置號(hào)m為1;接著再設(shè)定定時(shí)器 的定時(shí)時(shí)間為T(mén)秒，T為滿足香農(nóng)采樣定理要求的采樣周期，并開(kāi)放系統(tǒng)的定時(shí)中斷。
[0147] 圖2所示的定時(shí)中斷服務(wù)步驟中，首先保存中斷現(xiàn)場(chǎng)的各個(gè)寄存器的當(dāng)前值;接著 按采樣周期T對(duì)被測(cè)信號(hào)進(jìn)行采樣，獲得采樣數(shù)據(jù)us，如果采集的數(shù)據(jù)的序號(hào)k < 8，執(zhí)行算 法第一階段計(jì)算校正值增量數(shù)值并進(jìn)行校正處理;如果k的值大于8執(zhí)行算法第二階段計(jì)算 預(yù)估值、校正值增量并進(jìn)行校正處理;也就是算法第一階段和算法第二階段，執(zhí)行除數(shù)據(jù)輸 出（即第一階段除S108,第二階段除S206)外的所有步驟。再恢復(fù)中斷現(xiàn)場(chǎng)的各個(gè)寄存器的 值;然后中斷返回到主步驟執(zhí)行。
[0148] 主步驟循環(huán)讀取中斷返回時(shí)經(jīng)校正處理的數(shù)據(jù)，執(zhí)行數(shù)據(jù)計(jì)算和輸出（即執(zhí)行第 一階段S108或第二階段S206,這兩步功能相同）、數(shù)據(jù)分析，以及關(guān)中斷、程序結(jié)束等程序控 制功能。
[0149] 在線分析實(shí)施方式通過(guò)定時(shí)器循環(huán)產(chǎn)生定時(shí)中