JP3004853B2 - Phase demodulator - Google Patents
Phase demodulatorInfo
- Publication number
- JP3004853B2 JP3004853B2 JP5285084A JP28508493A JP3004853B2 JP 3004853 B2 JP3004853 B2 JP 3004853B2 JP 5285084 A JP5285084 A JP 5285084A JP 28508493 A JP28508493 A JP 28508493A JP 3004853 B2 JP3004853 B2 JP 3004853B2
- Authority
- JP
- Japan
- Prior art keywords
- harmonic component
- order harmonic
- output
- odd
- multiplier
- 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.)
- Expired - Fee Related
Links
Landscapes
- Optical Communication System (AREA)
Description
【0001】[0001]
【産業上の利用分野】本発明は、干渉型光ファイバセン
サなどのように、音波などが印加されたセンサファイバ
を通過する伝搬光とリファレンスファイバを通過する参
照光とによる干渉光から音波信号などを復調する位相復
調器に関するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an interference type optical fiber sensor and the like, and a sound wave signal derived from interference light generated by a propagation light passing through a sensor fiber to which a sound wave is applied and a reference light passing through a reference fiber. And a phase demodulator that demodulates the signal.
【0002】[0002]
【従来の技術】従来、このような分野の技術としては、
例えば、次のような文献に記載されるものがあった。 文献;アイイイイ ジャーナル オブ クォンタム エ
レクトロニクス(IEEEJOURNAL OF QUANTUM ELECTRONICS
)、QE−18[10](1982−10)IEEE
(米) A.Dandridge,A.B.Tveten,and T.G.Giallorenzi "Homody
ne Demodulation Scheme for Fiber Optic Sensors Usi
ng Phase Generated Carrier"P.1647−1653 図2は、前記文献に記載された従来の干渉型光ファイバ
センサに用いられる位相復調器の一構成例を示す図であ
る。この干渉型光ファイバセンサ位相復調器は、周波数
がω0 の正弦波信号を出力する信号発生器1を有し、そ
の出力側には、基本周波数ω0 /2πで周波数変調され
たFM変調光を出力する光源2及び位相復調器10が接
続されている。光源2の出力側には、音波信号が印加さ
れるセンシングファイバ3a及びリファレンスファイバ
3bを有するセンサ部3が接続されている。センサ部3
の出力側には、そのセンサ部3の出力の干渉光を電気信
号に変換する光/電気信号変換器4(以下、O/E変換
器と呼ぶ)が接続され、さらに、該O/E変換器4の出
力側には、位相復調器10が接続されている。2. Description of the Related Art Conventionally, techniques in such a field include:
For example, there is one described in the following literature. Literature; Iii Journal of Quantum
Electronics (IEEEJOURNAL OF QUANTUM ELECTRONICS
), QE-18 [10] (1982-10) IEEE
(US) A.Dandridge, ABTveten, and TGGiallorenzi "Homody
ne Demodulation Scheme for Fiber Optic Sensors Usi
ng Phase Generated Carrier "P. 1647-1654 Fig. 2 is a diagram showing a configuration example of a phase demodulator used in the conventional interference type optical fiber sensor described in the above-mentioned document. the demodulator has a signal generator 1 frequency to output a sine wave signal of omega 0, the output side, the light source 2 and the phase and outputs the FM-modulated light frequency modulated at a fundamental frequency omega 0/2 [pi A demodulator 10 is connected to a sensor unit 3 having a sensing fiber 3a to which a sound signal is applied and a reference fiber 3b to the output side of the light source 2.
Is connected to an optical / electrical signal converter 4 (hereinafter, referred to as an O / E converter) for converting the interference light output from the sensor unit 3 into an electric signal. The output side of the device 4 is connected to a phase demodulator 10.
【0003】位相復調器10は、O/E変換器4の干渉
出力の基本周波数ω0 /2πに対する第1次高調波成分
を抽出する第1次高調波成分抽出手段11及び第2次高
調波成分を抽出する第2次高調波成分抽出手段12を有
し、その第1次高調波成分抽出手段11の出力側には、
微分器14及び乗算器17が接続され、その第2次高調
波成分抽出手段12の出力側には、微分器15及び乗算
器16が接続されている。第1次高調波成分抽出手段1
1は、乗算器11a及び低域通過フィルタ11b(以
下、LPFと呼ぶ)を有し、その乗算器11aが、O/
E変換器4及び信号発生器1に接続され、LPF11b
が、乗算器11aに接続されている。第2次高調波成分
抽出手段12は、乗算器12a及びLPF12bを有
し、その乗算器12aが、O/E変換器4及び信号発生
器1に同期し周波数が2ω0 の正弦波信号を出力する同
期信号発生器13に接続され、LPF12bが、乗算器
12aに接続されている。同期信号発生器13は、信号
発生器1に接続されている。第1の微分器14の出力側
には、第1の乗算器16が接続され、さらに、その第1
の乗算器16の出力側に減算器18が接続されている。
第2の微分器15の出力側に第2の乗算器17が接続さ
れ、さらに、その第2乗算器の17の出力側に減算器1
8が接続されている。減算器18の出力側には、積分器
19が接続されている。The phase demodulator 10 includes a first harmonic component extracting means 11 for extracting a first harmonic component of the interference output of the O / E converter 4 with respect to a fundamental frequency ω 0 / 2π, and a second harmonic component. It has a second harmonic component extracting means 12 for extracting the components, and an output side of the first harmonic component extracting means 11 has:
A differentiator 14 and a multiplier 17 are connected, and a differentiator 15 and a multiplier 16 are connected to an output side of the second harmonic component extraction means 12. First harmonic component extraction means 1
1 has a multiplier 11a and a low-pass filter 11b (hereinafter referred to as LPF), and the multiplier 11a
LPF 11b connected to E converter 4 and signal generator 1
Are connected to the multiplier 11a. The second harmonic component extraction unit 12 has a multiplier 12a and LPF 12b, the multiplier 12a is synchronized with the frequency in the O / E converter 4 and the signal generator 1 outputs a sine wave signal of 2 [omega 0 The LPF 12b is connected to the multiplier 12a. The synchronization signal generator 13 is connected to the signal generator 1. A first multiplier 16 is connected to an output side of the first differentiator 14, and a first multiplier 16 is connected to the first multiplier 16.
The subtractor 18 is connected to the output side of the multiplier 16.
A second multiplier 17 is connected to an output side of the second differentiator 15, and a subtracter 1 is connected to an output side of the second multiplier 17.
8 are connected. An integrator 19 is connected to the output side of the subtractor 18.
【0004】以下、この干渉型光ファイバセンサ位相復
調器の動作を説明する。信号発生器1から出力された周
波数ω0 の正弦波により光源2が電流駆動されて、変調
周波数ω0 /2πのFM変調光をセンサ部3に出力す
る。センサ部3に設けられたセンシングファイバ3aに
音波が印加されるとセンシングファイバ3aの屈折率及
びファイバ長が変化し、センシングファイバ3aに入力
された光(以下、伝搬光と呼ぶ)の位相が変化する。一
方、リファレンスファイバ32に入力された光(以下、
参照光と呼ぶ)の位相は変化しないので、伝搬光と参照
光とにより干渉光が発生する。この干渉光がO/E変換
器4により電気信号に変換され、その干渉出力Iは、次
式(1)のようになる。 I=A+Bcosθcos(Ccosω0 t+φ(t)) =A+BJ0 (C)cosθcosφ(t) +2BJ1 (C)cosθcos(φ(t)+π/2)cosω0 t +2BJ2 (C)cosθcos(φ(t)+π)cos2ω0 t ・ +2BJk (C)cosθcos(φ(t)+kπ/2)coskω0 t ・ ・ ・ ・・・(1) 但し、A、B;入力光量に比例する定数 C;FM変調信号の最大周波数偏移及びセンシングファ
イバ3aとリファレンスファイバ3bとの間の光路差の
関数となる位相変調度 θ;伝搬光と参照光の偏光角 φ(t);音波信号 Jk (C);ベッセル関数(k=0,1,2,・・) 乗算器11aで干渉出力Iと、信号発生器1の角周波数
ω0 の正弦波信号とを乗算し、LPF11bを通すこと
により、次式(2)が出力される。 BJ1 (C)cosθsinφ(t) ・・・(2) 乗算器12aで出力電流Iと、同期信号発生器13で発
生した角周波数2ω0の正弦波信号とを乗算し、LPF
12bを通すことにより、次式(3)が出力される。 BJ2 (C)cosθcosφ(t) ・・・(3) 第1の微分器14,第2の微分器15で微分し、第1の
乗算器16,第2の乗算器17でクロス乗算すると、次
式(4),(5)が出力される。 B2 J1 (C)J2 (C)( dφ(t)/dt)cos 2 φ(t)cos 2 θ・・・(4) − B2 J1 (C)J2 (C)( dφ(t)/dt)sin 2 φ(t)cos 2 θ・・・(5) 減算器18で減算することにより、次式(6)が出力される。 B2 J1 (C)J2 (C)(dφ(t)/dt)cos2 θ・・・(6) 積分器19で積分することにより、次式(7)が出力され、音波信号φ(t) が復調される。 B2 J1 (C)J2 (C)φ(t)cos2 θ ・・・(7) [0004] The operation of the interference type optical fiber sensor phase demodulator will be described below. The light source 2 is current-driven by the sine wave of the frequency ω 0 output from the signal generator 1, and outputs the FM modulated light of the modulation frequency ω 0 / 2π to the sensor unit 3. When a sound wave is applied to the sensing fiber 3a provided in the sensor unit 3, the refractive index and the fiber length of the sensing fiber 3a change, and the phase of light (hereinafter, referred to as propagation light) input to the sensing fiber 3a changes. I do. On the other hand, the light input to the reference fiber 32 (hereinafter, referred to as
Since the phase of the reference light does not change, interference light is generated by the propagation light and the reference light. The interference light is converted into an electric signal by the O / E converter 4, and the interference output I is represented by the following equation (1). I = A + Bcosθ cos (Ccosω 0 t + φ (t)) = A + BJ 0 (C) cosθcosφ (t) + 2BJ 1 (C) cosθcos (φ (t) + π / 2) cosω 0 t + 2BJ 2 (C) cosθcos (φ ( t) + π) cos2ω 0 t · + 2BJ k (C) cosθcos (φ (t) + kπ / 2) coskω 0 t · · · ··· (1) where, a, B; constant C is proportional to the input light intensity; FM Phase modulation degree θ that is a function of the maximum frequency shift of the modulation signal and the optical path difference between the sensing fiber 3a and the reference fiber 3b; the polarization angle of the propagating light and the reference light φ (t); the sound wave signal J k (C) A Bessel function (k = 0, 1, 2,...) A multiplier 11a multiplies the interference output I by a sine wave signal having an angular frequency ω 0 of the signal generator 1 and passes through the LPF 11b to obtain the following equation. (2) is output. BJ 1 (C) cos θ sin φ (t) (2) The output current I is multiplied by the multiplier 12 a and the sine wave signal of the angular frequency 2ω 0 generated by the synchronization signal generator 13, and the LPF
By passing through 12b, the following equation (3) is output. BJ 2 (C) cosθcosφ (t) (3) Differentiation by the first differentiator 14 and the second differentiator 15 and cross multiplication by the first multiplier 16 and the second multiplier 17 give The following equations (4) and (5) are output. B 2 J 1 (C) J 2 (C) (dφ (t) / dt) cos 2 φ (t) cos 2 θ ··· (4) - B 2 J 1 (C) J 2 (C) (dφ (T) / dt) sin 2 φ (t) cos 2 θ (5) By subtraction by the subtractor 18, the following equation (6) is output. B 2 J 1 (C) J 2 (C) (dφ (t) / dt) cos 2 θ (6) By integrating with the integrator 19, the following equation (7) is output, and the sound wave signal φ (T) is demodulated. B 2 J 1 (C) J 2 (C) φ (t) cos 2 θ (7)
【0005】[0005]
【発明が解決しようとする課題】しかしながら、従来の
位相変調器10においては、次のような課題があった。
伝搬光と参照光の偏光角θに変動があると、位相変調器
10の出力は、式(7)に示すようにcos2 θで変動
することになる。そのため、光ファイバセンサ3の感度
劣化および変動となり問題であった。特に、光ファイバ
センサをアレイ状に並べた光ファイバセンサアレイのよ
うに、複数チャネルの場合では、偏光角の変動はチャネ
ル毎に異なるためチャネル間の感度のばらつきとなる。
本発明は前記従来技術が持っていた課題として、伝搬光
と参照光の偏光角によって光ファイバセンサの感度が劣
化あるいは変動する、特に複数チャネルの場合には偏光
角の変動によりチャネル間の感度がばらつくという点に
ついて解決をした位相変調器を提供するものである。However, the conventional phase modulator 10 has the following problems.
If the polarization angle θ between the propagating light and the reference light fluctuates, the output of the phase modulator 10 fluctuates by cos 2 θ as shown in Expression (7). Therefore, the sensitivity of the optical fiber sensor 3 is deteriorated and fluctuated, which is a problem. In particular, in the case of a plurality of channels, such as an optical fiber sensor array in which optical fiber sensors are arranged in an array, fluctuations in the polarization angle differ from channel to channel, resulting in variations in sensitivity between channels.
The present invention has a problem that the sensitivity of the optical fiber sensor deteriorates or fluctuates depending on the polarization angles of the propagating light and the reference light. Particularly, in the case of a plurality of channels, the sensitivity between the channels due to the fluctuation of the polarization angle is increased. An object of the present invention is to provide a phase modulator which solves the problem of variation.
【0006】[0006]
【課題を解決するための手段】第1の発明は、前記課題
を解決するために、FM変調された伝搬光と参照光との
光路差により生じた干渉光を電気信号に変換する光/電
気信号変換器と、前記FM変調の変調周波数を基本周波
数として、前記電気信号の奇数次の高調波成分を抽出す
る奇数次高調波成分抽出手段と、前記電気信号の偶数次
の高調波成分を抽出する偶数次高調波成分抽出器と、前
記奇数次高調波成分抽出器により抽出された奇数次高調
波成分を微分する第1の微分器と、前記偶数次高調波成
分と前記第1の微分器の出力とを乗算する第1の乗算器
と、前記偶数次高調波成分抽出手段により抽出された偶
数次高調波成分を微分する第2の微分器と、前記奇数次
高調波成分と前記第2の微分器の出力とを乗算する第2
の乗算器と、前記第1の乗算器と第2の乗算器の出力を
減算する減算器と、前記減算器の出力を積分する積分器
とを備えた位相変調器において、次のような回路を設け
ている。すなわち、前記奇数次高調波成分の補正係数及
び前記偶数次高調波成分の補正係数を出力する補正係数
生成手段と、前記奇数次高調波成分の補正係数と奇数次
高調波成分とを乗算する第1の補正乗算器と、前記第1
の補正乗算器の出力を自乗する第1の自乗器と、前記偶
数次高調波成分の補正係数と偶数次高調波成分とを乗算
する第2の補正乗算器と、前記第2の補正乗算器の出力
を自乗する第2の自乗器と、前記第1と第2の自乗器の
出力を加算する加算器と、前記積分器の出力を前記加算
器の出力との間で除算する除算器とを、設けている。第
2の発明では、第1の発明の位相変調器における補正係
数生成手段は、前記奇数次高調波成分抽出手段による奇
数次高調波成分又は前記偶数次高調波成分抽出手段によ
る偶数次高調波成分と奇数次高調波成分対又は偶数次高
調波成分対を形成するように高調波成分を抽出する高調
波成分抽出手段と、前記奇数次高調波成分対間又は偶数
次高調波成分間で除算する除算器と、前記除算器の出力
から前記奇数次高調波成分の補正係数及び前記偶数次高
調波成分の補正係数を出力する補正係数算出器とを備え
ている。第3の発明では、第1の発明の位相変調器にお
いて、加算器の出力の高周波成分を除去する低域通過フ
ィルタを設けている。According to a first aspect of the present invention, there is provided an optical / electrical device for converting an interference light generated by an optical path difference between an FM-modulated propagation light and a reference light into an electric signal. A signal converter, an odd-order harmonic component extracting unit that extracts an odd-order harmonic component of the electric signal using a modulation frequency of the FM modulation as a fundamental frequency, and an even-order harmonic component of the electric signal. An even-order harmonic component extractor, a first differentiator for differentiating the odd-order harmonic component extracted by the odd-order harmonic component extractor, the even-order harmonic component and the first differentiator A second multiplier for multiplying the output of the second harmonic with the output of the second harmonic, a second differentiator for differentiating the even harmonic component extracted by the even harmonic component extracting means, the odd harmonic component and the second harmonic component. Multiply by the output of the differentiator
And a subtractor for subtracting the outputs of the first and second multipliers, and an integrator for integrating the outputs of the subtractors. Is provided. That is, a correction coefficient generation unit that outputs a correction coefficient for the odd-order harmonic component and a correction coefficient for the even-order harmonic component, and a second multiplier that multiplies the odd-order harmonic component by the correction coefficient for the odd-order harmonic component. 1 correction multiplier, and the first
A first squarer for squaring the output of the correction multiplier, a second correction multiplier for multiplying the correction coefficient of the even-order harmonic component and the even-order harmonic component, and the second correction multiplier. A second squarer for squaring the output of the adder, an adder for adding the outputs of the first and second squarers, and a divider for dividing the output of the integrator between the output of the adder. Is provided. In the second invention, the correction coefficient generation means in the phase modulator according to the first invention is an odd-order harmonic component by the odd-order harmonic component extraction means or an even-order harmonic component by the even-order harmonic component extraction means. And harmonic component extraction means for extracting a harmonic component so as to form an odd harmonic component pair or an even harmonic component pair, and dividing the odd harmonic component pair or even harmonic component. A divider; and a correction coefficient calculator that outputs a correction coefficient of the odd-order harmonic component and a correction coefficient of the even-order harmonic component from an output of the divider. According to a third aspect, in the phase modulator of the first aspect, a low-pass filter for removing a high-frequency component of the output of the adder is provided.
【0007】[0007]
【作用】第1の発明によれば、以上のように位相変調器
を構成したので、補正係数生成手段は、奇数次高調波成
分抽出手段により抽出された奇数次高調波成分の補正係
数(J2m(C1 ))及び偶数次高調波成分抽出手段によ
り抽出された偶数次高調波成分の補正係数(J2k+1(C
1 ))を求める。第1の補正乗算器は、この奇数次高調
波成分の補正係数と奇数次高調波成分とを乗算し、例え
ば、BcosθJ2k+1(C)J2m(C1 )sinφ
(t)を出力する。一方、第2の補正乗算器は、例え
ば、BcosθJ2k+1(C1 )J2m(C)cosφ
(t)を出力する。但し、C1 はCの推定値である。第
1の自乗器は、第1の補正乗算器の出力を自乗し、例え
ば、B2 cos2 θJ2k+1 2 (C)J2m 2 (C1 )si
n2 φ(t)を出力する。第2の自乗器は、第2の補正
乗算器の出力を自乗し、例えば、B2 cos2 θJ2k+1
2 (C1 )J2m 2 (C)cos2 φ(t)を出力する。
C1 はCの値に近似されるので、加算器は、例えば、B
2 cos2 θJ2k+1 2 (C)J2m 2 (C)の値にほぼ等
しい値を出力する。積分器からは、例えば、B2 cos
2 θJ2k+1(C)J2m(C)が出力されるので、除算器
により、積分器の出力と加算器の出力との間で除算する
と、偏光角θを含む項がとり除かれる。第2の発明によ
れば、高調波成分抽出手段は、高調波成分を抽出して、
奇数次高調波成分抽出手段による奇数次高調波成分又は
偶数次高調波成分抽出手段による偶数次高調波成分と、
奇数次高調波成分対又は偶数次高調波成分対を形成す
る。除算器は、奇数次高調波成分対間又は偶数次高調波
成分間で除算して、例えば、J2k+1(C)/J
2n+1(C)又はJ2k(C)/J2n(C)を出力する。補
正係数算出器は、除算器からの出力に基づいて、例え
ば、J2k+1(C1 )及びJ2m(C1 )を求める。第3の
発明では、低域通過フィルタは、C1 がCよりずれた場
合の誤差による高周波成分を除去する。従って、前記課
題を解決できるのである。According to the first aspect of the present invention, since the phase modulator is configured as described above, the correction coefficient generation means can correct the odd-order harmonic component (J) extracted by the odd-order harmonic component extraction means. 2m (C 1 )) and the correction coefficient of the even-order harmonic component extracted by the even-order harmonic component extraction means (J 2k + 1 (C
1 )) The first correction multiplier multiplies the odd-order harmonic component correction coefficient by the odd-order harmonic component, and outputs, for example, Bcos θJ 2k + 1 (C) J 2m (C 1 ) sin φ.
(T) is output. On the other hand, the second correction multiplier is, for example, BcosθJ 2k + 1 (C 1 ) J 2m (C) cos φ
(T) is output. Here, C 1 is an estimated value of C. The first squarer squares the output of the first correction multiplier, for example, B 2 cos 2 θJ 2k + 1 2 (C) J 2m 2 (C 1 ) si
Output n 2 φ (t). The second squarer squares the output of the second correction multiplier, for example, B 2 cos 2 θJ 2k + 1
2 (C 1 ) J 2m 2 (C) cos 2 φ (t) is output.
Since C 1 is approximated by the value of C, the adder, for example,
2 cos 2 θJ 2k + 1 2 (C) Outputs a value substantially equal to the value of J 2m 2 (C). From the integrator, for example, B 2 cos
Since 2 θJ 2k + 1 (C) J 2m (C) is output, when the divider divides the output of the integrator and the output of the adder, the term including the polarization angle θ is removed. According to the second aspect, the harmonic component extracting means extracts the harmonic component,
An odd harmonic component by the odd harmonic component extraction means or an even harmonic component by the even harmonic component extraction means,
An odd harmonic component pair or an even harmonic component pair is formed. The divider divides between odd-order harmonic component pairs or even-order harmonic components to obtain, for example, J 2k + 1 (C) / J
2n + 1 (C) or J 2k (C) / J 2n (C) is output. The correction coefficient calculator obtains, for example, J 2k + 1 (C 1 ) and J 2m (C 1 ) based on the output from the divider. In the third aspect of the invention, the low-pass filter removes high-frequency components due to error when C 1 is shifted from C. Therefore, the above problem can be solved.
【0008】[0008]
【実施例】図1は、本発明の実施例を示す干渉型光ファ
イバセンサに用いられる位相変調器の構成図であり、従
来の図2中の要素と共通の要素には共通の符号が付され
ている。この干渉型光ファイバセンサ位相変調器が、従
来のものと異なる点はAGC回路20(Automatic gain
control)を新たに設けたことである。AGC回路20
は、信号発生器1に同期し周波数が3ω0 の正弦波信号
を出力する同期信号発生器21を有し、その出力側に3
次高調波成分を抽出する第3次高調波成分抽出手段22
が接続されている。第3次高調波成分抽出手段22は、
乗算器22a及びLPF22bを有し、その乗算器22
aが、O/E変換器4及び同期信号発生器21に接続さ
れ、さらにその出力側にLPF22bが接続されてい
る。LPF22bの出力側には、第1次高調波成分及び
第2次周波数成分の補正係数を発生する補正係数発生器
23が接続されている。補正係数発生器23は、除算器
23a及び補正係数算出器23bを有し、その除算器2
3aが、LPF22bに接続され、さらに補正係数算出
器23bが除算器23aに接続されている。補正係数算
出器23bは、図示しないアナログ/ディジタル変換器
(以下、A/D変換器と呼ぶ)、メモリ及びディジタル
/アナログ変換器(以下D/A変換器)などで構成さ
れ、その出力側に、第1の補正乗算器24及び第2の補
正乗算器25が接続されている。第1の補正乗算器24
は、位相復調器10のLPF11bにも接続され、その
出力側に第1の自乗器26が接続されている。第2の補
正乗算器25は、位相復調器10のLPF12bにも接
続され、その出力側に第2の自乗器27が接続されてい
る。第1の自乗器26及び第2の自乗器27の出力側に
は、加算器28が接続され、さらに、その加算器28の
出力側には、LPF29が接続されている。LPF29
の出力側には、除算器30が接続されている。この除算
器30は、位相復調器10の積分器19にも接続されて
いる。FIG. 1 is a block diagram of a phase modulator used in an interference type optical fiber sensor showing an embodiment of the present invention. Elements common to those in FIG. Have been. The point that this interferometric optical fiber sensor phase modulator differs from the conventional one is that the AGC circuit 20 (Automatic gain
control). AGC circuit 20
Has a synchronizing signal generator 21 which is synchronized with the frequency to the signal generator 1 outputs a sine wave signal of 3 [omega] 0, 3 on the output side
Third harmonic component extracting means 22 for extracting the second harmonic component
Is connected. The third harmonic component extraction means 22 includes:
A multiplier 22a and an LPF 22b;
a is connected to the O / E converter 4 and the synchronizing signal generator 21, and the output side thereof is connected to the LPF 22b. The output side of the LPF 22b is connected to a correction coefficient generator 23 that generates correction coefficients for the first harmonic component and the second frequency component. The correction coefficient generator 23 has a divider 23a and a correction coefficient calculator 23b.
3a is connected to the LPF 22b, and the correction coefficient calculator 23b is connected to the divider 23a. The correction coefficient calculator 23b includes an analog / digital converter (hereinafter, referred to as an A / D converter), a memory, a digital / analog converter (hereinafter, a D / A converter), and the like. , A first correction multiplier 24 and a second correction multiplier 25 are connected. First correction multiplier 24
Is also connected to the LPF 11b of the phase demodulator 10, and a first squarer 26 is connected to the output side. The second correction multiplier 25 is also connected to the LPF 12b of the phase demodulator 10, and a second squarer 27 is connected to the output side. An output side of the first squarer 26 and the second squarer 27 is connected to an adder 28, and an output side of the adder 28 is connected to an LPF 29. LPF29
Is connected to a divider 30. The divider 30 is also connected to the integrator 19 of the phase demodulator 10.
【0009】次に、図1に示す干渉型光ファイバセンサ
位相変調器の動作を説明する。信号発生器1から出力さ
れた周波数ω0 の正弦波により光源2が電流駆動され
て、変調周波数ω0 /2πのFM変調光をセンサ部3に
出力する。センサ部3に設けられたセンシングファイバ
3aに音波が印加されると、センシングファイバ3aの
屈折率及びファイバ長が変化し、センシングファイバ3
aに入力された伝搬光の位相が変化する。一方、リファ
レンスファイバ3bに入力された参照光の位相は変化せ
ず、伝搬光と参照光とにより干渉光が発生する。O/E
変換器4は、この干渉光を電気信号に変換し、(1)式
で示される干渉出力Iを出力する。乗算器11aで干渉
出力Iと、信号発生器1の角周波数ω0 の正弦波信号と
を乗算し、LPF11bを通すことにより、(2)式で
示される第1次高調波成分を得る。乗算器12aで干渉
出力Iと、同期信号発生器13で発生した角周波数2ω
0 の正弦波信号とを乗算し、LPF12bを通すことに
より、(3)式で示される第2次高調波成分を得る。第
1の微分器14及び第2の微分器15で微分し、第1の
乗算器16及び第2の乗算器17でクロス乗算すると、
(4)式,(5)式を得る。減算器18で、引き算をし
て、(6)式を得て、積分器19で積分することによ
り、(7)式を得る。乗算器22aで、干渉出力Iと同
期信号発生器21の角周波数3ω0 の正弦波信号とを乗
算し、LPF22bを通すことにより、(8)式によっ
て表される第3次高調波成分が出力される。Next, the operation of the interferometric optical fiber sensor phase modulator shown in FIG. 1 will be described. The light source 2 is current-driven by the sine wave of the frequency ω 0 output from the signal generator 1, and outputs the FM modulated light of the modulation frequency ω 0 / 2π to the sensor unit 3. When a sound wave is applied to the sensing fiber 3a provided in the sensor unit 3, the refractive index and the fiber length of the sensing fiber 3a change, and the sensing fiber 3a
The phase of the propagating light input to a changes. On the other hand, the phase of the reference light input to the reference fiber 3b does not change, and the propagation light and the reference light generate interference light. O / E
The converter 4 converts the interference light into an electric signal and outputs an interference output I represented by the equation (1). An interference output I in multiplier 11a, multiplies the sine wave signal of the angular frequency omega 0 of the signal generator 1, by passing the LPF 11b, to obtain a first-order harmonic components represented by the equation (2). The interference output I in the multiplier 12a and the angular frequency 2ω generated in the synchronization signal generator 13
By multiplying by the sine wave signal of 0 and passing through the LPF 12b, the second harmonic component represented by the equation (3) is obtained. Differentiation by the first differentiator 14 and the second differentiator 15 and cross multiplication by the first multiplier 16 and the second multiplier 17 give:
Equations (4) and (5) are obtained. Subtraction is performed by the subtractor 18 to obtain the expression (6), and integration is performed by the integrator 19 to obtain the expression (7). A multiplier 22a, multiplies the sine wave signal of an angular frequency 3 [omega] 0 of the interference output I and the synchronization signal generator 21, by passing the LPF 22b, output third-order harmonic component represented by the formula (8) Is done.
【0010】 BJ3 (C)sinφ(t)cosθ ・・・(8) 除算器23aで、(2)式で示されるLPE11bの出
力を、(8)式で示されるLPF22bの出力で割り算
をしてJ1 (C)/J3 (C)を得る。補正係数算出器
23bのA/D変換器で、除算器23aの出力J
1 (C)/J3 (C)をディジタル信号に変換し、J1
(C)/J3 (C)のディジタル値によりメモリに格納
される換算テーブルを参照し、J1 (C)とJ2 (C)
の推定値J1 (C1 )とJ2 (C1 )を求めて、この値
をD/A変換器でアナログ信号に変換して、第1の補正
乗算器24に第1次高調波成分の補正係数としてJ
2 (C1 )を、第2の補正乗算器25に第2次高調波成
分の補正係数としてJ1 (C1 )を出力する。第1の補
正乗算器24では、(2)式で示されるLPF11bの
出力と補正係数算出器23bの出力(J2 (C1 ))と
を乗算して、次式(9)が出力される。 BJ1 (C)J2 (C1 )sinφ(t)cosθ ・・・(9) 第1の自乗器26で、自乗して、次式(10)が出力さ
れる。 B2 J1 2 (C)J2 2 (C1 )sin2 φ(t)cos2 θ・・(10) 第2の補正乗算器25では、(3)式で示されるLPF
12bの出力と補正係数算出器23bの出力(J1 (C
1 ))とを乗算して、次式(11)が得られる。 BJ1 (C1 )J2 (C)cosφ(t)cosθ ・・・(11) 第2の自乗器27で、自乗して、次式(12)を得る。 B2 J1 2 (C1 )J2 2 (C)cos2 φ(t)cos2 θ・・(12) 加算器28で、加算すると、次式(13)が出力され
る。 B 2 J 1 2 (C)J2 2 (C)cos2 θ+ψ ・・・(13) ψは、C1 の値がCの値よりずれた場合の誤差であっ
て、次式(14)のように表される。 ψ= ((J 1 2 (C)J2 2 (C1 )−J 1 2 (C)J2 2 (C))B 2 sin 2 φ(t)cos2 θ +((J 1 2 (C1 )J2 2 (C) −J 1 2 (C)J2 2 (C))B 2 cos 2 φ(t)cos2 θ ・・・(14) C1 の値がCの値にほぼ等しいときは、加算器28で
の、加算結果は、次式( 15)のようになる。 B2 J1 2 (C)J2 2 (C)cos2 θ ・・・(15) さらにLPF29では、高周波成分として現れる誤差ψ
を取り除き、より(15)式に近づける。除算器30で
は、(7)式で示される積分器19の出力を、(15)
式で示されるLPF29の出力で割り算をし、φ(t)
/J1 (C)J2 (C)を得る。このように、除算器3
0の出力から伝搬光と参照光との偏光角θによる偏光角
θを含む係数を取り除くことができる。BJ 3 (C) sinφ (t) cos θ (8) The divider 23a divides the output of the LPE 11b expressed by the equation (2) by the output of the LPF 22b expressed by the equation (8). To obtain J 1 (C) / J 3 (C). The A / D converter of the correction coefficient calculator 23b and the output J of the divider 23a
1 (C) / J 3 (C) is converted to a digital signal, and J 1
Referring to the conversion table stored in the memory based on the digital value of (C) / J 3 (C), J 1 (C) and J 2 (C)
Seeking the estimate J 1 (C 1) and J 2 (C 1), into an analog signal the value in the D / A converter, the first harmonic component to the first correction multiplier 24 J as the correction coefficient of
2 (C 1 ) is output to the second correction multiplier 25 as J 1 (C 1 ) as a correction coefficient of the second harmonic component. The first correction multiplier 24 multiplies the output of the LPF 11b represented by the equation (2) by the output (J 2 (C 1 )) of the correction coefficient calculator 23b, and outputs the following equation (9). . BJ 1 (C) J 2 (C 1 ) sin φ (t) cos θ (9) The first squarer 26 squares and outputs the following equation (10). B 2 J 1 2 (C) J 2 2 (C 1 ) sin 2 φ (t) cos 2 θ (10) In the second correction multiplier 25, the LPF expressed by the equation (3) is used.
12b and the output of the correction coefficient calculator 23b (J 1 (C
1 )) to obtain the following equation (11). BJ 1 (C 1 ) J 2 (C) cos φ (t) cos θ (11) The second squarer 27 squares to obtain the following equation (12). B 2 J 1 2 (C 1 ) J 2 2 (C) cos 2 φ (t) cos 2 θ (12) When the addition is performed by the adder 28, the following equation (13) is output. B 2 J 1 2 (C) J 2 2 (C) cos 2 θ + ψ (13) where ψ is the error when the value of C 1 deviates from the value of C. Is represented as ψ = ((J 1 2 (C) J 2 2 (C 1 ) −J 1 2 (C) J 2 2 (C)) B 2 sin 2 φ (t) cos 2 θ + ((J 1 2 (C 1) J 2 2 (C) -J 1 2 (C) J 2 2 (C)) B 2 cos 2 φ (t) cos 2 θ ··· (14) the value of C 1 is approximately equal to the value of C case, the adder 28, the addition result is expressed by the following equation (15). B 2 J 1 2 (C) J 2 2 (C) cos 2 θ ··· (15) further LP F 29 Then, the error と し て
Is removed, and the expression (15) is brought closer. In the divider 30, the output of the integrator 19 expressed by the equation (7) is calculated by the following equation (15).
Divide by the output of the LPF 29 shown in the equation, φ (t)
/ J 1 (C) and J 2 (C). Thus, the divider 3
A coefficient including the polarization angle θ based on the polarization angle θ between the propagation light and the reference light can be removed from the output of 0.
【0011】以上のように、本実施例においては、次の
ような利点がある。AGC回路20は、位相復調回路1
0の出力から伝搬光と参照光との偏光角θによる偏光角
θを含む係数を取り除くことができる。従って、AGC
回路20からの最終出力結果は安定し、光ファイバセン
サとしての感度劣化・変動を抑制することが出来る。特
に、多重化の大きな光ファイバセンサアレイに用いると
チャネル間の感度ばらつきを抑制することが出来る。な
お、本発明は上記実施例に限定されず、種々の変形が可
能である。その変形例としては、例えば次のようなもの
がある。 (1) 第3次高調波成分の代わりに他の奇数次高調波
成分を同期信号発生器21及び乗算器22aによって抽
出し、この奇数次高調波成分と第1次高調波成分とを割
り算して、補正係数算出器23bによって第1次高調波
成分と第2次高調波成分の補正係数を求めることができ
る。また、第2次以外の偶数高調波成分、例えば第4次
高調波成分を抽出し、この偶数高調波成分と第2次高調
波成分とを割り算して、補正係数算出器23bによって
第1次高調波成分と第2次高調波成分の補正係数を求め
ることができる。 (2) 補正係数算出器23bでは、J1 (C)/J3
(C)の値に基づいて、メモリ内の換算テーブルを用い
ずに、CPUによる計算によりJ1 (C1 )及びJ
2 (C1 )を求めることができる。As described above, this embodiment has the following advantages. The AGC circuit 20 includes the phase demodulation circuit 1
A coefficient including the polarization angle θ based on the polarization angle θ between the propagation light and the reference light can be removed from the output of 0. Therefore, AGC
The final output result from the circuit 20 is stable, and the sensitivity deterioration and fluctuation as an optical fiber sensor can be suppressed. In particular, when used in an optical fiber sensor array with large multiplexing, sensitivity variations between channels can be suppressed. Note that the present invention is not limited to the above embodiment, and various modifications are possible. For example, there are the following modifications. (1) Instead of the third harmonic component, another odd harmonic component is extracted by the synchronization signal generator 21 and the multiplier 22a, and the odd harmonic component and the first harmonic component are divided. Thus, the correction coefficient of the first harmonic component and the second harmonic component can be obtained by the correction coefficient calculator 23b. Further, an even harmonic component other than the second harmonic component, for example, a fourth harmonic component is extracted, the even harmonic component is divided by the second harmonic component, and the first harmonic component is calculated by the correction coefficient calculator 23b. The correction coefficients for the harmonic component and the second harmonic component can be obtained. (2) In the correction coefficient calculator 23b, J 1 (C) / J 3
Based on the value of (C), J 1 (C 1 ) and J 1 are calculated by the CPU without using the conversion table in the memory.
2 (C 1 ) can be obtained.
【0012】[0012]
【発明の効果】以上詳細に説明したように、第1の発明
によれば、位相復調器が、奇数次高調波成分の補正係数
及び偶数次高調波成分の補正係数を出力する補正係数生
成手段と、第1の補正乗算器と、第1の自乗器と、第2
の補正乗算器と、第2の自乗器と、第1と第2の自乗器
の出力の加算器と、除算器とを備えている。そのため、
位相復調器の出力から伝搬光と参照光との偏光角θによ
る成分を除去することができ、最終出力を安定させるこ
とができる。第2の発明によれば、補正係数生成手段
が、高調波成分抽出手段と、除算器と、奇数次高調波成
分の補正係数及び前記偶数次高調波成分の補正係数を出
力する補正係数算出器とを備えている。そのため、第1
の補正乗算器には、奇数次高調波成分の補正係数、第2
の補正乗算器には、偶数次高調波成分の補正係数をそれ
ぞれ出力されるので、位相復調器の出力から伝搬光と参
照光との偏光角θによる成分を除去することができる。
第3の発明によれば、位相復調器が、低域フィルタを備
えている。そのため、奇数次高調波成分の補正係数及び
偶数次高調波成分の補正係数の誤差を取り除くことがで
きる。従って、位相復調器の出力から伝搬光と参照光と
の偏光角θによる成分を一層確実に除去することができ
る。As described above in detail, according to the first aspect, the phase demodulator outputs the correction coefficient of the odd-order harmonic component and the correction coefficient of the even-order harmonic component. , A first correction multiplier, a first squarer, and a second
, A second squarer, an adder for the outputs of the first and second squarers, and a divider . for that reason,
The component due to the polarization angle θ between the propagation light and the reference light can be removed from the output of the phase demodulator, and the final output can be stabilized. According to the second aspect, the correction coefficient generation unit includes a harmonic component extraction unit, a divider, and a correction coefficient calculator that outputs a correction coefficient of an odd-order harmonic component and a correction coefficient of the even-order harmonic component. It is equipped with a door. Therefore, the first
The correction multiplier of the odd harmonic component, the second
Since the correction multiplier of each of the above-mentioned correction multipliers outputs the correction coefficient of the even-order harmonic component, the component due to the polarization angle θ between the propagation light and the reference light can be removed from the output of the phase demodulator.
According to the third aspect, the phase demodulator includes the low-pass filter . Therefore, it is possible to remove an error between the correction coefficient of the odd-order harmonic component and the correction coefficient of the even-order harmonic component. Therefore, the component due to the polarization angle θ between the propagation light and the reference light can be more reliably removed from the output of the phase demodulator.
【図1】本発明の実施例を示す干渉型光ファイバセンサ
位相変調器の構成図である。FIG. 1 is a configuration diagram of an interferometric optical fiber sensor phase modulator showing an embodiment of the present invention.
【図2】従来の干渉型光ファイバセンサ位相変調器の構
成図である。FIG. 2 is a configuration diagram of a conventional interference type optical fiber sensor phase modulator.
11 第1次高調波成分抽出手段 12 第2次高調波成分抽出手段 14 第1の微分器 15 第2の微分器 16 第1の乗算器 17 第2の乗算器 18 減算器 19 積分器 22 第3次高調波成分抽出手段 23a 除算器 23b 補正係数算出器 24 第1の補正乗算器 25 第2の補正乗算器 29 LPF 26 第1の自乗器 27 第2の自乗器 28 加算器 29 LPF 30 除算器 11 First Harmonic Component Extraction Means 12 Second Harmonic Component Extraction Means 14 First Differentiator 15 Second Differentiator 16 First Multiplier 17 Second Multiplier 18 Subtractor 19 Integrator 22 Second Third harmonic component extracting means 23a Divider 23b Correction coefficient calculator 24 First correction multiplier 25 Second correction multiplier 29 LPF 26 First squarer 27 Second squarer 28 Adder 29 LPF 30 Division vessel
Claims (3)
差により生じた干渉光を電気信号に変換する光電気信号
変換器と、前記FM 変調の変調周波数を基本周波数として、前記電
気信号の奇数次の高調波成分を抽出する奇数次高調波成
分抽出手段と、 前記電気信号の偶数次の高調波成分を抽出する偶数次高
調波成分抽出手段と、 前記奇数次高調波成分抽出手段により抽出された奇数次
高調波成分を微分する第1の微分器と、 前記偶数次高調波成分と前記第1の微分器の出力とを乗
算する第1の乗算器と、 前記偶数次高周波成分抽出手段により抽出された偶数次
高調波成分を微分する第2の微分器と、 前記奇数次高調波成分と前記第2の微分器の出力とを乗
算する第2の乗算器と、 前記第1の乗算器の出力と第2の乗算器の出力とを減算
する減算器と、 前記減算器の出力を積分する積分器とを備えた位相復調
器において、 前記奇数次高調波成分の補正係数及び前記偶数次高調波
成分の補正係数を出力する補正係数生成手段と、 前記奇数次高調波成分の補正係数と奇数次高調波成分と
を乗算する第1の補正乗算器と、 前記第1の補正乗算器の出力を自乗する第1の自乗器
と、 前記偶数次高調波成分の補正係数と偶数次高調波成分と
を乗算する第2の補正乗算器と、 前記第2の補正乗算器の出力を自乗する第2の自乗器
と、 前記第1と第2の自乗器の出力を加算する加算器と、 前記積分器の出力と前記加算器の出力との間で除算する
除算器とを、 設けたことを特徴とする位相復調器。And 1. A light electrical signal converter for converting an electrical signal interference light generated by the optical path difference between the FM-modulated transmitted light and reference light, the modulation frequency of the FM modulated as a fundamental frequency, said collector
An odd-order harmonic component extraction unit for extracting an odd-order harmonic component of the air signal; an even-order harmonic component extraction unit for extracting an even-order harmonic component of the electric signal; and the odd-order harmonic component extraction. A first differentiator for differentiating the odd-order harmonic component extracted by the means; a first multiplier for multiplying the even-order harmonic component by an output of the first differentiator; A second differentiator for differentiating the even-order harmonic component extracted by the component extracting means; a second multiplier for multiplying the odd-order harmonic component by an output of the second differentiator; 1. A phase demodulator comprising: a subtractor for subtracting an output of a first multiplier and an output of a second multiplier; and an integrator for integrating an output of the subtractor, wherein a correction coefficient for the odd-order harmonic component is provided. And a correction coefficient generator for outputting a correction coefficient for the even-order harmonic component. Means, a first correction multiplier for multiplying the odd-order harmonic component correction coefficient by the odd-order harmonic component, a first squarer for squaring the output of the first correction multiplier, A second correction multiplier for multiplying the correction coefficient of the even-order harmonic component by the even-order harmonic component, a second squarer for squaring the output of the second correction multiplier, the first and the second 2. A phase demodulator comprising: an adder for adding an output of a squarer of 2; and a divider for dividing an output of the integrator and an output of the adder.
又は前記偶数次高調波成分抽出手段による偶数次高調波
成分で奇数次高調波成分対又は偶数次高調波成分対を形
成するように高調波成分を抽出する高調波成分抽出手段
と、 前記奇数次高調波成分対間又は偶数次高調波成分間で除
算する除算器と、 前記除算器の出力から前記奇数次高調波成分の補正係数
及び前記偶数次高調波成分の補正係数を出力する補正係
数算出器とを、 備えたことを特徴とする請求項1記載の位相復調器。2. The method according to claim 1, wherein the correction coefficient generation unit includes an odd-order harmonic component pair by the odd-order harmonic component extraction unit or an odd-order harmonic component pair or an even-order harmonic component by the even-order harmonic component extraction unit by the even-order harmonic component extraction unit. Harmonic component extracting means for extracting a harmonic component so as to form a higher harmonic component pair; a divider for dividing between the odd harmonic component pairs or between the even harmonic components; and an output of the divider. 2. The phase demodulator according to claim 1, further comprising: a correction coefficient calculator that outputs a correction coefficient of the odd-order harmonic component and a correction coefficient of the even-order harmonic component from the second component.
る低域通過フィルタを設けたことを特徴とする請求項1
記載の位相復調器。3. A filter according to claim 1, further comprising a low-pass filter for removing a high-frequency component of the output of said adder.
A phase demodulator as described.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5285084A JP3004853B2 (en) | 1993-11-15 | 1993-11-15 | Phase demodulator |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5285084A JP3004853B2 (en) | 1993-11-15 | 1993-11-15 | Phase demodulator |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH07143073A JPH07143073A (en) | 1995-06-02 |
| JP3004853B2 true JP3004853B2 (en) | 2000-01-31 |
Family
ID=17686935
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5285084A Expired - Fee Related JP3004853B2 (en) | 1993-11-15 | 1993-11-15 | Phase demodulator |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3004853B2 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116907627B (en) * | 2023-09-13 | 2023-12-19 | 之江实验室 | Optical path difference auxiliary-based large dynamic range distributed phase sensing method and device |
| CN118961648A (en) * | 2024-09-09 | 2024-11-15 | 国网江苏省电力有限公司电力科学研究院 | Detection device and method for dissolved acetylene in oil based on optical fiber photothermal heterodyne method |
-
1993
- 1993-11-15 JP JP5285084A patent/JP3004853B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH07143073A (en) | 1995-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104964735B (en) | A kind of detecting system and demodulation method of laser phase carrier doppler vibration signal | |
| CN110411486B (en) | PGC-DCDM demodulation method insensitive to phase delay and modulation depth | |
| US5883548A (en) | Demodulation system and method for recovering a signal of interest from an undersampled, modulated carrier | |
| JPH02300623A (en) | Signal processor for optical fiber gyro | |
| CN114964001A (en) | Linear PGC demodulation signal processing method based on hybrid phase modulation | |
| JP3004853B2 (en) | Phase demodulator | |
| JP2011214921A (en) | Interference type optical fiber sensor system and calculator | |
| JP2007163941A (en) | Four-phase phase modulation circuit | |
| CN116184815B (en) | Closed-loop control method and system for eliminating reset error of fiber optic gyroscope | |
| CN111579047A (en) | A Signal Demodulation Method of Fiber-optic Vector Hydrophone | |
| JP3459399B2 (en) | Gas concentration measurement processing equipment | |
| JPH09229718A (en) | Modulation and demodulation method for interference type optical fiber sensor | |
| CN113834561A (en) | A system and method for extracting and compensating phase modulation depth in PGC phase demodulation | |
| KR100581158B1 (en) | How to measure the phase difference of light waves propagated through two paths | |
| JP2975514B2 (en) | Phase demodulator for optical fiber sensor | |
| CN1752836A (en) | Demodulation method of optical fiber interference signal | |
| JP3124177B2 (en) | Phase demodulator | |
| JPH11108669A (en) | Fiber optic gyro | |
| JP4109757B2 (en) | Fiber optic gyro | |
| CN110581737B (en) | Control method and device of IQ modulator | |
| JP3568284B2 (en) | Demodulation method and demodulation device | |
| JP2005003530A (en) | Phase detector | |
| JPH07140043A (en) | Optical fiber sensor unit | |
| US9369321B2 (en) | Increasing sensor data carrying capability of phase generated carriers | |
| JP2598157B2 (en) | Signal processing device for optical gyro |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19991102 |
|
| LAPS | Cancellation because of no payment of annual fees |