This application is based upon and claims the benefit of priority from the prior Japanese Patent Application No. 2005-187029, filed on Jun. 27, 2005, the entire contents of which are incorporated herein by reference.
1. Field of the Invention
The present invention relates to a head position control method for a disk device for suppressing the displacement of the head by disturbance and the disk device, and more particularly to a head position control method for suppressing the displacement of the head by external vibration and the disk device.
2. Description of the Related Art
In the case of a disk device, which reads and writes data on a rotating storage medium (e.g. disk medium) by a head, the head is positioned on a desired track of the storage medium and data is read from/written to the track by the head. For such a disk device, such as a magnetic disk device or optical disk device, it is extremely important to accurately position the head on the target track for improving the recording density.
One factor that interrupts this positioning is eccentricity, which is caused by the center of the servo signals disposed concentrically on the disk that is different from the rotation center of the motor. A technology to correct this eccentricity is a control method using an observer (e.g. Japanese Patent No. 3,460,795).
Eccentricity generates sinusoidal positional fluctuation synchronizing with an integral multiple of the rotation frequency. The above-mentioned observer control method can suppress such sinusoidal positional fluctuation so as to position the head accurately on the target track. However for this eccentricity correction, the frequency to be corrected must be known in advance. For example, the frequency must be an integral multiple, multiplied once or twice, of the rotation frequency.
The second factor that interrupts this positioning is vibration, which is applied to the disk device from the outside. This vibration has various waveforms, but here a method for handling sinusoidal vibration will be described. By applying the above-mentioned eccentricity correction control, it becomes possible to handle a frequency which is not an integral multiple of the rotation frequency.
In the case of the above-mentioned conventional configuration, it is based on the assumption that the frequency of the disturbance is known in advance. The external vibration to be applied, however, is unknown at the point of designing the control system, so it is impossible to know the frequency in advance. Therefore some means of detecting this unknown frequency is necessary, and only if the frequency can be detected, the positional fluctuation by external vibration can be suppressed using such a control method as the above-mentioned Japanese Patent No. 3,460,795.
A frequency estimation unit (a) estimation) 106 for estimating disturbance, and a compensator (Cd) 104 for suppressing disturbance of a predetermined frequency by adaptive control are added to this controller 102. The sum of the output Un of the controller 102 (Cn) and the output Ud of the compensator 104 (Cd) that is U, is supplied to the control target 103 (P). This frequency estimation unit 106 estimates the angular frequency ω (=2πf) of the disturbance based on the position error ‘e’, and introduces it to the transfer function of the disturbance frequency suppression of the compensator 104. The compensator 104 calculates the recurrence formula (adaptive control formula) of the sine wave from the position error ‘e’ and the estimated angular frequency ‘ω’, and calculates the compensating current output Ud.
In this way, the frequency of disturbance is detected and the unknown frequency is suppressed so that conventional eccentricity correction control can handle the disturbance of an unknown frequency in a certain range (e.g. Proceedings of the 40th IEEE Conference on Decision and Control, pp. 4909-4914, issued December 2001). As a method of estimating an unknown frequency and suppressing the disturbance of an unknown frequency, assuming the recurrence formula of the sine wave (e.g. The Institute of Electrical Engineers of Japan, Technical Meeting on Industrial Instrumentation and Control, IIC-04-70, “Frequency following type peak filter”, issued Sep. 10, 2004), or correcting the drive amount of the control target by introducing an adaptive rule based on the above-mentioned error signal (e.g. Proceedings of the 40th IEEE Conference on Decision and Control, pp. 4909-4914, issued December 2001), have been proposed. Another proposal is estimating an unknown frequency by the error signal, generating a disturbance suppression signal in the position level, correcting the error signal and inputting it to the controller (e.g. U.S. Pat. No. 6,762,902).
For example, as an example of the open loop characteristics of the system in
Recently such a disk device, in particular a hard disk drive (HDD), is installed in mobile equipment, such as a portable personal computer, portable terminal, portable telephone and portable AV equipment. When the disk device is used in this environment the disturbance frequency is unknown and multiple, so the disk device must adapt to a wide range of disturbance frequencies.
Therefore if a conventional compensator for suppressing external vibration is attached to a conventional controller, a coefficient, according to the frequency of the external vibration, must be appropriately set in the suppression signal generation formula. In the case of the system in
The above-mentioned various conventional disturbance suppression control methods determine such coefficients as the gain and the phase in advance based on experience or on estimate using the design theory of the feedback system, and the gain and the phase are fixed values in a generation formula. In such a setting of coefficients, disturbance suppression by a disturbance frequency in a certain range (e.g. estimated frequency range of the servo system) is possible.
However if the disk device is used in the above-mentioned mobile environment, the range of the unknown disturbance frequency changes depending on the operating conditions, and cannot be easily estimated. In the case of the adaptive control of a disturbance in prior art, the frequency range that can be followed up is limited, so in the above-mentioned operation environment, an adjusted gain and phase may be inappropriate, and may make the control system unstable. The level of disturbance cannot be specified either, which makes the control system unstable, and may amplify errors or may even diverge errors.
With the foregoing in view, it is an object of the present invention to provide a head position control method, head position control device and disk device for implementing disturbance suppression control adapting to disturbance frequencies in a wider range.
It is another object of the present invention to provide a head position control method, head position control device and disk device for implementing disturbance suppression control even when used in a mobile environment.
It is still another object of the present invention to provide a head position control method, head position control device and disk device for implementing disturbance suppression control adapting to disturbance frequencies and the disturbance amplitudes in a wider range.
It is still another object of the present invention to provide a head position control method, head position control device and disk device for implementing disturbance suppression control adapting to disturbance frequencies in a wider range at high precision.
To achieve these objects, the present invention is a head positioning control method for controlling the positioning of a head onto a desired track of a disk by an actuator according to a position signal, which the head has read from the disk. The control method has: a step of obtaining a control value of the actuator according to a position error between a target position and a current position acquired from the position signal; a step of estimating sequentially an external vibration frequency from a signal that is based on the position error according to an adaptive rule; a step of estimating a rotation vector of the cyclic disturbance from the signal that is based on the position error and obtaining an external vibration suppression control value using a constant for external vibration suppression; a step of sequentially correcting the constant for external vibration suppression according to each of the sequentially estimated external vibration frequencies; and a step of generating a drive value of the actuator by adding the control value and the external vibration suppression control value.
Also the present invention is a head positioning control device for controlling the positioning of a head onto a desired track of a disk by an actuator according to a position signal, which the head has read from the disk. The device has: a control unit for obtaining a control value of the actuator according to a position error between a target position and a current position acquired from the position signal; an adaptive control unit for estimating sequentially an external vibration frequency from a signal that is based on the position error according to an adaptive rule, estimating a rotation vector of the cyclic disturbance from the signal based on the position error, obtaining an external vibration suppression control value using a constant for external vibration suppression, and sequentially correcting the constant for external vibration suppression according to each of the sequentially estimated external vibration frequencies; and an addition unit for adding the control value and the external vibration suppression control value, and generating a drive value of the actuator.
The disk device of the present invention has a head for reading information from a disk, an actuator for moving the head in track crossing direction of the disk, and a control unit for obtaining a control value of the actuator according to a position error between a target position and a current position acquired by a position signal read by the head, so as to position the head onto a desired track of a disk according to the position signal. The control unit further has an adaptive control block for estimating sequentially an external vibration frequencies from a signal that is based on the position error according to an adaptive rule, estimating a rotation vector of the cyclic disturbance from the signal that is based on the position error, obtaining an external vibration suppression control value using a constant for external vibration suppression, and sequentially correcting the constant for external vibration suppression according to each of the sequentially estimated external vibration frequencies, and an addition block for adding the control value and the external vibration suppression control value, and generating a drive value of the actuator.
In the present invention, it is preferable that the step of estimating further has a step of estimating an angular frequency of the sequential disturbance from the signal based on the position error according to the adaptive rule, the step of obtaining the external vibration suppression value further has a step of converting the estimated rotation vector of cyclic disturbance into a current value as the external vibration suppression control value using an output gain as the constant for external vibration suppression, and the step of correcting further has a step of sequentially correcting the output gain according to the sequentially estimated angular frequency.
In the present invention, it is also preferable that the step of correcting further has a step of referring to a table for storing the output gain for each of a plurality of angular frequency values by the estimated angular frequency, fetching the corresponding output gain, and updating the output gain as the constant for external vibration suppression.
In the present invention, it is also preferable that the step of correcting further has a step of referring to a table for storing the output gain for each of a plurality of angular frequency values by the estimated angular frequency, fetching the corresponding two output gains, and a step of obtaining the output gain of the estimated angular frequency by the interpolation of the two output gains and updating the output gain as the constant for external vibration suppression.
It is also preferable that the present invention further has a step of delaying the external vibration suppression control value by one sample, and the step of generating further has a step of generating a drive value of the actuator by adding the control value in the current sample and the one sample-delayed external vibration suppression control value.
In the present invention, it is also preferable that the step of estimating further has a step of limiting the estimated angular frequency by an upper limit value and a lower limit value.
In the present invention, it is also preferable that the step of estimating, the step of correcting and the step of obtaining the external vibration suppression value and the step of generating are executed for every one sample, and the step of obtaining the control value is executed a plurality of times during one sample.
It is also preferable that the present invention further has a step of executing the step of estimating, the step of correcting and the step of obtaining the external vibration suppression value for each of the plurality of external vibration frequencies for each of the external vibration frequencies, and a step of adding the plurality of external vibration suppression values.
In the present invention it is also preferable that the step of executing further has a step of correcting one of the estimated external vibration frequencies so that the estimated external vibration frequencies do not overlap.
In the present invention, an accurate external vibration frequency, which is integral-compensated, can be estimated, since the external vibration frequency is estimated from a signal based on the position error according to an adaptive rule, and a highly accurate follow up control can be performed for a wide range of external vibration frequencies, since the constant of the compensator is corrected using this estimated external vibration frequency. In this way, the operation of the compensator can always be maintained at an optimum according to the value of the estimated external vibration frequency, and follow up control for a wide range of external vibration frequencies becomes possible by correcting the constant of the compensator.
The present invention will now be described in the sequence of the disk device, first embodiment of the positioning control system, second embodiment of the positioning control system, third embodiment of the positioning control system, fourth embodiment of the positioning control system, fifth embodiment of the positioning control system, sixth embodiment of the positioning control system, seventh embodiment of the positioning control system, example, and other embodiments. In the present invention, the magnetic disk device (hard disk drive) will be described as an example. However the technology described in the present invention can also be applied to other disk devices, for example such optical disk devices as CD-ROM and DVD-ROM, and such magneto-optical disk devices as MO and MD.
Disk Device
The actuator 1 is comprised of a voice coil motor (VCM), which rotates around the rotation axis. In
The magnetic head 3 is comprised of a read element and a write element. In the magnetic head 3, the read elements including the magneto-resistance (MR) element are layered on the slider, and the write elements including the write coil are layered thereon.
The position detection circuit 7 converts the position signals (analog signals) read by the magnetic head 3 into digital signals. The read/write (R/W) circuit 10 controls the read and write of the magnetic head 3. The spindle motor (SPM) drive circuit 8 drives the spindle motor 5. The voice coil motor (VCM) drive circuit 6 supplies drive current to the voice coil motor (VCM) 1 and drives the VCH 1.
The microcontroller (MCU) 14 detects (demodulates) the current position from the digital position signals from the position detection circuit 7, and computes the VCM drive command value according to the error between the detected current position and the target position. In other words, the position demodulation and the servo control are performed. The read only memory (ROM) 13 stores the control programs of the MCU 14, for example. The random access memory (RAM) 12 stores the data for processing of the MCU 14, for example.
The hard disk controller (HDC) 11 judges the position on one rotation of the disk based on the sector number of the servo signal, and reads and writes data. The random access memory (RAM) for buffer 15 temporarily stores the read data and write data. The HDC 11 communicates with the host using such an interface IF as USB, ATA and SCSI. The bus 9 connects these composing elements.
As
For example, the sector number, when the index signal is detected, is set to No. 0, and the sector number is incremented each time a servo signal is detected so as to acquire the sector number of each sector of the track. The sector number of the servo signal becomes a reference when the data is read and written. One index signal exists in each track. A sector number may be assigned instead of the index number.
In other words, in seek control, the head can be moved to the target position in the transition of coarse control, settlement control and following control. Coarse control is basically velocity control, and settlement control and following control are basically position control, and for both controls, the current position of the head must be detected.
To confirm such a position, servo signals, as shown in
The frequency estimation unit 24 estimates the angular frequency ω of the external vibration using the position error ‘e’ and the internal variable of the disturbance suppression compensator 23 (Cd). The compensation table 25 stores the constants of the compensator for external vibration suppression 23 (Cd) corresponding to each frequency ω. The disturbance suppression compensator 23 (Cd) corrects the internal constant by the constant read from the compensation table 25 based on the angular frequency XD of the frequency estimation unit 24, and calculates the disturbance suppression control amount Ud from the position error ‘e’. The addition block 26 adds the control amount Un and the disturbance suppression control amount Ud and outputs it to the plant 22.
In this way, the internal constant of the compensator Cd is corrected according to the estimated value of the angular frequency (disturbance frequency) ω, so the operation of the compensator Cd can always be maintained at an optimum in a wide frequency range.
First the control block (controller) 21 will be described. In the design of the positioning control system of the magnetic disk device, the control target plant, that is the actuator, is often expressed by Expression (1).
In Expression (1), ‘u’ is the current to be supplied to the actuator, ‘y’ is the observed position, B1 is the force constant, ‘m’ is an equivalent mass and ‘s’ is a Laplace operator. The magnetic disk device is generally a rotary actuator, but when designing a controller, it is converted into an equivalent linear actuator, as shown in Expression (1).
A typical control of the magnetic disk device is observer control. The transfer function of the actuator is expressed by the equation of state shown in Expression (2).
‘x’ and ‘y’ are positions, ‘v’ is velocity, and ‘u’ is the drive current. Z-transformation is performed on this analog model at sample period T, and the units are converted into the units for digital control of the magnetic disk device. Specifically, position is in track units, velocity is in track/sample units, and current is in units normalized such that the maximum current becomes 1. Then the equation of state of the actuator becomes Expression (3).
Ka is the gain, including B1/m.
Using this expression, observer control is implemented. The observer control used for the digital control of the magnetic disk device has a configuration called the current observer in the following Expression (4).
As Expression (4) shows, the state variables x_hat [k] and v_hat [k] are corrected by the estimated position error of the observer, that is the difference between the observed position y [k] and the estimated position x_bar [k] calculated one sample before. Then the two state variables x_hat [k] and v_hat [k] are multiplied by the feedback gains Fx and Fv, and the current u [k] is output. After outputting the current, the state variable values of the next sample, x_bar [k+1] and v bar [k+1] are estimated.
Using this Expression (4), the controller 21 (Cn) executes observer control.
The above-mentioned observer corresponds to the PD control in analog control. Integral compensation is also required in the magnetic disk device. A method for adding control equivalent to integrated compensation is bias compensation. A stationary bias ‘b’ given by the following Expression (5) is added to the analog plant model.
s·b=0 (5)
Then the current observer is configured by using the expression (3) and dicretized expression (5), and the following Expression (6) is acquired.
In the above description, a known observer controller was used as the controller, but as another control method, PID control, having the characteristics shown in the following Expression (7), may be used.
In Expression (7), Ki is the integration gain, Kp is the proportion gain, Kd is the differential gain, and s is a Laplace operator.
This control is also possible with a controller where PI control and a LeadLag filter are connected in series. The transfer function is shown in the following Expression (8).
In this Expression (8), Ki is the integration gain, Kp is the proportion gain, a1 and b1 are the characteristics of the Lead Lag filter, and s is a Laplace operator.
In recent magnetic disk devices, the actuator is controlled not by analog control but by digital control. The controller 21 described in Expressions (7) and (8) is based on analog control expressions, but these expressions are converted for use in digital control, and are installed in the device.
The configuration of the controller 21 is described by the configuration of PID control, PI control x Lead Lag filter and observer control. The present invention can be applied not only to the above-mentioned controller 21, but also to a linear controller 21. For example, a filter designed by H-infinity control theory may be used for the controller.
Also the actuator of the magnetic disk device has resonance characteristics, so the vibration of the actuator, due to resonance, influences the resonance position accuracy. Therefore in many cases, a plurality of notch filters based on the following Expression (9) is cascade-connected to suppress resonance. Actually Expression (9) is discretized and used.
The compensator 23, table 25 and ω estimator 24 are described in
Now the configuration of the compensator 23 (Cd) of digital control will be described. When the disturbance frequency is known, the compensator 23 of digital control, which follows up cyclic disturbance, is designed by the equation of state in Expression (10).
In Expression (10), (x1, x2) is the state variable of the rotation vector of disturbance, ‘e’ is the position error, ‘L’ is the input gain, Udist is the disturbance suppression output current Ud, F is the output gain, ω is the angular frequency of a known disturbance, and T is the sample period of the control system. For example, the output gain F is designed to be the optimum at the center of frequency of the fluctuation range of the disturbance frequency to be assumed, and once designed, this constant F is the constant. For example, if the original characteristic of the controller, before attaching the compensator for disturbance suppression, is C and the characteristic of the plant is P, then the transfer function of the following Expression (11) is used.
P/(1+CP) (11)
The output gain F is adjusted considering the phase characteristic with respect to the angular frequency ω in the transfer function of this Expression (11) and the convergence response time with respect to a known disturbance.
In Expression (10), an element of the input gain L is “0” or “1”, but another fixed value may be used. The constant design of the gain F is decided using simulation. However a calculation delay, the resonance of the actuator, the insertion of various filters represented by a notch filter for suppressing the resonance, and the frequency characteristics of the amplifier for driving VCM, for example, have actual influence. Therefore the final value of F is checked and adjusted while actually operating the device.
In the description here, a position error is used as the input. The estimated position error of the observer described in
In Expression (10), the angular frequency ω of the disturbance is a fixed value. In this invention, the angular frequency ω of the disturbance is sequentially corrected to adapt to an unknown disturbance frequency. The adaptive rule for this will now be described. The adaptive rule is derived by transforming the equation of state of the compensator 23 of the cyclic disturbance of the discrete system in Expression (10), and the continuous system of Expression (12) is used.
In Expression (12), x1 is the sine component of disturbance, x2 is the cosine component of disturbance, L1 and L2 are the input gain, ‘e’ is the position error, and s is a Laplace operator.
For the angular frequency ω of the rotation vector, the rotation vector is the tangent of the sine component x1 and the cosine component x2, and the angular frequency ω is a differential value thereof, so the angular frequency ω of this rotation vector is given by the following Expression (13).
The equation of state of the analog compensator in Expression (12) is expanded, and the differential value of x1 and x2 is determined, and the differential value of x1 and x2 is substituted in Expression (13), then the following Expression (14) is acquired.
If the angular frequency of the estimated unknown disturbance is correct, the compensator 23 (Cd) can appropriately suppress the disturbance. As a result, the position error ‘e’ or the estimated position error of the observer becomes zero. In other words, in Expression (14), the conditions when the target angular frequency ω of the compensator 23 (Cd) and the angular frequency ω of the estimated disturbance are such that the term of the position error ‘e’ in the right side number in Expression (13) is zero. In other words, Expression (15) is established.
In other words, the angular frequency ω is sequentially corrected so as to satisfy Expression (15). For this, Expression (15) is expressed by the time differential form and the adaptive rule (integral compensation rule) of the following Expression (16).
Using this expression, the value of ω is sequentially corrected. Expression (16) is corrected to an integral form, and is expressed by a digital control expression, then the following Expression (17) is acquired.
Ka is the adaptive gain. In Expression (17), the adaptive rule is used in the addition format. This can be expressed in multiplication format, as shown in the following Expression (18).
In Expression (18), the adaptive gain Ka2 is a gain that is a different value from the above-mentioned adaptive gain Ka.
In the following description, the adaptive rule in the addition format shown in Expression (17) will be used, but the adaptive rule in the multiplication format shown in Expression (18) can also be used just the same.
First the target frequency of the compensator, with respect to the cyclic disturbance, is corrected using the adaptive rule, then the angular frequency changes along with time k. Therefore the compensator that follows up the cyclic disturbance, with respect to a known disturbance frequency shown in Expression (10), is configured such that the parameters are updated according to the changes of the angular frequency. In other words, the input gain L for adjusting the phase and the magnitude of the input, the matrix G for rotating the vector, and the output gain (matrix) F for adjusting the phase and the magnitude of the output are also updated according to the sequential correction. By this, the accuracy of the sequential correction improves. In other words, the adaptive controller of Expression (19) is designed by matching Expression (10) with the adaptive rule.
L, F [k] and ω[k] however are given by Expression (20) or Expression (21) according to the matrix L.
The matrix L is constant in Expression (19) to Expression (21), and as the adaptive rule ω[k] is updated, the matrixes G and F are changed to G[k] and F[k]. In other words, the angular frequency changes as time elapses, so the cosine (ω[k]T) and the sine (ω[k]T) of G and the value of the matrix F are also updated accordingly. Here Mag and Phs are the gain (magnitude) and the phase with respect to the angular frequency ω in the transfer function (shown in Expression (10)) from the current disturbance to the position error. And by changing to the gain Mag (ω[k]) and the phase Phs (ω[k]) according to ω[k], the phase and the gain of the output are adjusted by the output gain F[k].
For the input gain L, L1 and L2 are set to sine and cosine respectively as shown in Expression (20) and Expression (21), and in the adaptive control Expressions (20) and (21), one of the two variables of L is set to “1” and the other to “0”. Either one of the input side L and the output side F is sufficient to adjust the gain and the phase, and the above-mentioned adaptive control must not only suppress disturbance but also estimate the angular frequency ω of the disturbance.
The internal state variables x1 and x2 are used for estimation. To simplify the estimation formula of ω, one of L1 and L2 is set to “1” and the other is set to “0”, and are substituted in Expression (17), then the adaptive rule is transformed to the adaptive rule (ω[k]) of Expression (20) if L1=0 and L2=1, and is transformed to the adaptive rule (ω[k]) of Expression (21) if L1=1 and L2=0. By this, the operation time can be decreased.
For certain both values of the two variables in each matrix L and matrix F may be changed according to ω[k]. For example, L=(L1(ω[k]) and L2(ω[k]) can be set. The output side gain F may be fixed and the input side gain L may be changed according to ω[k]. However it takes time to compute ω[k] compared with Expression (20) and (21).
x1 and x2 are shifted 90 degrees, so by setting the matrix L, the matrix of sine and the matrix of the cosine of the terms of the phase of F[k] become opposite thereof as shown in Expression (20) and Expression (21).
Returning to
The adder 26 adds the control current value Un[k] of the controller 21 and the external disturbance suppression current value Ud[k], outputs the result to the plant 22 (VCM1), and drives it.
This will be described in more detail with reference to
Table 25, on the other hand, stores the values G and F according to the value of each estimated ω, as shown in
Also as the output gain F, F1 (1·ωr)−F1(n·ωr) and F2(1·ωr)−F2(n·ωr) of ω[k]=1·ωr−n ωr are stored. In Expression (20), F1 is=−Mag (ω[k])·cos(Phs (ω[k]), and F2 is=−Mag (ω[k])·sin(Phs(ω[k]). This is the same for the case of Expression (21).
The compensator 23 computes x1[k+1], x2[k+1] and Udist [k] in Expression (18). In other words, the compensator 23 comprises the multiplication block 23-1 for multiplying the position error e[k] by “L”, the multiplication block 23-4 for multiplying x1[k] and x2[k] by G[k] from the table 25, the addition block 23-2 for adding the outputs of both multiplication blocks 23-1 and 23-4 and outputting x1[k+1] and x2[k+1], the delay block 23-3 for delaying the output of the addition block 23-2 by one sample and outputting x1[k] and x2[k], and the multiplication block 23-5 for multiplying the output x1[k] and x2[k] of the delay block 23-3 by the output gains F1[k] and F2[k] from the table 25.
In
SIN=(1−COS2)1/2 (22)
It must be remembered that in Expression (22) that the SIN component, and not the COS component, is deleted. In this case, values up to the Nyquist frequency must be represented. In other words, as
Therefore to decrease the table size, it is preferable to have the COS values remain. In the case when such a consideration is unnecessary, the COS components may be deleted and the SIN components remain.
As
In
As the prior art in
In the simulation result when the present invention is applied, as shown in
Therefore by applying the present invention, disturbance can be accurately and quickly suppressed, even if an unknown disturbance frequency is provided. For example, as the simulation result in
The frequency estimation unit 24 estimates the angular frequency o of the external vibration using the position error ‘e’ and the internal variable of the disturbance suppression compensator 23 (Cd). The compensation table 25 stores the constants of the compensator 23 (Cd) for the extended vibration suppression corresponding to each frequency ω. The disturbance suppression comparator 23 (Cd) corrects the internal constant by the constant read from the compensation table 25 using the angular frequency ω of the frequency estimation unit 24, and calculates the disturbance suppression control amount Ud from the position error ‘e’. The addition block 26 adds the control amount Un and the disturbance suppression control amount Ud, and outputs it to the plant 22. The delay unit 27 delays the output Ud of the compensator 23 by one sample and outputs it to the addition block 26.
If one sample is delayed, a shift in the phase characteristics (phase lag) is generated. Therefore the set value of the output gain F is different from the values in
The disturbance suppression output is intentionally delayed by one sample, as shown in
Expressions (17) to (21) for adaptive control of the present invention, on the other hand, require extra calculation processing time compared to prior art. In other words, as
Then the table 25 is referred to by this estimated angular frequency ω, and G and F of the compensator 23 are updated. And the output Ud (Udist) of the compensator (Cd) 23 is calculated by Expression (19). Finally the output U[k] of the plant 22 is calculated by Un[k]+Ud[k]. This calculated value is output to the plant 22, and changes the output current. In other words, time from the acquisition of the observed position y to the output of the drive current U becomes long and the response of the control system drops (in particular a phase lag is generated). This delay is certainly not guaranteed if a high-speed MCU 14 is installed. However to provide a low cost disk device, installing a high-speed MCU 14 is not the best solution.
Therefore the value ω[k] of the current sample is not used for the current for disturbance suppression, but a value ω[k−1], which estimated one sample before, is used, and the correction current Ud is calculated. In other words, as
Then the adaptive rule of the frequency estimation unit 24 is calculated based on Expression (17) or Expression (18), and the estimated angular frequency ω[k] is updated. Then the table 25 is referred to by this estimated angular frequency ω[k], and G and F of the compensator 23 are updated. And by Expression (19), the output Ud (Udist) [k] of the compensator (Cd) 23 is calculated by Expression (19). In this way, the time from the acquisition of the observed position ‘y’ to the output of the drive current U decreases, and the response of the control system improves (phase lag is prevented).
Also the estimated values ω[k−1] and Ud [k−1] with the previous sample may be used, as mentioned above, since the response of the control system (phase accuracy) becomes quicker, which improves disturbance frequency estimation speed.
(S10) First the MCU 14 acquires the observed position y from the head 3, and calculates the position error ‘e’.
(S12) Then the MCU 14 calculates the output Un[k] of the controller (Cn) 21 by Expression (6).
(S14) Then the MCU 14 calculates the output U[k] of the plant 22 by Un[k]+Ud[k−1] using the Ud[k−1] calculated with the previous sample. This calculated value U is output to the plant 22, and the output current changes.
(S16) The MCU 14 then executes calculation of adaptive control. First the update calculation of ω is executed. In other words, the estimated angular frequency ω[k] is updated by calculating the adaptive rule of the frequency estimation unit (24) in Expression (17) or Expression (18).
(S18) Then the MCU 14 refers to tables 25, and updates the values G and F.
(S20) Then the MCU 14 performs the update calculation of the state variables x1[k +1] and x2[k +1] in Expression (19), and at the same time calculates the correction current Ud[k] of the next sample.
In
The table 25 stores the values G and F at each predetermined frequency of the angular frequency ω, as shown in
The interpolation-computing unit 28 reads the two G and F closest to the estimated angular frequency ω[k] output from the estimation unit 24 from the table 25, performs linear interpolation, and calculates G[k] and F[k] corresponding to the estimated angular frequency ω[k].
The comparator 23 computes x1[k+1], x2[k+1] and Udist[k] of Expression (18). In other words, the compensator 23 comprises a multiplication block 23-1 for multiplying the position error e[k] by “L”, a multiplication block 23-4 for multiplying x1[k] and x2[k] by G[k] from the interpolation computing unit 28, an addition block 23-2 for adding the outputs of both the multiplication blocks 23-1 and 23-4 and outputting x1[k+1] and x2[k+1], a delay block 23-3 for delaying the output of the addition block 23-2 by one sample and outputting x1[k] and x2[k], and the multiplication block 23-5 for multiplying the outputs x1[k] and x2[k] of the delay block 23-3 by the output gains F1[k] and F2[k] from the interpolation computing unit 25.
The configuration of the table 25 may be the one in
By this configuration, a disturbance suppression control the same as the first embodiment can be implemented with decreasing the size of the table 25, which is effective for installation.
In other words, the value ω does not become “0” or a negative value. And it is unnecessary to estimate a frequency higher than the Nyquist frequency, which is an angular frequency corresponding to a frequency that is half of the sample frequency, since the present invention is implemented by digital control. Also because of such an influence as the calculation delay, a disturbance around the Nyquist frequency is difficult to handle. Also an extremely low frequency disturbance, such as 1 Hz, can be handled by the integral compensation of the controller 21, so this need not be dealt with. Therefore as
In
The table 25 stores the values G and F at each predetermined frequency of the angular frequency ω, as shown in
The interpolation computing unit 28 reads the two G and F values which are closest to the estimated angular frequency ω[k] from the ω estimation unit 24 in the table 25, performs linear interpolation, and calculates G[k] and F[k] corresponding to the estimated angular frequency ω[k].
The compensator 23 computes x1[k +1], x2[k +1] and Udist[k] of Expression (19). In other words, the compensator 23 comprises the multiplication block 23-1 for multiplying the position error e[k] by “L”, a multiplication block 23-4 for multiplying x1[k] and x2[k] by G[k] from the interpolation computing unit 28, an addition block 23-2 for adding the output of both multiplication blocks 23-1 and 23-4 and outputting x1[k+1] and x2[k+1], a delay block 23-3 for delaying the output of the addition block 23-2 by one sample and outputting x1[k] and x2[k], and a multiplication block 23-5 for multiplying the outputs x1[k] and x2[k] of the delay block 23-3 by the output gains F1[k] and F2[k] from the interpolation computing unit 28.
The configuration of the table 25 may be the one in
By this configuration, the estimated angular frequency is limited by the upper limit ω max (e.g. half of the sampling frequency=Nyquist frequency), and is also limited by the lower limit ω min (e.g. rotation frequency of the disk). Therefore excessive disturbance suppression control can be prevented and stable disturbance suppression control can be implemented. Also the upper limit may be a lower limit frequency of the fluctuation range due to the individual difference in the resonant frequency of the actuator and the temperature change.
(S30) First the MCU 14 acquires the observed position ‘y’ from the head 3, and calculates the position error ‘e’.
(S32) Then the MCU 14 calculates the output Un[k] of the controller (Cn) 21 by Expression (6).
(S34) Then the MCU 14 calculates the output U[k] of the plant 22 by Un[k]+Ud[k−1] using the Ud[k−1] calculated with the previous sample. This calculation value U is output to the plant 22, and the output current changes.
(S36) The MCU 14 then executes calculation of adaptive control. First the update calculation of ω is performed. In other words, the adaptive rule of the frequency estimation unit 24 in Expression (17) or Expression (18) is performed, and the estimated angular frequency ω[k] is updated.
(S40) Then the MCU 14 judges whether the estimated angular frequency ω[k] is the lower limit a) min or less. If the estimated angular frequency ω[k] is the lower limit ω min or less, the estimated angular frequency ω[k] is ω[k]=ω min.
(S42) Then the MCU 14 refers to the table 25 and executes the above-mentioned linear interpolation to update the values G and F.
(S44) Then the MCU 14 performs the update calculation of the state variables x1[k+1] and x2[k+1] of Expression (19), and at the same time calculates the correction current Ud[k] of the next sample.
In other words, in the description of the above-mentioned embodiments, the number of frequencies of unknown disturbance is one. However depending on the operating format of the disk device, two different frequencies may be causing an unknown disturbance. For example, when the disk device is installed in a portable device, it is installed via an elastic element. In such an installation example, the disk device not only receives disturbance directly but also receives disturbance via the elastic element that is the disk device receives two different frequency disturbances even if there is one disturbance. Therefore suppression control is performed for the disturbance of two different unknown frequencies.
In
Two disturbance suppression control mechanisms are installed. In other words, the first disturbance suppression control mechanisms 23A, 24A and 25 and the second disturbance suppression control mechanisms 23B, 24B and 25 are installed. The frequency estimation units 24A and 24B estimate the angular frequencies ω1 and ω2 of the external vibration, using the position error ‘e’ and the internal variables of the disturbance suppression compensators 23A and 23B (Cd). The compensation table 25 stores the constants of the compensators 23A and 23B (Cd) for the external vibration suppression corresponding to each frequency ω1 and ω2. The disturbance suppression compensators 23A and 23B (Cd) correct the internal constants by the constants read from the compensation table 25 by the angular frequencies ω1 and ω2 of the frequency estimation units 24A and 24B, and calculate the disturbance suppression control amounts Ud1 and Ud2 from the position error ‘e’.
The addition block 29 adds the two disturbance suppression control amounts Ud1 and Ud2. The addition block 26 adds the control amount Un and the disturbance suppression control amount Ud (=Ud1+Ud2), and outputs it to the plant 22.
In other words, two correction values, Ud1 and Ud2, are added and output. Here the angular frequency ω is corrected between the first disturbance suppression control mechanisms 23A, 24A and 25, and the second disturbance suppression control mechanisms 23B, 24B and 25. In other words, if adaptive control is performed, the estimated angular frequencies ω1 and ω2 may become the same. Then only one disturbance can be supported. Therefore, as mentioned later, the values of the estimated angular frequencies ω1 and ω2 are adjusted as not to overlap.
The table 25 stores the values G and F at each predetermined frequency of the angular frequency ω, as shown in
The interpolation computing unit 28A/28B reads the two G and F values which are closest to the estimated angular frequency ω1[k] or ω2[k] from the ω estimation unit 24A/24B in the table 25, performs linear interpolation, and calculates G1[k], G2[k], FA[k]=F1[k], F2[k], FB[k]=F1[k], F2[k] corresponding to the estimated angular frequency ω1[k] or ω2[k].
The compensator 23A/23B computes x1[k+1], x2[k+1], Udist1[k] and Udist2[k] in Expression (19). In other words, the compensator 23A/23B comprises a multiplication block 23-1A/23-1B for multiplying the position error e[k] by [LA=L1, L2], [LB=L1, L2], a multiplication block 23-4A/23-4B for multiplying x1[k] and x2[k] by G1[k] and G2[k] from the interpolation computing unit 28A/28B, an addition block 23-2A/23-2B for adding the outputs of both multiplication blocks 23-1A/23-B and 23-4A/23-4B, and outputting x1[k+1] and x2[k+1], a delay block 23-3A/23-3B for delaying the output of the addition block 23-2A/23-2B by one sample and outputting x1[k] and x2[k], and a multiplication block 23-5A/23-5B for multiplying the outputs x1[k] and x2[k] of the delay block 23-3A/23-3B by the output gains FA[k] and FB[k] from the interpolation computing unit 28A/28B.
The addition block 29 adds the outputs Ud1 and Ud2 of the two multiplication blocks 23-5A and 23-5B, and outputs the disturbance suppression control value Ud[k]. The configuration of the table 25 may be the one in
Also in the second ω estimation section 24B, a ω2 correction block 24-5 for correcting the estimated ω2 by the output ω1 of the delay block 24-2A of the first ω estimation block 24A is installed. In other words, when the adaptive control configuration for two unknown frequencies like this is used, it must be remembered that there is the case when ω1 and ω2 always take a same value, and that virtually only one disturbance can be handled while two disturbance suppression control mechanisms are provided.
To solve this problem, a means of adjusting ω1 and ω2, so that the values ω1 and ω2 do not overlap, is installed. In other words, the ω2 correction block 24-5 compares the updated ω2 with ω1, and corrects ω2 if these values are similar. Details will be described in the processing flow in
(S50) First the MCU 14 acquires the observed position y from the head 3, and calculates the position error ‘e’.
(S52) Then the MCU 14 calculates the output Un[k] of the controller (Cn) 21 by Expression (6).
(S54) Then the MCU 14 calculates the output U[k] of the plant 22 by Un[k]+Ud[k−1] using the Ud[k−1] calculated with the previous sample. This calculation value U is output to the plant 22, and the output current changes.
(S56) The MCU 14 then executes the calculation of adaptive control. First the update calculation of ω1 is performed. In other words, the adaptive rule of the frequency estimation unit 24A in Expression (17) or Expression (18) is performed, and the estimated angular frequency ω1[k] is updated.
(S58) Then the MCU 14 refers to the table 25, executes the above-mentioned linear interpolation, and updates the values G1 and FA.
(S60) Then the MCU 14 performs the update calculation of state variables x1[k+1] and x2[k+1] of Expression (19), and at the same time calculates the correction current Ud1[k] of the next sample.
(S62) Then the MCU 14 performs the update calculation of the estimated angular frequency ω2[k]. First ω[k−1] of the previous sample is set for the variable ω2_OLD. Then the MCU 14 executes the adaptive rule of the frequency estimation unit 24B of Expression (17) or Expression (18), and updates the estimated angular frequency ω2[k].
(S64) The MCU 14 compares the updated ω2[k] with ω1[k], and if both are close, ω2[k] is corrected. In other words, Ax for defining a predetermined frequency range is set for the first estimated angular frequency ω1, and it is judged whether the second estimated angular frequency ω2[k] is smaller than (ω1+Δω) and greater than (ω1−Δω). This range from (ω1+Δω) to (ω1−Δω) is the range to judge whether ω2 is close to ω1. If the second estimated angular frequency ω2[k] is not smaller than (ω1+Δω) or is not greater than (ω1−Δω), then ω2 and ω1 are not close, so it is not necessary to correct ω2, and the processing advances to step S68.
(S66) If the second estimated angular frequency ω2[k] is smaller than (ω1+Δω) and greater than (ω1−Δω), it is judged that ω2 is close to ω1, and ω2[k] is corrected. That is, in order to judge whether ω2 is smaller than or greater than ω1, it is judged whether ω1 is ω_OLD[k], where ω2[k−1) of the previous sample is set, or less. If ω1 is greater than ω_OLD[k] where ω2[k−1] is set, ω2 is smaller than ω1, so ω2[k] is corrected to (ω1−Δω), that is the lower limit of the tolerance. If ω1 is not greater than ω_OLD[k], where ω2[k−1] is set, ω2 is greater than ω1, so ω2[k] is corrected to ω1+Δω), that is to the upper limit of the tolerance.
(S68) Now the MCU 14 refers to the table 25 by ω2[k], executes the above-mentioned linear interpolation, and updates the values G2 and FB.
(S70) Then the MCU 14 performs the update calculation of the state variables x1[k+1] and x2[k+1] in Expression (19), and at the same time calculates the correction current Ud2[k] of the next sample. And the MCU 14 adds Ud1[k] in step S60 and Ud2[k], and determines Ud[k].
In
In this way the value Δω is set. This value is the minimum value of the difference between ω2 and ω1. If the difference is smaller than Δω, the value of ω2 is shifted. This shift has direction. The shift in the positive direction and the shift in the negative direction must be distinguished.
Generally an improvement of positioning accuracy can be expected if the sample frequency is increased. In the magnetic disk device, however, increasing the sample frequency increases the recording area of servo signals, which decreases the area where data is recorded. To improve this, multi-rate control may be used for the magnetic disk device. Multi-rate control can also be implemented by using a filter. In the magnetic disk device, a method of implementing multi-rate control using observer control is generally used.
In the present embodiment, the configuration of a controller of the observer control using multi-rate control is used. In other words, multi-rate control is a method for implementing the output of a frequency higher than the sample frequency of the input. In this example, a double rate of multi-rate control is constructed. In other words, current U is output at double the frequency of the detected sample frequency at the observed position ‘y’, which is input.
The controllers 21A and 21B of observer control using the multi-rate control in
As Expression (23) shows, basically Expression (6) is executed twice for the detected sample at the observed position y. u[k] and u[k+0.5] are calculated, and current u[k] and u[k+0.5] are output with double frequency.
In other words, as
Specifically a clock (multi-rate sample clock) 30 with a double frequency of the detected sample frequency at the observed position ‘y’ is supplied to the controllers 21A and 21B, and Un1 (U[k]) and Un2 (U[k+0.5]) are output as Un synchronizing with this clock.
For the blocks 23, 24 and 25 of adaptive control, the disturbance suppression control value Ud is calculated, as mentioned above, using the estimated position error of the observer 21A. In this multi-rate control as well, the blocks 23, 24 and 25 of adaptive control are operated once per one sample, and the disturbance suppression control value Ud is calculated. This is because calculation processing becomes heavy and the processing time becomes long, and also suppression performance is not much different between operating once per one sample and twice per one sample. The addition block 26 adds Un and Ud, and outputs it to the plant 22.
In other words, as shown in
(S80) First the MCU 14 acquires the observed position ‘y’ from the head 3, and calculates the position error ‘e’.
(S82) Then the MCU 14 calculates the output Un1[k] of the controller (Cn1) 21A by Expression (33).
(S84) Then the MCU 14 calculates the output U1[k] of the plant 22 by Un1[k]+Ud[k−1] using Ud[k−1] calculated with the previous sample. This calculated value U1 is output to the plant 22, and the output current changes.
(S86) Then the MCU 14 calculates the output Un2[k] (=Un[k+0.5]) of the controller (Cn2) by Expression (33).
(S88) Then the MCU 14 calculates the output U2[k] of the plant 22 by Un2[k]+Ud[k−1] using Ud[k−1] calculated with the previous sample. This calculated value U2 is output to the plant 22 at a specified time, and the output current changes.
(S90) The MCU 14 then executes the calculation of adaptive control. First the update calculation of a is performed. In other words, the adaptive rule of the frequency estimation unit 24 of Expression (17) or Expression (18) is calculated, and the estimated angular frequency ω[k] is updated.
(S92) Then the MCU 14 refers to the table 25 and updates the values of G and F.
(S94) Then the MCU 14 performs the update calculation of the state variables x1[k+1] and x2[k+1] in Expression (19), and at the same time calculates the correction current Ud[k] of the next sample.
One or plurality of cascade-connected notch filters 31 suppresses the resonance of the actuator 1. In other words, an actuator has resonance characteristics. Therefore by using a filter represented by a notch filter, the controller output is passed through the filter and the filter output is supplied to the actuator.
(S100) First the MCU 14 acquires the observed position ‘y’ from the head 3, and calculates the position error ‘e’.
(S102) Then the MCU 14 calculates the output Un1[k] of the controller (Cn1) 21A by Expression (33).
(S104) Then the MCU 14 calculates the output U1[k] of the plant 22 by Un1[k]+Ud[k−1] using Ud[k−1] calculated with the previous sample.
(S106) The MCU 14 performs filter calculation processing for passing the output U1[k] through the notch filter 31 in Expression (9), and this calculated value U1 is output to the plant 22, and the output current changes.
(S108) Then the MCU 14 calculates the output Un2[k] (=Un[k+0.5]) of the controller (Cn2) 21A by Expression (33).
(S110) Then the MCU 14 calculates the output U2[k] of the plant by Un2[k]+Ud[k−1] using Ud[k−1] calculated with the previous sample.
(S112) The MCU 14 performs filter calculation processing for passing the output U2[k] through the notch filter 31 in Expression (9), and this calculated value U2 is output to the plant 22, and the output current changes. (S114) The MCU 14 then executes the calculation of adaptive control. First the update calculation of o is performed. In other words, the adaptive rule of the frequency estimation unit 24 of Expression (17) or Expression (18) is calculated, and the estimated angular frequency ω[k] is updated.
(S116) Then the MCU 14 refers to the table 25, and updates the values of G and F.
(S118) Then the MCU 14 performs the update calculation of the state variables x1[k+1] and x2[k+1] in Expression (19), and at the same time calculates the correction current Ud[k] of the next sample.
In the above description, adaptive control is started in a state where external vibration is applied when following up a track, so as to suppress the disturbance. In an actual magnetic disk, however, a same track is not always followed up, but seek control is also executed.
The operation timing of the adaptive control of the present invention when control is switched in this way will be described.
The angular frequency estimation unit 24 estimates the angular frequency ω of external vibration, using the position error e and the internal variable of the disturbance suppression compensator 23 (Cd). The compensation table 25 stores the constants of the compensator 23 (Cd) for external vibration suppression corresponding to each frequency ω. The disturbance suppression compensator 23 (Cd) corrects the internal constants using the constants read from the compensation table 25 by the angular frequency ω of the frequency estimation unit 24, and calculates the disturbance suppression control amount Ud from the position error ‘e’. The addition block 26 adds the control amount Un and the disturbance suppression control amount Ud, and outputs it to the plant 22.
This controller 21 is a block for performing track following. This positioning control system further comprises a known seek control block 40, settlement control block 41 and switching block 42 for switching the seek control block 40, settlement control block 41 and controller (track following control block) 21 according to the position error, so as to perform seek control and settlement control, as shown in
In this configuration, when the head moving control in
Therefore the load of processing on the MCU 14 can be decreased, and the track following control of the head during read/writing can be executed with suppressing disturbance.
In other words, as the configuration of the control system in
When it is judged that the positioning accuracy has dropped due to a disturbance in the above-mentioned monitoring, the positioning accuracy monitoring block 43 sends a startup instruction to both switches 44 and 45, and turns the switches 44 and 45 ON. By this, the position error ‘e’ is input to the compensator 23 and the ω estimation block 24, and as mentioned above, the disturbance suppression control value Ud is output from the compensator 23 to the addition block 26 via the switch 45. Therefore the output U to the plant 22 (1, 3) during track following changes from the output Un of the controller 21 to (Un+Ud), and disturbance suppression control is performed.
In this case, if the positioning accuracy monitoring block 43 judged that positioning accuracy improved, the switches 44 and 45 can be turned OFF.
By this, the load of processing on the MCU 14 can be further decreased, and when disturbance drops the positioning accuracy, the track following control of the head during a read/write operation can be executed with suppressing disturbance.
The output Ud of the adaptive control 23 is always supplied, but the input of the adaptive controls 23 and 24 is supplied during track following. By this configuration, adaptive control can be executed without being subject to the influence of a major position error generated during seek response.
In other words, as
On the other hand, the addition block 26 is installed between the switching block 42 and the plant 22 (1, 3). Therefore the disturbance suppression control value Ud is added to the seek control value and the settlement control value even when track following is not performed (when seek or settlement control is performed), and is output to the plant 22 (1, 3). For the output U to the plant 22 (1, 3) during track following, on the other hand, the position error ‘e’ is input, and Ud, according to the position error ‘e’, is added to the output Un of the controller 21 calculated by the adaptive control blocks 23, and 24, so as to perform disturbance suppression control.
When track following is not performed (when seek or settlement control is performed), the value “0” (position error e=0) is input to the adaptive control blocks 23 and 24 by the switch 44, so the angular frequency X remains as the estimated value in track following, as shown in Expression (20). Therefore the seek operation can be started at track following, without being subject to the influence of a major position error (number of seek tracks) during seeking, even if disturbance cannot be completely suppressed.
In other words, the disturbance suppression control value Ud is added to a seek control value and settlement control value and is output to the plant 22 (1, 3). In this configuration, the positioning control system can be used continuously when heads are switched.
Examples of the present invention will now be described.
The top of
Now examples using an actual magnetic disk device will be described.
For the magnetic disk device, a 2.5-inch device of which rotation frequency is 70 Hz (4200 rpm) was used. The graphs show the response when pseudo-vibration is generated by applying a sinusoidal disturbance at the position y by a program of the MCU 14 of this magnetic disk device. The disturbance frequencies are shifted from an integral multiple of the rotation frequency. The result is expanded as numeric data on the memory of the circuit, and after operation, the values are fetched and shown as a graphic.
In
In the same way, in
This convergence time can be adjusted by the adaptive gain Ka and F in Expression (20) and Expression (21). As the adaptive gain Ka or the value F increases, the frequency follows up more quickly, and the convergence time, to match the disturbance frequency, decreases.
In all the above described examples, disturbance is suppressed by starting adaptive control in a status where an external vibration is applied during track following. However in an actual magnetic device, not only is a same track always followed, but seek control as well is executed to move to another track or to switch heads.
In the above embodiments, the disk device was described using a magnetic disk device, but the present invention can also be applied to other disk device, such as an optical disk deice and a magneto-optical disk device, and the shape of the disk is not limited to circular, but may be square, for example. Program control by the MPU was described, but individual hardware may be used.
The present invention was described above using embodiments, but the present invention can be modified in various ways within the scope of the essential character thereof, and these shall not be excluded from the scope of the present invention.
Since the external vibration frequency is estimated by an adaptive rule using signals based on the position error, the accurate external vibration frequency can be estimated based on integral compensation, and since the constants of the compensator are sequentially corrected using this external vibration frequency, highly accurate follow up control is possible for a wide range of external vibration frequencies. In the same way, the operation of the compensator can always be kept at an optimum according to the value of the estimated external vibration frequency, and in particular the present invention contributes to guaranteeing normal operation of the disk drive, even in an environment where external vibration is easily applied.
Number | Date | Country | Kind |
---|---|---|---|
2005-187029 | Jun 2005 | JP | national |