JP2021113792A - Effective value calculator - Google Patents
Effective value calculator Download PDFInfo
- Publication number
- JP2021113792A JP2021113792A JP2020007774A JP2020007774A JP2021113792A JP 2021113792 A JP2021113792 A JP 2021113792A JP 2020007774 A JP2020007774 A JP 2020007774A JP 2020007774 A JP2020007774 A JP 2020007774A JP 2021113792 A JP2021113792 A JP 2021113792A
- Authority
- JP
- Japan
- Prior art keywords
- sample
- signal
- effective value
- zero
- outputs
- 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.)
- Granted
Links
- 238000005070 sampling Methods 0.000 claims abstract description 65
- 238000001514 detection method Methods 0.000 claims abstract description 31
- 238000012360 testing method Methods 0.000 description 38
- 230000006870 function Effects 0.000 description 16
- 238000000605 extraction Methods 0.000 description 11
- 230000000052 comparative effect Effects 0.000 description 5
- 230000001934 delay Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
Landscapes
- Measurement Of Current Or Voltage (AREA)
Abstract
Description
本開示は、実効値算出装置に関する。 The present disclosure relates to an effective value calculation device.
従来、入力信号を所定のサンプリング間隔でサンプルし、入力信号の半波毎の実効値としてサンプル値の二乗平均平方根を算出する実効値演算回路が知られている(例えば、特許文献1参照)。 Conventionally, there is known an effective value calculation circuit that samples an input signal at a predetermined sampling interval and calculates the root mean square of the sample value as an effective value for each half wave of the input signal (see, for example, Patent Document 1).
実効値の算出精度の向上が求められる。 It is required to improve the calculation accuracy of the effective value.
本開示は、上述の点に鑑みてなされたものであり、実効値の算出精度を向上できる実効値算出装置を提供することを目的とする。 The present disclosure has been made in view of the above points, and an object of the present disclosure is to provide an effective value calculation device capable of improving the calculation accuracy of the effective value.
幾つかの実施形態に係る実効値算出装置は、入力信号を所定のサンプリング間隔でサンプルしたサンプル値を含むサンプル信号を出力するAD変換器と、前記サンプル信号の高周波数成分を減衰させてゼロクロス検出用信号として出力するゼロクロスフィルタと、前記ゼロクロス検出用信号に基づいて検出された連続する2つのゼロクロス点の間のカウント期間に含まれる前記サンプル値の数をカウントし、サンプリングポイント数として出力するサンプル数カウンタと、前記サンプル信号に含まれる前記サンプル値のうち、前記サンプル数カウンタから取得した前記サンプリングポイント数だけ連続する前記サンプル値の二乗平均平方根を算出し、前記入力信号の実効値として出力する演算器とを備える。このようにすることで、ゼロクロス点の検出精度が高められる。その結果、実効値算出装置は、実効値の算出精度を向上できる。 The effective value calculation device according to some embodiments is an AD converter that outputs a sample signal including a sample value obtained by sampling an input signal at a predetermined sampling interval, and zero-cross detection by attenuating a high frequency component of the sample signal. A sample that counts the number of sample values included in the count period between the zero cross filter output as a signal for zero cross and two consecutive zero cross points detected based on the signal for zero cross detection, and outputs it as the number of sampling points. Of the number counter and the sample value included in the sample signal, the squared average square root of the sample value continuous with the number of sampling points acquired from the sample number counter is calculated and output as an effective value of the input signal. It is equipped with a calculator. By doing so, the detection accuracy of the zero cross point is improved. As a result, the effective value calculation device can improve the calculation accuracy of the effective value.
一実施形態に係る実効値算出装置において、前記サンプル数カウンタは、少なくとも2つの前記カウント期間それぞれに含まれる前記サンプル値の数をカウントし、前記各カウント期間に含まれる前記サンプル値の数の平均値を前記サンプリングポイント数として出力してよい。このようにすることで、実効値を算出するための二乗平均の演算における除数となるサンプリングポイント数のカウントの精度が高められる。その結果、実効値算出装置は、実効値の算出精度を向上できる。 In the effective value calculation device according to the embodiment, the sample number counter counts the number of the sample values included in each of at least two counting periods, and averages the number of the sample values included in each counting period. The value may be output as the number of sampling points. By doing so, the accuracy of counting the number of sampling points, which is a divisor in the calculation of the root mean square for calculating the effective value, is improved. As a result, the effective value calculation device can improve the calculation accuracy of the effective value.
一実施形態に係る実効値算出装置において、前記サンプル数カウンタは、前記サンプリングポイント数を逐次更新する更新モード、及び、前記サンプリングポイント数を保持するホールドモードのうち一方のモードで動作してよい。このようにすることで、サンプリング値の二乗平均を算出する際の除数となるサンプリングポイント数に対して、入力信号の振幅の変化が及ぼす影響が低減され得る。その結果、実効値算出装置は、実効値の算出精度を向上できる。 In the effective value calculation device according to the embodiment, the sample number counter may operate in one of an update mode in which the number of sampling points is sequentially updated and a hold mode in which the number of sampling points is held. By doing so, the influence of the change in the amplitude of the input signal on the number of sampling points, which is the divisor when calculating the root mean square of the sampling values, can be reduced. As a result, the effective value calculation device can improve the calculation accuracy of the effective value.
本開示によれば、実効値の算出精度を向上できる実効値算出装置が提供される。 According to the present disclosure, an effective value calculation device capable of improving the calculation accuracy of the effective value is provided.
本開示に係る実施形態が、比較例と対比しながら説明される。 An embodiment according to the present disclosure will be described in comparison with a comparative example.
(比較例)
図1に示されるように、比較例に係る実効値算出装置900は、AD変換器910と、演算器960と、ゼロクロス検出器970と、サンプル数カウンタ980とを備える。演算器960は、2乗回路920と、総和回路930と、除算回路940と、開平回路950とを備える。
(Comparison example)
As shown in FIG. 1, the effective
AD変換器910は、入力されるアナログ信号をサンプリングクロックに基づいてデジタル信号に変換し、ゼロクロス検出器970と演算器960とに出力する。デジタル信号は、サンプリングしたタイミングとそのときのアナログ信号の値であるサンプル値とを対応付けたデータのセットとして表される。
The
ゼロクロス検出器970は、取得したデジタル信号のゼロクロス点を検出し、検出したゼロクロス点を演算器960とサンプル数カウンタ980とに出力する。ゼロクロス点は、デジタル信号の符号が反転するタイミングに対応する。2つの連続するゼロクロス点は、1つの半波の開始点と終了点に対応する。
The zero-
サンプル数カウンタ980は、2つのゼロクロス点の間、つまり1つの半波の期間に入力されるサンプル信号の数をカウントし、演算器960に出力する。
The
演算器960の各構成部は、以下説明するように機能する。2乗回路920は、サンプル値を2乗し、総和回路930に出力する。総和回路930は、ゼロクロス検出器970からゼロクロス点を取得し、ゼロクロス点を半波の開始点又は終了点として1つの半波の期間を特定し、その期間内において2乗回路920の出力値を累積加算する。つまり、総和回路930は、アナログ信号の半波毎に2乗回路920の出力値を累積加算する。総和回路930は、累積加算した値を除算回路940に出力する。除算回路940は、サンプル数カウンタ980から1つの半波の期間に入力されるサンプル信号の数のカウント値を取得し、総和回路930の出力値を、取得したカウント値で除算し、開平回路950に出力する。開平回路950は、除算回路940の出力値を開平する。つまり、開平回路950は、除算回路940の出力値の平方根を算出する。開平回路950は、算出した値をアナログ信号の実効値として出力する。
Each component of the
ここで、アナログ信号が入力される場合における実効値算出装置900の動作が説明される。図2に例示される波形を有するアナログ信号が実効値算出装置900に入力されるとする。アナログ信号が正の半波となっている期間は、P1及びP2として表されている。P1の期間において、アナログ信号は、ノイズを含まない。その結果、ゼロクロス点は、期間の最初と最後で検出される。したがって、D1で表されている、P1において半波として検出される長さは、P1の期間の長さに一致する。
Here, the operation of the effective
一方で、P2の期間において、アナログ信号は、二点鎖線で示される真の波形に対して小さい値となっているノイズ波形を含むとする。ノイズ波形は、符号が反転する点を含むとする。この場合、ゼロクロス点は、アナログ信号が正の半波となる期間(P2)の最初と最後だけでなく途中のノイズ波形の部分でも検出される。P2において半波として検出される長さは、D2で表されており、P2の期間より短くなっている。 On the other hand, during the period of P2, it is assumed that the analog signal includes a noise waveform whose value is smaller than the true waveform indicated by the alternate long and short dash line. It is assumed that the noise waveform includes a point where the sign is inverted. In this case, the zero cross point is detected not only at the beginning and the end of the period (P2) when the analog signal becomes a positive half wave, but also at the part of the noise waveform in the middle. The length detected as a half wave in P2 is represented by D2, which is shorter than the period of P2.
ゼロクロス検出器970は、検出したゼロクロス点を、総和回路930に出力する。総和回路930は、1つのゼロクロス点から次のゼロクロス点までの間に含まれるサンプル値の2乗値を累積加算する。一方で、ゼロクロス検出器970は、検出したゼロクロス点を、サンプル数カウンタ980にも出力する。サンプル数カウンタ980は、1つのゼロクロス点から次のゼロクロス点までの間に含まれるサンプル信号の数をカウントし、除算回路940に出力する。除算回路940は、総和回路930の出力値を、サンプル数カウンタ980が出力したカウント値で除算する。
The zero-
P1の期間において、実効値算出装置900がゼロクロス点を検出する間隔は、アナログ信号の半波の長さに一致する。これによって、実効値算出装置900は、入力されたアナログ信号の実効値を本来の半波の長さで算出できる。一方で、P2の期間において、実効値算出装置900がゼロクロス点を検出する間隔は、アナログ信号の実際の半波の長さと異なる。これによって、実効値算出装置900は、入力されたアナログ信号の実効値を本来の半波の長さと異なる長さで算出する。その結果、実効値算出装置900によるアナログ信号の実効値の算出精度が低下する。
In the period of P1, the interval at which the effective
以上述べてきたように、比較例に係る実効値算出装置900は、ノイズ波形を含むアナログ信号が入力された場合にゼロクロス点を誤って検出する可能性が高く、アナログ信号の実効値の算出精度が低下するという課題を有している。
As described above, the effective
そこで、本開示は、ゼロクロス点の検出誤差を低減できる実効値算出装置1(図3等参照)を説明する。 Therefore, the present disclosure describes an effective value calculation device 1 (see FIG. 3 and the like) that can reduce the detection error of the zero cross point.
(本開示の一実施形態)
図3に示されるように、一実施形態に係る実効値算出装置1は、AD変換器10と、ゼロクロスフィルタ20と、ゼロクロス検出器25と、サンプル数カウンタ30と、ディレイ回路50と、演算器60とを備える。図4に示されるように、演算器60は、2乗回路62と、総和回路64と、除算回路66と、開平回路68とを備える。実効値算出装置1は、必須ではないが、ラインフィルタ40を更に備える。実効値算出装置1は、入力される信号の半波の長さを検出して、その信号の実効値を算出する。
(One Embodiment of the present disclosure)
As shown in FIG. 3, the effective value calculation device 1 according to the embodiment includes an
実効値算出装置1は、必須ではないが、制御部70を更に備える。制御部70は、信号の実効値の算出結果を取得してよい。制御部70は、取得した算出結果を外部装置へ出力してもよい。制御部70は、後述するように、実効値算出装置1の各構成部の動作を定めるパラメータを設定してもよい。
The effective value calculation device 1 is not essential, but further includes a
制御部70は、例えばCPU(Central Processing Unit)等のプロセッサを含んで構成されてよい。制御部70は、プロセッサに所定のプログラムを実行させることによって所定の機能を実現してもよい。制御部70は、記憶部を備えてもよい。記憶部は、制御部70の動作に用いられる各種情報、又は、制御部70の機能を実現するためのプログラム等を格納してよい。記憶部は、制御部70のワークメモリとして機能してよい。記憶部は、例えば半導体メモリ等で構成されてよい。記憶部は、制御部70に含まれてもよいし、制御部70と別体として構成されてもよい。
The
実効値算出装置1は、必須ではないが、操作部72を更に備える。操作部72は、ユーザの操作入力を受け付けることができるように構成される。操作部72は、例えば、マウス等のポインティングデバイス、物理キー、又はタッチパネル等の入力デバイスを含んでよい。制御部70は、操作部72によってユーザからの操作入力を受け付け、ユーザの操作に基づいて、実効値算出装置1の各構成部のパラメータを設定してもよい。
The effective value calculation device 1 is not essential, but further includes an
実効値算出装置1は、必須ではないが、表示部74を更に備える。表示部74は、ユーザに報知する情報を表示するディスプレイ又は発光素子等の表示デバイスを含んでよい。制御部70は、ユーザに実効値の算出結果を知らせるために、信号の実効値の算出結果を表示部74に表示させてもよい。
The effective value calculation device 1 is not essential, but further includes a
(実効値算出装置1の動作例)
AD変換器10は、入力されるアナログ信号をサンプリングクロックに基づいてデジタル信号に変換する。AD変換器10は、サンプリングクロックを外部のクロック発生回路から取得してよい。サンプリングクロックの周波数は、入力されるアナログ信号の周波数よりも十分に高い周波数であるとする。実効値算出装置1に入力されるアナログ信号は、所定のピーク周波数を有する。実効値算出装置1に入力されるアナログ信号は、単一の周波数を有する正弦波であってよい。
(Operation example of effective value calculation device 1)
The
AD変換器10に入力されるアナログ信号は、入力信号とも称される。AD変換器10は、所定のサンプリング間隔で入力信号をサンプルする。サンプリング間隔は、サンプリングクロックの周波数に基づいて定まる。AD変換器10が入力信号をサンプルして得られる値は、サンプル値とも称される。AD変換器10は、サンプル値に対応する信号を順次出力することによって、サンプル値を含むデジタル信号を出力する。サンプル値を含むデジタル信号は、サンプル信号とも称される。
The analog signal input to the
AD変換器10は、サンプル信号を、ゼロクロスフィルタ20と、ラインフィルタ40とに出力する。実効値算出装置1がラインフィルタ40を備えない場合、AD変換器10は、サンプル信号を、ディレイ回路50に出力する。サンプル信号は、サンプリング番号(i)と、各サンプリング番号に対応するサンプル値とを含むデータとして表される。サンプリング番号(i)に対応するサンプル値は、u(i)と表されるとする。
The
ゼロクロスフィルタ20は、デジタル信号としてのサンプル信号をフィルタリングするデジタルフィルタとして構成される。ゼロクロスフィルタ20は、第1カットオフ周波数以上の高周波数成分を減衰させるローパスフィルタとして機能する。ゼロクロスフィルタ20は、ローパスフィルタとして機能することによって、サンプル信号から高周波数のノイズ成分を除去できる。ゼロクロスフィルタ20は、サンプル信号の第1カットオフ周波数以上の周波数成分を所定値以上の減衰率で減衰させる一方で、サンプル信号の第1カットオフ周波数未満の周波数成分を所定値未満の減衰率でしか減衰させない。所定値は、各周波数成分のゲインが−3dB(デシベル)となる減衰率に対応してよい。第1カットオフ周波数は、デジタルフィルタの可変パラメータであるとする。制御部70は、第1カットオフ周波数を設定できるとする。
The zero
ゼロクロスフィルタ20は、ローパスフィルタとして機能する際に、入力されるサンプル信号と、フィルタリングしたサンプル信号との間に所定の位相遅延を生じさせる。ゼロクロスフィルタ20で生じる位相遅延は、周波数ごとに異なり得る。本実施形態において、ゼロクロスフィルタ20で生じる位相遅延は、フィルタリングされる信号のピーク周波数における位相遅延によって代表され、φZで表されるとする。
When the zero-
ゼロクロスフィルタ20は、フィルタリングしたサンプル信号をゼロクロス検出器25に出力する。フィルタリングしたサンプル信号は、ゼロクロス検出用信号とも称される。
The zero-
図5に示されるように、ゼロクロス検出用信号は、サンプル信号に含まれるノイズを減衰させたノイズ除去後波形を有している。ノイズ除去後波形は、真の波形に近づけられている。ゼロクロス検出用信号は、サンプル信号に対してφZだけ遅延して出力される。 As shown in FIG. 5, the zero-cross detection signal has a noise-removed waveform obtained by attenuating the noise contained in the sample signal. The waveform after noise removal is closer to the true waveform. The zero-cross detection signal is output with a delay of φZ with respect to the sample signal.
ゼロクロス検出器25は、ゼロクロス検出用信号の符号が反転するタイミングをゼロクロス点として検出する。ゼロクロス検出器25は、ゼロクロス点を検出したタイミングを演算器60とサンプル数カウンタ30とに出力する。
The zero-
サンプル数カウンタ30は、ゼロクロス検出器25が検出した1つのゼロクロス点から次のゼロクロス点までの間に含まれるサンプル信号の数をカウントする。ノイズを含まないサンプル信号において、1つのゼロクロス点から次のゼロクロス点までの長さは、サンプル信号の1つの半波の長さに対応する。サンプル信号の1つの半波の長さに含まれるサンプル信号の数は、サンプリングポイント数とも称され、Nで表される。サンプル数カウンタ30は、サンプリングポイント数(N)を演算器60に出力する。サンプル数カウンタ30は、サンプリングポイント数(N)を制御部70に出力してもよい。言い換えれば、サンプル数カウンタ30は、連続する2つのゼロクロス点の間の期間に含まれるサンプル値の数をカウントする。連続する2つのゼロクロス点の間の期間は、カウント期間とも称される。
The sample number counter 30 counts the number of sample signals included between one zero cross point detected by the zero
ここで、ゼロクロス検出器25は、サンプル信号が正の半波となっている期間において、ゼロクロス検出用信号からゼロクロス点を検出する。サンプル信号が正の半波となっている期間は、P3で表されている。1つのゼロクロス点から次のゼロクロス点までの間が1つの半波とみなされる。1つのゼロクロス点から次のゼロクロス点までの長さは、D3で表される。したがって、P3で表される期間において1つの半波として検出される長さ(D3)は、P3の期間の長さに一致する。
Here, the zero-
一方、仮に、ゼロクロス検出器25がサンプル信号からゼロクロス点を検出する場合、1つのゼロクロス点から次のゼロクロス点までの長さは、D3’で表される。この場合、P3で表される期間において1つの半波として検出される長さ(D3’)は、P3の期間の長さと異なる。
On the other hand, if the zero-
したがって、本実施形態に係る実効値算出装置1は、ゼロクロス検出器25にゼロクロス検出用信号からゼロクロス点を検出させることによって、サンプル信号からゼロクロス点を検出させる場合よりも、1つの半波の長さの検出精度を向上できる。
Therefore, the effective value calculation device 1 according to the present embodiment has a length of one half wave as compared with the case where the zero
サンプル数カウンタ30は、ゼロクロス検出器25で検出されたゼロクロス点に基づいて定まる1つの半波についてサンプリングポイント数をカウントする。サンプル数カウンタ30は、正の半波又は負の半波にかかわらず、半波のサンプリングポイント数をカウントする。
The sample number counter 30 counts the number of sampling points for one half wave determined based on the zero cross point detected by the zero
演算器60の各構成部は、以下説明するように機能する。2乗回路62は、サンプル値を2乗し、総和回路64に出力する。総和回路64は、ゼロクロス検出器25からゼロクロス点を取得する。総和回路64は、ゼロクロス点を半波の開始点又は終了点として、サンプル信号の1つの半波の期間を特定する。総和回路64は、1つの半波の期間として特定した期間内において2乗回路62の出力値を累積加算する。つまり、総和回路64は、サンプル信号の半波毎に2乗回路62の出力値を累積加算した値を得る。総和回路64は、累積加算した値を除算回路66に出力する。除算回路66は、サンプル数カウンタ30からサンプリングポイント数(N)を取得する。除算回路66は、総和回路64の出力値を、サンプリングポイント数(N)で除算し、開平回路68に出力する。開平回路68は、除算回路66の出力値を開平する。つまり、開平回路68は、除算回路66の出力値の平方根を算出する。開平回路68は、算出した平方根を制御部70に出力する。
Each component of the
以上述べてきたように、演算器60は、サンプル信号の1つの半波の期間に含まれるサンプル値の二乗平均平方根をサンプル信号の実効値として算出する。ある1つの半波の期間において、サンプル信号は、1からNまでのサンプリング番号と、各サンプリング番号に対応するu(i)で表されるサンプル値(i=1〜N)とによって特定されるとする。演算器60は、各サンプル値を2乗して累積加算し、サンプル数カウンタ30で算出されたサンプリングポイント数(N)で除算し、開平する。以上の演算の結果は、サンプル信号の実効値に対応し、Uで表される。サンプル信号の実効値(U)は、以下の式(1)によって表される。
演算器60は、サンプル信号の正の半波又は負の半波にかかわらず、半波の実効値を算出する。演算器60は、算出したサンプル信号の実効値を制御部70に出力する。制御部70は、算出した実効値を、例えば表示部74に表示してよい。
The
以上述べてきたように、本実施形態に係る実効値算出装置1は、ゼロクロスフィルタ20によってサンプル信号に含まれるノイズを低減することによって、ゼロクロス検出器25で検出するゼロクロス点を、真の波形におけるゼロクロス点に近づけることができる。つまり、本実施形態に係る実効値算出装置1は、ゼロクロス点の検出精度を向上できる。その結果、本実施形態に係る実効値算出装置1は、サンプル信号の実効値の算出精度を向上できる。
As described above, the effective value calculation device 1 according to the present embodiment sets the zero cross point detected by the zero
入力信号の周波数は、ほぼ一定であるとする。この場合、入力信号の半波の長さは、ほぼ一定である。したがって、サンプル数カウンタ30で検出されるゼロクロス点の間隔はほぼ一定である。つまり、1つの半波の長さの検出値は、ほぼ一定である。サンプル数カウンタ30は、複数の半波の長さに含まれるサンプリングポイント数の平均値(M)を算出し、平均値(M)を演算器60に出力してよい。このようにすることで、1つの半波の始点と終点のゼロクロス点の検出誤差、又は、1つの半波の長さの検出誤差が平均化によって低減され得る。その結果、ゼロクロス点の検出誤差又は1つの半波の長さの検出誤差に起因するサンプリングポイント数のカウント誤差が平均化によって低減され得る。サンプル数カウンタ30は、正の半波に含まれるサンプリングポイント数と負の半波に含まれるサンプリングポイント数とを区別せずに平均値を算出してよい。サンプル数カウンタ30は、サンプリングポイント数の平均値(M)を制御部70に出力してもよい。
It is assumed that the frequency of the input signal is almost constant. In this case, the length of the half wave of the input signal is almost constant. Therefore, the interval between the zero cross points detected by the
演算器60は、各サンプル値を2乗して累積加算した値を、サンプル数カウンタ30で算出されたサンプリングポイント数の平均値(M)で除算し、開平してもよい。この場合、サンプル信号の実効値(U)は、以下の式(2)によって表される。
演算器60が二乗平均の際にサンプリングポイント数の平均値(M)で除算することによって、ゼロクロス点の検出誤差が低減され得る。このようにすることで、サンプリングポイント数の検出精度が向上する。その結果、サンプル信号の実効値の算出精度が向上する。
By dividing by the average value (M) of the number of sampling points when the
式(2)において、2乗して累積したサンプル値の数(N)と、除数となる平均値(M)とが異なることがある。ここで、2乗したサンプル値のうち、ゼロクロス点に近い点におけるサンプル値の2乗は、無視できる程度に小さくなる。したがって、二乗平均平方根(U)の値の精度に対して除数となる平均値(M)の精度が及ぼす影響は、サンプル値の数(N)が及ぼす影響よりも大きくなる。その結果、サンプル値の数(N)と平均値(M)とが異なる値となっても、平均値(M)が高精度で算出されることによって、二乗平均平方根(U)が高精度で算出される。 In the formula (2), the number of sample values (N) obtained by squaring and accumulating may be different from the average value (M) which is a divisor. Here, among the squared sample values, the square of the sample value at a point close to the zero cross point becomes negligibly small. Therefore, the effect of the accuracy of the mean value (M), which is a divisor, on the accuracy of the root mean square (U) value is greater than the effect of the number of sample values (N). As a result, even if the number of sample values (N) and the average value (M) are different, the root mean square (U) can be calculated with high accuracy by calculating the average value (M) with high accuracy. Calculated.
ラインフィルタ40は、デジタル信号としてのサンプル信号をフィルタリングするデジタルフィルタとして構成される。ラインフィルタ40は、第2カットオフ周波数以上の高周波数成分を減衰させるローパスフィルタとして機能する。ラインフィルタ40は、ローパスフィルタとして機能することによって、サンプル信号から高周波数のノイズ成分を除去できる。ラインフィルタ40の第2カットオフ周波数は、ゼロクロスフィルタ20の第1カットオフ周波数よりも高いとする。このようにすることで、ラインフィルタ40でフィルタリングされて演算器60に入力されるサンプル信号の変形が少なくなる。その結果、演算器60で算出されるサンプル値の二乗平均平方根、すなわちサンプル信号の実効値に対するラインフィルタ40の影響が低減され得る。ラインフィルタ40は、ローパスフィルタとして機能する際に、入力されるサンプル信号と、フィルタリングしたサンプル信号との間に所定の位相遅延を生じさせる。ラインフィルタ40で生じる位相遅延は、フィルタリングされる信号のピーク周波数における位相遅延によって代表され、φLで表されるとする。第2カットオフ周波数は、デジタルフィルタの可変パラメータであるとする。制御部70は、第2カットオフ周波数を設定できるとする。
The
ディレイ回路50は、ラインフィルタ40の出力をサンプリングクロック単位で遅延させる。言い換えれば、ディレイ回路50は、ラインフィルタ40から入力されたフィルタリング後のサンプル信号を所定の遅延時間だけ遅延させて演算器60に出力する。サンプル信号を遅延させた信号は、ディレイ信号とも称される。図5に示されるように、ディレイ信号の位相は、ゼロクロス検出用信号の位相に合わせられてよい。
The
遅延時間は、例えば、サンプリングクロックの周期の自然数倍の時間に設定されてよい。遅延時間は、ゼロクロスフィルタ20の位相遅延と、ラインフィルタ40の位相遅延とに基づいて設定される。位相遅延が周波数成分によって異なる場合、遅延時間は、入力信号のピーク周波数における位相遅延に基づいて設定される。上述のとおり、入力信号のピーク周波数における、ゼロクロスフィルタ20及びラインフィルタ40の位相遅延は、それぞれφZ及びφLで表される。位相遅延の単位は、「°」又は「度」(degree)で表されるとする。入力信号のピーク周波数は、fpで表されるとする。遅延時間は、Tdで表されるとする。この場合、遅延時間(Td)を算出する式は、Td=(1/fp)×(φZ−φL)/360と表される。位相遅延の単位は、「ラジアン」(rad)等の他の単位で表されてもよい。
The delay time may be set to, for example, a time that is a natural number multiple of the period of the sampling clock. The delay time is set based on the phase delay of the zero
遅延時間は、サンプリングクロックの周期の倍数に対応する遅延量に変換され得る。遅延量は、Dで表されるとする。サンプリングクロックの周期は、Tsで表されるとする。この場合、D=Td/Tsが成立する。制御部70は、ディレイ回路50のパラメータとして遅延時間(Td)を設定してもよいし、遅延量(D)を設定してもよい。
The delay time can be converted into a delay amount corresponding to a multiple of the sampling clock period. The amount of delay is represented by D. It is assumed that the period of the sampling clock is represented by Ts. In this case, D = Td / Ts holds. The
演算器60は、ディレイ信号のサンプル値と、ゼロクロス検出器25から取得したゼロクロス点の検出結果とに基づいて、1つの半波の期間における実効値を算出してよい。
The
制御部70は、操作部72からの操作入力によって、ゼロクロスフィルタ20の第1カットオフ周波数、又は、ラインフィルタ40の第2カットオフ周波数を設定してよい。制御部70は、操作部72からの操作入力によって入力信号のピーク周波数(fp)を取得してよい。操作部72は、入力信号のピーク周波数(fp)として、50Hz又は60Hz等の所定値を選択して入力できるように構成されてよいし、ピーク周波数(fp)の数値を入力できるように構成されてもよい。制御部70は、入力されたピーク周波数(fp)に基づいてピーク周波数における位相遅延を算出し、遅延時間(Td)又は遅延量(D)を設定してもよい。
The
制御部70は、サンプリングポイント数の平均値(M)から、入力信号のピーク周波数(fp)を算出してもよい。ピーク周波数(fp)を算出する式は、fp=1/(M×2×Ts)と表される。制御部70は、算出したピーク周波数(fp)に基づいてピーク周波数における位相遅延を算出し、遅延時間(Td)又は遅延量(D)を設定してもよい。
The
実効値算出装置1の各構成部は、1又は複数の回路として実現され得る。回路は、電気回路又は電子回路を含む。例えば、サンプル数カウンタ30は、ゼロクロスを検出する回路と、サンプル信号の数をカウントする回路とをそれぞれ備えてよい。例えば、演算器60は、2乗回路62と、総和回路64と、除算回路66と、開平回路68とをそれぞれ備えてもよいし、これらの回路をまとめた1つの回路を備えてもよい。実効値算出装置1に含まれる複数の構成部が1つの回路として実現されてもよい。例えば、ディレイ回路50と演算器60とは、1つの回路として実現されてもよい。例えば、演算器60と制御部70とは、1つの回路として実現されてもよい。回路の少なくとも一部は、集積回路で置き換えられてもよい。
Each component of the effective value calculation device 1 can be realized as one or a plurality of circuits. The circuit includes an electric circuit or an electronic circuit. For example, the
実効値算出装置1の各構成部の機能は、1又は複数のプロセッサにプログラムを実行させることによって実現されてもよい。例えば、制御部70として機能するプロセッサが、演算器60等の他の構成部の機能を実現してもよい。実効値算出装置1全体としての機能が1つの電気回路又は電子回路によって実現されてもよい。例えば、1つのプロセッサが実効値算出装置1全体としての機能を実現してもよい。
The function of each component of the effective value calculation device 1 may be realized by having one or a plurality of processors execute the program. For example, the processor that functions as the
(ホールド機能)
図6に示されるように、サンプル信号の振幅が急激に変化することがある。ここで、サンプル信号の振幅がQで表されるゼロクロス点の前後で変化すると仮定する。上述したように、ゼロクロス検出用信号の位相は、ゼロクロスフィルタ20によってφZで表される分だけ遅れる。Qで表されるタイミングから、φZで表される分だけ遅れたタイミングは、Q’として表されている。
(Hold function)
As shown in FIG. 6, the amplitude of the sample signal may change abruptly. Here, it is assumed that the amplitude of the sample signal changes before and after the zero cross point represented by Q. As described above, the phase of the zero-cross detection signal is delayed by the amount represented by φZ by the zero-
ここで、ゼロクロスフィルタ20は、サンプル信号の高周波成分を減衰させた信号をゼロクロス検出用信号として出力する。サンプル信号の振幅が変化するタイミングにおいて、ゼロクロス検出用信号におけるゼロクロス点は、サンプル信号の周期に対してずれる。その結果、ゼロクロス検出用信号におけるゼロクロス点と、Q’で表されるタイミングとの間に、φerrで表される誤差が生じる。ゼロクロス点の検出誤差は、サンプリングポイント数の誤差を引き起こす。
Here, the zero-
実効値算出装置1は、二乗和をサンプリングポイント数で除算することによって、サンプル値の二乗平均を算出する。したがって、サンプリングポイント数の誤差は、二乗平均の算出精度を低下させるとともに、サンプル信号の実効値の算出精度を低下させる。 The effective value calculation device 1 calculates the root mean square of the sample values by dividing the sum of squares by the number of sampling points. Therefore, the error in the number of sampling points lowers the calculation accuracy of the root mean square and also lowers the calculation accuracy of the effective value of the sample signal.
一方で、入力信号のピーク周波数(fp)がほぼ一定であることに基づけば、入力信号の振幅の変化にかかわらず、サンプリングポイント数の平均値(M)がそのまま二乗平均を算出するための除数として用いられ得る。つまり、入力信号の振幅が変化した場合にも、入力信号の振幅が変化していない期間において検出したサンプリングポイント数の平均値(M)がそのまま二乗平均を算出するための除数として用いられてよい。このようにすることで、振幅の変化によって生じるゼロクロス点の検出誤差の影響が回避され得る。 On the other hand, based on the fact that the peak frequency (fp) of the input signal is almost constant, the mean value (M) of the number of sampling points is a divisor for calculating the root mean square as it is regardless of the change in the amplitude of the input signal. Can be used as. That is, even when the amplitude of the input signal changes, the average value (M) of the number of sampling points detected during the period when the amplitude of the input signal does not change may be used as it is as a divisor for calculating the root mean square. .. By doing so, the influence of the detection error of the zero cross point caused by the change in amplitude can be avoided.
サンプル数カウンタ30は、サンプリングポイント数の平均値(M)を更新する更新モード、及び、平均値(M)を更新せずに既に算出済みの値を保持するホールドモードのうち一方のモードで動作し得る。制御部70は、サンプル数カウンタ30がどちらのモードで動作するか設定してよい。
The
図7に示されるように、本実施形態に係る実効値算出装置1は、電気機器又は電子機器等の規格試験に用いられる測定装置80に含まれてよい。測定装置80は、電気機器又は電子機器の規格試験において、試験対象となる供試装置82(EUT:Equipment Under Test)に入力される電圧を測定する。実効値算出装置1は、供試装置82に入力される電圧の実効値を算出するために用いられてよい。
As shown in FIG. 7, the effective value calculation device 1 according to the present embodiment may be included in the measuring
供試装置82は、単相2線式の配線で電源84に接続されているが、三相4線式の配線で電源84に接続されてもよい。三相4線式の配線の場合、測定装置80には、実効値算出装置1を3個使用する。単相2線式の配線は、L相の配線とN相の配線とを含むとする。N相の配線が接地されてもよい。供試装置82は、インピーダンスネットワーク86を介して電源84に接続されてもよい。インピーダンスネットワーク86があることで、供試装置82の消費電流の変動が、電圧変動に変換される。インピーダンスネットワーク86は、例えば、L相に接続される抵抗(RA)及びリアクタンス(jXA)、並びに、N相に接続される抵抗(RN)及びリアクタンス(jXN)を有する。抵抗及びリアクタンスの値は、電源84の仕様に基づいて適宜定められ得る。
The
電源84の出力電圧は、Usで表されている。インピーダンスネットワーク86を介して供試装置82に印加される電圧は、Ueで表されている。供試装置82に流れる電流は、Ieで表されている。Ueは、Ieの大きさに応じて変動する。実効値算出装置1は、Ueで表されている電圧信号の実効値を算出できる。
The output voltage of the
測定装置80は、実効値算出装置1で算出された電圧信号の実効値に基づいて、供試装置82に起因する電圧変動のパラメータを算出してよい。電圧変動のパラメータは、例えば、相対定常電圧変化、又は、最大相対電圧変化を含んでよい。電圧信号は、定常状態と変動状態とを含む。定常状態は、電圧信号の1つの半波の期間の実効値が1秒以上にわたって安定している状態に対応する。変動状態は、定常状態となっていない状態に対応する。
The measuring
相対定常電圧変化は、1回の電圧変動に挟まれた前後の2つの定常状態の電圧の差を、定格電圧で割った値を%で表したものである。相対定常電圧変化は、dcとも表される。例えば、定格電圧が230Vの電源84において、変動前の定常状態の電圧が229Vであり、変動後の定常状態の電圧が228Vであった場合、相対定常電圧変化は、|(228−229)/230|×100=0.43%となる。測定期間中に電圧変動が一度も発生しない場合、dcの値はゼロになるとする。測定期間中に電圧信号が一度も定常状態にならない場合、測定装置80は、dcの算出結果としてゼロを出力してもよいし、未定義の結果である旨を出力してもよい。
The relative steady-state voltage change is the value obtained by dividing the difference between the two steady-state voltages before and after being sandwiched by one voltage fluctuation by the rated voltage and expressing it in%. The relative steady-state voltage change is also expressed as dc. For example, in a
最大相対電圧変化は、1回の電圧変動(2つの定常状態の間の状態)での最大値と最小値の差を定格電圧で割った値を%で表したものである。あるいは、最大相対電圧変化は、1回の電圧変動(2つの定常状態の間の状態)での最大値と最小値の直前の定常状態との差を絶対値で比較して、大きい方の値を定格電圧で割った値を%で表したものである。 The maximum relative voltage change is the value obtained by dividing the difference between the maximum value and the minimum value in one voltage fluctuation (state between two steady states) by the rated voltage and expressing it in%. Alternatively, the maximum relative voltage change is the larger value by comparing the difference between the maximum value in one voltage fluctuation (the state between two steady states) and the steady state immediately before the minimum value by an absolute value. Is divided by the rated voltage and expressed in%.
測定装置80は、電圧変動のパラメータを、実効値算出装置1に算出させてよい。測定装置80は、電圧変動のパラメータを算出するパラメータ算出装置を備えてもよい。パラメータ算出装置は、実効値算出装置1で算出された電圧信号の実効値に基づいて、電圧変動のパラメータを算出してよい。
The measuring
実効値算出装置1に入力される電圧信号(Ue)の周波数は、供試装置82に電流が流れていない場合と電流が流れている場合とで同じである。一方で、電圧信号(Ue)の振幅は、供試装置82に電流が流れていない場合と電流が流れている場合とで異なる。つまり、入力信号の振幅は、供試装置82に電流が流れ始めるタイミングで変動する。入力信号の振幅の変動は、図6に例示される、ゼロクロス点の検出誤差を生じさせる。
The frequency of the voltage signal (Ue) input to the effective value calculation device 1 is the same when no current is flowing through the
図8に、供試装置82(EUT)の状態と、サンプル信号の波形との関係が示されている。供試装置82に電流が流れる状態は、オン状態に対応する。供試装置82に電流が流れていない状態は、オフ状態に対応する。供試装置82がオン状態である場合とオフ状態である場合とで、サンプル信号の振幅が変化する。
FIG. 8 shows the relationship between the state of the test device 82 (EUT) and the waveform of the sample signal. The state in which the current flows through the
制御部70は、供試装置82がオフ状態である期間にサンプル数カウンタ30の動作モードを更新モードに設定することによって、サンプリングポイント数の平均値(M)を算出する。供試装置82がオフ状態である期間において、サンプル信号の振幅は、ほぼ一定である。
The
入力信号の振幅は、供試装置82がオフ状態からオン状態に変化するタイミングで変化する。入力信号の振幅の変化は、ゼロクロス検出用信号におけるゼロクロス点の誤差を生じさせる。したがって、供試装置82がオフ状態からオン状態に変化するタイミングでサンプル数カウンタ30が検出するゼロクロス点の誤差が大きくなり得る。
The amplitude of the input signal changes at the timing when the
制御部70は、供試装置82がオフ状態からオン状態に変化する前に、サンプル数カウンタ30の動作モードをホールドモードに変更する。制御部70は、供試装置82から、供試装置82の状態がオフ状態からオン状態に変化することを予告する信号を取得し、その信号に基づいて、サンプル数カウンタ30の動作モードをホールドモードに変更してよい。制御部70は、供試装置82に対して、サンプル数カウンタ30の動作モードをホールドモードに変更した後、供試装置82をオフ状態からオン状態に遷移させる信号を出力してもよい。
The
ホールドモードで動作するサンプル数カウンタ30は、供試装置82の状態の変化の有無にかかわらず、供試装置82がオフ状態であったときのサンプリングポイント数の平均値(M)を保持して演算器60に出力する。このようにすることで、供試装置82の状態の変化が演算器60で算出されるサンプル信号の実効値の算出結果に対して影響を及ぼさなくなる。その結果、サンプル信号の実効値の算出精度が向上する。
The
制御部70は、供試装置82がオン状態からオフ状態に変化した後に、サンプル数カウンタ30の動作モードを更新モードに変更してよい。制御部70は、供試装置82から、供試装置82の状態がオン状態からオフ状態に変化したことを表す信号を取得し、その信号に基づいて、サンプル数カウンタ30の動作モードを更新モードに変更してよい。制御部70は、供試装置82に対して、供試装置82をオン状態からオフ状態に遷移させる信号を出力した後、サンプル数カウンタ30の動作モードを更新モードに変更してもよい。制御部70は、操作部72からサンプル数カウンタ30の動作モードを設定する操作入力を受け付け、入力に基づいてサンプル数カウンタ30の動作モードを更新モード又はホールドモードに変更してよい。
The
本実施形態に係る実効値算出装置1は、ホールドモードで動作することによって、入力信号の振幅の変化にかかわらず、サンプリングポイント数の誤差を低減できる。その結果、実効値算出装置1は、二乗平均の算出精度を向上できるとともに、サンプル信号の実効値の算出精度を向上できる。 By operating the effective value calculation device 1 according to the present embodiment in the hold mode, it is possible to reduce the error in the number of sampling points regardless of the change in the amplitude of the input signal. As a result, the effective value calculation device 1 can improve the calculation accuracy of the root mean square and the calculation accuracy of the effective value of the sample signal.
本開示に係る実施形態について、諸図面及び実施例に基づき説明してきたが、当業者であれば本開示に基づき種々の変形又は修正を行うことが容易であることに注意されたい。従って、これらの変形又は修正は本開示の範囲に含まれることに留意されたい。例えば、各構成部に含まれる機能などは論理的に矛盾しないように再配置可能であり、複数の構成部を1つに組み合わせたり、或いは分割したりすることが可能である。 Although the embodiments according to the present disclosure have been described based on the drawings and examples, it should be noted that those skilled in the art can easily make various modifications or modifications based on the present disclosure. It should be noted, therefore, that these modifications or modifications are within the scope of this disclosure. For example, the functions included in each component can be rearranged so as not to be logically inconsistent, and a plurality of components can be combined or divided into one.
1 実効値算出装置
10 AD変換器
20 ゼロクロスフィルタ
25 ゼロクロス検出器
30 サンプル数カウンタ
40 ラインフィルタ
50 ディレイ回路
60 演算器(62:2乗回路、64:総和回路、66:除算回路、68:開平回路)
70 制御部(72:操作部、74:表示部)
80 測定装置
82 供試装置(EUT)
84 電源
86 インピーダンスネットワーク
1 Efficient
70 Control unit (72: operation unit, 74: display unit)
80
84
Claims (3)
前記サンプル信号の高周波数成分を減衰させてゼロクロス検出用信号として出力するゼロクロスフィルタと、
前記ゼロクロス検出用信号に基づいて検出された連続する2つのゼロクロス点の間のカウント期間に含まれる前記サンプル値の数をカウントし、サンプリングポイント数として出力するサンプル数カウンタと、
前記サンプル信号に含まれる前記サンプル値のうち、前記サンプル数カウンタから取得した前記サンプリングポイント数だけ連続する前記サンプル値の二乗平均平方根を算出し、前記入力信号の実効値として出力する演算器と
を備える実効値算出装置。 An AD converter that outputs a sample signal including a sample value that samples the input signal at a predetermined sampling interval, and
A zero-cross filter that attenuates the high frequency component of the sample signal and outputs it as a zero-cross detection signal.
A sample number counter that counts the number of sample values included in the count period between two consecutive zero cross points detected based on the zero cross detection signal and outputs it as the number of sampling points.
Among the sample values included in the sample signal, an arithmetic unit that calculates the root mean square of the sample values that are continuous by the number of sampling points acquired from the sample number counter and outputs the effective value of the input signal. Equipped with an effective value calculation device.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020007774A JP7441052B2 (en) | 2020-01-21 | 2020-01-21 | Effective value calculation device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020007774A JP7441052B2 (en) | 2020-01-21 | 2020-01-21 | Effective value calculation device |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021113792A true JP2021113792A (en) | 2021-08-05 |
JP7441052B2 JP7441052B2 (en) | 2024-02-29 |
Family
ID=77077570
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020007774A Active JP7441052B2 (en) | 2020-01-21 | 2020-01-21 | Effective value calculation device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7441052B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113866488A (en) * | 2021-10-21 | 2021-12-31 | 重庆华虹仪表有限公司 | Method and device for calculating true effective value |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9766278B2 (en) | 2012-01-12 | 2017-09-19 | Cirrus Logic, Inc. | Three phase power quality measurement using asynchronous, isolated single phase circuits |
-
2020
- 2020-01-21 JP JP2020007774A patent/JP7441052B2/en active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113866488A (en) * | 2021-10-21 | 2021-12-31 | 重庆华虹仪表有限公司 | Method and device for calculating true effective value |
Also Published As
Publication number | Publication date |
---|---|
JP7441052B2 (en) | 2024-02-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Karimi-Ghartemani et al. | Measurement of harmonics/inter-harmonics of time-varying frequencies | |
CA2695373C (en) | Circuit device and method of measuring clock jitter | |
JP2006250934A (en) | Method and device for measuring duty cycle | |
KR101487788B1 (en) | Clock detection method and device | |
JP2007519005A (en) | Method and apparatus for measuring jitter | |
JP7441052B2 (en) | Effective value calculation device | |
JPS62110168A (en) | Digital trouble point locating apparatus | |
CN114252837A (en) | High-precision equal sampling system and method for harmonic measurement of digital power meter | |
US20020008548A1 (en) | Device for detecting abnormality of clock signal | |
JP4955196B2 (en) | AC signal measuring device | |
WO2020155068A1 (en) | Current measuring device, method, and apparatus | |
EP3388846A1 (en) | Frequency multiplying device | |
KR20110025680A (en) | Capacitance measurement circuit | |
CN214375934U (en) | Wide-range signal generator device | |
JP2005214932A (en) | Signal processor, and voltage measuring instrument and current measuring instrument using signal processor | |
US11188114B2 (en) | Digital circuit to detect presence and quality of an external timing device | |
JP4572536B2 (en) | Sampling type measuring device | |
JP2007198763A (en) | Harmonic analyzer | |
US9797936B2 (en) | Counter enhancements for improved performance and ease-of-use | |
Pardhu et al. | Design and simulation of digital frequency meter using VHDL | |
JPH0980090A (en) | Frequency detector | |
KR20220165658A (en) | Abnormality detection method and abnormality detection apparatus | |
JP2002368700A (en) | Noise eliminator | |
CN114280467A (en) | Time-frequency parameter measuring device based on FPGA | |
CN111092600A (en) | FPGA frequency doubling method based on phase superposition method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20221213 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230719 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230801 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230830 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20231024 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240115 |
|
A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20240125 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20240213 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240216 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7441052 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |