[發明專利]擾動源追溯方法

有效

• 申請號: 201611032532.9
• 申請日: 2016-11-15
• 公開（公告）號: CN107871197B
• 公開（公告）日: 2022-02-01
• 發明（設計）人: 李俊賢;鄭儀誠;徐振凱;陳俊彥
• 申請（專利權）人: 財團法人工業技術研究院
• 主分類號: G06Q10/06
• 分類號: G06Q10/06;G06F17/16
• 代理公司: 中科專利商標代理有限責任公司 11021
• 代理人: 任巖

1.一種擾動源追溯方法，適于由計算裝置從包括多個候選節點的系統中追溯擾動源，該方法包括下列步驟：

對該系統的多個節點進行一頻譜特征分析；

依據該頻譜特征分析的一分析結果，篩選該頻譜特征相近的所述節點作為所述候選節點，其中所述候選節點之間的該頻譜特征的差值小于預設值；

獲得由多種流程分析方法，分別分析該系統所得的拓撲信息，該拓撲信息記錄該系統中所述候選節點之間的因果關系；

將各所述拓撲信息中的所述因果關系導入機率分布演算法，計算各所述候選節點的靜態機率分布；

融合針對各所述流程分析方法所計算的所述靜態機率分布，計算各所述候選節點為該擾動源的機率分布；

排序該機率分布所對應的各所述候選節點；以及

判定排序最高的該候選節點為該系統的該擾動源；

其中將各所述拓撲信息中的所述因果關系導入該機率分布演算法，計算各所述候選節點的該靜態機率分布的步驟包括：

建立記錄所述候選節點之間的所述因果關系的鄰接矩陣；

正規化該鄰接矩陣為轉移矩陣；以及

以迭代方式計算該轉移矩陣在穩定狀態的數值解，以作為各所述候選節點的該靜態機率分布。

2.如權利要求1所述的擾動源追溯方法，其中以迭代方式計算該轉移矩陣在該穩定狀態的該數值解，以作為各所述候選節點的該靜態機率分布的步驟包括：

計算各所述轉移矩陣在特征值為1時的一特征向量，以作為所述候選節點的該靜態機率分布。

3.如權利要求1所述的擾動源追溯方法，其中融合針對各所述候選節點所計算的所述靜態機率分布，計算各所述候選節點為該擾動源的機率分布的步驟包括：

計算針對各所述候選節點所計算的所述靜態機率分布的總和，以作為該候選節點為該擾動源的機率分布。

4.如權利要求1所述的擾動源追溯方法，其中所述流程分析方法包括數據驅動方法(data-driven method)及基于模型的方法(model-based method)。

5.如權利要求4所述的擾動源追溯方法，其中該數據驅動方法包括格蘭杰因果關系檢驗法(Granger causality test)，并且該基于模型的方法包括管道儀表流程圖法(Pipingand Instrument Diagram，PID)。

6.如權利要求4所述的擾動源追溯方法，其中該數據驅動方法包括格蘭杰因果關系檢驗法、熵轉移法(Transfer entropy)、貝氏網路法(Bayesian network)以及互相關法(Cross-correlation)的至少其中之一。

7.如權利要求4所述的擾動源追溯方法，其中該機率分布演算法基于馬可夫鏈(Markovchain)。