Abstract:
An NMR apparatus provided with a lock device having improved monotonicity without reducing a control range. In a lock device including a lock transmitter, a lock receiver, a lock detector, a lock corrector, a lock power supply and a lock coil, the lock detector performs complex detection of a detection signal and outputs to the lock corrector an absorption signal SA which is a real number component of the detection signal and a dispersion signal SD which is an imaginary number component of the detection signal; the lock corrector obtains a value SC by multiplying together SD/(SD 2 +SA 2 ) and a value G, and the lock power supply performs adjustment so that a current flowing through the lock coil is proportional to SC. Different constants G are selectively used according to the value of SD/(SD 2 +SA 2 ).

Description:
BACKGROUND OF THE INVENTION  
       [0001]     The present invention relates to a magnetic field stabilizer in a nuclear magnetic resonance spectrograph (hereinafter referred to as “NMR spectrograph”) and a method of stabilizing an NMR magnetic field with the stabilizer.  
         [0002]     NMR spectrographs are being used to measure various specimens including organic materials, solid materials and biological macromolecules. As described in a non-patent document 1 (Gerhard Wider, Technical aspect of NMR spectroscopy with biological macromolecules and studies of hydration in solution, Progress in Nuclear Magnetic Resonance Spectroscopy, vol. 32, p 225, 1998”), an apparatus of this kind measures a resonance frequency of a specimen placed in a static magnetic field and it is required that this static magnetic field be stable and have a rate of change of 1 ppb/hr or less with respect to time. To produce a stable static magnetic field, most of the existing NMR spectrographs are equipped with a superconducting magnet. Even in a superconducting magnet, however, it is virtually impossible to achieve a stability of 1 ppb/hr. A lock device for stabilizing a static magnetic field by feedback control is important in enabling NMR spectrographs to deliver the desired performance. With the development of studies on proteins requiring long-time measurement processes, the importance of such a lock device is increasing.  
         [0003]     A non-patent document 2 (Benedict W. Bangerter, Field/frequency lock monitor for signal averaging with high resolution NMR sepectrometers, Review of Scientific Instruments, vol. 46, p 617, 1975) discloses a technique used for the lock device. An alternating-current magnetic field having a radio frequency is radiated from a lock transmitter to a specimen placed in a magnet for use in an NMR apparatus to excite a magnetization vector in the specimen. The motion of the excited magnetization vector is received through an antenna for detection and is detected by using a receiving circuit. A component called a dispersion signal (SD) is obtained by adjusting the phase of a reference signal used at the time of detection.  
         [0004]     The dispersion signal SD is multiplied by a predetermined proportion constant G and the result of this multiplication is delivered to a lock power supply. The lock power supply causes a current proportional to G×SD (correction signal SC) to flow through a coil provided around the specimen (hereinafter referred to as “lock coil”) to form a correction magnetic field proportional to the correction signal SC around the specimen. The proportion constant G is determined so that the dispersion SD obtained after adding the correction magnetic field is close to zero, thereby stabilizing the static magnetic field.  
         [0005]     As described above, the technique disclosed in non-patent document 2 is characterized by stabilizing the static magnetic field by controlling the current flowing through the lock coil, and the lock corrector can be easily implemented, for example, by using only an analog circuit. The above-described lock detector is also capable of obtaining an absorption signal SA having a 90° phase difference from the dispersion signal SD, but this signal is not used for computation of the correction signal SC.  FIG. 5  shows an example of the ideal dispersion signal SD and the absorption signal SA. A detailed description of the dispersion signal SD and the absorption signal SA will be made below.  
         [0006]     A patent document 1 (JP-B-2504666) discloses a technique for the lock device, which differs from that disclosed in non-patent document 2 in that the absorption signal SA is used together with the dispersion signal SD to compute the correction signal SC output from the lock corrector. The correction signal SC in patent document 1 is a weighted combination of SD/SA and (1/SA)(dSD/dt).  
         [0007]     According to patent document 1, SD/SA is perfectly proportional to a frequency offset DF=ω-ω 0  due to a disturbance in the static magnetic field when the frequency of the disturbance is low, and the proportion constant is 1/T 2 . T 2  is a time period relating to the rate of attenuation of the detection signal VD and called a lateral relaxation time. According to patent document 1, T 2  can also be obtained in the control loop and DF=(1/T 2 )(SD/SA) can therefore be computed and used as the correction signal SC in the lock control section.  
         [0008]     According to patent document 1, the frequency offset DF is also equal to (1/T 2 )×(SD/SA)+(1/SA)×(dSD/dt). This equation is accurate and has the advantage of being unrelated to any condition and the advantage of including a noise component. Therefore, the technique disclosed in patent document 1 is suitable for use in a case where the disturbance amplitude is comparatively large.  
         [0009]     The technique disclosed in patent document 1 has an advantage over the technique disclosed in non-patent document 2 in that the magnetic field can be returned to a stable value by a smaller number of control loops in comparison with the technique because of the proportional relationship between the frequency offset DF and the correction signal SC and good monotonicity.  
         [0010]     The technique disclosed in non-patent document 2 still has a problem in terms of monotonicity of the correction signal. As is apparent from the graph of  FIG. 5  showing the dispersion signal SD, the dispersion signal SD is divided into three regions having boundaries corresponding to the maximum and minimum values of the dispersion signal SD, and the sign of the change rate dSC/dDF is reversed between adjacent two of the regions. Accordingly, a plurality of frequency offsets DF are associated with one dispersion signal SD. Since the correction signal SC in non-patent document 2 is G×SD, the correction signal apparently lacks monotonicity. Therefore, it is difficult to obtain the desired constant G and a number of control loops are required to return the static magnetic field to a stable value, resulting in a reduction in control speed.  
         [0011]     The technique disclosed in patent document 1 increases the control speed by solving the monotonicity problem with the technique disclosed in non-patent document 2. However, since division by the absorption signal SA is introduced, there is a possibility of the correction signal SC becoming unstable when the frequency offset is increased and when the absorption signal becomes smaller. This is because when the magnitude of the absorption signal SA becomes approximately equal to the magnitude of noise, the sign of the correction signal SC can be reversed by the influence of noise. Also, the extent of variation appeared in the correction signal SC due to an error a in receiving phase is larger than that in the case of the technique disclosed in non-patent document 2. Accordingly, the range of frequency offset DF controllable by using the technique disclosed in non-patent document 1 is narrower than that in the case of the technique disclosed in non-patent document 2. The technique disclosed in patent document 1 is advantageous in terms of monotonicity but has a problem that noise and instability with respect to an error in receiving phase are increased and a problem that the controllable range is reduced.  
         [0012]     It is an object of the present invention to provide a high-performance NMR apparatus by implementing a lock device having improved monotonicity without reducing the control range.  
       SUMMARY OF THE INVENTION  
       [0013]     The present invention provides, as a means for achieving the object of the present invention, a lock device for an NMR apparatus constituted by a lock transmitter, a lock receiver, a lock detector, a lock corrector, a lock power supply and a lock coil, wherein the lock detector performs complex detection of a detection signal, and outputs to the lock corrector an absorption signal SA which is a real number component of the detection signal and a dispersion signal SD which is an imaginary number component of the detection signal; the lock corrector obtains a value SC by multiplying together SD/(SD 2 +SA 2 ) and a predetermined constant G; and the lock power supply performs adjustment so that a current flowing through the lock coil is proportional to SC.  
         [0014]     A lock device having improved monotonicity without reducing the control range is thereby implemented to enable provision of a high-performance NMR apparatus.  
         [0015]     Other objects, features and advantages of the invention will become apparent from the following description of the embodiments of the invention taken in conjunction with the accompanying drawings. 
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0016]      FIG. 1  is a block diagram showing the entire configuration of a lock device in an embodiment of the present invention;  
         [0017]      FIG. 2  is a block diagram showing the configuration of a lock corrector;  
         [0018]      FIG. 3  is a flowchart showing the operation of the lock corrector;  
         [0019]      FIG. 4  is a timing chart showing the operation of the lock device;  
         [0020]      FIG. 5  is a graph showing an ideal dispersion signal SD and absorption signal SA for explanation of the present invention;  
         [0021]      FIG. 6  is a graph of a correction signal obtained by using the present invention and setting G=1;  
         [0022]      FIG. 7  is a graph of a correction signal obtained by using the technique disclosed in non-patent document 2 and setting G=1;  
         [0023]      FIG. 8  is a graph of a correction signal obtained by using the technique disclosed in patent document 1 and setting G=1; and  
         [0024]      FIG. 9  is a flowchart showing the operation of a lock corrector in another embodiment of the present invention. 
     
    
     DESCRIPTION OF THE INVENTION  
     EMBODIMENT 1  
       [0025]      FIG. 1  is a block diagram showing an embodiment of a lock device in accordance with the present invention. A lock transmitter  10  transmits an alternating current having a frequency f to a switch  12 . The frequency f is, for example, 92 MHz. The switch  12  is an active or passive electric circuit having the function of connecting/shutting off two circuits: a circuit for connection for transmission from the lock transmitter  10  to a probe  14  and a circuit for connection for transmission from the probe  14  to the lock receiver  16 .  
         [0026]     When the lock transmitter  10  is operating, the current flows through the circuit for transmission from the lock transmitter  10  to the probe  14 . The probe  14  is inserted in a magnet  18  for use in an NMR apparatus. A specimen (not shown) on which a measurement is performed is placed inside the magnet  18 . The alternating-current signal that has entered the probe  14  via the switch  12  produces an alternating-current magnetic field having the frequency f through a radiation antenna (not shown) provided in the probe  14 . The magnetic field is radiated to the specimen placed in the probe  14 .  
         [0027]     The frequency f is a resonance frequency which a predetermined reference atomic nucleus carries at a target static magnetic field intensity. As the reference atomic nucleus, heavy hydrogen is ordinarily used. In some case, an atomic nucleus of, for example, fluorine other than heavy hydrogen is used. The spin of the reference atomic nucleus in the specimen is excited by the alternating-current magnetic field to produce a motion. The microscopic spin motion appears as the motion of a macroscopic magnetization vector and induces a voltage of the frequency f (hereinafter referred to as “detection signal VD”) in a detection antenna provided in the probe  14 .  
         [0028]     The radiation antenna may also be used as the detection antenna. The detection signal VD is sent to the lock receiver from the detection antenna via the switch  12 . The lock receiver  16  performs one or both of frequency conversion and amplitude conversion on the detection signal VD and transmits the converted signal (hereinafter referred to as the reception signal VR) to a lock detector  20 .  
         [0029]     The lock detector  20  performs complex detection by multiplying the reception signal VR by a reference signal VF (hereinafter referred to as VF) having a predetermined frequency, phase and amplitude, thereby forms two signals called a dispersion signal SD and an absorption signal (SA) and delivers these signals to a lock corrector  22 .  
         [0030]     The lock corrector  22  multiplies SD 2 /(SD 2 +SA 2 ) by a predetermined constant G to obtain a value as a correction signal SC and delivers the correction signal SC to the lock power supply  24 . The lock power supply  24  causes a current proportional to the correction signal SC to flow through a lock coil  26  provided around the specimen, thereby applying a correction magnetic field proportional to the correction signal SC around the specimen. Accordingly, the magnitude of the current flowing through the lock coil  26  is adjusted as a function of the dispersion signal SD and the absorption signal SA. The operations of the lock transmitter  10 , the switch  12 , the lock receiver  16 , the lock detector  20  and the lock corrector  22  are controlled by a lock controller  27 .  
         [0031]      FIG. 2  is a block diagram showing a concrete configuration of the lock corrector  22 . The lock corrector  22  is constituted by a computation section  28  in which the correction signal SC is obtained by using the dispersion signal SD and the absorption signal SA supplied from the lock detector  20 , and a correction section  30  from which the correction signal SC obtained in the computation section  28  is output to the lock power supply  24 .  
         [0032]      FIG. 3  is a flowchart showing the operation of the computation section  28 . The operation shown in  FIG. 3  is started when the computation section  28  receives a control start instruction from the lock controller  27 . When control is started, a register in the computation section  28  is initialized and a control loop is started. Subsequently, a loop gain G and a control range discrimination threshold value J are updated on the basis of a value supplied from the lock controller  27 . Default values for the loop gain G and the control range discrimination threshold value J are set at the time of register initialization. If there is no input supplied from the lock controller  27 , the initial values remain as they are.  
         [0033]     Subsequently, the dispersion signal SD and the absorption signal SA are updated to values supplied from the lock detector  20 , and weight averaging is performed by using the updated dispersion signal SD and absorption signal SA. The order of weight averaging and weight coefficients are set by a user through the lock controller  27 . The order may be set to 1 to make weight averaging ineffective. If the S/N of the lock signal is high, weight averaging may be made ineffective to improve the response of the lock device.  
         [0034]     When a weight-averaged dispersion signal SD′ and a weight-averaged absorption signal SA′ are obtained, A 2 =SD′ 2 +SA′ 2  and SC=G×SD′/A 2  are successively obtained. A 2  is compared with the control range discrimination threshold value J. If A 2  is smaller than the control range discrimination threshold value J, it is determined that a variation is out of the control range, and transition to control loop termination processing is made. If A 2  is equal to or larger than the control range discrimination threshold value J, the correction signal SC is output to the lock power supply  24  through the corrector  30 . Comparison between A 2  and J may be made immediately after obtaining A 2 .  
         [0035]     Subsequently, a check is made as to whether or not a control termination instruction has been provided from the lock controller  27 . The control termination instruction issued during control is stored by using, for example, a queue memory so as not to be lost during control and is executed after output of the correction signal SC. If the control termination instruction has not been provided, the process returns to the initial point in the control loop to repeat the control loop. If the control termination instruction has been provided, transition to loop termination processing is made and control is terminated when the termination processing is completed.  
         [0036]      FIG. 4  is a timing chart used for implementation of the present invention. Transmission (TX) of the alternating current signal from the lock transmitter  10 , receiving operation (RX) of the lock receiver  16 , transmission/reception change (TR-SW) by the switch  12 , operation (ADC) of an analog-digital converter which belongs to the lock detector  20  and timing (LCC) of operation of the lock power supply  24  are shown in this order from the top.  
         [0037]     The abscissa of  FIG. 4  represents time and the ordinate of  FIG. 4  represents the voltage of a trigger signal used to designate timing of each operation. Each of the abscissa and the ordinate is of an arbitrary scale. In  FIG. 4 , LCC is performed at intervals of a TX-RX period. However, a timing setting can may be easily selected such that LCC is performed one time in several TX-RX periods.  
         [0038]     A method of computing the correction signal SC in the present invention will be described. The atomic nucleus spin of the specimen placed in the magnetic field filed is oriented to form a macroscopic magnetization vector M. It is well known that, under a static magnetic field (0, 0, B 0 ) applied in the z-axis direction and a high-frequency magnetic field (2B 1  cos(ft), 0, 0) applied in the x-axis direction, the magnetization vector M=(Mx, My, Mz) exhibits a motion expressed by equations shown below.  
             [     Expression   ⁢           ⁢   1     ]                                     ⅆ     M   x         ⅆ   t       =         -     M   x       ⁢     ω   1     ⁢   sin   ⁢           ⁢     (     ω   ⁢           ⁢   t     )       -       M   y     ⁢           ⁢     ω   1     ⁢   cos   ⁢           ⁢     (     ω   ⁢           ⁢   t     )       -         M   z     -     M   0         T   1                         ⅆ     M   x         ⅆ   t       =         M   y     ⁢     ω   0       +       M   x     ⁢           ⁢     ω   1     ⁢   sin   ⁢           ⁢     (     ω   ⁢           ⁢   t     )       -       M   x       T   2                         ⅆ     M   y         ⅆ   t       =         M   z     ⁢     ω   1     ⁢   cos   ⁢           ⁢     (     ω   ⁢           ⁢   t     )       -       M   x     ⁢           ⁢     ω   0       -       M   y       T   2                       ω   =     2   ⁢   π   ⁢           ⁢   f       ,           ⁢       ω   0     =     γ   ×     B   0         ,           ⁢       ω   1     =     γ   ×     B   1                       (   1   )             
 
         [0039]     In equations (1), γ is an intrinsic physical property value of the reference atomic nucleus called a gyromagnetic ratio; T 2  is the above-mentioned lateral relaxation time; and T 1  is the longitudinal relaxation time relating to the rate of attenuation of the detection signal VD like T 2 . Magnetic field intensities B 0  and B 1  are expressed by ω 0  and ω 1  as shown in equation (1) at the final line because they appear together with the gyromagnetic ratio γ at all times.  
         [0040]     It is well known from an electromagnetic law that the detection signal VD induced in the detection coil by the motion of the magnetization vector M is VD=−C×dMx/dt when the detection coil is placed on the x-axis. C is a proportion constant determined by the characteristics of the detection coil, the positional relationship between the detection coil and the specimen and other factors. If equations (1) can be solved, the detection signal VD can be obtained by time differentiation of the Mx solution.  
         [0041]     Since the object of the present invention is to provide a lock device used to stabilize the static magnetic field intensity B 0 , it is desirable that equations (1) be solved under a condition in which ω 0  changes randomly with time. However, it is impossible to solve equations (1) under a condition in which ω 0  changes randomly with time. According to the present invention, therefore, an approximation is made such that the magnetization vector M is in a steady state. To obtain a solution, equations (1) are rewritten by using a rotational coordinate system having a rotational speed of ω. Mx and My which are the x- and y-components of the magnetization vector M in the fixed coordinate system have relations shown by equations (2) below with u and v which are the x- and y-components of M in the above-descrived rotational coordinate system.  
         [0000]     [Expression 2]
 
 M   x   =u  cos ω t−v  sin ω t  
 
 M   y   =−u  sin ω t−v  cos ω t   (2) 
 
         [0042]     Equations (1) are rewritten by using the rotational coordinate system to obtain equations (3)  
             [     Expression   ⁢           ⁢   3     ]                                     ⅆ   u       ⅆ   t       =         -     (       ω   0     -   ω     )       ⁢           ⁢   v     -     u     T   2                         ⅆ   v       ⅆ   t       =         -     ω   1       ⁢     M   z       +       (       ω   0     -   ω     )     ⁢           ⁢   u     -     v     T   2                         ⅆ     M   z         ⅆ   t       =         ω   1     ⁢   v     -         M   z     -     M   0         T   1                       (   3   )             
 
         [0043]     If an approximation is made such that the magnetization vector M is in a steady state, each of the differential terms on the left side of equations (3) is zero. Equations (4) are the results of solving equations (3) under a condition du/dt=dv/dt=dMz/dt=0.  
             [     Expression   ⁢           ⁢   4     ]                                 u   =       M   0     ⁢         (     ω   -     ω   0       )     ⁢           ⁢     ω   1     ⁢     T   2   2         1   +         (     ω   -     ω   0       )     2     ⁢     T   2   2       +       ω   1   2     ⁢     T   1     ⁢     T   2                         v   =       M   0     ⁢         -     ω   1       ⁢     T   2         1   +         (     ω   -     ω   0       )     2     ⁢     T   2   2       +       ω   1   2     ⁢     T   1     ⁢     T   2                           M   z     =       M   0     ⁢       1   +         (     ω   -     ω   0       )     2     ⁢           ⁢     T   2   2           1   +         (     ω   -     ω   0       )     2     ⁢     T   2   2       +       ω   1   2     ⁢     T   1     ⁢     T   2                           (   4   )             
 
         [0044]     Equations (2) and the law of electromagnetic induction are used to obtain an ideal detection signal VD 0  induced in the detection antenna placed in the x-axis direction, as shown below.  
             [     Expression   ⁢           ⁢   5     ]                                   VD   0     =       -   ω     ⁢           ⁢     C   ⁢           [       u   ⁢           ⁢   sin   ⁢           ⁢     (     ω   ⁢           ⁢   t     )       +     v   ⁢           ⁢   cos   ⁢           ⁢     (     ω   ⁢           ⁢   t     )         ]                   =     ω   ⁢           ⁢     CA   0     ⁢   cos   ⁢           ⁢     (       ω   ⁢           ⁢   t     +   ϕ     )                     (   5   )             
 
         [0045]     In equation (5), C is a proportion constant determined by the characteristics of the detection coil, the positional relationship between the detection coil and the specimen and other factors, A 0   2 =u 2 +v 2 , and Tan(φ)=−u/v. In the actual detection signal VD, noise η and a phase deviation φ in the circuit from the detection antenna to the lock detector are added to VD 0 .  
         [0046]     When the lock detector  20  performs detection, the phase of the reference signal VF is adjusted so that the phase deviation φ is cancelled out. In theory it is possible to completely cancel out the phase deviation φ. In actuality, however, it is impossible to do so and a small phase error α is allowed to remain. Therefore, the dispersion signal SD and the absorption signal SA output from the lock detector  20  are accompanied by the phase error α and noise η as shown in equation (7).  
             [     Expression   ⁢           ⁢   6     ]                             (         SD           SA         )     =       ω   ⁢           ⁢     C   ′     ⁢     A   0     ⁢           ⁢     (           sin   ⁢           ⁢     (     ϕ   +   α     )                 cos   ⁢           ⁢     (     ϕ   +   α     )             )       +     η   ⁢           ⁢     (         1           1         )                 (   6   )             
 
         [0047]     In equation (6), C′ is a value obtained by multiplying C in equation (5) by a change in signal amplitude due to the gain of the receiving circuit. To obtain the desired expression of the correction signal SC, equation (6) is rewritten into equation (7) shown below by substituting u and v in equation (6) and using a matrix.  
             [     Expression   ⁢           ⁢   7     ]                                   (         SD           SA         )     =       ⁢     ω   ⁢           ⁢     C   ′     ⁢     M   0     ⁢               ⁢       ω   1     ⁢     T   2           1   +         (     ω   -     ω   0       )     2     ⁢     T   2   2       +       ω   1   2     ⁢     T   1     ⁢     T   2                             ⁢         (           cos   ⁢           ⁢     (   α   )             sin   ⁢           ⁢     (   α   )                   -   sin     ⁢           ⁢     (   α   )             cos   ⁢           ⁢     (   α   )             )     ⁢           ⁢     (             (     ω   -     ω   0       )     ⁢           ⁢     T   2               1         )       +     η   ⁢           ⁢     (         1           1         )                       =       ⁢     ω   ⁢           ⁢     C   ′     ⁢     M   0     ⁢         ω   1     ⁢     T   2         f   ⁢           ⁢     (     ω   0     )     ⁢       1   +         (     ω   -     ω   0       )     2     ⁢     T   2   2               ⁢     (           cos   ⁢           ⁢     (   α   )             sin   ⁢           ⁢     (   α   )                   -   sin     ⁢           ⁢     (   α   )             cos   ⁢           ⁢     (   α   )             )         ⁢                           ⁢         (             u   0     ⁡     (     ω   0     )                   v   0     ⁡     (     ω   0     )             )     +     η   ⁢           ⁢     (         1           1         )         ,                     f   ⁢           ⁢     (     ω   0     )       =       ⁢     1   +         ω   1   2     ⁢     T   1     ⁢     T   2         1   +         (     ω   -     ω   0       )     2     ⁢     T   2   2               ,                     u   0     ⁡     (     ω   0     )       =       ⁢         (     ω   -     ω   0       )     ⁢           ⁢     T   2           1   +         (     ω   -     ω   0       )     2     ⁢     T   2   2               ,           ⁢         v   0     ⁡     (     ω   0     )       =     1       1   +         (     ω   -     ω   0       )     2     ⁢     T   2   2                             (   7   )             
 
         [0048]     Equation (7) is considered under a condition w 1   2 T 1 T 2 &lt;&lt;1. Then f(ω 0 )=1 in equation (7). The dependence of the dispersion signal SD and the absorption signal SA on ω 0  is only u 0 (ω 0 ) and v 0 (ω 0 ). For example, if no phase error and no noise exist, (ω-ω 0 )=ωC′M 0 ω 1 ×SD/(SD 2 +SA 2 ) and it is possible to immediately achieve correction of the static magnetic field by equalizing the loop gain G to 1/(ωC′M 0 ω 1 ) and obtaining the correction signal SC by the method of the present invention.  FIG. 5  is a graph showing the ideal dispersion signal SD and absorption signal SA.  FIG. 6  shows an example of the correction signal SC according to the present invention. SC 0  on the ordinate is the value of the correction signal SC when the loop gain G is 1.  
         [0049]      FIGS. 7 and 8  show the values of the correction signal SC 0  obtained by using the techniques disclosed in non-patent document 2 and patent document 1. From comparison between  FIG. 6  and  FIGS. 7 and 8 , it is apparent that the correction signal SC 0  obtained by the present invention has improved monotonicity through a wide ω-ω 0  range. That is, according to the present invention, the correction signal SC changes linearly through a wide range.  
         [0050]     When the phase error α is not zero, the correction signal SC obtained by the method of the present invention is as shown below.  
             [     Expression   ⁢           ⁢   8     ]                                   ⁢           SC   =       ⁢     G   ×     SD   /     (       SD   2     +     SA   2       )                     =       ⁢     G   ×       1       ω   1     ⁢     T   2         ⁡     [       cos   ⁢           ⁢     (   α   )     ⁢           ⁢       T   2     ⁡     (     ω   -     ω   0       )         +     sin   ⁢           ⁢     (   α   )         ]                     =       ⁢       G   ×   cos   ⁢           ⁢     (   α   )     ⁢           ⁢       (     ω   -     ω   0       )     /     ω   1         +     G   ×   sin   ⁢           ⁢       (   α   )     /     (       ω   1     ⁢     T   2       )                       ≈       ⁢       G   ×       (     ω   -     ω   0       )     /     ω   1         +     G   ×     α   /     (       ω   1     ⁢     T   2       )                           (   8   )             
 
         [0051]     At the final line of equation (8), the correction signal SC is expressed by assuming that the phase error α is sufficiently small and by using Cosine and Sine to the primary term expanded by Taylor expansion with respect to the phase error α. As shown by using equation (8), the method of computing the correction signal SC in accordance with the present invention ensures stable control even when the receiving phase error α is not zero.  
       EMBODIMENT 2  
       [0052]     Another embodiment of the lock device of the present invention will be described with reference to the flowchart shown in  FIG. 9 . In the lock device having the configuration shown in  FIGS. 1 and 2 , the computation section  28  receives a control start instruction from the lock controller  27  and starts control, and starts the control loop after initializing the register in the computation section  28 .  
         [0053]     Subsequently, a table GA and a table J for the loop gain are updated on the basis of a value supplied from the lock controller  27 . If there is no input supplied from the lock controller  27 , default values for the table GA and the table J set at the time of register initialization are used for computation of the correction signal SC subsequently performed. The dispersion signal SD and the absorption signal SA are then updated to values supplied from the lock detector  20 , and weight averaging is performed by using the values of the updated dispersion signal SD and absorption signal SA.  
         [0054]     The order of weight averaging and weight coefficients are set by a user through the lock controller  27 . The order may be set to 1 to make weight averaging ineffective. If the S/N of the lock signal is high, weight averaging may be made ineffective to improve the response of the lock device. When a weight-averaged dispersion signal SD′ and a weight-averaged absorption signal SA′ are obtained, A 2 =SD′ 2 +SA′ 2  and SC 0 =SD′/A 2  are successively obtained. After obtaining SC 0 , a suitable loop gain value G is selected from the table GA for the loop gain according to the value of SC 0 .  
         [0055]     The correspondence relationship between the loop gain G and the correction signal SC 0  is, for example, G=X−X/(SC 0   2 +Y). X and Y are constants designated by the user. A suitable value may be selected as the loop gain G from GA according to a correspondence relationship such as that mentioned above, or the loop main G may be directly computed without using GA. After obtaining the loop gain G, the correction signal SC=G×SC 0  is obtained. Subsequently, A 2  is compared with the control range discrimination threshold value J. If A 2  is smaller than J, it is determined that a variation is out of the control range, and transition to control loop termination processing is made. If A 2  is equal to or larger than J, the correction signal SC is output to the lock power supply  24  through the corrector  30 . Comparison between A 2  and J may be made at any stage between the step of obtaining A 2  and the step of outputting the correction signal SC.  
         [0056]     Subsequently, a check is made as to whether or not a control termination instruction has been provided from the lock controller  27 . The control termination instruction issued during control is stored by using, for example, a queue memory so as not to be lost during control and is executed after output of the correction signal SC. If the control termination instruction has not been provided, the process returns to the initial point in the control loop to repeat the control loop. If the control termination instruction has been provided, transition to loop termination processing is made and control is terminated when the termination processing is completed. If the method of computing the correction signal SC in this embodiment is used, the correction signal SC can be obtained with high accuracy regardless of the condition w 1   2 T 1 T 2 &lt;&lt;1.  
         [0057]     It should be further understood by those skilled in the art that although the foregoing description has been made on embodiments of the invention, the invention is not limited thereto and various changes and modifications may be made without departing from the spirit of the invention and the scope of the appended claims.