JP2010127698A - Apparatus and method for estimation of radiated electromagnetic wave frequency - Google Patents

Apparatus and method for estimation of radiated electromagnetic wave frequency Download PDF

Info

Publication number
JP2010127698A
JP2010127698A JP2008301061A JP2008301061A JP2010127698A JP 2010127698 A JP2010127698 A JP 2010127698A JP 2008301061 A JP2008301061 A JP 2008301061A JP 2008301061 A JP2008301061 A JP 2008301061A JP 2010127698 A JP2010127698 A JP 2010127698A
Authority
JP
Japan
Prior art keywords
electromagnetic wave
time
analysis
signal
radiated
Prior art date
Legal status (The legal status 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 status listed.)
Pending
Application number
JP2008301061A
Other languages
Japanese (ja)
Inventor
Fumihiko Ishiyama
文彦 石山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2008301061A priority Critical patent/JP2010127698A/en
Publication of JP2010127698A publication Critical patent/JP2010127698A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a radiated electromagnetic wave frequency estimation apparatus that estimates the radiation of an unnecessary electromagnetic wave with a frequency to be addressed from a change in a current passing through a wire. <P>SOLUTION: The apparatus includes an A/D conversion unit 15 which obtains an analog signal indicating a current passing through the wire 12 of an electrical circuit and converts it into a digital time-series signal and an analyzer 13 which calculates and estimates the frequency of an electromagnetic wave radiated in accordance with a temporal change in the current passing through the wire from the temporal change in the current on the basis of the time-series signal from the A/D conversion unit. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

本発明は、プリント基板等に形成された電気回路から放射される電磁波の周波数を解析する放射電磁波周波数推定装置および放射電磁波周波数推定方法に関する。   The present invention relates to a radiated electromagnetic wave frequency estimation device and a radiated electromagnetic wave frequency estimation method for analyzing the frequency of an electromagnetic wave radiated from an electric circuit formed on a printed circuit board or the like.

従来、プリント基板等に形成された電気回路から放射される電磁波の周波数を解析するために種々の技術が知られている。例えば、特許文献1は、測定回数を減らし周波数の変化によるノイズ強度分布解析を短時間かつ容易にできるプリント基板のノイズ源探査装置を開示している。   Conventionally, various techniques are known for analyzing the frequency of an electromagnetic wave radiated from an electric circuit formed on a printed circuit board or the like. For example, Patent Document 1 discloses a noise source search device for a printed circuit board that can reduce the number of measurements and easily perform noise intensity distribution analysis by changing the frequency in a short time.

この特許文献1に開示されたノイズ源探査装置では、被測定プリント基板の画像データをパーソナルコンピュータに取り込み、スキャニング部で微小アンテナプローブを被測定プリント基板の2次元平面上に沿って走査させて放射ノイズを検出し、放射ノイズをスペクトラムアナライザで周波数分析を行い、その強度レベルと各測定区画の位置情報をパーソナルコンピュータに取り込んでメモリに格納し、マップデータを作成し、画像データにマップデータを重ね合わせてディスプレイに表示する。
特開2000−19204号公報
In the noise source search apparatus disclosed in Patent Document 1, image data of a printed circuit board to be measured is taken into a personal computer, and a scanning unit scans a minute antenna probe along a two-dimensional plane of the printed circuit board to be radiated. Detect noise, perform frequency analysis of radiated noise with a spectrum analyzer, capture the intensity level and position information of each measurement section into a personal computer, store it in memory, create map data, and overlay map data on image data Also display on the display.
JP 2000-19204 A

しかしながら、上述した従来の技術では、放射される電磁波の時間平均情報しか得ることができず、配線上の電流変化の、どのタイミングで対策すべき周波数の不要電磁波が放射されるのかに関する情報や、どの配線上の電流変化が不要電磁波を放射しているのか等に関する情報を得ることができず、不便であった。   However, in the conventional technology described above, it is possible to obtain only the time average information of the radiated electromagnetic wave, and information regarding the timing of the unnecessary electromagnetic wave of the frequency to be dealt with at which timing of the current change on the wiring, It was inconvenient because it was impossible to obtain information on which wiring current change radiated unnecessary electromagnetic waves.

本発明は、上述した不便を解消するためになされたものであり、その課題は、対策すべき周波数の不要電磁波の放射を、配線上を流れる電流変化から推定して、プリント基板等に形成された電気回路の構成を改善するための指針を与えることができる放射電磁波周波数推定装置および放射電磁波周波数推定方法を提供することにある。   The present invention has been made to solve the above-mentioned inconvenience, and the problem is that an unnecessary electromagnetic wave having a frequency to be counteracted is estimated from a change in current flowing on the wiring and formed on a printed circuit board or the like. Another object of the present invention is to provide a radiated electromagnetic wave frequency estimation device and a radiated electromagnetic wave frequency estimation method capable of providing a guideline for improving the configuration of an electric circuit.

上記課題を解決するために、請求項1の放射電磁波周波数推定装置は、電気回路の配線を流れる電流を示すアナログ信号を取り込んでデジタルの時系列信号に変換するA/D変換部と、A/D変換部からの時系列信号に基づき、配線を流れる電流の時間変化に伴って放射される電磁波の周波数を、電流の時間変化から算出して推定する解析装置を備えたことを特徴とする。   In order to solve the above-mentioned problem, the radiated electromagnetic wave frequency estimation apparatus according to claim 1 includes an A / D converter that takes in an analog signal indicating a current flowing through a wiring of an electric circuit and converts it into a digital time-series signal; An analysis device is provided that calculates and estimates the frequency of an electromagnetic wave radiated with the time change of the current flowing through the wiring based on the time series signal from the D conversion unit from the time change of the current.

請求項2の放射電磁波周波数推定装置は、請求項1の放射電磁波周波数推定装置において、解析装置は、A/D変換部からの時系列信号を用いて線形予測解析を行う信号解析部と、信号解析部における線形予測解析により得られた予測係数を用いて配線を流れる電流の時間変化に伴って放射される電磁波の周波数を推定する周波数推定部を備えたことを特徴とする。   The radiated electromagnetic wave frequency estimation device according to claim 2 is the radiated electromagnetic wave frequency estimation device according to claim 1, wherein the analysis device uses a time series signal from the A / D conversion unit to perform linear prediction analysis, and a signal A frequency estimation unit is provided that estimates a frequency of an electromagnetic wave radiated with time change of a current flowing through the wiring using a prediction coefficient obtained by linear prediction analysis in the analysis unit.

請求項3の放射電磁波周波数推定装置は、請求項2の放射電磁波周波数推定装置において、信号解析部は、時系列信号のパワーが時間的に定常である場合の自己相関行列を用いて線形予測解析を行うことを特徴とする。   The radiated electromagnetic wave frequency estimation apparatus according to claim 3 is the radiated electromagnetic wave frequency estimation apparatus according to claim 2, wherein the signal analysis unit performs linear prediction analysis using an autocorrelation matrix when the power of the time series signal is stationary in time. It is characterized by performing.

請求項4の放射電磁波周波数推定装置は、請求項2の放射電磁波周波数推定装置において、時系列信号のパワーが時間的に定常でない場合の自己相関行列を用いて線形予測解析を行うことを特徴とする。   The radiated electromagnetic wave frequency estimation apparatus according to claim 4 is characterized in that, in the radiated electromagnetic wave frequency estimation apparatus according to claim 2, linear prediction analysis is performed using an autocorrelation matrix when the power of the time series signal is not temporally stationary. To do.

請求項5の放射電磁波周波数推定方法は、電気回路の配線を流れる電流を示すアナログ信号を取り込んでデジタルの時系列信号に変換するステップと、時系列信号に基づき、配線を流れる電流の時間変化に伴って放射される電磁波の周波数を、電流の時間変化から算出して推定する解析ステップを備えたことを特徴とする。   The method of estimating a frequency of radiated electromagnetic wave according to claim 5 includes a step of taking an analog signal indicating a current flowing through the wiring of the electric circuit and converting it into a digital time series signal, and a time change of the current flowing through the wiring based on the time series signal. An analysis step of calculating and estimating the frequency of the electromagnetic wave radiated from the time change of the current is provided.

請求項6の放射電磁波周波数推定方法は、請求項5の放射電磁波周波数推定方法において、解析ステップは、時系列信号を用いて線形予測解析を行う信号解析ステップと、線形予測解析により得られた予測係数を用いて配線を流れる電流の時間変化に伴って放射される電磁波の周波数を推定する周波数推定ステップを備えたことを特徴とする。   The radiated electromagnetic wave frequency estimation method according to claim 6 is the radiated electromagnetic wave frequency estimation method according to claim 5, wherein the analysis step includes a signal analysis step of performing linear prediction analysis using a time series signal, and a prediction obtained by linear prediction analysis. A frequency estimation step for estimating the frequency of an electromagnetic wave radiated with a time change of the current flowing through the wiring using a coefficient is provided.

請求項7の放射電磁波周波数推定方法は、請求項6の放射電磁波周波数推定方法において、信号解析ステップでは、時系列信号のパワーが時間的に定常である場合の自己相関行列を用いて線形予測解析が行われることを特徴とする。   The radiated electromagnetic wave frequency estimation method according to claim 7 is the radiated electromagnetic wave frequency estimation method according to claim 6, wherein in the signal analysis step, linear prediction analysis is performed using an autocorrelation matrix when the power of the time series signal is stationary in time. Is performed.

請求項8の放射電磁波周波数推定方法は、請求項6の放射電磁波周波数推定方法において、信号解析ステップでは、時系列信号のパワーが時間的に定常でない場合の自己相関行列を用いて線形予測解析が行われることを特徴とする。   The radiated electromagnetic wave frequency estimation method according to claim 8 is the radiated electromagnetic wave frequency estimation method according to claim 6, wherein in the signal analysis step, linear prediction analysis is performed using an autocorrelation matrix when the power of the time series signal is not temporally stationary. It is performed.

本発明によれば、プリント基板等に形成された電気回路上の配線を流れる電流の時間変化に伴い放射される電磁波の周波数を、電流の時間変化から推定することができる。その結果、プリント基板等に形成された電気回路の構成を改善するための指針を与えることができる。   ADVANTAGE OF THE INVENTION According to this invention, the frequency of the electromagnetic waves radiated | emitted with the time change of the electric current which flows through the wiring on the electric circuit formed in the printed circuit board etc. can be estimated from the time change of an electric current. As a result, it is possible to provide a guideline for improving the configuration of an electric circuit formed on a printed circuit board or the like.

以下、本発明の放射電磁波周波数推定装置および放射電磁波周波数推定方法の実施の形態を図面を参照しながら詳細に説明する。   Hereinafter, embodiments of a radiated electromagnetic wave frequency estimation apparatus and a radiated electromagnetic wave frequency estimation method of the present invention will be described in detail with reference to the drawings.

図1は本発明の実施例1に係る放射電磁波周波数推定装置の構成を示すブロック図である。この放射電磁波周波数推定装置は、プリント基板11と解析装置13とから構成されている。プリント基板11には、電気回路を構成する配線12が設けられており、この配線12を電流が流れる。   FIG. 1 is a block diagram showing a configuration of a radiated electromagnetic wave frequency estimation apparatus according to Embodiment 1 of the present invention. This radiated electromagnetic wave frequency estimation device is composed of a printed circuit board 11 and an analysis device 13. The printed circuit board 11 is provided with a wiring 12 constituting an electric circuit, and a current flows through the wiring 12.

解析装置13は、プローブ14、A/D変換部15、信号解析部16、周波数推定部17および表示部18を備えている。   The analysis device 13 includes a probe 14, an A / D conversion unit 15, a signal analysis unit 16, a frequency estimation unit 17, and a display unit 18.

プローブ14は、プリント基板11に設けられた配線12に接触することにより、配線12を流れる電流の時間変化を測定する。具体的には、プローブ14は、配線12に流れる電流の大きさを示すアナログ信号を継続的に取り込んでA/D変換部15に送る。   The probe 14 measures the time change of the current flowing through the wiring 12 by contacting the wiring 12 provided on the printed circuit board 11. Specifically, the probe 14 continuously takes in an analog signal indicating the magnitude of the current flowing through the wiring 12 and sends it to the A / D converter 15.

A/D変換部15は、プローブ14から送られてくるアナログ信号をデジタル信号に変換する。A/D変換部15における変換によって得られたデジタル信号は、時系列信号Sとして信号解析部16および表示部18に送られる。 The A / D converter 15 converts the analog signal sent from the probe 14 into a digital signal. Digital signal obtained by conversion in the A / D converter 15 is sent as a time-series signal S n to the signal analyzing unit 16 and a display unit 18.

信号解析部16は、例えばCPU(Central Processing Unit)で実行されるプログラムにより構成されており、A/D変換部15から送られてくる時系列信号Sを用いて線形予測解析を行う。この線形予測解析にあたり、線形予測式が含むパラメタ(予測係数)の数をM(Mの値は、必要に応じ、3〜10程度を目安に選ぶのが好ましい)とし、M個の予測係数をaとして、線形予測式を下記(1)式で表すものとする。

Figure 2010127698
Signal analyzing section 16 performs linear prediction analysis using the time-series signal S n, for example, a CPU is configured by a program executed (Central Processing Unit), sent from the A / D converter 15. In this linear prediction analysis, the number of parameters (prediction coefficients) included in the linear prediction formula is M (the value of M is preferably selected from about 3 to 10 if necessary), and M prediction coefficients are set. as a m, denote the linear prediction equation in equation (1) below.
Figure 2010127698

この(1)式に含まれる予測係数aの値を求めるために、下記(2)式で示される予測残差が最小になる条件を与えるものとする。

Figure 2010127698
To determine the value of the prediction coefficients a m included in the equation (1), shall be given the condition that prediction residual represented by the following equation (2) is minimized.
Figure 2010127698

ここで、N+1は1回の解析に用いられる時系列信号Sの数であり、必要に応じて、Mより大きな任意の値で与えられる。推定周波数f(詳細は後述)は、この1回の解析区間ごとに、M個算出される。 Here, N + 1 is the number of time-series signal S n which is used for one analysis, if necessary, be given a large arbitrary value from M. M estimated frequencies f m (details will be described later) are calculated for each analysis section.

予測残差を最小とするために、下記(3)式に示すm’に関するM個の連立方程式を満足するような予測係数aの値を求める。

Figure 2010127698
In order to minimize the prediction residual to obtain the value of the prediction coefficients a m that satisfies the M simultaneous equations relating m 'shown in the following equation (3).
Figure 2010127698

この連立方程式を解くにあたり、下記(4)式に示す、時系列信号のパワーが時間的に定常である場合の自己相関係数(係数の値が定数nの取り方によらないとする仮定。これは、周波数推定部17で副次的に算出されるλをゼロとみなすことに相当する。)を用いる。

Figure 2010127698
In solving this simultaneous equation, the autocorrelation coefficient (the value of the coefficient does not depend on how to take the constant n 0 ) when the power of the time-series signal is temporally steady as shown in the following equation (4). This is equivalent to assuming that λ m that is secondarily calculated by the frequency estimation unit 17 is zero.
Figure 2010127698

この仮定に基づき、改めて自己相関係数を表記すると下記(5)式になる。

Figure 2010127698
Based on this assumption, the autocorrelation coefficient is newly expressed as the following equation (5).
Figure 2010127698

この(5)式に示す自己相関係数を用いた、時系列信号のパワーが時間的に定常である場合の自己相関行列を下記(6)式に示す。この(6)式に示す自己相関行列を対角化することにより、簡便に予測係数aが求められる。この求められた予測係数aは、周波数推定部17に送られる。

Figure 2010127698
The autocorrelation matrix when the power of the time series signal is stationary in time using the autocorrelation coefficient shown in the equation (5) is shown in the following equation (6). By diagonalized autocorrelation matrix shown in equation (6), conveniently prediction coefficient a m is obtained. The prediction coefficients a m obtained is sent to the frequency estimation unit 17.
Figure 2010127698

周波数推定部17は、例えばCPUで実行されるプログラムにより構成されており、信号解析部16から送られてくる予測係数aを用いて、配線12を流れる電流の時間変化に伴って放射される電磁波の周波数を推定する。具体的には、予測係数aを用いたXに関する下記(7)式がXに関して因数分解される。

Figure 2010127698
Frequency estimation unit 17, for example, is constituted by a program executed by a CPU, a using the prediction coefficients a m sent from the signal analyzer section 16, is radiated with the time variation of the current flowing through the wiring 12 Estimate the frequency of electromagnetic waves. Specifically, the following equation (7) is factored with respect to X in the X using prediction coefficients a m.
Figure 2010127698

この因数分解により、下記(8)式に示すM個の係数bが得られる。

Figure 2010127698
By this factorization, M coefficients b m shown in the following equation (8) are obtained.
Figure 2010127698

この因数分解により得られM個の係数bを用い、これをXについて解くことにより(9)式に示すfとλが得られる。この解は、一般に、複素数になるため、λが現れるが、信号解析部16における自己相関行列の算出にあたり、λをゼロに仮定する近似が暗黙に用いられているので、λは使用されない。この(9)式に示すfが推定周波数であり、表示部18に送られる。

Figure 2010127698
ここで、△TはA/D変換部15におけるサンプリング間隔である。 By using M coefficients b m obtained by this factorization and solving for X, f m and λ m shown in equation (9) are obtained. Since this solution is generally a complex number, λ m appears. However, in calculating the autocorrelation matrix in the signal analysis unit 16, since an approximation assuming λ m to be zero is implicitly used, λ m is used. Not. The (9) and f m is the estimated frequency in the expression, are sent to the display unit 18.
Figure 2010127698
Here, ΔT is a sampling interval in the A / D converter 15.

表示部18は、周波数推定部17から送られてくるM個の推定周波数fを、図2(b)に示すように、複数区間に渡ってプロットすると同時に、A/D変換部15からの時系列信号Sを、図2(a)に示すように、並べてプロットして表示する。これにより、1区間N+1個の時系列信号Sのプロット毎に、M個一組の推定周波数fがプロットされることになる。 Display unit 18, the M estimated frequency f m which is sent from the frequency estimation unit 17, as shown in FIG. 2 (b), at the same time is plotted over a plurality of sections, from the A / D converter 15 the time series signal S n, as shown in FIG. 2 (a), and displays the alongside plotted. Thus, for each plot of one section (N + 1) time-series signal S n, so that the estimated frequency f m of the M set is plotted.

図3は本発明の実施例1に係る放射電磁波周波数推定装置の動作を放射電磁波周波数推定処理を中心に示すフローチャートである。   FIG. 3 is a flowchart showing the operation of the radiated electromagnetic wave frequency estimating apparatus according to the first embodiment of the present invention, centering on the radiated electromagnetic wave frequency estimating process.

この放射電磁波周波数推定処理では、まず、配線を流れる電流が取り込まれる(ステップS11)。すなわち、プローブ14は、配線12に流れる電流の大きさを示すアナログ信号を継続的に取り込んでA/D変換部15に送る。   In this radiation electromagnetic wave frequency estimation process, first, a current flowing through the wiring is captured (step S11). That is, the probe 14 continuously takes in an analog signal indicating the magnitude of the current flowing through the wiring 12 and sends it to the A / D converter 15.

次いで、A/D変換が行われる(ステップS12)。すなわち、A/D変換部15は、プローブ14から送られてくるアナログ信号をデジタル信号に変換し、時系列信号Sとして信号解析部16および表示部18に送る。 Next, A / D conversion is performed (step S12). That, A / D converter 15, an analog signal sent from the probe 14 into a digital signal, when sent to the signal analyzer 16 and a display unit 18 as a sequence signal S n.

次いで、線形予測解析が行われる(ステップS13)。すなわち、信号解析部16は、A/D変換部15から送られてくる時系列信号Sに対し、時系列信号のパワーが時間的に定常である場合の自己相関行列を用いて線形予測解析を実行し、線形予測解析により得られた予測係数aを周波数推定部17に送る。 Next, linear prediction analysis is performed (step S13). That is, the signal analyzing unit 16, a linear predictive analysis using autocorrelation matrix when to time-series signal S n sent from the A / D converter 15, the power of the time series signal is temporally constant It is executed, and sends the prediction coefficients a m obtained by linear prediction analysis to the frequency estimation unit 17.

次いで、周波数推定が行われる(ステップS14)。すなわち、周波数推定部17は、信号解析部16から送られてくる予測係数aを用いて、配線12を流れる電流の時間変化に伴って放射される電磁波の周波数を推定する。この周波数推定部17で推定された推定周波数fは表示部18に送られる。 Next, frequency estimation is performed (step S14). That is, the frequency estimation unit 17 uses the prediction coefficients a m sent from the signal analyzer section 16, it estimates the frequency of the electromagnetic waves radiated with time changes in the current flowing through the wiring 12. Estimated frequency f m estimated by the frequency estimation unit 17 is sent to the display unit 18.

次いで、表示が行われる(ステップS15)。すなわち、表示部18は、図2に示すように、周波数推定部17から送られてくるM個の推定周波数fを複数区間に渡ってプロットすると同時に、A/D変換部15からの時系列信号Sを並べてプロットして表示する。その後、放射電磁波周波数推定処理は終了する。 Next, display is performed (step S15). That is, the display unit 18, as shown in FIG. 2, at the same time is plotted across the M estimated frequency f m which is sent from the frequency estimation unit 17 in plurality of sections, when the A / D converter 15 sequence It plotted to display side-by-side signal S n. Thereafter, the radiated electromagnetic wave frequency estimation process ends.

以上説明したように、本発明の実施例1に係る放射電磁波周波数推定装置によれば、特許文献1に記載された方法等によって別途スペクトルアナライザを用いて取得されたパワースペクトルから、対策すべき不用電磁波の周波数情報を取得し、推定周波数fが、その周波数に一致する時刻をプロット上で確認する。その時刻前後における時系列信号Sの波形(電気回路上の電流変化)により、不要電磁波が放射される。このタイミングの波形を、電気回路の構成を変更することによって改善し、不用電磁波の放射を抑制することができる。 As described above, according to the radiated electromagnetic wave frequency estimating apparatus according to the first embodiment of the present invention, it is unnecessary to take measures from the power spectrum separately obtained by using the spectrum analyzer by the method described in Patent Document 1. to obtain frequency information of the electromagnetic wave, the estimated frequency f m is, to check the time that matches the frequency on the plot. The time series signal S n of the waveform (current change on the electric circuit) at the time before and after, unnecessary electromagnetic waves are radiated. The waveform of this timing can be improved by changing the configuration of the electric circuit, and the emission of unnecessary electromagnetic waves can be suppressed.

図4は本発明の実施例2に係る放射電磁波周波数推定装置の構成を示すブロック図である。この放射電磁波周波数推定装置は、実施例1に係る放射電磁波周波数推定装置の解析装置13に含まれる信号解析部16の機能が変更され、他の信号解析部21に置き換えられて構成されている。   FIG. 4 is a block diagram showing the configuration of the radiated electromagnetic wave frequency estimation apparatus according to Embodiment 2 of the present invention. This radiated electromagnetic wave frequency estimation device is configured by changing the function of the signal analysis unit 16 included in the analysis device 13 of the radiated electromagnetic wave frequency estimation device according to the first embodiment and replacing it with another signal analysis unit 21.

以下においては、実施例1に係る放射電磁波周波数推定装置と同じ構成要素には、実施例1で使用した符号と同一の符号を付して説明を省略または簡略化し、異なる部分を中心に説明する。   In the following, the same components as those in the radiated electromagnetic wave frequency estimation apparatus according to the first embodiment are denoted by the same reference numerals as those used in the first embodiment, and the description thereof will be omitted or simplified, and different portions will be mainly described. .

信号解析部21は、例えばCPUで実行されるプログラムにより構成されており、A/D変換部15から送られてくる時系列信号Sを用いて線形予測解析を実行する。この線形予測解析にあたり、線形予測式が含むパラメタ(予測係数)の数をM(Mの値は、必要に応じ、3〜10程度を目安に選ぶのが好ましい)とし、M個の予測係数をaとして、線形予測式を上記(1)式で表すものとする。 Signal analyzer 21, for example, is constituted by a program executed by the CPU, and performs linear prediction analysis using the time-series signal S n sent from the A / D converter 15. In this linear prediction analysis, the number of parameters (prediction coefficients) included in the linear prediction formula is M (the value of M is preferably selected from about 3 to 10 if necessary), and M prediction coefficients are set. as a m, denote the linear prediction equation in equation (1).

この(1)式に含まれる予測係数aの値を求めるために、上記(2)式で示される予測残差が最小になる条件を与えるものとする。(2)式で用いられるN+1は1回の解析に用いる時系列信号Sの数であり、必要に応じて、Mより大きな任意の値で与える。推定周波数fは、この1回の解析区間ごとに、M個算出される。 To determine the value of the prediction coefficients a m included in the equation (1), shall be given the condition that prediction residual represented by the formula (2) is minimized. (2) the number of time-series signal S n to be used for N + 1 is one of analysis used in the formula, if necessary, give a large arbitrary value from M. M estimated frequencies f m are calculated for each analysis section.

予測残差を最小とするために、上記(3)式に示すm’に関するM個の連立方程式を満足するような予測係数aの値を求める。この連立方程式を解くにあたり、下記(10)式に示す時系列信号のパワーが時間的に定常でない(発散や減衰を伴う)場合の自己相関係数(以下、「時間依存自己相関係数」という)を用いた自己相関行列、つまり、下記(11)式に示す、時系列信号のパワーが時間的に定常でない場合の自己相関行列(以下、「時間依存自己相関行列」という)が使用される。 In order to minimize the prediction residual to obtain the value of the prediction coefficients a m that satisfies the M simultaneous equations relating m 'shown in the above (3). In solving this simultaneous equation, the autocorrelation coefficient (hereinafter referred to as “time-dependent autocorrelation coefficient”) when the power of the time series signal shown in the following equation (10) is not temporally steady (with divergence and attenuation) ), That is, an autocorrelation matrix (hereinafter referred to as “time-dependent autocorrelation matrix”) when the power of the time-series signal is not temporally steady, as shown in the following equation (11). .

この時間依存自己相関行列を対角化することにより、予測係数aが求められる。この時間依存自己相関行列を用いることにより、より正確に周波数fの推定を行うことができる。以上のようにして信号解析部21で求められた予測係数aは、周波数推定部17に送られる。

Figure 2010127698
By diagonalizing the time-dependent autocorrelation matrix, the prediction coefficient a m is obtained. By using the time-dependent autocorrelation matrix, it is possible to more accurately estimate the frequency f m. Prediction coefficients a m obtained by the signal analyzer 21 as described above is sent to the frequency estimation unit 17.
Figure 2010127698

周波数推定部17は、実施例1で説明したように、信号解析部16から送られてくる予測係数aを用いて、配線12を流れる電流の時間変化に伴って放射される電磁波の周波数を推定する。具体的には、予測係数aを用いたXに関する上記(7)式をXに関して因数分解する。この因数分解により、上記(8)式に示すM個の係数bが得られる。 Frequency estimation unit 17, as described in Example 1, using the prediction coefficients a m sent from the signal analyzer section 16, the frequency of the electromagnetic waves radiated with time changes in the current flowing through the wiring 12 presume. Specifically, factoring in the X above (7) about the X using prediction coefficients a m. By this factorization, M coefficients b m shown in the above equation (8) are obtained.

この因数分解により得られM個の係数bを用い、これをXについて解くことにより上記(9)式に示すfとλが得られる。この解は、一般に、複素数になるため、λが現れる。λは、放射波の減衰率に対応し、減衰率が高いほど、ブロードなスペクトルを形成する。この(9)式に示すfが推定周波数であり、表示部18に送られる。 Using M coefficients b m obtained by this factorization and solving for X, f m and λ m shown in the above equation (9) are obtained. Since this solution is generally complex, λ m appears. λ m corresponds to the attenuation rate of the radiated wave, and the higher the attenuation rate, the broader the spectrum. The (9) and f m is the estimated frequency in the expression, are sent to the display unit 18.

図5は本発明の実施例2に係る放射電磁波周波数推定装置の動作を、放射電磁波周波数推定処理を中心に、フローチャートで示したものである。以下では、図3のフローチャートに示した実施例1に係る放射電磁波周波数推定装置の放射電磁波周波数推定処理と同じ処理を行うステップには、図3で使用した符号と同一の符号を付して説明を簡略化し、異なる部分を中心に説明する。   FIG. 5 is a flowchart showing the operation of the radiated electromagnetic wave frequency estimation apparatus according to the second embodiment of the present invention, focusing on the radiated electromagnetic wave frequency estimation process. In the following, the same reference numerals as those used in FIG. 3 are attached to the steps for performing the same processing as the radiated electromagnetic wave frequency estimation process of the radiated electromagnetic wave frequency estimation apparatus according to the first embodiment shown in the flowchart of FIG. The explanation will be focused on the different parts.

この放射電磁波周波数推定処理では、まず、配線を流れる電流が取り込まれる(ステップS11)。次いで、A/D変換が行われる(ステップS12)。次いで、線形予測解析が行われる(ステップS21)。すなわち、信号解析部16は、A/D変換部15から送られてくる時系列信号Sに対し、時系列信号のパワーが時間的に定常でない場合の自己相関行列(時間依存自己相関行列)を用いて線形予測解析を実行し、この線形予測解析により得られた予測係数aを周波数推定部17に送る。 In this radiation electromagnetic wave frequency estimation process, first, a current flowing through the wiring is captured (step S11). Next, A / D conversion is performed (step S12). Next, linear prediction analysis is performed (step S21). That is, the signal analyzing unit 16, with respect to the time-series signal S n sent from the A / D converter 15, the autocorrelation matrix when the power of the time-series signal is not stationary in time (time-dependent autocorrelation matrix) performs linear predictive analysis using sends the prediction coefficients a m obtained by the linear predictive analysis to the frequency estimation unit 17.

次いで、周波数推定が行われる(ステップS14)。次いで、表示が行われる(ステップS15)。その後、放射電磁波周波数推定処理は終了する。   Next, frequency estimation is performed (step S14). Next, display is performed (step S15). Thereafter, the radiated electromagnetic wave frequency estimation process ends.

以上説明したように、本発明の実施例2に係る放射電磁波周波数推定装置によれば、特許文献1に記載された方法等によって別途スペクトルアナライザを用いて取得されたパワースペクトルから、対策すべき不用電磁波の周波数情報を取得し、推定周波数fが、その周波数に一致する時刻をプロット上で確認する。その時刻前後における時系列信号Sの波形(電気回路上の電流変化)により、不要電磁波が放射される。このタイミングの波形を、電気回路の構成を変更することによって改善し、不用電磁波の放射を抑制することができる。 As described above, according to the radiated electromagnetic wave frequency estimation apparatus according to the second embodiment of the present invention, it is unnecessary to take countermeasures from the power spectrum separately obtained by using the spectrum analyzer by the method described in Patent Document 1. to obtain frequency information of the electromagnetic wave, the estimated frequency f m is, to check the time that matches the frequency on the plot. The time series signal S n of the waveform (current change on the electric circuit) at the time before and after, unnecessary electromagnetic waves are radiated. The waveform of this timing can be improved by changing the configuration of the electric circuit, and the emission of unnecessary electromagnetic waves can be suppressed.

本発明は、電気回路から放射される電磁波の抑圧が要求される種々の装置に利用可能である。   The present invention is applicable to various devices that are required to suppress electromagnetic waves radiated from an electric circuit.

本発明の実施例1に係る放射電磁波周波数推定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radiation electromagnetic wave frequency estimation apparatus which concerns on Example 1 of this invention. 本発明の実施例1に係る放射電磁波周波数推定装置の表示部に表示される時系列信号および推定周波数の表示例を示す図である。It is a figure which shows the example of a display of the time series signal and estimated frequency which are displayed on the display part of the radiation electromagnetic wave frequency estimation apparatus which concerns on Example 1 of this invention. 本発明の実施例1に係る放射電磁波周波数推定装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the radiation electromagnetic wave frequency estimation apparatus which concerns on Example 1 of this invention. 本発明の実施例2に係る放射電磁波周波数推定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radiation electromagnetic wave frequency estimation apparatus which concerns on Example 2 of this invention. 本発明の実施例2に係る放射電磁波周波数推定装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the radiation electromagnetic wave frequency estimation apparatus which concerns on Example 2 of this invention.

符号の説明Explanation of symbols

11 プリント基板
12 配線
13 解析装置
14 プローブ
15 A/D変換部
16、21 信号解析部
17 周波数推定部
18 表示部
DESCRIPTION OF SYMBOLS 11 Printed circuit board 12 Wiring 13 Analysis apparatus 14 Probe 15 A / D conversion part 16, 21 Signal analysis part 17 Frequency estimation part 18 Display part

Claims (8)

電気回路の配線を流れる電流を示すアナログ信号を取り込んでデジタルの時系列信号に変換するA/D変換部と、
前記A/D変換部からの時系列信号に基づき、前記配線を流れる電流の時間変化に伴って放射される電磁波の周波数を、電流の時間変化から算出して推定する解析装置と、
を備えたことを特徴とする放射電磁波周波数推定装置。
An A / D converter that takes in an analog signal indicating the current flowing through the wiring of the electric circuit and converts it into a digital time-series signal;
Based on the time-series signal from the A / D conversion unit, an analysis device that calculates and estimates the frequency of the electromagnetic wave radiated with time change of the current flowing through the wiring from the time change of current,
A radiated electromagnetic wave frequency estimation device comprising:
前記解析装置は、
前記A/D変換部からの時系列信号を用いて線形予測解析を行う信号解析部と、
前記信号解析部における線形予測解析により得られた予測係数を用いて前記配線を流れる電流の時間変化に伴って放射される電磁波の周波数を推定する周波数推定部と、
を備えたことを特徴とする請求項1記載の放射電磁波周波数推定装置。
The analysis device includes:
A signal analysis unit that performs linear prediction analysis using a time-series signal from the A / D conversion unit;
A frequency estimation unit that estimates a frequency of an electromagnetic wave radiated with a time change of a current flowing through the wiring using a prediction coefficient obtained by linear prediction analysis in the signal analysis unit;
The radiated electromagnetic wave frequency estimation apparatus according to claim 1, comprising:
前記信号解析部は、前記時系列信号のパワーが時間的に定常である場合の自己相関行列を用いて線形予測解析を行うことを特徴とする請求項2記載の放射電磁波周波数推定装置。   The radiated electromagnetic wave frequency estimation apparatus according to claim 2, wherein the signal analysis unit performs linear prediction analysis using an autocorrelation matrix when the power of the time series signal is stationary in time. 前記信号解析部は、前記時系列信号のパワーが時間的に定常でない場合の自己相関行列を用いて線形予測解析を行うことを特徴とする請求項2記載の放射電磁波周波数推定装置。   The radiated electromagnetic wave frequency estimation apparatus according to claim 2, wherein the signal analysis unit performs linear prediction analysis using an autocorrelation matrix when the power of the time series signal is not stationary in time. 電気回路の配線を流れる電流を示すアナログ信号を取り込んでデジタルの時系列信号に変換するステップと、
前記時系列信号に基づき、前記配線を流れる電流の時間変化に伴って放射される電磁波の周波数を、電流の時間変化から算出して推定する解析ステップ
と、
を備えたことを特徴とする放射電磁波周波数推定方法。
Taking an analog signal indicating the current flowing through the wiring of the electric circuit and converting it into a digital time-series signal;
Based on the time series signal, an analysis step for calculating and estimating the frequency of the electromagnetic wave radiated with the time change of the current flowing through the wiring from the time change of the current;
A method for estimating the frequency of radiated electromagnetic waves, comprising:
前記解析ステップは、
前記時系列信号を用いて線形予測解析を行う信号解析ステップと、
前記線形予測解析により得られた予測係数を用いて前記配線を流れる電流の時間変化に伴って放射される電磁波の周波数を推定する周波数推定ステップと、
を備えたことを特徴とする請求項5記載の放射電磁波周波数推定方法。
The analysis step includes
A signal analysis step of performing linear prediction analysis using the time series signal;
A frequency estimation step for estimating a frequency of an electromagnetic wave radiated with a time change of a current flowing through the wiring using a prediction coefficient obtained by the linear prediction analysis;
The radiated electromagnetic wave frequency estimation method according to claim 5, further comprising:
前記信号解析ステップでは、前記時系列信号のパワーが時間的に定常である場合の自己相関行列を用いて線形予測解析が行われることを特徴とする請求項6記載の放射電磁波周波数推定方法。   7. The method of estimating a radiated electromagnetic wave frequency according to claim 6, wherein in the signal analysis step, linear prediction analysis is performed using an autocorrelation matrix when the power of the time series signal is stationary in time. 前記信号解析ステップでは、前記時系列信号のパワーが時間的に定常でない場合の自己相関行列を用いて線形予測解析が行われることを特徴とする請求項6記載の放射電磁波周波数推定方法。   7. The method of estimating a radiated electromagnetic wave frequency according to claim 6, wherein, in the signal analysis step, linear prediction analysis is performed using an autocorrelation matrix when the power of the time series signal is not temporally stationary.
JP2008301061A 2008-11-26 2008-11-26 Apparatus and method for estimation of radiated electromagnetic wave frequency Pending JP2010127698A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008301061A JP2010127698A (en) 2008-11-26 2008-11-26 Apparatus and method for estimation of radiated electromagnetic wave frequency

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008301061A JP2010127698A (en) 2008-11-26 2008-11-26 Apparatus and method for estimation of radiated electromagnetic wave frequency

Publications (1)

Publication Number Publication Date
JP2010127698A true JP2010127698A (en) 2010-06-10

Family

ID=42328217

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008301061A Pending JP2010127698A (en) 2008-11-26 2008-11-26 Apparatus and method for estimation of radiated electromagnetic wave frequency

Country Status (1)

Country Link
JP (1) JP2010127698A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012098255A (en) * 2010-11-05 2012-05-24 Nippon Telegr & Teleph Corp <Ntt> Superimposed mitigation signal analysis apparatus, superimposed mitigation signal analysis method and superimposed mitigation signal analysis program
JP2013120253A (en) * 2011-12-07 2013-06-17 Nippon Telegr & Teleph Corp <Ntt> Spectrum vector storage device, spectrum vector storage method and spectrum vector storage program

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09222442A (en) * 1995-04-05 1997-08-26 Nippon Telegr & Teleph Corp <Ntt> Specifying method of electronmagnetic interference entering route, electromagnetic interference search device, and non-contact voltage probe device
JPH10185973A (en) * 1996-12-25 1998-07-14 Matsushita Electric Works Ltd Method and device for measuring electromagnetic interference of circuit substrate
JPH10307159A (en) * 1997-05-08 1998-11-17 Nec Corp Method and device for evaluating radiation
JPH11190752A (en) * 1997-12-26 1999-07-13 Matsushita Electric Ind Co Ltd Method and apparatus for evaluating electromagnetic radiation from electric apparatus and countermeasure therefor
JP2005249967A (en) * 2004-03-02 2005-09-15 Nippon Telegr & Teleph Corp <Ntt> Method and apparatus for frequency analysis
JP2008197873A (en) * 2007-02-13 2008-08-28 Nippon Telegr & Teleph Corp <Ntt> Spectral distribution and statistical distribution analysis method, spectral distribution and statistical distribution analysis device, and spectral distribution and statistical distribution analysis program

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09222442A (en) * 1995-04-05 1997-08-26 Nippon Telegr & Teleph Corp <Ntt> Specifying method of electronmagnetic interference entering route, electromagnetic interference search device, and non-contact voltage probe device
JPH10185973A (en) * 1996-12-25 1998-07-14 Matsushita Electric Works Ltd Method and device for measuring electromagnetic interference of circuit substrate
JPH10307159A (en) * 1997-05-08 1998-11-17 Nec Corp Method and device for evaluating radiation
JPH11190752A (en) * 1997-12-26 1999-07-13 Matsushita Electric Ind Co Ltd Method and apparatus for evaluating electromagnetic radiation from electric apparatus and countermeasure therefor
JP2005249967A (en) * 2004-03-02 2005-09-15 Nippon Telegr & Teleph Corp <Ntt> Method and apparatus for frequency analysis
JP2008197873A (en) * 2007-02-13 2008-08-28 Nippon Telegr & Teleph Corp <Ntt> Spectral distribution and statistical distribution analysis method, spectral distribution and statistical distribution analysis device, and spectral distribution and statistical distribution analysis program

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JPN6012058849; 村川 外3名: 'ノイズ波形に対する相関処理及び積分操作の組合せによる線形予測法の周波数抽出精度の改善' 電子情報通信学会論文誌B Vol.J87-B,No.2, 200402, pp.235-245 *
JPN6012058850; 村川 外2名: '複素周波数を用いた線形予測法によるノイズ識別法について' 電子情報通信学会通信ソサイエティ大会講演論文集 , 2003, p.331 *
JPN6012058851; 石山文彦: '広帯域複素周波数解析装置の構築とインバータ系ノイズの測定' 電子情報通信学会総合大会講演論文集 通信1 , 2006, p.309 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012098255A (en) * 2010-11-05 2012-05-24 Nippon Telegr & Teleph Corp <Ntt> Superimposed mitigation signal analysis apparatus, superimposed mitigation signal analysis method and superimposed mitigation signal analysis program
JP2013120253A (en) * 2011-12-07 2013-06-17 Nippon Telegr & Teleph Corp <Ntt> Spectrum vector storage device, spectrum vector storage method and spectrum vector storage program

Similar Documents

Publication Publication Date Title
US20160238460A1 (en) Temperature measuring system and temperature measuring method
CA2821462C (en) An active sonar apparatus with a noise elimination unit
JP2006266907A (en) Radar device and its radar signal processing method
US11275005B2 (en) Fatigue limit stress specification system, fatigue limit stress specification device, and fatigue limit stress specification method
RU2473917C2 (en) Method and device to measure electromagnet wave
JP2010127698A (en) Apparatus and method for estimation of radiated electromagnetic wave frequency
JP6288945B2 (en) Noise source position estimation apparatus and noise source position estimation program
JP2003323466A (en) Design support system
CN115047022A (en) Time domain reconstruction method and system for thermal diffusion process
JP2006242672A (en) Electromagnetic wave measuring instrument
JP2013080787A (en) Life estimation device, life estimation method and program
JP5756826B2 (en) Server device and radioactivity measurement system
JP2010243459A (en) Fluorescence measuring apparatus
CN112991260A (en) Infrared nondestructive testing system with light and ultrasonic composite excitation
JP5023029B2 (en) Radar equipment
JP5150386B2 (en) Electromagnetic noise diagnostic device, electromagnetic noise diagnostic system, and electromagnetic noise diagnostic method
WO2018062323A1 (en) Electronic device, control method, and program
JP2010206244A (en) Image degradation evaluation system, image degradation evaluation method, and image degradation evaluation program
JP2006084413A (en) Source identifying method for electromagnetic interference signal
JP2020041823A (en) Inspection method, inspection apparatus, inspection system, and computer program
EP4321881A1 (en) Insulation degradation diagnosis model creation device, insulation degradation diagnosis device, and insulation degradation diagnosis method
KR20190011167A (en) System and method for predicting thickness of damage of target object
KR101294405B1 (en) Method for voice activity detection using phase shifted noise signal and apparatus for thereof
JP2009002680A (en) Method of measuring electric field or magnetic field
WO2012102321A1 (en) System for analyzing behavior of ions in insulation film

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20101216

RD03 Notification of appointment of power of attorney

Effective date: 20120509

Free format text: JAPANESE INTERMEDIATE CODE: A7423

RD04 Notification of resignation of power of attorney

Effective date: 20120511

Free format text: JAPANESE INTERMEDIATE CODE: A7424

A977 Report on retrieval

Effective date: 20121009

Free format text: JAPANESE INTERMEDIATE CODE: A971007

A131 Notification of reasons for refusal

Effective date: 20121113

Free format text: JAPANESE INTERMEDIATE CODE: A131

A02 Decision of refusal

Effective date: 20130312

Free format text: JAPANESE INTERMEDIATE CODE: A02