JP2003337062A - 重量測定用ノイズ除去装置および重量測定用ノイズ除去方法 - Google Patents
重量測定用ノイズ除去装置および重量測定用ノイズ除去方法Info
- Publication number
- JP2003337062A JP2003337062A JP2002143110A JP2002143110A JP2003337062A JP 2003337062 A JP2003337062 A JP 2003337062A JP 2002143110 A JP2002143110 A JP 2002143110A JP 2002143110 A JP2002143110 A JP 2002143110A JP 2003337062 A JP2003337062 A JP 2003337062A
- Authority
- JP
- Japan
- Prior art keywords
- noise
- filter
- frequency
- notches
- weight measurement
- 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
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0248—Filters characterised by a particular frequency response or filtering method
- H03H17/025—Notch filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G23/00—Auxiliary devices for weighing apparatus
- G01G23/06—Means for damping oscillations, e.g. of weigh beams
- G01G23/10—Means for damping oscillations, e.g. of weigh beams by electric or magnetic means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G3/00—Weighing apparatus characterised by the use of elastically-deformable members, e.g. spring balances
- G01G3/12—Weighing apparatus characterised by the use of elastically-deformable members, e.g. spring balances wherein the weighing element is in the form of a solid body stressed by pressure or tension during weighing
- G01G3/14—Weighing apparatus characterised by the use of elastically-deformable members, e.g. spring balances wherein the weighing element is in the form of a solid body stressed by pressure or tension during weighing measuring variations of electrical resistance
- G01G3/142—Circuits specially adapted therefor
- G01G3/147—Circuits specially adapted therefor involving digital counting
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0294—Variable filters; Programmable filters
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
(57)【要約】
【課題】 定周期振動ノイズの周波数変動に対してフィ
ルタ特性を追従、変更させて、周波数変動する定周期振
動ノイズを高精度に減衰させること。 【解決手段】 ノイズ除去装置１は、計量装置３から出
力された測定信号ｍｓを増幅するアンプ１０、アナログ
・フィルタ１１、ＡＤＣ１２、デジタル・フィルタ１
４、係数制御部１３およびデータ入力部１５を備える。
デジタル・フィルタ１４は、ＡＤＣ１２から出力された
計量信号Ｄ（ｎ）の周波数成分を複数のノッチ（阻止
域）で減衰させるノッチ・フィルタである。このデジタ
ル・フィルタ１４の伝達関数は、複数のノッチに対応す
る根をもつ複数の因子の積で表現される。係数制御部１
３は、前記複数のノッチの中から指定されたノッチに対
応する根の周波数成分を変更した後、再構成されたデジ
タル・フィルタ１４の伝達関数を多項式に展開し、その
多項式に基づいてフィルタ係数の新たな系列を算出す
る。
ルタ特性を追従、変更させて、周波数変動する定周期振
動ノイズを高精度に減衰させること。 【解決手段】 ノイズ除去装置１は、計量装置３から出
力された測定信号ｍｓを増幅するアンプ１０、アナログ
・フィルタ１１、ＡＤＣ１２、デジタル・フィルタ１
４、係数制御部１３およびデータ入力部１５を備える。
デジタル・フィルタ１４は、ＡＤＣ１２から出力された
計量信号Ｄ（ｎ）の周波数成分を複数のノッチ（阻止
域）で減衰させるノッチ・フィルタである。このデジタ
ル・フィルタ１４の伝達関数は、複数のノッチに対応す
る根をもつ複数の因子の積で表現される。係数制御部１
３は、前記複数のノッチの中から指定されたノッチに対
応する根の周波数成分を変更した後、再構成されたデジ
タル・フィルタ１４の伝達関数を多項式に展開し、その
多項式に基づいてフィルタ係数の新たな系列を算出す
る。
Description
【０００１】
【発明の属する技術分野】本発明は、デジタル・フィル
タを用いた重量測定用ノイズ除去装置および重量測定用
ノイズ除去方法に関する。
タを用いた重量測定用ノイズ除去装置および重量測定用
ノイズ除去方法に関する。
【０００２】
【従来の技術】歪みゲージ式のロードセルやフォースバ
ランスなどの荷重検出器から出力される計量信号には、
計量系に起因する多種多様なノイズが含まれている。搬
送系の荷重測定を行う場合には、例えば、ベルトコンベ
ヤーを駆動するモーターや搬送ローラなどの回転系振動
に起因するノイズや、商用電源などに起因する電気的な
ノイズなどが、定周期性の振動ノイズとして計量信号に
重畳しているのが一般的である。また、構造物を含む計
量系は当該構造物の固有振動数をもつ。この固有振動数
を含む帯域の外部振動が計量系に作用したとき、その外
部振動はその固有振動数付近で増幅され、振動ノイズと
なって計量信号中に出現する。以下、定周期性の振動ノ
イズと固有振動数に起因する振動ノイズとを合わせて
「定周期振動ノイズ」と呼ぶ。
ランスなどの荷重検出器から出力される計量信号には、
計量系に起因する多種多様なノイズが含まれている。搬
送系の荷重測定を行う場合には、例えば、ベルトコンベ
ヤーを駆動するモーターや搬送ローラなどの回転系振動
に起因するノイズや、商用電源などに起因する電気的な
ノイズなどが、定周期性の振動ノイズとして計量信号に
重畳しているのが一般的である。また、構造物を含む計
量系は当該構造物の固有振動数をもつ。この固有振動数
を含む帯域の外部振動が計量系に作用したとき、その外
部振動はその固有振動数付近で増幅され、振動ノイズと
なって計量信号中に出現する。以下、定周期性の振動ノ
イズと固有振動数に起因する振動ノイズとを合わせて
「定周期振動ノイズ」と呼ぶ。
【０００３】近年、その定周期振動ノイズの除去のため
に、経年変化や周辺環境などの影響を受け難く、フィル
タ特性を比較的容易に変更できるデジタル・フィルタが
採用されており、デジタル・フィルタの中でも、特に安
定性に優れた有限インパルス応答（ＦＩＲ）型フィルタ
が好まれて使用される。
に、経年変化や周辺環境などの影響を受け難く、フィル
タ特性を比較的容易に変更できるデジタル・フィルタが
採用されており、デジタル・フィルタの中でも、特に安
定性に優れた有限インパルス応答（ＦＩＲ）型フィルタ
が好まれて使用される。
【０００４】また、従来は、計量信号から定周期振動ノ
イズを除去するため、その定周期振動ノイズに対応する
周波数帯域において十分な減衰量を得られるローパスフ
ィルタが採用されていた。この種のローパスフィルタを
用いて、定周期振動ノイズを比較的短時間で減衰させる
従来技術としては、特公昭６３−０５２６８４号公報記
載の「定周期振動波除去方法」が知られている。この定
周期振動波除去方法は、定周期振動ノイズの略整数周期
分に相当するサンプリング数をタップ数とした移動平均
フィルタを採用している。しかしながら、そのタップ数
は、ノイズの周期（周波数）に応じて変動し得るため、
その移動平均フィルタの遅れ時間が変動し得るという問
題がある。また、定周期振動ノイズの整数周期分に相当
するサンプリング数にタップ数が一致しない場合は、そ
の不一致に相当する誤差が出現するという問題もある。
イズを除去するため、その定周期振動ノイズに対応する
周波数帯域において十分な減衰量を得られるローパスフ
ィルタが採用されていた。この種のローパスフィルタを
用いて、定周期振動ノイズを比較的短時間で減衰させる
従来技術としては、特公昭６３−０５２６８４号公報記
載の「定周期振動波除去方法」が知られている。この定
周期振動波除去方法は、定周期振動ノイズの略整数周期
分に相当するサンプリング数をタップ数とした移動平均
フィルタを採用している。しかしながら、そのタップ数
は、ノイズの周期（周波数）に応じて変動し得るため、
その移動平均フィルタの遅れ時間が変動し得るという問
題がある。また、定周期振動ノイズの整数周期分に相当
するサンプリング数にタップ数が一致しない場合は、そ
の不一致に相当する誤差が出現するという問題もある。
【０００５】また、その定周期振動波除去方法では、周
波数が相対的に高い電源ノイズや、それが相対的に低い
機械振動ノイズのように、互いに周波数の異なる複数の
定周期振動ノイズに対処するために、各ノイズに対応し
た移動平均フィルタを多段で用意しなければならない。
１段のフィルタで複数のノイズ振動を減衰させるには、
当該フィルタが、それらノイズ振動に対応した複数の阻
止域（ノッチ）をもつ必要がある。
波数が相対的に高い電源ノイズや、それが相対的に低い
機械振動ノイズのように、互いに周波数の異なる複数の
定周期振動ノイズに対処するために、各ノイズに対応し
た移動平均フィルタを多段で用意しなければならない。
１段のフィルタで複数のノイズ振動を減衰させるには、
当該フィルタが、それらノイズ振動に対応した複数の阻
止域（ノッチ）をもつ必要がある。
【０００６】しかしながら、搬送計量系の場合であっ
て、被計量物の搬送速度が可変であるがために、その搬
送計量系を構成するモーターやローラの振動周波数が搬
送速度に合わせて変化する場合では、定周期振動ノイズ
の周波数変動に対応するために複数段のフィルタを用意
しなければならなかった。また、例えばモーターやロー
ラに起因するノイズのようにその周波数が変動する場合
と、電源周波数に起因するノイズのようにその周波数が
殆ど不変の場合とにそれぞれ対応した複数種のフィルタ
特性を用意する必要があるのが実状であった。
て、被計量物の搬送速度が可変であるがために、その搬
送計量系を構成するモーターやローラの振動周波数が搬
送速度に合わせて変化する場合では、定周期振動ノイズ
の周波数変動に対応するために複数段のフィルタを用意
しなければならなかった。また、例えばモーターやロー
ラに起因するノイズのようにその周波数が変動する場合
と、電源周波数に起因するノイズのようにその周波数が
殆ど不変の場合とにそれぞれ対応した複数種のフィルタ
特性を用意する必要があるのが実状であった。
【０００７】
【発明が解決しようとする課題】以上の問題や実状など
に鑑みて本発明が目的とするところは、定周期振動ノイ
ズの周波数変動に対してフィルタ特性を追従して変更さ
せ得て、周波数変動する定周期振動ノイズを高精度に減
衰させ得るノイズ除去装置およびノイズ除去方法を提供
する点にある。
に鑑みて本発明が目的とするところは、定周期振動ノイ
ズの周波数変動に対してフィルタ特性を追従して変更さ
せ得て、周波数変動する定周期振動ノイズを高精度に減
衰させ得るノイズ除去装置およびノイズ除去方法を提供
する点にある。
【０００８】
【課題を解決するための手段】上記課題を解決するた
め、請求項１に係る発明は、被計量物の重量を検出して
得た信号をデジタル化した計量信号に含まれる定周期振
動ノイズを除去する重量測定用ノイズ除去装置であっ
て、前記計量信号の周波数成分を複数のノッチで減衰さ
せるフィルタリングを、可変のフィルタ係数を用いて実
行する帯域通過フィルタと、前記フィルタ係数の系列を
算出して前記帯域通過フィルタに供給する係数制御部
と、を備え、前記帯域通過フィルタは、前記複数のノッ
チを含む周波数特性を有すると共に、当該帯域通過フィ
ルタの伝達関数は、前記複数のノッチに対応する根をも
つ複数の因子の積を含んで構成されており、前記係数制
御部は、前記複数のノッチの中から指定されたノッチに
対応する根を含む因子の当該根の周波数成分を変更する
周波数変更手段と、前記帯域通過フィルタの伝達関数を
構成する複数の因子の積を多項式に展開する展開手段
と、前記多項式に基づいて前記フィルタ係数の系列を算
出する係数算出手段と、を備えることを特徴としたもの
である。
め、請求項１に係る発明は、被計量物の重量を検出して
得た信号をデジタル化した計量信号に含まれる定周期振
動ノイズを除去する重量測定用ノイズ除去装置であっ
て、前記計量信号の周波数成分を複数のノッチで減衰さ
せるフィルタリングを、可変のフィルタ係数を用いて実
行する帯域通過フィルタと、前記フィルタ係数の系列を
算出して前記帯域通過フィルタに供給する係数制御部
と、を備え、前記帯域通過フィルタは、前記複数のノッ
チを含む周波数特性を有すると共に、当該帯域通過フィ
ルタの伝達関数は、前記複数のノッチに対応する根をも
つ複数の因子の積を含んで構成されており、前記係数制
御部は、前記複数のノッチの中から指定されたノッチに
対応する根を含む因子の当該根の周波数成分を変更する
周波数変更手段と、前記帯域通過フィルタの伝達関数を
構成する複数の因子の積を多項式に展開する展開手段
と、前記多項式に基づいて前記フィルタ係数の系列を算
出する係数算出手段と、を備えることを特徴としたもの
である。
【０００９】請求項２に係る発明は、請求項１記載の重
量測定用ノイズ除去装置であって、前記帯域通過フィル
タは有限インパルス応答（ＦＩＲ）型のデジタル・フィ
ルタである。
量測定用ノイズ除去装置であって、前記帯域通過フィル
タは有限インパルス応答（ＦＩＲ）型のデジタル・フィ
ルタである。
【００１０】請求項３に係る発明は、請求項１または２
記載の重量測定用ノイズ除去装置であって、前記周波数
変更手段で変更される周波数成分の変化量を手動で入力
せしめるデータ入力部を更に備えたものである。
記載の重量測定用ノイズ除去装置であって、前記周波数
変更手段で変更される周波数成分の変化量を手動で入力
せしめるデータ入力部を更に備えたものである。
【００１１】請求項４に係る発明は、請求項１〜３の何
れか１項に記載の重量測定用ノイズ除去装置であって、
前記計量信号の周波数特性を算出するフーリエ変換部
と、前記周波数特性を解析して当該計量信号に含まれる
定周期振動ノイズの周波数帯域の変化量を検出する変動
検出部と、前記変動検出部で検出された当該変化量に基
づいて前記周波数変更手段で変更される当該周波数成分
の変化量を指定する制御部と、を更に備えたものであ
る。
れか１項に記載の重量測定用ノイズ除去装置であって、
前記計量信号の周波数特性を算出するフーリエ変換部
と、前記周波数特性を解析して当該計量信号に含まれる
定周期振動ノイズの周波数帯域の変化量を検出する変動
検出部と、前記変動検出部で検出された当該変化量に基
づいて前記周波数変更手段で変更される当該周波数成分
の変化量を指定する制御部と、を更に備えたものであ
る。
【００１２】次に、請求項５に係る発明は、被計量物の
重量を検出して得た信号をデジタル化した計量信号に含
まれる定周期振動ノイズを除去する重量測定用ノイズ除
去方法であって、（ａ）帯域通過フィルタ用のフィルタ
係数の系列を算出する工程と、（ｂ）前記計量信号の周
波数成分を複数のノッチで減衰させるフィルタリング
を、前記フィルタ係数の系列をもつ帯域通過フィルタに
よって実行する工程と、を備え、前記帯域通過フィルタ
は、前記複数のノッチを含む周波数特性を有すると共
に、当該帯域通過フィルタの伝達関数は、前記複数のノ
ッチに対応する根をもつ複数の因子の積を含んで構成さ
れており、前記工程（ａ）は、（ａ−１）前記複数のノ
ッチの中から指定されたノッチに対応する根を含む因子
の当該根の周波数成分を変更する工程と、（ａ−２）前
記帯域通過フィルタの伝達関数を構成する複数の因子の
積を多項式に展開する工程と、（ａ−３）前記多項式に
基づいて前記フィルタ係数の系列を算出する工程と、を
備える、ことを特徴としたものである。
重量を検出して得た信号をデジタル化した計量信号に含
まれる定周期振動ノイズを除去する重量測定用ノイズ除
去方法であって、（ａ）帯域通過フィルタ用のフィルタ
係数の系列を算出する工程と、（ｂ）前記計量信号の周
波数成分を複数のノッチで減衰させるフィルタリング
を、前記フィルタ係数の系列をもつ帯域通過フィルタに
よって実行する工程と、を備え、前記帯域通過フィルタ
は、前記複数のノッチを含む周波数特性を有すると共
に、当該帯域通過フィルタの伝達関数は、前記複数のノ
ッチに対応する根をもつ複数の因子の積を含んで構成さ
れており、前記工程（ａ）は、（ａ−１）前記複数のノ
ッチの中から指定されたノッチに対応する根を含む因子
の当該根の周波数成分を変更する工程と、（ａ−２）前
記帯域通過フィルタの伝達関数を構成する複数の因子の
積を多項式に展開する工程と、（ａ−３）前記多項式に
基づいて前記フィルタ係数の系列を算出する工程と、を
備える、ことを特徴としたものである。
【００１３】請求項６に係る発明は、請求項５記載の重
量測定用ノイズ除去方法であって、前記帯域通過フィル
タは有限インパルス応答（ＦＩＲ）型のデジタル・フィ
ルタである。
量測定用ノイズ除去方法であって、前記帯域通過フィル
タは有限インパルス応答（ＦＩＲ）型のデジタル・フィ
ルタである。
【００１４】請求項７に係る発明は、請求項５または６
記載の重量測定用ノイズ除去方法であって、（ｃ）前記
工程（ａ−１）で変更される周波数成分の変化量を手動
で入力せしめる工程、を更に備えたものである。
記載の重量測定用ノイズ除去方法であって、（ｃ）前記
工程（ａ−１）で変更される周波数成分の変化量を手動
で入力せしめる工程、を更に備えたものである。
【００１５】そして、請求項８に係る発明は、請求項５
〜７の何れか１項に記載の重量測定用ノイズ除去方法で
あって、（ｅ）前記計量信号の周波数特性を算出するフ
ーリエ変換工程と、（ｆ）前記周波数特性を解析して当
該計量信号に含まれる定周期振動ノイズの周波数帯域の
変化量を検出する工程と、（ｇ）前記工程（ｆ）で検出
された当該変化量に基づいて前記工程（ａ−１）で変更
される当該周波数成分の変化量を指定する工程と、を更
に備えたものである。
〜７の何れか１項に記載の重量測定用ノイズ除去方法で
あって、（ｅ）前記計量信号の周波数特性を算出するフ
ーリエ変換工程と、（ｆ）前記周波数特性を解析して当
該計量信号に含まれる定周期振動ノイズの周波数帯域の
変化量を検出する工程と、（ｇ）前記工程（ｆ）で検出
された当該変化量に基づいて前記工程（ａ−１）で変更
される当該周波数成分の変化量を指定する工程と、を更
に備えたものである。
【００１６】
【発明の実施の形態】図１は、本発明の実施の形態に係
る重量測定用のノイズ除去装置１を示す機能ブロック図
である。
る重量測定用のノイズ除去装置１を示す機能ブロック図
である。
【００１７】このノイズ除去装置１は、計量装置３から
出力された測定信号ｍｓを増幅するアンプ１０、アナロ
グ・フィルタ１１、Ａ／Ｄ変換器（ＡＤＣ）１２、デジ
タル・フィルタ１４、係数制御部１３およびデータ入力
部１５を備えて構成されている。また、計量装置（ロー
ドセル）３は、計量系２から受けた重量を測定して得た
測定信号ｍｓをノイズ除去装置１に出力する。
出力された測定信号ｍｓを増幅するアンプ１０、アナロ
グ・フィルタ１１、Ａ／Ｄ変換器（ＡＤＣ）１２、デジ
タル・フィルタ１４、係数制御部１３およびデータ入力
部１５を備えて構成されている。また、計量装置（ロー
ドセル）３は、計量系２から受けた重量を測定して得た
測定信号ｍｓをノイズ除去装置１に出力する。
【００１８】前記アンプ１０は、計量装置３から入力し
た測定信号ｍｓを増幅して得た増幅信号ＭＳをアナログ
・フィルタ１１に出力する。アナログ・フィルタ１１
は、その増幅信号ＭＳから不要な高域成分を除去したア
ナログ信号Ａｓを出力する。Ａ／Ｄ変換器１２は、アナ
ログ・フィルタ１１から入力するアナログ信号Ａｓに対
して、所定のサンプリング周期でサンプリングし、所定
の量子化ビット数で量子化したデジタル信号（以下、計
量信号と呼ぶ。）Ｄ（ｎ）（ｎ：整数）をデジタル・フ
ィルタ１４と係数制御部１３とにそれぞれ出力する。
た測定信号ｍｓを増幅して得た増幅信号ＭＳをアナログ
・フィルタ１１に出力する。アナログ・フィルタ１１
は、その増幅信号ＭＳから不要な高域成分を除去したア
ナログ信号Ａｓを出力する。Ａ／Ｄ変換器１２は、アナ
ログ・フィルタ１１から入力するアナログ信号Ａｓに対
して、所定のサンプリング周期でサンプリングし、所定
の量子化ビット数で量子化したデジタル信号（以下、計
量信号と呼ぶ。）Ｄ（ｎ）（ｎ：整数）をデジタル・フ
ィルタ１４と係数制御部１３とにそれぞれ出力する。
【００１９】図２は、その計量信号Ｄ（ｎ）をフィルタ
リングする有限インパルス応答（ＦＩＲ）型のデジタル
・フィルタ１４を示す概略構成図である。
リングする有限インパルス応答（ＦＩＲ）型のデジタル
・フィルタ１４を示す概略構成図である。
【００２０】このデジタル・フィルタ１４は、係数制御
部１３から供給されたフィルタ係数の系列｛ａi｝（ｉ
＝０，１，２，…，Ｋ−１）を用いたフィルタリング機
能を有する。デジタル・フィルタ１４は、Ｋタップ
（Ｋ：２以上の整数）の帯域通過フィルタ（ローパスフ
ィルタ）であり、入力する一連の計量信号Ｄ（ｎ）を１
周期遅延させる遅延素子１７1，１７2，１７3，…，１
７K-1と、係数乗算器１８0，１８1，…，１８K-1と、加
算器１９とで構成されている。（Ｋ−１）個の遅延素子
１７1〜１７K-1は、Ｄフリップ・フロップなどで構成さ
れ、（Ｋ−１）段で直列接続されている。これら遅延素
子１７1〜１７K-1は、サンプリング周期Ｔでパルス（図
示せず）が供給される度に、前段から入力する計量信号
Ｄ（ｓ）（ｓ＝ｎ〜ｎ−Ｋ＋１）を取り込み保持して出
力する。また、係数乗算器１８0〜１８K-1は、それぞ
れ、係数制御部１３から供給されたフィルタ係数ａ0〜
ａK-1を用いて係数乗算処理を行う。各係数乗算器１８m
（ｍ＝０〜Ｋ−１）は、入力信号Ｄ（ｎ−ｍ）にフィル
タ係数ａmを乗算して加算器１９に出力する。そして、
加算器１９は、係数乗算器１８0〜１８K-1から伝達した
全ての信号を加算して得た積和信号Ｘ（ｎ）を出力す
る。
部１３から供給されたフィルタ係数の系列｛ａi｝（ｉ
＝０，１，２，…，Ｋ−１）を用いたフィルタリング機
能を有する。デジタル・フィルタ１４は、Ｋタップ
（Ｋ：２以上の整数）の帯域通過フィルタ（ローパスフ
ィルタ）であり、入力する一連の計量信号Ｄ（ｎ）を１
周期遅延させる遅延素子１７1，１７2，１７3，…，１
７K-1と、係数乗算器１８0，１８1，…，１８K-1と、加
算器１９とで構成されている。（Ｋ−１）個の遅延素子
１７1〜１７K-1は、Ｄフリップ・フロップなどで構成さ
れ、（Ｋ−１）段で直列接続されている。これら遅延素
子１７1〜１７K-1は、サンプリング周期Ｔでパルス（図
示せず）が供給される度に、前段から入力する計量信号
Ｄ（ｓ）（ｓ＝ｎ〜ｎ−Ｋ＋１）を取り込み保持して出
力する。また、係数乗算器１８0〜１８K-1は、それぞ
れ、係数制御部１３から供給されたフィルタ係数ａ0〜
ａK-1を用いて係数乗算処理を行う。各係数乗算器１８m
（ｍ＝０〜Ｋ−１）は、入力信号Ｄ（ｎ−ｍ）にフィル
タ係数ａmを乗算して加算器１９に出力する。そして、
加算器１９は、係数乗算器１８0〜１８K-1から伝達した
全ての信号を加算して得た積和信号Ｘ（ｎ）を出力す
る。
【００２１】このデジタル・フィルタ１４のｚ変換形式
の伝達関数Ｈ（ｚ）は次式（１）の通りである。
の伝達関数Ｈ（ｚ）は次式（１）の通りである。
【００２２】
【数１】
【００２３】前記デジタル・フィルタ１４は、複数のノ
ッチ（阻止域）を含む周波数特性を有するノッチ・フィ
ルタである。一般に、デジタル・フィルタの通過域の振
幅特性は略一定になることが望ましく、そのノッチは等
リップル特性をもつことが望ましい。等リップル特性を
もつデジタル・フィルタの設計法としてはＲｅｍｅｚア
ルゴリズムが知られている。図３は、Ｒｅｍｅｚアルゴ
リズムによって設計されたデジタル・フィルタ１４の振
幅特性を示すグラフである。ここで、同グラフの縦軸は
信号の減衰量（ゲイン）を示し、その横軸は規格化周波
数（ωＴ／π）を示している。その振幅特性は、次式
（２）に示す伝達関数Ｈ（ｅjωT）の絶対値｜Ｈ（ｅ
jωT）｜に従って算出される。
ッチ（阻止域）を含む周波数特性を有するノッチ・フィ
ルタである。一般に、デジタル・フィルタの通過域の振
幅特性は略一定になることが望ましく、そのノッチは等
リップル特性をもつことが望ましい。等リップル特性を
もつデジタル・フィルタの設計法としてはＲｅｍｅｚア
ルゴリズムが知られている。図３は、Ｒｅｍｅｚアルゴ
リズムによって設計されたデジタル・フィルタ１４の振
幅特性を示すグラフである。ここで、同グラフの縦軸は
信号の減衰量（ゲイン）を示し、その横軸は規格化周波
数（ωＴ／π）を示している。その振幅特性は、次式
（２）に示す伝達関数Ｈ（ｅjωT）の絶対値｜Ｈ（ｅ
jωT）｜に従って算出される。
【００２４】
【数２】
【００２５】上式（２）中、ｊ2＝−１、ωは角周波
数、Ｔはサンプリング周期を表している。
数、Ｔはサンプリング周期を表している。
【００２６】図３に示す通り、その振幅特性は、減衰量
が大きい部分にノッチ１６，１６，…を形成する。デジ
タル・フィルタ１４の伝達関数Ｈ（ｚ）は、次式（３）
に示すように各ノッチ１６，…の一部または全部に対応
するゼロ点の複素根λ1，λ2，…，λM（Ｍ：根の数）
を有している。ここで、「ゼロ点」とは、伝達関数Ｈ
（ｚ）がゼロ値をとるときのｚ平面の点を意味してお
り、計量時の基準点を意味するものではない。
が大きい部分にノッチ１６，１６，…を形成する。デジ
タル・フィルタ１４の伝達関数Ｈ（ｚ）は、次式（３）
に示すように各ノッチ１６，…の一部または全部に対応
するゼロ点の複素根λ1，λ2，…，λM（Ｍ：根の数）
を有している。ここで、「ゼロ点」とは、伝達関数Ｈ
（ｚ）がゼロ値をとるときのｚ平面の点を意味してお
り、計量時の基準点を意味するものではない。
【００２７】
【数３】
【００２８】また、前記伝達関数Ｈ（ｚ）は、複素根λ
1，λ2，…，λMを含む因子（ｚ-1−λ1），（ｚ-1−λ
2），…，（ｚ-1−λM）の積からなり、因数分解された
形をもつ。図４は、上式（３）で表現される伝達関数Ｈ
（ｚ）の複素根の配置を示すグラフである。同グラフの
横軸は複素数ｚの実数部を示し、その縦軸は複素数ｚの
虚数部を示している。また、全ての複素根λ1，λ2，
…，λMが単位円２４上に配置されている。
1，λ2，…，λMを含む因子（ｚ-1−λ1），（ｚ-1−λ
2），…，（ｚ-1−λM）の積からなり、因数分解された
形をもつ。図４は、上式（３）で表現される伝達関数Ｈ
（ｚ）の複素根の配置を示すグラフである。同グラフの
横軸は複素数ｚの実数部を示し、その縦軸は複素数ｚの
虚数部を示している。また、全ての複素根λ1，λ2，
…，λMが単位円２４上に配置されている。
【００２９】尚、本実施の形態では、前記伝達関数Ｈ
（ｚ）は、複数の因子（ｚ-1−λ1），（ｚ-1−λ2），
…，（ｚ-1−λM）の積だけで構成されるが、本発明で
はこれに限らず、その積を含むものであればよい。例え
ば、伝達関数Ｈ（ｚ）がその分母に極を含む形を有して
いてもよい。
（ｚ）は、複数の因子（ｚ-1−λ1），（ｚ-1−λ2），
…，（ｚ-1−λM）の積だけで構成されるが、本発明で
はこれに限らず、その積を含むものであればよい。例え
ば、伝達関数Ｈ（ｚ）がその分母に極を含む形を有して
いてもよい。
【００３０】前記複素根λk（ｋ＝１，２，…，Ｍ）
は、極座標で表現されると、次式（４）のようになる。
は、極座標で表現されると、次式（４）のようになる。
【００３１】
【数４】
【００３２】上式（４）中、ｒkは複素根λkの振幅（絶
対値）、Re(λk)は複素根λkの実数部、Im(λk)は複素
根λkの虚数部を表している。また、上式（４）中、θk
は、図４に示すグラフの原点回りの角度を意味し、本実
施の形態では、周波数成分ωとサンプリング周期Ｔとの
積ω×Ｔに等しい。デジタル・フィルタ１４は、その積
ω×Ｔを−π〜πの区画で規格化した規格化周波数（ω
×Ｔ／π）の−１〜＋１の区画でゼロ点を有するが、図
４に示すように、ゼロ点は、実軸２４Ｒに対して軸対称
の位置に存在するため、図３に示すグラフでは、規格化
周波数の０〜＋１の区画に限定した振幅特性が示されて
いる。
対値）、Re(λk)は複素根λkの実数部、Im(λk)は複素
根λkの虚数部を表している。また、上式（４）中、θk
は、図４に示すグラフの原点回りの角度を意味し、本実
施の形態では、周波数成分ωとサンプリング周期Ｔとの
積ω×Ｔに等しい。デジタル・フィルタ１４は、その積
ω×Ｔを−π〜πの区画で規格化した規格化周波数（ω
×Ｔ／π）の−１〜＋１の区画でゼロ点を有するが、図
４に示すように、ゼロ点は、実軸２４Ｒに対して軸対称
の位置に存在するため、図３に示すグラフでは、規格化
周波数の０〜＋１の区画に限定した振幅特性が示されて
いる。
【００３３】従来、ノッチ１６に対応する周波数（以
下、ノッチ周波数と呼ぶ。）を変える方法としては、デ
ジタル・フィルタのタップ数を変更する方法、或いは、
サンプリング周期Ｔを変更する方法の何れかが採用され
ていた。しかしながら、これら従来法では、周波数変動
する定周期振動ノイズに対して実際上対応することが困
難であった。また、サンプリング周期Ｔを変更する従来
法を採用して、サンプリング周期Ｔを、定周期振動ノイ
ズの変動に合わせて連続的に変化させることは現実的で
はない。たとえサンプリング周期Ｔの連続的な変化が実
現可能であるとしても、その変化によって、例えば商用
電源に起因するノイズのように、周波数変動しないノイ
ズに対応したノッチ周波数が一緒に変化してしまう。以
下に説明するように、本発明では、複数のノッチ１６，
１６，…の中から任意のノッチ１６を選択し、選択した
当該ノッチ１６のみを変化させることが可能である。
下、ノッチ周波数と呼ぶ。）を変える方法としては、デ
ジタル・フィルタのタップ数を変更する方法、或いは、
サンプリング周期Ｔを変更する方法の何れかが採用され
ていた。しかしながら、これら従来法では、周波数変動
する定周期振動ノイズに対して実際上対応することが困
難であった。また、サンプリング周期Ｔを変更する従来
法を採用して、サンプリング周期Ｔを、定周期振動ノイ
ズの変動に合わせて連続的に変化させることは現実的で
はない。たとえサンプリング周期Ｔの連続的な変化が実
現可能であるとしても、その変化によって、例えば商用
電源に起因するノイズのように、周波数変動しないノイ
ズに対応したノッチ周波数が一緒に変化してしまう。以
下に説明するように、本発明では、複数のノッチ１６，
１６，…の中から任意のノッチ１６を選択し、選択した
当該ノッチ１６のみを変化させることが可能である。
【００３４】図５は、図１に示した係数制御部１３の概
略構成を示す機能ブロック図である。この係数制御部１
３は、ノッチ選択部２０、周波数変更手段２１、展開手
段２２および係数算出手段２３を備えて構成されてい
る。周波数変更手段２１は、ノッチ選択部２０もしくは
データ入力部１５から、ノッチ１６とこのノッチ１６に
対応するノッチ周波数の変化量δωとを指定される。周
波数変更手段２１は、複数のノッチ１６，１６，…の中
から指定されたノッチ１６に対応する複素根λkの周波
数成分ωkを（ωk＋δω）に変更する。また、周波数変
更手段２１は、変更後の周波数成分（ωk＋δω）を用
いて、変更後の複素根λk’を生成して次段の展開手段
２２に出力する。その変更後の複素根λk’は、次式
（５）のようになる。
略構成を示す機能ブロック図である。この係数制御部１
３は、ノッチ選択部２０、周波数変更手段２１、展開手
段２２および係数算出手段２３を備えて構成されてい
る。周波数変更手段２１は、ノッチ選択部２０もしくは
データ入力部１５から、ノッチ１６とこのノッチ１６に
対応するノッチ周波数の変化量δωとを指定される。周
波数変更手段２１は、複数のノッチ１６，１６，…の中
から指定されたノッチ１６に対応する複素根λkの周波
数成分ωkを（ωk＋δω）に変更する。また、周波数変
更手段２１は、変更後の周波数成分（ωk＋δω）を用
いて、変更後の複素根λk’を生成して次段の展開手段
２２に出力する。その変更後の複素根λk’は、次式
（５）のようになる。
【００３５】
【数５】
【００３６】次に、展開手段２２は、その変更後の複素
根λk’を用いて再構成されるデジタル・フィルタ１４
の伝達関数Ｈ（ｚ）を多項式に展開する。再構成された
伝達関数Ｈ（ｚ）は次式（６Ａ）のようになり、この式
（６Ａ）が多項式に展開されると次式（６Ｂ）のように
なる。
根λk’を用いて再構成されるデジタル・フィルタ１４
の伝達関数Ｈ（ｚ）を多項式に展開する。再構成された
伝達関数Ｈ（ｚ）は次式（６Ａ）のようになり、この式
（６Ａ）が多項式に展開されると次式（６Ｂ）のように
なる。
【００３７】
【数６】
【００３８】次に、係数算出手段２３は、前記展開手段
２２から出力された上式（６Ｂ）の多項式の展開情報に
基づいてフィルタ係数の新たな系列｛ａi’｝を得て出
力する。従って、変更前のフィルタ係数の系列｛ａi｝
から、新たな系列｛ａi’｝が生成され、デジタル・フ
ィルタ１４に出力される。そして、デジタル・フィルタ
１４においては、図２に示した係数乗算器１８0〜１８
K-1は、係数制御部１３から伝達した更新後のフィルタ
係数ａ0’〜ａK-1’を保持し、これらフィルタ係数
ａ0’〜ａK-1’と入力信号Ｄ（ｎ）〜Ｄ（ｎ−Ｋ＋１）
とを乗算して加算器１９に出力する。加算器１９は、係
数乗算器１８0〜１８K-1から出力された全信号を加算
し、上記定周期振動ノイズが除去された計量信号Ｘ
（ｎ）を出力する。
２２から出力された上式（６Ｂ）の多項式の展開情報に
基づいてフィルタ係数の新たな系列｛ａi’｝を得て出
力する。従って、変更前のフィルタ係数の系列｛ａi｝
から、新たな系列｛ａi’｝が生成され、デジタル・フ
ィルタ１４に出力される。そして、デジタル・フィルタ
１４においては、図２に示した係数乗算器１８0〜１８
K-1は、係数制御部１３から伝達した更新後のフィルタ
係数ａ0’〜ａK-1’を保持し、これらフィルタ係数
ａ0’〜ａK-1’と入力信号Ｄ（ｎ）〜Ｄ（ｎ−Ｋ＋１）
とを乗算して加算器１９に出力する。加算器１９は、係
数乗算器１８0〜１８K-1から出力された全信号を加算
し、上記定周期振動ノイズが除去された計量信号Ｘ
（ｎ）を出力する。
【００３９】このように上記ノイズ除去装置１によれ
ば、周波数変動する定周期振動ノイズに対して、デジタ
ル・フィルタ１４の複数のノッチ１６，１６，…の中か
ら指定された任意のノッチ１６に対応する複素根の値を
変化させて、当該ノッチ周波数をその定周期振動ノイズ
の変動に合わせて変更できる。従来法のように、サンプ
リング周期Ｔを変更したりタップ数を変更したりする必
要が無いことから、デジタル・フィルタ１４の応答が速
く、複数種類の定周期振動ノイズを効率良く且つ迅速に
除去し得るノイズ除去装置１を実現できる。更に、一組
のフィルタ係数の系列｛ａi｝から新たな系列｛ａi’｝
を算出できるため、フィルタ係数を格納するメモリ容量
が小さくて済むという利点もある。
ば、周波数変動する定周期振動ノイズに対して、デジタ
ル・フィルタ１４の複数のノッチ１６，１６，…の中か
ら指定された任意のノッチ１６に対応する複素根の値を
変化させて、当該ノッチ周波数をその定周期振動ノイズ
の変動に合わせて変更できる。従来法のように、サンプ
リング周期Ｔを変更したりタップ数を変更したりする必
要が無いことから、デジタル・フィルタ１４の応答が速
く、複数種類の定周期振動ノイズを効率良く且つ迅速に
除去し得るノイズ除去装置１を実現できる。更に、一組
のフィルタ係数の系列｛ａi｝から新たな系列｛ａi’｝
を算出できるため、フィルタ係数を格納するメモリ容量
が小さくて済むという利点もある。
【００４０】また、上記ノイズ除去装置１は、周波数変
更手段２１で変更される周波数成分ωの変化量δωを手
動で入力せしめるデータ入力部（ユーザー・インターフ
ェース）１５を備えている。このデータ入力部１５は、
ユーザーのデータ入力によって指定された変更対象のノ
ッチ１６とノッチ周波数の変化量δωとを周波数変更手
段２１に供給する。これにより、ユーザーは、データ入
力部１５を操作することで、定周期振動ノイズの変動に
合わせてノッチ周波数を手動で変化させることが可能で
ある。また、ユーザーは、ノッチ１６を指定し変化量δ
ωを指定するだけで、フィルタ係数ａ0〜ａK-1を簡便且
つ正確に微調整できる。
更手段２１で変更される周波数成分ωの変化量δωを手
動で入力せしめるデータ入力部（ユーザー・インターフ
ェース）１５を備えている。このデータ入力部１５は、
ユーザーのデータ入力によって指定された変更対象のノ
ッチ１６とノッチ周波数の変化量δωとを周波数変更手
段２１に供給する。これにより、ユーザーは、データ入
力部１５を操作することで、定周期振動ノイズの変動に
合わせてノッチ周波数を手動で変化させることが可能で
ある。また、ユーザーは、ノッチ１６を指定し変化量δ
ωを指定するだけで、フィルタ係数ａ0〜ａK-1を簡便且
つ正確に微調整できる。
【００４１】また、周波数変動するノッチ１６とそのノ
ッチ周波数の変化量δωとが予測可能な場合には、計量
系の状況に応じて変更対象となるノッチ１６や変化量δ
ωなどを予約登録しておくのが望ましい。前記データ入
力部１５は、その予約登録内容を保持し一覧で表示する
ことできる。ユーザーは、一覧で表示された予約登録内
容を確認して手動で選択することが可能である。
ッチ周波数の変化量δωとが予測可能な場合には、計量
系の状況に応じて変更対象となるノッチ１６や変化量δ
ωなどを予約登録しておくのが望ましい。前記データ入
力部１５は、その予約登録内容を保持し一覧で表示する
ことできる。ユーザーは、一覧で表示された予約登録内
容を確認して手動で選択することが可能である。
【００４２】また、上記ノイズ除去装置１は、計量信号
Ｄ（ｎ）に含まれる定周期振動ノイズの周波数帯域の変
動に合わせてノッチ１６，１６，…を自動的に変化させ
追従させる機能をもつ。この機能は、図５に示すノッチ
選択部２０によって実現される。図６は、そのノッチ選
択部２０の概略構成を示す機能ブロック図である。ノッ
チ選択部２０は、ＦＦＴ（高速フーリエ変換）を実行す
るＦＦＴ演算器３０、変動検出部３１および制御部３４
を備えて構成される。ＦＦＴ演算器３０は、Ａ／Ｄ変換
器１２から入力する計量信号Ｄ（ｎ）に対してＦＦＴを
実行し、その計量信号Ｄ（ｎ）の周波数特性ｄ（ｎ）を
算出して変動検出部３１に出力する。
Ｄ（ｎ）に含まれる定周期振動ノイズの周波数帯域の変
動に合わせてノッチ１６，１６，…を自動的に変化させ
追従させる機能をもつ。この機能は、図５に示すノッチ
選択部２０によって実現される。図６は、そのノッチ選
択部２０の概略構成を示す機能ブロック図である。ノッ
チ選択部２０は、ＦＦＴ（高速フーリエ変換）を実行す
るＦＦＴ演算器３０、変動検出部３１および制御部３４
を備えて構成される。ＦＦＴ演算器３０は、Ａ／Ｄ変換
器１２から入力する計量信号Ｄ（ｎ）に対してＦＦＴを
実行し、その計量信号Ｄ（ｎ）の周波数特性ｄ（ｎ）を
算出して変動検出部３１に出力する。
【００４３】前記変動検出部３１は、その計量信号Ｄ
（ｎ）に含まれる定周期振動ノイズを検出するノイズ検
出手段３２と、単数または複数の定周期振動ノイズの周
波数帯域の変化量を検出する変化量検出手段３３とで構
成されている。
（ｎ）に含まれる定周期振動ノイズを検出するノイズ検
出手段３２と、単数または複数の定周期振動ノイズの周
波数帯域の変化量を検出する変化量検出手段３３とで構
成されている。
【００４４】前記ノイズ検出手段３２は、ＦＦＴ演算器
３０から入力する周波数特性ｄ（ｎ）を解析して計量信
号Ｄ（ｎ）に含まれる定周期振動ノイズを検出し、定周
期振動ノイズの中心周波数などの検出情報を変化量検出
手段３３に供給する。図７に、その周波数特性ｄ（ｎ）
のスペクトルの例を示す。実線で示したスペクトルのノ
イズ成分４０，４１，４２，４３，４４は、それぞれ、
互いに異種の定周期振動ノイズの周波数成分を示すもの
である。これら定周期振動ノイズを除去するデジタル・
フィルタ１４の振幅特性を図８に示す。図８に示すノッ
チ（阻止域）１６1，１６2，１６3，１６4，１６5は、
それぞれ、図７に示すノイズ成分４０，４１，４２，４
３，４４を減衰させる特性を示している。このデジタル
・フィルタ１４の伝達関数Ｈ（ｚ）は、各ノッチ１６1
〜１６5に対応したゼロ点（複素根）を有する。
３０から入力する周波数特性ｄ（ｎ）を解析して計量信
号Ｄ（ｎ）に含まれる定周期振動ノイズを検出し、定周
期振動ノイズの中心周波数などの検出情報を変化量検出
手段３３に供給する。図７に、その周波数特性ｄ（ｎ）
のスペクトルの例を示す。実線で示したスペクトルのノ
イズ成分４０，４１，４２，４３，４４は、それぞれ、
互いに異種の定周期振動ノイズの周波数成分を示すもの
である。これら定周期振動ノイズを除去するデジタル・
フィルタ１４の振幅特性を図８に示す。図８に示すノッ
チ（阻止域）１６1，１６2，１６3，１６4，１６5は、
それぞれ、図７に示すノイズ成分４０，４１，４２，４
３，４４を減衰させる特性を示している。このデジタル
・フィルタ１４の伝達関数Ｈ（ｚ）は、各ノッチ１６1
〜１６5に対応したゼロ点（複素根）を有する。
【００４５】前記変化量検出手段３３は、ノイズ検出手
段３２から入力する検出情報に基づいて各定周期振動ノ
イズのスペクトル上の変化量を個別に検出し、検出した
変化量を制御部３４に出力する。その変化量としては、
各ノイズ成分のピークにおける中心周波数などが挙げら
れる。制御部３４は、前記変化量に基づいて周波数変更
手段２１で変更される周波数成分ωkの変化量δωを決
定し、それを次段の周波数変更手段２１に供給する。こ
れにより、ノイズ成分の変動に合わせてノッチ１６，…
を自動的に変化させることができる。
段３２から入力する検出情報に基づいて各定周期振動ノ
イズのスペクトル上の変化量を個別に検出し、検出した
変化量を制御部３４に出力する。その変化量としては、
各ノイズ成分のピークにおける中心周波数などが挙げら
れる。制御部３４は、前記変化量に基づいて周波数変更
手段２１で変更される周波数成分ωkの変化量δωを決
定し、それを次段の周波数変更手段２１に供給する。こ
れにより、ノイズ成分の変動に合わせてノッチ１６，…
を自動的に変化させることができる。
【００４６】図７に示したスペクトルは搬送計量系から
得たものであり、同図中、点線で示したスペクトルは、
その搬送計量系の搬送速度を増した後のスペクトルを示
している。その搬送速度の変更によってノイズ成分４
１，４２の分布は高域へシフトし、ノイズ成分４１’，
４２’へ変化している。その他のノイズ成分４０，４
３，４４に変化は見られない。図８に示した点線は、本
実施の形態に係るノイズ除去装置１によって、ノイズ成
分４１，４２の変動に合わせてノッチ１６2，１６3を高
域へシフトさせたときの振幅特性を示している。この図
８に示すように、ノイズ除去装置１は、多数のノッチ１
６1〜１６5の中からノッチ１６2，１６3を指定し、当該
ノッチ１６2，１６3を指定した変化量だけシフトさせる
ことが可能である。
得たものであり、同図中、点線で示したスペクトルは、
その搬送計量系の搬送速度を増した後のスペクトルを示
している。その搬送速度の変更によってノイズ成分４
１，４２の分布は高域へシフトし、ノイズ成分４１’，
４２’へ変化している。その他のノイズ成分４０，４
３，４４に変化は見られない。図８に示した点線は、本
実施の形態に係るノイズ除去装置１によって、ノイズ成
分４１，４２の変動に合わせてノッチ１６2，１６3を高
域へシフトさせたときの振幅特性を示している。この図
８に示すように、ノイズ除去装置１は、多数のノッチ１
６1〜１６5の中からノッチ１６2，１６3を指定し、当該
ノッチ１６2，１６3を指定した変化量だけシフトさせる
ことが可能である。
【００４７】以上に述べたように、本実施の形態に係る
ノイズ除去装置１では、計量信号Ｄ（ｎ）に含まれる複
数種の定周期振動ノイズの周波数変動を監視しつつ、各
周波数変動に合わせてフィルタ係数ａ0〜ａK-1を変化さ
せ、デジタル・フィルタ１４のノッチ１６，１６，…を
自動的に変化させてその周波数変動に追従させることが
できる。また、定周期振動ノイズの小さく速い変動に対
しても高精度且つ迅速にノッチ１６，１６，…を追従さ
せることが可能である。
ノイズ除去装置１では、計量信号Ｄ（ｎ）に含まれる複
数種の定周期振動ノイズの周波数変動を監視しつつ、各
周波数変動に合わせてフィルタ係数ａ0〜ａK-1を変化さ
せ、デジタル・フィルタ１４のノッチ１６，１６，…を
自動的に変化させてその周波数変動に追従させることが
できる。また、定周期振動ノイズの小さく速い変動に対
しても高精度且つ迅速にノッチ１６，１６，…を追従さ
せることが可能である。
【００４８】
【発明の効果】以上の如く、本発明の請求項１に係るノ
イズ除去装置および請求項５に係るノイズ除去方法によ
れば、周波数変動する定周期振動ノイズに対し、帯域通
過フィルタの複数のノッチのうち指定されたノッチに対
応する根の値のみを変化させて、当該ノッチをその定周
期振動ノイズの変動に追従させることができるため、定
周期振動ノイズを高精度且つ迅速に除去することが可能
である。更に、一組のフィルタ係数の系列から、定周期
振動ノイズの変動に対応した新たな系列を算出できるた
め、フィルタ係数を格納するメモリ容量が小さくて済
む。従って、定周期振動ノイズを効率良く除去し得るノ
イズ除去装置が実現可能である。
イズ除去装置および請求項５に係るノイズ除去方法によ
れば、周波数変動する定周期振動ノイズに対し、帯域通
過フィルタの複数のノッチのうち指定されたノッチに対
応する根の値のみを変化させて、当該ノッチをその定周
期振動ノイズの変動に追従させることができるため、定
周期振動ノイズを高精度且つ迅速に除去することが可能
である。更に、一組のフィルタ係数の系列から、定周期
振動ノイズの変動に対応した新たな系列を算出できるた
め、フィルタ係数を格納するメモリ容量が小さくて済
む。従って、定周期振動ノイズを効率良く除去し得るノ
イズ除去装置が実現可能である。
【００４９】請求項２および請求項６によれば、ＦＩＲ
型のデジタル・フィルタを用いることで動作安定性に優
れた重量測定用ノイズ除去装置を実現できる。
型のデジタル・フィルタを用いることで動作安定性に優
れた重量測定用ノイズ除去装置を実現できる。
【００５０】請求項３および請求項７によれば、この重
量測定用ノイズ除去装置のユーザーは、データ入力部を
介して、定周期振動ノイズの変動に合わせて手動でノッ
チを変化させ追従させることが可能である。また、ユー
ザーは、ノッチを指定し周波数成分の変化量を指定する
だけで、フィルタ係数を簡便且つ正確に微調整すること
ができる。
量測定用ノイズ除去装置のユーザーは、データ入力部を
介して、定周期振動ノイズの変動に合わせて手動でノッ
チを変化させ追従させることが可能である。また、ユー
ザーは、ノッチを指定し周波数成分の変化量を指定する
だけで、フィルタ係数を簡便且つ正確に微調整すること
ができる。
【００５１】請求項４および請求項８によれば、前記定
周期振動ノイズの周波数帯域の変動を監視しつつ、その
変動に合わせてフィルタ係数を変化させ、ノッチをノイ
ズ変動に自動的に追従させられる。また、定周期振動ノ
イズの小さく速い変動に対しても高精度且つ迅速にフィ
ルタのノッチを追従させることが可能となる。
周期振動ノイズの周波数帯域の変動を監視しつつ、その
変動に合わせてフィルタ係数を変化させ、ノッチをノイ
ズ変動に自動的に追従させられる。また、定周期振動ノ
イズの小さく速い変動に対しても高精度且つ迅速にフィ
ルタのノッチを追従させることが可能となる。
【図１】本発明の実施の形態に係る重量測定用のノイズ
除去装置を示す機能ブロック図である。
除去装置を示す機能ブロック図である。
【図２】ＦＩＲ型のデジタル・フィルタを示す概略構成
図である。
図である。
【図３】デジタル・フィルタの振幅特性を示すグラフで
ある。
ある。
【図４】デジタル・フィルタの伝達関数Ｈ（ｚ）の複素
根の配置を示すグラフである。
根の配置を示すグラフである。
【図５】本実施の形態に係るノイズ除去装置を構成する
係数制御部の概略構成を示す機能ブロック図である。
係数制御部の概略構成を示す機能ブロック図である。
【図６】本実施の形態に係るノイズ除去装置を構成する
ノッチ選択部の概略構成を示す機能ブロック図である。
ノッチ選択部の概略構成を示す機能ブロック図である。
【図７】計量信号の周波数スペクトルの例を示すグラフ
である。
である。
【図８】デジタル・フィルタの振幅特性を示すグラフで
ある。
ある。
１ ノイズ除去装置
２ 計量系
３ 計量装置
１０ アンプ
１１ アナログ・フィルタ
１２ Ａ／Ｄ変換器
１３ 係数制御部
１４ デジタル・フィルタ
１５ データ入力部
フロントページの続き
(72)発明者 村田 修二
滋賀県栗東市下鈎959番地１ 株式会社イ
シダ滋賀事業所内
(72)発明者 若狭 由喜夫
滋賀県栗東市下鈎959番地１ 株式会社イ
シダ滋賀事業所内
(72)発明者 栖原 一浩
滋賀県栗東市下鈎959番地１ 株式会社イ
シダ滋賀事業所内
Claims (8)
- 【請求項１】 被計量物の重量を検出して得た信号をデ
ジタル化した計量信号に含まれる定周期振動ノイズを除
去する重量測定用ノイズ除去装置であって、 前記計量信号の周波数成分を複数のノッチで減衰させる
フィルタリングを、可変のフィルタ係数を用いて実行す
る帯域通過フィルタと、 前記フィルタ係数の系列を算出して前記帯域通過フィル
タに供給する係数制御部と、を備え、 前記帯域通過フィルタは、前記複数のノッチを含む周波
数特性を有すると共に、当該帯域通過フィルタの伝達関
数は、前記複数のノッチに対応する根をもつ複数の因子
の積を含んで構成されており、 前記係数制御部は、 前記複数のノッチの中から指定されたノッチに対応する
根を含む因子の当該根の周波数成分を変更する周波数変
更手段と、 前記帯域通過フィルタの伝達関数を構成する複数の因子
の積を多項式に展開する展開手段と、 前記多項式に基づいて前記フィルタ係数の系列を算出す
る係数算出手段と、を備える、ことを特徴とする重量測
定用ノイズ除去装置。 - 【請求項２】 請求項１記載の重量測定用ノイズ除去装
置であって、前記帯域通過フィルタは有限インパルス応
答（ＦＩＲ）型のデジタル・フィルタである、重量測定
用ノイズ除去装置。 - 【請求項３】 請求項１または２記載の重量測定用ノイ
ズ除去装置であって、前記周波数変更手段で変更される
周波数成分の変化量を手動で入力せしめるデータ入力部
を更に備える重量測定用ノイズ除去装置。 - 【請求項４】 請求項１〜３の何れか１項に記載の重量
測定用ノイズ除去装置であって、 前記計量信号の周波数特性を算出するフーリエ変換部
と、 前記周波数特性を解析して当該計量信号に含まれる定周
期振動ノイズの周波数帯域の変化量を検出する変動検出
部と、 前記変動検出部で検出された当該変化量に基づいて前記
周波数変更手段で変更される当該周波数成分の変化量を
指定する制御部と、を更に備える重量測定用ノイズ除去
装置。 - 【請求項５】 被計量物の重量を検出して得た信号をデ
ジタル化した計量信号に含まれる定周期振動ノイズを除
去する重量測定用ノイズ除去方法であって、 （ａ）帯域通過フィルタ用のフィルタ係数の系列を算出
する工程と、 （ｂ）前記計量信号の周波数成分を複数のノッチで減衰
させるフィルタリングを、前記フィルタ係数の系列をも
つ帯域通過フィルタによって実行する工程と、を備え、 前記帯域通過フィルタは、前記複数のノッチを含む周波
数特性を有すると共に、当該帯域通過フィルタの伝達関
数は、前記複数のノッチに対応する根をもつ複数の因子
の積を含んで構成されており、 前記工程（ａ）は、 （ａ−１）前記複数のノッチの中から指定されたノッチ
に対応する根を含む因子の当該根の周波数成分を変更す
る工程と、 （ａ−２）前記帯域通過フィルタの伝達関数を構成する
複数の因子の積を多項式に展開する工程と、 （ａ−３）前記多項式に基づいて前記フィルタ係数の系
列を算出する工程と、を備える、ことを特徴とする重量
測定用ノイズ除去方法。 - 【請求項６】 請求項５記載の重量測定用ノイズ除去方
法であって、前記帯域通過フィルタは有限インパルス応
答（ＦＩＲ）型のデジタル・フィルタである、重量測定
用ノイズ除去方法。 - 【請求項７】 請求項５または６記載の重量測定用ノイ
ズ除去方法であって、 （ｃ）前記工程（ａ−１）で変更される周波数成分の変
化量を手動で入力せしめる工程、を更に備える重量測定
用ノイズ除去方法。 - 【請求項８】 請求項５〜７の何れか１項に記載の重量
測定用ノイズ除去方法であって、 （ｅ）前記計量信号の周波数特性を算出するフーリエ変
換工程と、 （ｆ）前記周波数特性を解析して当該計量信号に含まれ
る定周期振動ノイズの周波数帯域の変化量を検出する工
程と、 （ｇ）前記工程（ｆ）で検出された当該変化量に基づい
て前記工程（ａ−１）で変更される当該周波数成分の変
化量を指定する工程と、を更に備える重量測定用ノイズ
除去方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002143110A JP3771195B2 (ja) | 2002-05-17 | 2002-05-17 | 重量測定用ノイズ除去装置および重量測定用ノイズ除去方法 |
| EP03009836A EP1363112A1 (en) | 2002-05-17 | 2003-05-13 | Noise cancelling apparatus and method for weight measurement |
| US10/436,075 US6907128B2 (en) | 2002-05-17 | 2003-05-13 | Noise canceling apparatus for weight measurement, and noise canceling method for weight measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002143110A JP3771195B2 (ja) | 2002-05-17 | 2002-05-17 | 重量測定用ノイズ除去装置および重量測定用ノイズ除去方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003337062A true JP2003337062A (ja) | 2003-11-28 |
| JP3771195B2 JP3771195B2 (ja) | 2006-04-26 |
Family
ID=29267837
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002143110A Expired - Fee Related JP3771195B2 (ja) | 2002-05-17 | 2002-05-17 | 重量測定用ノイズ除去装置および重量測定用ノイズ除去方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US6907128B2 (ja) |
| EP (1) | EP1363112A1 (ja) |
| JP (1) | JP3771195B2 (ja) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005233897A (ja) * | 2004-02-23 | 2005-09-02 | Ishida Co Ltd | 重量測定装置 |
| JP2010118798A (ja) * | 2008-11-11 | 2010-05-27 | Gifu Univ | ノイズ除去装置、ノイズ除去方法および動作補正装置 |
| JP2010203972A (ja) * | 2009-03-04 | 2010-09-16 | Yamato Scale Co Ltd | デジタル秤用デジタルフィルタ、それを備えたデジタル秤及び濾波処理方法 |
| JP2012199692A (ja) * | 2011-03-18 | 2012-10-18 | Ishida Co Ltd | ノイズ除去装置、重量測定装置、ノイズ除去方法、及びディジタルフィルタの設計方法 |
| JP2013076669A (ja) * | 2011-09-30 | 2013-04-25 | Minebea Co Ltd | フィルタ装置、指示計および計測装置 |
| JP2014096613A (ja) * | 2012-11-07 | 2014-05-22 | Asahi Kasei Electronics Co Ltd | 慣性センサ用のノイズ除去装置及びそのノイズ除去方法 |
| JP2015073254A (ja) * | 2013-10-04 | 2015-04-16 | 日本電信電話株式会社 | 濾波器 |
| CN106104228A (zh) * | 2014-03-20 | 2016-11-09 | 株式会社石田 | 计量装置 |
| KR101694531B1 (ko) * | 2016-05-25 | 2017-01-11 | (주)나우시스템즈 | 진동잡음 추적이 가능한 중량 선별기 |
| JP2017134030A (ja) * | 2016-01-29 | 2017-08-03 | オムロン株式会社 | 信号処理装置、信号処理装置の制御方法、制御プログラム、および記録媒体 |
| JP2017134031A (ja) * | 2016-01-29 | 2017-08-03 | オムロン株式会社 | 信号処理装置、信号処理装置の制御方法、制御プログラム、および記録媒体 |
| US9927283B2 (en) | 2010-09-08 | 2018-03-27 | Yamato Scale Co., Ltd. | Digital filter for digital weigher |
| US9970809B2 (en) | 2010-09-08 | 2018-05-15 | Yamato Scale Co., Ltd. | Digital filter for digital weigher, digital weigher including digital filter for digital weigher and wave filtering process method |
| WO2018110115A1 (ja) * | 2016-12-15 | 2018-06-21 | オムロン株式会社 | 自動フィルタリング方法及び装置 |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007129408A (ja) * | 2005-11-02 | 2007-05-24 | Ishida Co Ltd | ノイズ除去装置、重量測定装置、ノイズ除去方法、及びディジタルフィルタの設計方法 |
| US8156794B2 (en) * | 2006-06-20 | 2012-04-17 | Shimadzu Corporation | Indenting type material testing machine, testing method, and testing program product |
| JP2008182367A (ja) * | 2007-01-23 | 2008-08-07 | Ishida Co Ltd | ノイズ除去装置、重量測定装置、ノイズ除去方法、及びディジタルフィルタの設計方法 |
| DE102007013333A1 (de) | 2007-03-20 | 2008-09-25 | Bizerba Gmbh & Co. Kg | Verfahren zur Parametrierung und zum Betrieb von Waagen |
| DE102007040301B4 (de) | 2007-08-24 | 2013-07-18 | Sartorius Weighing Technology Gmbh | Verfahren und Vorrichtung zur dynamischen Kontrollwägung |
| TW201010311A (en) * | 2008-08-20 | 2010-03-01 | Ic Plus Corp | Echo canceller and echo cancellation method |
| JP2010271091A (ja) | 2009-05-20 | 2010-12-02 | Seiko Epson Corp | 周波数測定装置 |
| JP5440999B2 (ja) | 2009-05-22 | 2014-03-12 | セイコーエプソン株式会社 | 周波数測定装置 |
| JP5517033B2 (ja) | 2009-05-22 | 2014-06-11 | セイコーエプソン株式会社 | 周波数測定装置 |
| JP5582447B2 (ja) | 2009-08-27 | 2014-09-03 | セイコーエプソン株式会社 | 電気回路、同電気回路を備えたセンサーシステム、及び同電気回路を備えたセンサーデバイス |
| JP5815918B2 (ja) | 2009-10-06 | 2015-11-17 | セイコーエプソン株式会社 | 周波数測定方法、周波数測定装置及び周波数測定装置を備えた装置 |
| JP5876975B2 (ja) | 2009-10-08 | 2016-03-02 | セイコーエプソン株式会社 | 周波数測定装置及び周波数測定装置における変速分周信号の生成方法 |
| JP2011259293A (ja) * | 2010-06-10 | 2011-12-22 | Alps Electric Co Ltd | デジタルフィルタ |
| JP5883558B2 (ja) | 2010-08-31 | 2016-03-15 | セイコーエプソン株式会社 | 周波数測定装置及び電子機器 |
| US10197434B2 (en) * | 2014-05-15 | 2019-02-05 | Ncr Corporation | Adaptive filtering of a load cell of a point-of-sale (POS) device |
| US9319057B1 (en) * | 2015-03-19 | 2016-04-19 | Nxp B.V. | Device and method for providing filtering for an analog-to-digital converter (ADC) |
| JP6512123B2 (ja) * | 2016-01-29 | 2019-05-15 | オムロン株式会社 | 信号処理装置、信号処理装置の制御方法、制御プログラム、および記録媒体 |
| JP6471705B2 (ja) * | 2016-01-29 | 2019-02-20 | オムロン株式会社 | 信号処理装置、信号処理装置の制御方法、情報処理プログラム、および記録媒体 |
| US9812114B2 (en) * | 2016-03-02 | 2017-11-07 | Cirrus Logic, Inc. | Systems and methods for controlling adaptive noise control gain |
| WO2020010257A1 (en) * | 2018-07-06 | 2020-01-09 | University Of Massachusetts | Three-dimensional location estimation using multiplicative processing of sensor measurements |
| CN112414526B (zh) * | 2020-11-13 | 2022-07-01 | 迪比（重庆）智能科技研究院有限公司 | 散袋小包装中药的快速称量方法 |
| CN116488613B (zh) * | 2023-06-26 | 2024-05-28 | 深圳深浦电气有限公司 | 一种滤波方法、装置、终端设备及可读存储介质 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06105867B2 (ja) * | 1987-12-18 | 1994-12-21 | 松下電器産業株式会社 | フィルタ係数演算装置 |
| JP2678047B2 (ja) * | 1989-02-18 | 1997-11-17 | 大和製衡株式会社 | 計量機 |
| US5130938A (en) * | 1989-11-29 | 1992-07-14 | Yamato Scale Company, Limited | Device and method for filtering weight indicative signal from weighing device |
| US5717772A (en) * | 1995-08-07 | 1998-02-10 | Motorola, Inc. | Method and apparatus for suppressing acoustic feedback in an audio system |
| US5960091A (en) * | 1997-04-25 | 1999-09-28 | White; Stanley A. | Adaptive removal of resonance-induced noise |
-
2002
- 2002-05-17 JP JP2002143110A patent/JP3771195B2/ja not_active Expired - Fee Related
-
2003
- 2003-05-13 EP EP03009836A patent/EP1363112A1/en not_active Withdrawn
- 2003-05-13 US US10/436,075 patent/US6907128B2/en not_active Expired - Fee Related
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4693353B2 (ja) * | 2004-02-23 | 2011-06-01 | 株式会社イシダ | 重量測定装置 |
| JP2005233897A (ja) * | 2004-02-23 | 2005-09-02 | Ishida Co Ltd | 重量測定装置 |
| JP2010118798A (ja) * | 2008-11-11 | 2010-05-27 | Gifu Univ | ノイズ除去装置、ノイズ除去方法および動作補正装置 |
| JP2010203972A (ja) * | 2009-03-04 | 2010-09-16 | Yamato Scale Co Ltd | デジタル秤用デジタルフィルタ、それを備えたデジタル秤及び濾波処理方法 |
| US9927283B2 (en) | 2010-09-08 | 2018-03-27 | Yamato Scale Co., Ltd. | Digital filter for digital weigher |
| US9970809B2 (en) | 2010-09-08 | 2018-05-15 | Yamato Scale Co., Ltd. | Digital filter for digital weigher, digital weigher including digital filter for digital weigher and wave filtering process method |
| JP2012199692A (ja) * | 2011-03-18 | 2012-10-18 | Ishida Co Ltd | ノイズ除去装置、重量測定装置、ノイズ除去方法、及びディジタルフィルタの設計方法 |
| JP2013076669A (ja) * | 2011-09-30 | 2013-04-25 | Minebea Co Ltd | フィルタ装置、指示計および計測装置 |
| JP2014096613A (ja) * | 2012-11-07 | 2014-05-22 | Asahi Kasei Electronics Co Ltd | 慣性センサ用のノイズ除去装置及びそのノイズ除去方法 |
| JP2015073254A (ja) * | 2013-10-04 | 2015-04-16 | 日本電信電話株式会社 | 濾波器 |
| CN106104228A (zh) * | 2014-03-20 | 2016-11-09 | 株式会社石田 | 计量装置 |
| US10126162B2 (en) | 2014-03-20 | 2018-11-13 | Ishida Co., Ltd. | Weighing device using filtering to improve accuracy |
| JP2017134030A (ja) * | 2016-01-29 | 2017-08-03 | オムロン株式会社 | 信号処理装置、信号処理装置の制御方法、制御プログラム、および記録媒体 |
| JP2017134031A (ja) * | 2016-01-29 | 2017-08-03 | オムロン株式会社 | 信号処理装置、信号処理装置の制御方法、制御プログラム、および記録媒体 |
| KR101694531B1 (ko) * | 2016-05-25 | 2017-01-11 | (주)나우시스템즈 | 진동잡음 추적이 가능한 중량 선별기 |
| WO2018110115A1 (ja) * | 2016-12-15 | 2018-06-21 | オムロン株式会社 | 自動フィルタリング方法及び装置 |
| JPWO2018110115A1 (ja) * | 2016-12-15 | 2019-10-24 | オムロン株式会社 | 自動フィルタリング方法及び装置 |
| US10644676B2 (en) | 2016-12-15 | 2020-05-05 | Omron Corporation | Automatic filtering method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| US20030215100A1 (en) | 2003-11-20 |
| EP1363112A1 (en) | 2003-11-19 |
| US6907128B2 (en) | 2005-06-14 |
| JP3771195B2 (ja) | 2006-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2003337062A (ja) | 重量測定用ノイズ除去装置および重量測定用ノイズ除去方法 | |
| RU2155325C2 (ru) | Способ и устройство для измерения массового расхода материала в массовом расходомере кориолиса | |
| US7702702B2 (en) | Signal processing device for computing phase difference between alternating current signals | |
| US20140067289A1 (en) | Integrated vibration measurement and analysis system | |
| US20040215426A1 (en) | Cable detector with decimating filter and filtering method | |
| US9970809B2 (en) | Digital filter for digital weigher, digital weigher including digital filter for digital weigher and wave filtering process method | |
| JP2008182367A (ja) | ノイズ除去装置、重量測定装置、ノイズ除去方法、及びディジタルフィルタの設計方法 | |
| JP4546926B2 (ja) | 流量計フィルタ・システム及び方法 | |
| CN108168680B (zh) | 一种动态称重滤波方法及其系统 | |
| EP3855620A1 (en) | Method and apparatus for improving mems accelerometer frequency response | |
| JP2007129408A (ja) | ノイズ除去装置、重量測定装置、ノイズ除去方法、及びディジタルフィルタの設計方法 | |
| JP5901126B2 (ja) | ノイズ除去装置、重量測定装置、ノイズ除去方法、及びディジタルフィルタの設計方法 | |
| JP5325266B2 (ja) | 回転体のアンバランス量測定方法および装置 | |
| AU2021203535A1 (en) | A notch filter in a vibratory flow meter | |
| JP4113758B2 (ja) | 重量測定装置、ノイズ除去方法およびディジタルフィルタの設計方法 | |
| JP7009984B2 (ja) | ノイズ除去フィルタ装置及びノイズ除去方法 | |
| US9927283B2 (en) | Digital filter for digital weigher | |
| JP5570787B2 (ja) | 流量計フィルタ・システム及び方法 | |
| US20190145818A1 (en) | Processing apparatus, processing system, physical quantity measurement apparatus, and measurement method | |
| JP4572536B2 (ja) | サンプリング式測定装置 | |
| JP4693353B2 (ja) | 重量測定装置 | |
| JP6985700B2 (ja) | 複数の荷重検出ユニットを備える計量装置とその計量方法 | |
| JP3036985B2 (ja) | スペクトル演算装置 | |
| JP2021135187A (ja) | 計測震度算出装置、計測震度算出方法及び計測震度算出プログラム | |
| JP2021099256A (ja) | 巻線診断装置および巻線診断方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050310 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20051130 |
|
| 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: 20060207 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060208 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 3771195 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090217 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100217 Year of fee payment: 4 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100217 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130217 Year of fee payment: 7 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130217 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140217 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |