【正文】 (num)idb(num) = samib1(num) samib2(num)irb(num) = samib1(num) + samib2(num)idc(num) = samic1(num) samic2(num)irc(num) = samic1(num) + samic2(num)samia1(num + 48) = samia1(num) 39。采用暫存區(qū)擴展法保存采樣濾波數(shù)據(jù)samib1(num + 48) = samib1(num)samic1(num + 48) = samic1(num)samia2(num + 48) = samia2(num)samib2(num + 48) = samib2(num)samic2(num + 48) = samic2(num)ida(num + 48) = ida(num)ira(num + 48) = ira(num)idb(num + 48) = idb(num)irb(num + 48) = irb(num)idc(num + 48) = idc(num)irc(num + 48) = irc(num)If WAIT = 1 Then GoTo FLG 39。若等待標志為1，則退出定時中斷服務(wù)過程If cktype 0 Then GoTo tzck 39。有出口信號則直接轉(zhuǎn)跳閘出口If QDB = 1 Then GoTo guzhangchuli 39。若啟動標志為1，轉(zhuǎn)入故障處理子程序If qdyb = 1 Then 39。若相電流突變量啟動壓板投入則 J = num + 48 X = Abs(samia1(J) samia1(J 12)) 39。相電流突變量啟動判斷DI Y = Abs(samia1(J 12) samia1(J 24)) dia = Abs(X Y) X = Abs(samib1(J) samib1(J 12)) Y = Abs(samib1(J 12) samib1(J 24)) dib = Abs(X Y) X = Abs(samic1(J) samic1(J 12)) Y = Abs(samic1(J 12) samic1(J 24)) dic = Abs(X Y) If dia diqd Then di1tia = di1tia + 1 Else di1tia = 0 If dib diqd Then di1tib = di1tib + 1 Else di1tib = 0 If dic diqd Then di1tic = di1tic + 1 Else di1tic = 0 If di1tia = 4 Or di1tib = 4 Or di1tic = 4 Then 39。DI啟動則 QDB = 1 39。啟動標志置1 st = 0 39。故障啟動后計時指針清0 End IfElse QDB = 1 st = 0End IfGoTo FLGguzhangchuli: 39。故障處理子程序 = QBColor(12) 39。點亮保護啟動燈If st = 48 Then 39。故障啟動后48點以前更新故障后計時指針 st = st + 1End IfIf yb = 1 Then 39。差動保護壓板投入則進行差動判斷 GoSub idm_irm If idma idz_max Or idmb idz_max Or idmc idz_max Then 39。滿足差動速斷條件則差動速斷動作 rept1$ = 差動速斷動作 cktype = 3: tzlreturn = 0: tzlout = 0: tzhtime = 0: GoTo tzck End If If iam2 / iam1 k2 And ibm2 / ibm1 k2 And icm2 / icm1 k2 Then 39。無二次諧波制動則進入比率差動 If idma idz_min + (irma ir_min) * k1 Then rept1$ = 比率差動動作 cktype = 3: tzlreturn = 0: tzlout = 0: tzhtime = 0: GoTo tzck End If If idmb idz_min + (irmb ir_min) * k1 Then rept1$ = 比率差動動作 cktype = 3: tzlreturn = 0: tzlout = 0: tzhtime = 0: GoTo tzck End If If idmc idz_min + (irmc ir_min) * k1 Then rept1$ = 比率差動動作 cktype = 3: tzlreturn = 0: tzlout = 0: tzhtime = 0: GoTo tzck End If End IfEnd IfIf wzbdy = 1 Then 39。若斷路器位置不對應(yīng)(偷跳造成),則轉(zhuǎn)三跳后出口(跳閘令發(fā)出標志tzlout設(shè)為1） cktype = 3: tzlreturn = 0: tzlout = 1: tzhtime = 0: GoTo FLGElse 39。否則(說明故障已由非本線路保護切除),則整組復歸并返回 GoTo allreturnEnd Iftzck: 39。跳閘出口子程序REP = 1 39。有報告標志置1If cktype = 3 Then 39。若出口類型為3則執(zhí)行跳閘 If tzlreturn = 0 And tzlout = 0 Then 39。發(fā)三跳令 If = 正常 Then = 分 = QBColor(10) = 分 = QBColor(10) = 分 = QBColor(10) = 分 = QBColor(10) = 分 = QBColor(10) = 分 = QBColor(10) End If tzlout = 1 39。跳令發(fā)出標志置1 End If If tzhtime 3000 Then GoTo allreturn 39。延時5s到(3000*Ts=3000*5/3ms=5000ms)則整組復歸 If tzhtime 24 Then 39。延時40ms未到(24*Ts=24*5/3ms=40ms)則中斷返回 If = 正常 Then break = 1 39。已跳閘標志設(shè)為1（仿真跳閘后的無流過程） tzhtime = tzhtime + 1: GoTo FLG End If If tzlreturn = 1 Then tzhtime = tzhtime + 1: GoTo FLG 39。若跳閘令已收回則中斷返回 GoSub nocurrent If noi = 1 Then tzlreturn = 1: tzhtime = tzhtime + 1: GoTo FLG 39。若無流則收回跳閘令后中斷返回 If tzhtime 150 Then tzhtime = tzhtime + 1: GoTo FLG 39。延時250ms未到(150*5/3ms=250ms)則中斷返回 If tzagain = 1 Then rept4$ = 跳閘失敗: GoTo FLG 39。重跳標志為1則報告跳閘失敗并返回 tzagain = 1 39。重跳標志置1 cktype = 2: tzlreturn = 0: tzlout = 0: tzhtime = 0 39。重新永跳一次End IfGoTo FLGnocurrent: 39。無流判別子模塊(noi為無流標志，1無流，2有電流）noi = 2iw = 2 * iwl ^ 2If idam iw And idbm iw And idcm iw Then 39。A、B、C三相電流皆定值，則判為無流 noi = 1End IfReturnidm_irm: 39。電流幅值平方計算子模塊(idam,idbm,idcm,iram,irbm,ircm) pi = J = num + 48 M = 0: n = 0 39。富氏算法求ida幅值 For I = J 11 To J M = M + ida(I) * Cos(I * pi / 6) / 6 n = n + ida(I) * Sin(I * pi / 6) / 6 Next I idam = Sqr(M ^ 2 + n ^ 2) M = 0: n = 0 39。富氏算法求idb幅值 For I = J 11 To J M = M + idb(I) * Cos(I * pi / 6) / 6 n = n + idb(I) * Sin(I * pi / 6) / 6 Next I idbm = Sqr(M ^ 2 + n ^ 2) M = 0: n = 0 39。富氏算法求idc幅值 For I = J 11 To J M = M + idc(I) * Cos(I * pi / 6) / 6 n = n + idc(I) * Sin(I * pi / 6) / 6 Next I idcm = Sqr(M ^ 2 + n ^ 2) M = 0: n = 0 39。富氏算法求ira幅值 For I = J 11 To J M = M + ira(I) * Cos(I * pi / 6) / 6 n = n + ira(I) * Sin(I * pi / 6) / 6 Next I iram = Sqr(M ^ 2 + n ^ 2) M = 0: n = 0 39。富氏算法求irb幅值 For I = J 11 To J M = M + irb(I) * Cos(I * pi / 6) / 6 n = n + irb(I) * Sin(I * pi / 6) / 6 Next I irbm = Sqr(M ^ 2 + n ^ 2) M = 0: n = 0 39。富氏算法求irc幅值 For I = J 11 To J M = M + irc(I) * Cos(I * pi / 6) / 6 n = n + irc(I) * Sin(I * pi / 6) / 6 Next I ircm = Sqr(M ^ 2 + n ^ 2)Returnallreturn: 39。整組復歸QDB = 0: REP = 0: WAIT = 0: YUNXING = 0cktype = 0: tzlreturn = 0: tzlout = 0: tzagain = 0di1tia = 0: di1tib = 0: di1tic = 0st = 0: tzhtime = 0: break = 0rept1$ = = FLG:For I = 0 To 9 If numb 48 * I And numb = 6 + 48 * I Then = QBColor(7) 39。熄滅 ElseIf numb 6 + 48 * I And numb = 12 + 48 * I Then = QBColor(13) 39。點亮 ElseIf numb 12 + 48 * I And numb = 18 + 48 * I Then = QBColor(7) 39。熄滅 ElseIf numb 18 + 48 * I And numb = 24 + 48 * I Then = QBColor(13) 39。點亮 ElseIf numb 24 + 48 * I And numb = 30 + 48 * I Then = QBColor(7) 39。熄滅 ElseIf numb 30 + 48 * I And numb = 36 + 48 * I Then = QBColor(13) 39。點亮 ElseIf numb 36 + 48 * I And numb = 42 + 48 * I Then = QBColor(7) 39。熄滅 ElseIf numb 42 + 48 * I And numb = 48 + 48 * I Then = QBColor(13) 39。點亮 End IfNext IIf numb = 8