WO2022019020A1 - 異常検知装置及び異常検知プログラム - Google Patents

異常検知装置及び異常検知プログラム Download PDF

Info

Publication number
WO2022019020A1
WO2022019020A1 PCT/JP2021/023245 JP2021023245W WO2022019020A1 WO 2022019020 A1 WO2022019020 A1 WO 2022019020A1 JP 2021023245 W JP2021023245 W JP 2021023245W WO 2022019020 A1 WO2022019020 A1 WO 2022019020A1
Authority
WO
WIPO (PCT)
Prior art keywords
time
abnormality
time window
detection device
series data
Prior art date
Application number
PCT/JP2021/023245
Other languages
English (en)
French (fr)
Inventor
結人 久保
響輔 竹腰
英司 大島
Original Assignee
本田技研工業株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 本田技研工業株式会社 filed Critical 本田技研工業株式会社
Priority to CN202180049128.5A priority Critical patent/CN115989465A/zh
Priority to US18/005,897 priority patent/US20230322240A1/en
Priority to MX2023000909A priority patent/MX2023000909A/es
Publication of WO2022019020A1 publication Critical patent/WO2022019020A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • G06N3/0455Auto-encoder networks; Encoder-decoder networks
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W50/00Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
    • B60W50/02Ensuring safety in case of control system failures, e.g. by diagnosing, circumventing or fixing failures
    • B60W50/0205Diagnosing or detecting failures; Failure detection models
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M17/00Testing of vehicles
    • G01M17/007Wheeled or endless-tracked vehicles
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • G06N20/10Machine learning using kernel methods, e.g. support vector machines [SVM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W50/00Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
    • B60W2050/0001Details of the control system
    • B60W2050/0043Signal treatments, identification of variables or parameters, parameter estimation or state estimation
    • B60W2050/005Sampling
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W50/00Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
    • B60W50/02Ensuring safety in case of control system failures, e.g. by diagnosing, circumventing or fixing failures
    • B60W50/0205Diagnosing or detecting failures; Failure detection models
    • B60W2050/021Means for detecting failure or malfunction
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W2556/00Input parameters relating to data
    • B60W2556/10Historical data

Definitions

  • the present invention relates to a technique for detecting an abnormality in a vehicle signal.
  • Patent Documents 1 and 2 reference data indicating normality is collected in advance using another vehicle at the time of the durability test of the vehicle, and the data measured during the durability test using the reference data. The technology to detect anomalies is described.
  • the present invention has been made in view of the above points, and provides an anomaly detection device and an anomaly detection program capable of performing anomaly detection while suppressing the effects of individual differences in vehicles and wear during testing. Make it an issue.
  • the abnormality detection device of the present invention includes a time window processing unit that extracts data included in the time window by dividing a plurality of vehicle signals including time series data by a time window.
  • the learning unit that generates learned information by learning using the time-series data of the vehicle signal separated by the time window, and the time-series data of the vehicle signal separated by the time window are used as the learned information. It is characterized in that it is reconstructed based on the above, the error before and after the reconstruction is calculated, and the abnormality degree calculation unit for calculating the abnormality degree based on the calculated error and the learned information is provided.
  • the abnormality detection program of the present invention is characterized in that the computer functions as the abnormality detection device.
  • the present invention when detecting an abnormality in a vehicle, it is possible to suppress the influence of individual differences in the vehicle and wear during the test.
  • the abnormality detection system 1 is a system for detecting an abnormality in the vehicle 2, and includes a measuring device 10, an abnormality detecting device 20, a notification device 30, and the like. To prepare for.
  • the measuring device 10 measures various signals (vehicle signals) of the vehicle 2 during a simulated test run in the tabletop durability test, and transmits a plurality of measured vehicle signals to the abnormality detecting device 20.
  • the measured vehicle signal includes signals of different orders and / or units, and includes wheel rotation speed (wheel speed), engine rotation speed, engine rotation torque, engine throttle opening, and engine temperature (exhaust temperature). , Water temperature, oil temperature), output torque of the drive motor, current generated by the drive motor, failure code from the engine control unit, and the like.
  • the abnormality detection device 20 acquires a plurality of vehicle signals measured by the measuring device 10, and detects an abnormality for each of the acquired vehicle signals.
  • the abnormality detection device 20 is composed of, for example, a CPU (Central Processing Unit), a ROM (Read-Only Memory), a RAM (Random Access Memory), an input / output circuit, and the like.
  • the abnormality detection device 20 includes a sampling rate adjustment unit 21, a time window processing unit 22, an abnormality degree calculation unit 23, an abnormality determination unit 24, and a learning unit 25 as functional units.
  • the sampling rate adjusting unit 21 acquires a plurality of (n types) vehicle signals measured by the measuring device 10 and adjusts and matches the sampling rates of the acquired plurality of vehicle signals (step S2 in FIG. 2).
  • the plurality of vehicle signals include time series data for different sampling rates.
  • the signal included in each vehicle signal is indicated by a black circle on the axis indicating the time.
  • the sampling rate adjusting unit 21 is based on a predetermined sampling rate stored in advance, and among the data in the vehicle signal, the data closest to each of the timings (every time interval t) of the predetermined sampling rate (dotted circles in FIG. 4).
  • the data enclosed by the mark is extracted as the data at the relevant timing, and the other data is thinned out.
  • the sampling rate adjustment unit 21 outputs the vehicle signal whose sampling rate has been adjusted to the time window processing unit 22.
  • the predetermined sampling rate is at least the highest sampling rate among the sampling rates of the plurality of vehicle signals.
  • the sampling rate adjusting unit 21 can improve the correlation of each vehicle signal by such processing. Further, the sampling rate adjusting unit 21 can prevent overfitting of a vehicle signal having a high sampling rate.
  • the time window processing unit 22 is mounted on the volatile memory 20a.
  • the time window processing unit 22 acquires n types of vehicle signals whose sampling rates have been adjusted by the sampling rate adjustment unit 21, and from the vehicle signals whose sampling rates have been adjusted, a predetermined number (m pieces) included in the time window.
  • Data B n is extracted (step S3 in FIG. 2).
  • t is the slide width of the time window and the cycle of abnormality determination (time required to execute steps S2 to S9).
  • B n (or B n (0) ) is the data (m ⁇ 1 matrix) included in the time window of this time.
  • the time window processing unit 22 is a process of extracting m pieces of data included in the time window by sliding the time window for a predetermined period with the time t as the slide width. Generates B n.
  • the slide width t of the time window may be smaller than the predetermined period of the time window. Further, it is desirable that the slide width t of the time window is the same as the time interval t of the predetermined sampling rate.
  • the time window processing unit 22 outputs the generated B n to the abnormality degree calculation unit 23, and stores C n , which is a matrix in which B n is arranged, in the storage 20 b as storage information.
  • C n is the storage information is represented by the following formula.
  • B n (p) is data (m ⁇ 1 matrix) included in the time window p times before B n.
  • C n is stored information (m ⁇ p matrix) stored in the storage 20b.
  • p is the number of times the data Bn is extracted (in the present embodiment, a value equal to the number of times of abnormality determination in the abnormality determination unit 24), and increases as time passes from the start of operation of the abnormality detection device 20.
  • the time window processing unit 22 enables real-time learning and abnormality calculation by dividing the vehicle signal into each input number m of the self-encoder 23a (see FIG. 7) described later.
  • the abnormality degree calculation unit 23 acquires the learned information stored in the volatile memory 20a (step S1 in FIG. 2), and also acquires the B n generated by the time window processing unit 22. Then, the degree of abnormality of the vehicle signal is calculated based on the acquired Bn and the learned information.
  • the learned information will be described in detail later.
  • the abnormality degree calculation unit 23 includes the current B n (0) generated by the time window processing unit 22 and the previous C n, that is, B n (1) , B n (2) , ..., B.
  • the degree of abnormality of the vehicle information is determined. calculate.
  • the learning data L n used here is an m ⁇ (k + 1) matrix in which B n (0) and the data of columns 1 to k extracted from the storage data C n are combined, and is expressed by the following equation. Represented by.
  • the data in column 0 (first column) is B n (0) this time
  • the data in column k (k + 1st column) is B n (k) k times before.
  • k is determined by trial and error so as to satisfy the learning time for updating the learned information at an arbitrary timing and to secure the desired monitoring accuracy.
  • the abnormality degree calculation unit 23 scales (normalizes) the time window data B n based on the learned information, and calculates the scaled data B s, n (step S4 in FIG. 2).
  • Scaled data B s, n is the following formula, is calculated by calculating all the matrix elements of B S.
  • min n is the minimum value of the data in the nth signal in the learned information.
  • max n is the maximum value of the data in the nth signal in the learned information.
  • B s and n are data (m ⁇ 1 matrix) obtained by scaling B n.
  • [B n ] i is the i-th scalar quantity of B n.
  • [B s, n ] i is the i-th scalar quantity of B s, n. All the vector elements of the scaled data Bs , n fall within the range of the minimum value 0 to the maximum value 1.
  • the abnormality degree calculation unit 23 can handle a plurality of vehicle information having different orders in the same manner and calculate the abnormality degree by such a scaling process.
  • the anomaly degree calculation unit 23 reconstructs the scaled data Bs and n by self-coding using the self-encoder 23a (see FIG. 7) (step S5 in FIG. 2). ..
  • the self-encoder 23a includes an input layer 23a1 to which input data is input, an output layer 23a2 to which reconstruction data is output, and an intermediate layer (hidden layer) 23a3 provided between the input layer 23a1 and the output layer 23a2. , Equipped with.
  • the latent variable z and the reconstruction data x ⁇ in the intermediate layer 23a3 are calculated by the following equations, respectively.
  • x is a matrix input to the self-encoder 23a.
  • W x (1) , b x (1) , W x (2) , and b x (2) are self-encoder information corresponding to the input x.
  • z is a latent variable
  • x ⁇ ( ⁇ on x) is a matrix (reconstruction data) output from the self-encoder 23a.
  • the abnormality degree calculation unit 23 calculates the error (reconstruction error) e included in the reconstruction data (step S6 in FIG. 2).
  • the reconstruction error e is calculated by the following formula.
  • ⁇ x is a difference
  • e is the Euclidean norm of ⁇ x.
  • the abnormality degree calculation unit 23 can collectively process them to calculate the reconstruction error.
  • the numbers of the vehicle signals having a correlation are a and b
  • the input data and x and the reconstruction data x ⁇ are expressed by the following equations.
  • the anomaly degree calculation unit 23 calculates the reconstruction error e a from B s, a and B ⁇ s, a ( ⁇ on B), and B s, b and B ⁇ s, b (on B). ⁇ ) to calculate the reconstruction error e b from.
  • the self-encoder information is Wab (1) , bab (1) , Wab (2) , bab (2) .
  • the self-encoder 23a is not limited to one that uses the sigmoid function ⁇ as the activation function, and may use, for example, a hyperbolic tangent function or the like. Further, the intermediate layer 23a in the self-encoder is not limited to one layer, and a plurality of layers may be laminated. Further, the method for calculating the reconstruction error is not limited to the one using the Euclidean norm (L2 norm), and may be, for example, the one using the L1 norm or the like.
  • L2 norm Euclidean norm
  • the reconstruction method is not limited to the one using the self-encoder 23a, and may be, for example, the one using principal component analysis or the like. Since the vehicle information at the time of determination is almost only normal data, the abnormality degree calculation unit 23 cannot perform processing using another method using a neural network such as classification learning. On the other hand, the abnormality degree calculation unit 23 extracts latent variables using the self-encoder 23a that has learned only normal data, and reconstructs the data based on the extracted latent variables. The anomaly degree calculation unit 23 generates normal waveform data learned by reconstruction even when there is anomalous data deviating from the waveform data. Therefore, the abnormality degree calculation unit 23 calculates the reconstruction error e as information deviating from the normal waveform data.
  • error probability calculation unit 23 a learned information determination function f n (x), and, on the basis of the reconstruction error e n, scores of the n-th vehicle information (error probability) is calculated S n (Step S7 in FIG. 2), the calculation result is output to the abnormality determination unit 24.
  • the degree of abnormality Sn is calculated by the following formula.
  • the anomaly degree calculation unit 23 uses a one-class SVM (Support Vector Machine) as the determination function f n (x). Therefore, the score S n, determined function f n (x) takes a value greater than or equal to zero if the normal value is entered, the decision function f n minus if the abnormal value is input to the (x) Take a value.
  • SVM Serial Vector Machine
  • Abnormality determination unit 24 the abnormality score (error probability) calculated by the calculating unit 23 obtains the S n, based on the obtained scores, determines whether the vehicle signal is abnormal (in FIG. 2 Step S8).
  • the abnormality determination unit 24 determines that the corresponding vehicle signal is abnormal when a state in which the score Sn is negative occurs continuously for a plurality of time windows, and responds in other cases. It is determined that the vehicle signal is normal.
  • the abnormality determination unit 24 outputs the determination result to the notification device 30 (step S9 in FIG. 2).
  • the abnormality determination unit 24 can prevent the vehicle signal from being determined as an abnormality when pulse noise is generated, and can detect a true abnormality.
  • the learning unit 25 learns in real time in parallel with the operations of the sampling rate adjustment unit 21, the time window processing unit 22, the abnormality degree calculation unit 23, and the abnormality determination unit 24. And generate learned information.
  • the learning unit 25 extracts the maximum value max n and the minimum value min n included in the learning data L n (step S12 in FIG. 3). Subsequently, the learning unit 25 scales the learning data L n based on the extracted maximum value max n and the minimum value min n , and calculates the scaled learning data L s, n. Subsequently, the learning unit 25 learns the self-encoder 23a using the learned data Ls and n , and the self-encoder information W x (1) , b x (1) , W x (2). , B x (2) is generated (step S12 in FIG. 3).
  • a scaling conjugate gradient descent method can be mentioned.
  • a preset value is used as the hyperparameter.
  • the learning unit 25 calculates the reconstruction error based on the self-encoder information W x (1) , b x (1) , W x (2) , b x (2) , and recalculates the calculated error.
  • the decision function f n (x) is generated with the configuration error as normal data (step S12 in FIG. 3).
  • the learning unit 25 has a maximum value max n and a minimum value min n , self-encoder information W x (1) , b x (1) , W x (2) , b x (2) , and a determination function f n (x). ) Is stored in the storage 20b as learned information.
  • the learning unit 25 learns independently of the processing of the sampling rate adjustment unit 21, the time window processing unit 22, the abnormality degree calculation unit 23, and the abnormality determination unit 24, it is separate from the abnormality detection device 20 ( It can also be embodied by another processor, computer (PC), etc.).
  • the update unit 26 reads the updated learned information and stores it in the volatile memory 20a (step S13 in FIG. 3).
  • the notification device 30 is composed of a display capable of outputting an image, a speaker capable of outputting voice, and the like, and acquires a determination result by the abnormality determination unit 24 and notifies the user.
  • the abnormality detection device 20 has a time window processing unit 22 for extracting data included in the time window by dividing a plurality of vehicle signals including time series data by a time window, and the time window.
  • the learning unit 25 that generates learned information by learning using the time-series data of the vehicle signal separated by, and the time-series data of the vehicle signal separated by the time window are based on the learned information. It is characterized by including an abnormality degree calculation unit 23 that reconstructs, calculates an error before and after the reconstruction, and calculates an abnormality degree based on the calculated error and the learned information.
  • the abnormality detection device 20 can calculate the degree of abnormality of the vehicle information delimited by the time window by using the learned information using the vehicle information delimited by the time window before that. That is, the abnormality detection device 20 does not need to prepare the learned information in advance of the detection work, and can suitably detect the abnormality based on the learned information updated in real time.
  • the abnormality degree calculation unit 23 uses the time series data of the vehicle signal separated by the time window this time and the time series of the vehicle signal separated by the time window before the previous time. It is characterized in that the degree of abnormality is calculated based on the learned information using the data. Therefore, the abnormality detection device 20 can more preferably execute the generation of learned data and the calculation of the degree of abnormality in parallel.
  • the learned information includes the maximum value and the minimum value of the time series data of the vehicle signal separated by the time window before the previous time
  • the abnormality degree calculation unit includes the maximum value.
  • the time-series data of the vehicle signal separated by the time window of this time is scaled, and the scaled time-series data of the vehicle signal is self-encoded included in the learned information. It is characterized by being reconstructed by a vessel 23a. Therefore, the abnormality detection device 20 can handle a plurality of vehicle information having different orders in the same manner and appropriately calculate the degree of abnormality.
  • the abnormality detection device 20 includes a sampling rate adjusting unit 21 for matching the sampling rates of the plurality of vehicle signals, and the time window processing unit 22 has time-series data of the plurality of vehicle signals whose sampling rates have been adjusted. Is characterized by being separated by the time window. Therefore, the abnormality detection device 20 can improve the correlation of each vehicle signal by matching the sampling rates of the plurality of vehicle information, and can prevent over-learning of the vehicle signal having a large sampling rate. ..
  • the abnormality detection device 20 has an abnormality determination unit 24 that determines that the vehicle signal is abnormal when an abnormality occurs continuously in a predetermined number of the time windows based on the calculated abnormality degree. It is characterized by being prepared. Therefore, the abnormality detection device 20 can prevent the vehicle signal from being determined as an abnormality when pulse noise is generated, and can detect a true abnormality.
  • the present invention is not limited to the above-described embodiments and can be appropriately modified without departing from the gist of the present invention.
  • the abnormality detection device 20 may be configured to omit the abnormality determination unit 24.
  • the calculation result of the abnormality degree calculation unit 23 is output to the notification device 30 by voice or image.
  • the present invention can also be realized as an abnormality detection program that causes the computer to function as the abnormality detection device 20.
  • Anomaly detection system 10 Measuring device 20 Anomaly detection device 21 Sampling rate adjustment unit 22 Time window processing unit 23 Abnormality calculation unit 24 Abnormality determination unit 25 Learning unit

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • Mechanical Engineering (AREA)
  • Transportation (AREA)
  • Human Computer Interaction (AREA)
  • Biomedical Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Testing And Monitoring For Control Systems (AREA)
  • Traffic Control Systems (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

車両の個体差及び試験中の消耗による影響を抑えた異常検知を行うことが可能な異常検知装置を提供する。 異常検知装置(20)は、時系列データを含む複数の車両信号を時間窓で区切ることによって、当該時間窓に含まれるデータを抽出する時間窓処理部(22)と、時間窓で区切られた車両信号の時系列データを用いた学習によって、学習済み情報を生成する学習部(25)と、時間窓で区切られた車両信号の時系列データを学習済み情報に基づいて再構成し、再構成前後の誤差を算出するとともに、算出された誤差と学習済み情報とに基づいて異常度を算出する異常度算出部(23)と、を備える。

Description

異常検知装置及び異常検知プログラム
 本発明は、車両信号の異常を検知する技術に関する。
 特許文献1,2には、車両の耐久試験の際に、正常を示す基準データを、別の車両を用いて予め採取しておき、当該基準データを用いて耐久試験中に計測されたデータの異常を検知する技術が記載されている。
特開2015-170121号公報 特開2016-045861号公報
 特許文献1,2に記載の技術では、基準データが試験対象とは異なる車両から採取されているため、試験対象の車両の基準データに対する個体差が大きい場合には、正常であっても異常と検知される可能性がある。また、基準データが消耗の少ない車両から採取されている場合には、試験中の消耗による計測データの変動を異常と検知する可能性が高くなる。
 本発明は、前記の点に鑑みてなされたものであり、車両の個体差及び試験中の消耗による影響を抑えた異常検知を行うことが可能な異常検知装置及び異常検知プログラムを提供することを課題とする。
 前記の課題を解決するために、本発明の異常検知装置は、時系列データを含む複数の車両信号を時間窓で区切ることによって、当該時間窓に含まれるデータを抽出する時間窓処理部と、前記時間窓で区切られた前記車両信号の時系列データを用いた学習によって、学習済み情報を生成する学習部と、前記時間窓で区切られた前記車両信号の時系列データを前記学習済み情報に基づいて再構成し、再構成前後の誤差を算出するとともに、算出された前記誤差と前記学習済み情報とに基づいて異常度を算出する異常度算出部と、を備えることを特徴とする。
 また、本発明の異常検知プログラムは、コンピュータを前記異常検知装置として機能させることを特徴とする。
 本発明によると、車両の異常検知に際して、車両の個体差及び試験中の消耗による影響を抑えることができる。
本発明の実施形態に係る異常検知システムを模式的に示すブロック図である。 本発明の実施形態に係る異常検知装置の動作例を説明するためのフローチャートである。 本発明の実施形態に係る異常検知装置の動作例を説明するためのフローチャートである。 車両信号のサンプリングレートを一致させる動作例を説明するための模式図である。 サンプリングレートが調整された車両信号を時間窓で区切る動作例を説明するための模式図である。 サンプリングレートが調整された車両信号の波形と時間窓との関係を説明するための模式図である。 自己符号化器の一例を説明するための模式図である。 自己符号化器の入力データ及び出力データ並びに再構成誤差を説明するための模式図であり、車両信号が正常である場合を示す図である。 自己符号化器の入力データ及び出力データ並びに再構成誤差を説明するための模式図であり、車両信号が異常である場合を示す図である。
 次に、本発明の実施形態について、適宜図面を参照しながら詳細に説明する。なお、同様の部分には同一の符号を付し、重複する説明を省略する。
 図1に示すように、本発明の実施形態に係る異常検知システム1は、車両2の異常を検知するためのシステムであって、計測装置10と、異常検知装置20と、通知装置30と、を備える。
<計測装置>
 計測装置10は、台上耐久試験における模擬的な試験走行中の車両2の各種信号(車両信号)を計測し、計測された複数の車両信号を異常検知装置20へ送信する。計測される車両信号は、オーダ及び/又は単位が異なる信号を含み、車輪の回転速度(車輪速)、エンジンの回転速度、エンジンの回転トルク、エンジンのスロットル開度、エンジンの各部温度(排気温、水温、油温)、駆動モータの出力トルク、駆動モータの発生電流、エンジン制御ユニットからの故障コード等が例として挙げられる。
<異常検知装置>
 異常検知装置20は、計測装置10によって計測された複数の車両信号を取得し、取得された車両信号ごとに異常を検知する。異常検知装置20は、例えば、CPU(Central Processing Unit)、ROM(Read-Only Memory)、RAM(Random Access Memory)、入出力回路等によって構成されている。異常検知装置20は、機能部として、サンプリングレート調整部21と、時間窓処理部22と、異常度算出部23と、異常判定部24と、学習部25と、を備える。
≪サンプリングレート調整部≫
 サンプリングレート調整部21は、計測装置10によって計測された複数(n種類)の車両信号を取得し、取得された複数の車両信号のサンプリングレートを調整して一致させる(図2のステップS2)。図4に示すように、複数の車両信号は、それぞれ異なるサンプリングレートごとの時系列データを含む。図4において、各車両信号に含まれる信号は、時刻を示す軸上の黒丸によって示されている。サンプリングレート調整部21は、予め記憶された所定サンプリングレートに基づいて、車両信号内のデータのうち、所定サンプリングレートのタイミング(時間間隔tごと)のそれぞれに最も近いデータ(図4において、点線丸印で囲まれたデータ)を、当該タイミングにおけるデータとして抽出するとともに、それ以外のデータを間引く。サンプリングレート調整部21は、サンプリングレートが調整された車両信号を時間窓処理部22へ出力する。ここで、所定サンプリングレートは、複数の車両信号のサンプリングレートのうち、最も大きいもの以上であることが望ましい。
 サンプリングレート調整部21は、かかる処理によって、各車両信号の相関性を向上することができる。また、サンプリングレート調整部21は、サンプリングレートが大きい車両信号の過学習を防止することができる。
≪時間窓処理部≫
 時間窓処理部22は、揮発性メモリ20a上に実装されている。時間窓処理部22は、サンプリングレート調整部21によってサンプリングレートが調整されたn種類の車両信号を取得し、サンプリングレートが調整された車両信号から、時間窓に含まれる所定個数(m個)のデータBを抽出する(図2のステップS3)。ここで、tは、時間窓のスライド幅であるとともに、異常判定の周期(ステップS2~S9の実行に要する時間)である。また、B(又はB (0))は、今回の時間窓に含まれるデータ(m×1行列)である。
 図5及び図6に示すように、時間窓処理部22は、所定期間である時間窓を、時間tをスライド幅として滑走させることによって、当該時間窓に含まれるm個のデータを抽出する処理によって、Bを生成する。ここで、時間窓のスライド幅tは、時間窓の所定期間よりも小さくてもよい。また、時間窓のスライド幅tは、所定サンプリングレートの時間間隔tと同一であることが望ましい。
 時間窓処理部22は、生成されたBを異常度算出部23へ出力するとともに、Bを並べた行列であるCを、保存情報としてストレージ20b内に記憶させる。今回の時間窓データをB (0)と定義すると、保存情報であるCは、下記式によって表される。
Figure JPOXMLDOC01-appb-M000001
 ここで、B (p)は、Bのp回前の時間窓に含まれるデータ(m×1行列)である。また、Cは、ストレージ20bに記憶される保存情報(m×p行列)である。pは、データBの抽出回数(本実施形態では、異常判定部24における異常判定回数と等しい値)であり、異常検知装置20の動作開始から時間が経つほど増大する。
 時間窓処理部22は、車両信号を後記する自己符号化器23a(図7参照)の入力数mごとに区切ることによって、リアルタイムでの学習及び異常度算出を可能とする。
≪異常度算出部≫
 図1に示すように、異常度算出部23は、揮発性メモリ20aに記憶された学習済み情報を取得する(図2のステップS1)とともに、時間窓処理部22によって生成されたBを取得し、取得されたB及び学習済み情報に基づいて、車両信号の異常度を算出する。学習済み情報に関しては、後で詳細に説明する。
 本実施形態において、異常度算出部23は、時間窓処理部22によって生成された今回のB (0)と、前回のCすなわちB (1),B (2),…,B (p)を用いて学習部25によって学習された学習用データL(段落0046参照)を用いて学習された学習済み情報(段落0048参照)と、に基づいて、車両情報の異常度を算出する。ここで用いられる学習用データLは、B (0)と、ストレージデータであるCから抽出された列1~kのデータと、を組み合わせたm×(k+1)行列であり、下記式によって表される。
Figure JPOXMLDOC01-appb-M000002
 ここで、列0(1列目)のデータは、今回のB (0)であり、列k(k+1列目)のデータは、k回前のB (k)である。ここで、kが大きい場合には、学習済みモデルの精度が向上する一方、ステップS12における学習時間が長くなる。そのため、kは、学習済み情報を任意のタイミングで更新するための学習時間を満たすとともに、所望の監視精度を確保することができるよう、試行錯誤的に決定される。
 まず、異常度算出部23は、学習済み情報に基づいて、時間窓データBをスケーリング(正規化)し、スケーリングされたデータBs,nを算出する(図2のステップS4)。スケーリングされたデータBs,nは、下記式により、Bの行列要素を全て計算することによって算出される。
Figure JPOXMLDOC01-appb-M000003
 ここで、minは、学習済み情報におけるn番目の信号内のデータの最小値である。また、maxは、学習済み情報におけるn番目の信号内のデータの最大値である。また、Bs,nは、Bをスケーリングしたデータ(m×1行列)である。また、[Bは、Bのi番目のスカラー量である。また、[Bs,nは、Bs,nのi番目のスカラー量である。スケーリングされたデータBs,nのベクトル要素は、全て、最小値0から最大値1の範囲内に収まる。
 異常度算出部23は、かかるスケーリング処理によって、オーダが異なる複数の車両情報を同様に取り扱って異常度を算出することができるようにする。
 続いて、異常度算出部23は、スケーリングされたデータBs,nを、自己符号化器23a(図7参照)を用いて自己符号化することによって、再構成する(図2のステップS5)。自己符号化器23aは、入力データが入力される入力層23a1と、再構成データが出力される出力層23a2と、入力層23a1及び出力層23a2の間に設けられる中間層(隠れ層)23a3と、を備える。中間層23a3における潜在変数z及び再構成データx^は、それぞれ下記式によって算出される。
Figure JPOXMLDOC01-appb-M000004
 ここで、xは、自己符号化器23aに入力される行列である。また、W (1),b (1),W (2),b (2)は、入力xに対応する自己符号化器情報である。また、zは、潜在変数であり、x^(xの上に^)は、自己符号化器23aから出力される行列(再構成データ)である。
 続いて、異常度算出部23は、再構成データに含まれる誤差(再構成誤差)eを算出する(図2のステップS6)。再構成誤差eは、下記式によって算出される。
Figure JPOXMLDOC01-appb-M000005
 ここで、Δxは、差分であり、再構成誤差eは、Δxのユークリッドノルムである。
 異常度算出部23は、Bs,nを自己符号化器23aに入力することによって(x=Bs,n)、n番目の車両信号の再構成誤差eを算出する。図8A及び図8Bに示すように、自己符号化器23aに入力される車両情報が正常である場合(図8A)には、再構成誤差eは小さく、自己符号化器23aに入力される車両情報が異常である場合(図8B)には、再構成誤差eは大きい。
 また、異常度算出部23は、複数の車両信号が強い相関性を有する場合には、これらをまとめて処理して再構成誤差を算出することができる。ここで、相関性を有する車両信号の番号をa,bとすると、入力データ及びx及び再構成データx^は、下記式によって表される。
Figure JPOXMLDOC01-appb-M000006
 異常度算出部23は、Bs,a及びB^s,a(Bの上に^)から再構成誤差eを算出するとともに、Bs,b及びB^s,b(Bの上に^)から再構成誤差eを算出する。この場合には、自己符号化器情報は、Wab (1),bab (1),Wab (2),bab (2)である。
 自己符号化器23aは、活性化関数としてシグモイド関数σを用いるものに限定されず、例えば、双曲線正接関数等を用いるものであってもよい。また、自己符号化器における中間層23aは、1層に限定されず、複数層が積層されていてもよい。また、再構成誤差の算出手法は、ユークリッドノルム(L2ノルム)を用いるものに限定されず、例えば、L1ノルム等を用いるものであってもよい。
 また、再構成の手法は、自己符号化器23aを用いるものに限定されず、例えば、主成分分析等を用いるものであってもよい。判定を行う際の車両情報は、殆ど正常データのみであるため、異常度算出部23は、分類学習等といったニューラルネットワークを用いた他の手法を用いた処理を行うことができない。これに対し、異常度算出部23は、正常データのみを学習した自己符号化器23aを用いて潜在変数を抽出し、抽出された潜在変数に基づいてデータを再構成する。異常度算出部23は、波形データから逸脱した異常データが存在する場合でも、再構成によって学習された正常な波形データを生成する。したがって、異常度算出部23は、正常な波形データから逸脱した情報として、再構成誤差eを算出する。
 続いて、異常度算出部23は、学習済み情報である決定関数f(x)、及び、再構成誤差eに基づいて、n番目の車両情報のスコア(異常度)Sを算出し(図2のステップS7)、算出結果を異常判定部24へ出力する。異常度Sは、下記式によって算出される。
Figure JPOXMLDOC01-appb-M000007
 異常度算出部23は、決定関数f(x)として1クラスSVM(Support Vector Machine)を用いている。したがって、スコアSは、決定関数f(x)に正常値が入力された場合にはゼロ以上の値をとり、決定関数f(x)に異常値が入力された場合にはマイナスの値をとる。
≪異常判定部≫
 異常判定部24は、異常度算出部23によって算出されたスコア(異常度)Sを取得し、取得されたスコアに基づいて、車両信号が異常であるか否かを判定する(図2のステップS8)。異常判定部24は、スコアSがマイナスである状態が複数の時間窓に対して連続して発生した場合に、対応する車両信号は異常であると判定し、それ以外の場合には、対応する車両信号は正常であると判定する。異常判定部24は、異常であると判定した場合に、判定結果を通知装置30へ出力する(図2のステップS9)。
 異常判定部24は、車両信号にパルス的ノイズが発生した場合に異常と判定されることを防止し、真の異常を検知することができる。
≪学習部≫
 学習部25は、ストレージ20aに記憶された保存データに基づいて、サンプリングレート調整部21、時間窓処理部22、異常度算出部23及び異常判定部24の動作と並行して、リアルタイムで学習を行い、学習済み情報を生成する。学習部25は、ストレージ20bに記憶されたCのうち、学習に必要な時間の行列を読み出し、読み出されたデータを学習用データLとする(図3のステップS11)。
 続いて、学習部25は、学習用データLに含まれる最大値max及び最小値minを抽出する(図3のステップS12)。続いて、学習部25は、抽出された最大値max及び最小値minに基づいて学習用データLをスケーリングし、スケーリング済み学習用データLs,nを算出する。続いて、学習部25は、学習済みデータLs,nを用いて自己符号化器23aの学習を行い、自己符号化器情報W (1),b (1),W (2),b (2)を生成する(図3のステップS12)。学習部25による自己符号化器23aの学習の手法(アルゴリズム)としては、例えばスケーリング共役勾配降下法が挙げられる。この場合、ハイパーパラメータとして、予め選定された値が用いられる。続いて、学習部25は、自己符号化器情報W (1),b (1),W (2),b (2)に基づいて再構成誤差を算出し、算出された再構成誤差を正常データとして決定関数f(x)を生成する(図3のステップS12)。
 学習部25は、最大値max及び最小値min、自己符号化器情報W (1),b (1),W (2),b (2)並びに決定関数f(x)を学習済み情報としてストレージ20bに記憶させる。
 なお、学習部25は、サンプリングレート調整部21、時間窓処理部22、異常度算出部23及び異常判定部24の処理とは独立して学習を行うため、異常検知装置20とは別体(別のプロセッサ、計算機(PC)等)によって具現化することも可能である。
 更新部26は、ストレージ20bに記憶された学習済み情報が更新された場合に、更新された学習済み情報を読み出し、揮発性メモリ20aに記憶させる(図3のステップS13)。
<通知装置>
 通知装置30は、画像出力可能なディスプレイ、音声出力可能なスピーカ等によって構成されており、異常判定部24による判定結果を取得してユーザに通知する。
 本発明の実施形態に係る異常検知装置20は、時系列データを含む複数の車両信号を時間窓で区切ることによって、当該時間窓に含まれるデータを抽出する時間窓処理部22と、前記時間窓で区切られた前記車両信号の時系列データを用いた学習によって、学習済み情報を生成する学習部25と、前記時間窓で区切られた前記車両信号の時系列データを前記学習済み情報に基づいて再構成し、再構成前後の誤差を算出するとともに、算出された前記誤差と前記学習済み情報とに基づいて異常度を算出する異常度算出部23と、を備えることを特徴とする。
 したがって、異常検知装置20は、時間窓で区切られた車両情報の異常度を、それよりも前の時間窓で区切られた車両情報を用いた学習済み情報を用いて算出することができる。すなわち、異常検知装置20は、学習済み情報を検知作業よりも事前に準備する必要が無く、リアルタイムで更新された学習済み情報に基づいて、異常を好適に検知することができる。
 また、異常検知装置20は、前記異常度算出部23が、今回の前記時間窓によって区切られた前記車両信号の時系列データと、前回以前の前記時間窓によって区切られた前記車両信号の時系列データを用いた前記学習済み情報と、に基づいて、前記異常度を算出することを特徴とする。
 したがって、異常検知装置20は、学習済みデータの生成及び異常度の算出をより好適に並行して実行することができる。
 また、異常検知装置20は、前記学習済み情報が、前回以前の前記時間窓によって区切られた前記車両信号の時系列データの最大値及び最小値を含み、前記異常度算出部は、前記最大値及び前記最小値に基づいて、今回の前記時間窓によって区切られた前記車両信号の時系列データをスケーリングし、スケーリングされた前記車両信号の時系列データを、前記学習済み情報に含まれる自己符号化器23aによって再構成することを特徴とする。
 したがって、異常検知装置20は、オーダが異なる複数の車両情報を同様に取り扱って異常度を好適に算出することができる。
 また、異常検知装置20は、複数の前記車両信号のサンプリングレートを一致させるサンプリングレート調整部21を備え、前記時間窓処理部22は、サンプリングレートが調整された複数の前記車両信号の時系列データを、前記時間窓で区切ることを特徴とする。
 したがって、異常検知装置20は、複数の車両情報のサンプリングレートを一致させることによって、各車両信号の相関性を向上することができるとともに、サンプリングレートが大きい車両信号の過学習を防止することができる。
 また、異常検知装置20は、算出された前記異常度に基づいて、所定数の前記時間窓に連続して異常が発生した場合に、前記車両信号は異常であると判定する異常判定部24を備えることを特徴とする。
 したがって、異常検知装置20は、車両信号にパルス的ノイズが発生した場合に異常と判定されることを防止し、真の異常を検知することができる。
 以上、本発明の実施形態について説明したが、本発明は前記実施形態に限定されず、本発明の要旨を逸脱しない範囲で適宜変更可能である。例えば、異常検知装置20は、異常判定部24を省略する構成であってもよい。この場合、通知装置30には、異常度算出部23の算出結果が音声又は画像によって出力される。また、本発明は、コンピュータを異常検知装置20として機能させる異常検知プログラムとしても実現可能である。
 1   異常検知システム
 10  計測装置
 20  異常検知装置
 21  サンプリングレート調整部
 22  時間窓処理部
 23  異常度算出部
 24  異常判定部
 25  学習部

Claims (6)

  1.  時系列データを含む複数の車両信号を時間窓で区切ることによって、当該時間窓に含まれるデータを抽出する時間窓処理部と、
     前記時間窓で区切られた前記車両信号の時系列データを用いた学習によって、学習済み情報を生成する学習部と、
     前記時間窓で区切られた前記車両信号の時系列データを前記学習済み情報に基づいて再構成し、再構成前後の誤差を算出するとともに、算出された前記誤差と前記学習済み情報とに基づいて異常度を算出する異常度算出部と、
     を備えることを特徴とする異常検知装置。
  2.  前記異常度算出部は、今回の前記時間窓によって区切られた前記車両信号の時系列データと、前回以前の前記時間窓によって区切られた前記車両信号の時系列データを用いた前記学習済み情報と、に基づいて、前記異常度を算出する
     ことを特徴とする請求項1に記載の異常検知装置。
  3.  前記学習済み情報は、前回以前の前記時間窓によって区切られた前記車両信号の時系列データの最大値及び最小値を含み、
     前記異常度算出部は、前記最大値及び前記最小値に基づいて、今回の前記時間窓によって区切られた前記車両信号の時系列データをスケーリングし、スケーリングされた前記車両信号の時系列データを、前記学習済み情報に含まれる自己符号化器によって再構成する
     ことを特徴とする請求項2に記載の異常検知装置。
  4.  複数の前記車両信号のサンプリングレートを一致させるサンプリングレート調整部を備え、
     前記時間窓処理部は、サンプリングレートが調整された複数の前記車両信号の時系列データを、前記時間窓で区切る
     ことを特徴とする請求項1から請求項3のいずれか一項に記載の異常検知装置。
  5.  算出された前記異常度に基づいて、所定数の前記時間窓に連続して異常が発生した場合に、前記車両信号は異常であると判定する異常判定部を備える
     ことを特徴とする請求項1から請求項4のいずれか一項に記載の異常検知装置。
  6.  コンピュータを請求項1から請求項5のいずれか一項に記載の異常検知装置として機能させる異常検知プログラム。
PCT/JP2021/023245 2020-07-20 2021-06-18 異常検知装置及び異常検知プログラム WO2022019020A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202180049128.5A CN115989465A (zh) 2020-07-20 2021-06-18 异常检测装置及异常检测程序
US18/005,897 US20230322240A1 (en) 2020-07-20 2021-06-18 Abnormality detection device and abnormality detection program
MX2023000909A MX2023000909A (es) 2020-07-20 2021-06-18 Dispositivo de deteccion de anomalias y programa de deteccion de anomalias.

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020124062A JP7072611B2 (ja) 2020-07-20 2020-07-20 異常検知装置及び異常検知プログラム
JP2020-124062 2020-07-20

Publications (1)

Publication Number Publication Date
WO2022019020A1 true WO2022019020A1 (ja) 2022-01-27

Family

ID=79729448

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/023245 WO2022019020A1 (ja) 2020-07-20 2021-06-18 異常検知装置及び異常検知プログラム

Country Status (5)

Country Link
US (1) US20230322240A1 (ja)
JP (1) JP7072611B2 (ja)
CN (1) CN115989465A (ja)
MX (1) MX2023000909A (ja)
WO (1) WO2022019020A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115656847A (zh) * 2022-12-08 2023-01-31 东莞先知大数据有限公司 一种电动车电池异常确定方法、电子设备和存储介质

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210253116A1 (en) * 2018-06-10 2021-08-19 Osr Enterprises Ag System and method for enhancing sensor operation in a vehicle

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014234113A (ja) * 2013-06-04 2014-12-15 株式会社デンソー 車両用基準値生成装置
JP2019175462A (ja) * 2018-03-29 2019-10-10 コリア インスティテュート オブ オーシャン サイエンス テクノロジー 教師なし学習方法を用いた船舶の異常運航状態自動識別システムおよびその方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014234113A (ja) * 2013-06-04 2014-12-15 株式会社デンソー 車両用基準値生成装置
JP2019175462A (ja) * 2018-03-29 2019-10-10 コリア インスティテュート オブ オーシャン サイエンス テクノロジー 教師なし学習方法を用いた船舶の異常運航状態自動識別システムおよびその方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115656847A (zh) * 2022-12-08 2023-01-31 东莞先知大数据有限公司 一种电动车电池异常确定方法、电子设备和存储介质

Also Published As

Publication number Publication date
JP7072611B2 (ja) 2022-05-20
MX2023000909A (es) 2023-02-22
US20230322240A1 (en) 2023-10-12
CN115989465A (zh) 2023-04-18
JP2022020512A (ja) 2022-02-01

Similar Documents

Publication Publication Date Title
WO2022019020A1 (ja) 異常検知装置及び異常検知プログラム
JP6740247B2 (ja) 異常検出システム、異常検出方法、異常検出プログラム及び学習済モデル生成方法
US20240255942A1 (en) Automated analysis of non-stationary machine performance
JP2021528745A (ja) アプリケーション情報に関連する時系列ズデータに関する深層学習を使用した異常検出
CN102265227B (zh) 用于在机器状况监视中创建状态估计模型的方法和设备
US11544554B2 (en) Additional learning method for deterioration diagnosis system
EP3422517B1 (en) A method for recognizing contingencies in a power supply network
EP4068140A1 (en) Method and system for optimizing a simulation model using machine learning
WO2012164243A1 (en) Heterogeneous data fusion using gaussian processes
Choudhury Prognosticating autism spectrum disorder using artificial neural network: Levenberg-Marquardt algorithm
CN111079348B (zh) 一种缓变信号检测方法和装置
Deng et al. A novel sub-label learning mechanism for enhanced cross-domain fault diagnosis of rotating machinery
Zamba et al. A multivariate change-point model for change in mean vector and/or covariance structure
JP7238378B2 (ja) 異常検出装置、異常検出プログラム、及び、異常検出方法
AU2012264478A1 (en) Heterogeneous data fusion using gaussian processes
Ensley et al. Extrapolation of Mackey-Glass data using cascade correlation
CN116720073A (zh) 一种基于分类器的异常检测提取方法与系统
Yu et al. A Novel 1D‐CNN‐Based Diagnosis Method for a Rolling Bearing with Dual‐Sensor Vibration Data Fusion
JP2021071586A (ja) 音抽出システム及び音抽出方法
WO2020183781A1 (ja) 異常診断装置
CN114881110B (zh) 一种在轨航天器舱内总压变化模式的实时检测方法
Zhang et al. Bearing performance degradation assessment based on a combination of multi-scale entropy and K-medoids clustering
JP7281394B2 (ja) 異常診断装置およびプログラム
CN116821594B (zh) 基于频谱选择机制的图神经网络工业控制系统异常检测方法及装置
Vazquez et al. Fault Diagnosis Approach for Pedelec Drive Units Based on Support Vector Machines

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21846218

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21846218

Country of ref document: EP

Kind code of ref document: A1