System identification and mechanical resonance frequency suppression for servo control used in single gimbal control moment gyroscope.

Effective identification of the control model is one of the key aspects in improving the performance of the single gimbal control moment gyroscope (SGCMG) servo system. The accuracy and stability of the servo system can be improved by studying system identification and mechanical resonance frequency. In this study, firstly, the SGCMG gimbal servo system was simplified to a two-mass block model. The theoretical mathematical model of the system's transfer function and mechanical resonance frequency was derived. Secondly, this paper studied the effective suppression method for mechanical resonance. Thirdly, the mathematical model of the orthogonal correlation analysis method was deduced for system identification. Then, an experimental platform was investigated to obtain the frequency characteristic curve and the transfer function. Finally, the frequency characteristic curve obtained using the transfer function model was plotted and compared with the frequency characteristic curve obtained experimentally. Our results indicate that the orthogonal correlation analysis has a high identification accuracy.

1. Introduction

The single gimbal control moment gyroscope (SGCMG)has the advantages of a large output torque, simple structure, frameless locking, and simple control, because of which it is widely used in spacecraft. Its application in spacecraft attitude control has always been an active area of research [1-3]. The rotor used in the SGCMG has a constant angular momentum. Therefore, the control accuracy mainly depends on the accuracy of the gimbal servo control system. The need to improve the accuracy of servo control system has led to renewed interest in the development of servo control algorithms in recent studies, including the closed-loop I/f control scheme [4], the sliding mode control method using a new approach combined with an iterative learning controller and extended state observer [5], and torque ripple minimization by modulation of the phase current [6].

The motor selected for the SGCMG servo system is a permanent magnet synchronous motor (PMSM) [7,8], which is a complex nonlinear system with strong coupling and multivariable characteristics. Therefore, it is necessary to establish an accurate mathematical model before studying the servo control algorithms. There are two methods to establish a mathematical model, one is by theoretical analysis [9-11], the other is by system identification through experiment. System identification obtains the dynamic characteristics of the system through its input and outpu to establish the mathematical model of the system. One of the most common problems in servo motion control systems is mechanical resonance [12,13]. Mechanical resonance is caused by the flexible connection between two or more components during mechanical transmission. Among them, the most critical part is caused by the flexible connection between the motor and the load. Mechanical resonance affects the control accuracy of a precision servo system and the mechanical resonance frequency determines the upper limit of the system bandwidth. When the mechanical resonance frequency is low and the system bandwidth requirement is high, it may cause the entire servo system to resonate and make the system unable to operate normally. In operation, strong resonance can even lead to bearing fracture and damage to the entire servo system [14]. Therefore, it is of considerable importance to identify control system models that can predict the mechanical resonance frequency. Kang et al. [15] indicated mechanical resonance often degrades the performance of the servo system and have employed shifted discrete Fourier translations to determine the properties of the mechanical resonance. Chen et al. [16] used a binary multifrequency excitation test signal to perform system identification. Yang et al. [17] explained the cause of mechanical resonance and have presented a tracking scheme using the velocity error and bandpass filters to track the mechanical resonance frequency.

In this study, the mechanical resonance frequency of the servo control was identified using orthogonal correlation analysis and a method for suppressing the mechanical resonance frequency was developed. Owing to the relative complexity of the PMSM and the difficulty in establishing an accurate system model through theoretical derivation and simulation, identifying the precise system transfer function is a prerequisite for accurate system control. The literature includes a strong body of work on the identification of the PMSM system model. Previous studies [18-22] identified the parameters of PMSM and then constructed the system model through theoretical model derivations. However, some simplifications and trade-offs are often made in the theoretical derivation process, and the real system is subjected to many external sources of interference. Therefore, theoretical model derivations would always be prone to a certain degree of error. Consequently, the determination of an accurate method of system model identification is still an active research subject in various fields of industrial research. An identification method was devised using the iterative process of the linearized and weighted total least-squares method [23]. Wen et al. [24] used a pseudo-random binary sequence as the input signal, acquired the impulse response based on correlation analysis, and derived the state space model from the impulse response through singular value decomposition and the frequency-domain model, obtaining the explicit values of variable parameters by parameter fitting. Ishak et al. [25] studied system identification using a fractional-order model because real systems and dynamical processes display fractional-orderbehavior.

Previous studies on system identification have not considered the identification accuracy and suppression of the mechanical resonance frequency. In this study, we conducted a system identification experiment to establish an accurate system model. Simultaneously, we developed a method for suppressing the mechanical resonance frequency. Our experiment verified the effectiveness of the mechanical resonance frequency suppression method and the identification algorithm. We obtained the amplitude–frequency and phase–frequency characteristics of the system through the orthogonal correlation analysis method. Further, we obtained the transfer function and mechanical resonance frequency of the entire system by fitting the frequency characteristics.

The contribution of this study is the development of a system identification method that considers the mechanical resonance frequency and the system model. Moreover, methods for suppressing the mechanical resonance frequency were developed. In addition, we demonstrated that the orthogonal correlation analysis has a high identification accuracy.

The remainder of this paper is organized as follows. Section 2 describes the model of the SGCMG servo system and its mechanical resonance; it also details the method of mechanical resonance frequency suppression. Section 3 derives the orthogonal correlation analysis method for system identification. Section 4 describes the principle and the elements of the system identification experiment. Section 5 presents detailed results for the mechanical resonance frequency and the servo system model identification. Section 6 concludes the study.

2. Modeling and mechanical resonance frequency suppression of the SGCMG servo system

2.1 Modeling of the SGCMG servo system

The relationship between the motor assembly and the load of the SGCMG system is illustrated in Fig 1, where the shaft assembly, flywheel connector, and flywheel as a load are driven by the motor assembly.

Fig 1

SGCMG system structure diagram.

To facilitate the analysis of the flexible coupling relationship between the motor and the load, we recognize that it is, in principle, equivalent to a two-mass model. Fig 2 depicts the equivalent model diagram.

Fig 2

Servo system two-mass model.

The command current i receives the feedback current i through the current controller, and the feedback current is multiplied by the motor torque constant to obtain the electromagnetic torque T. The electromagnetic torque directly drives the PMSM to generate the motor acceleration A. The motor acceleration is integrated with respect to time to obtain the motor speed V, and the motor speed is integrated with respect to time to obtain the motor position θ. When the motor starts running, the two-mass spring deforms and tightens so that the load produces a torque proportional to the difference in position of the motor and the load. The load torque causes the load to accelerate, producing an acceleration A, the load acceleration is integrated to obtain the load speed V, and the load speed is integrated to obtain the load position θ. The torque transmitted by the motor causes the load to rotate, and the torque transmitted by the load further restricts the rotation of the motor.

The transmission principle block diagram is depicted in Fig 3. The stiffness torque K*(θ-θ) transmitted by the stiffness coefficient K applies a negative torque to the motor and a positive torque to the load based on the position difference. The transmission of a damping torque is based on the velocity difference, and the rest of the transmission principle is the same as that for the stiffness torque.

Fig 3

Servo system two-mass model transfer function block diagram.

According to the transfer function block diagram of the two-mass model of the servo system shown in Fig 3, the following matrix equation can be obtained: where T is the electromagnetic torque output by the motor, J is the moment of inertia of the motor, J is the moment of inertia of the load, K is the stiffness coefficient of the system, K is the damping coefficient of the system, θ is the angular position of the motor output, and θ is the angular position of the load output.

According to Eq (1), it is easy to deduce the transfer function between the electromagnetic torque T of the motor and the output angular velocity V of the motor as shown in Eq (2):

2.2 Mechanical resonance frequency suppression

In order to facilitate the introduction of the mechanical resonance frequency, the frequency characteristic curve containing the mechanical resonance frequency is shown in Fig 4. The gain in the amplitude–frequency characteristic of the valley frequency is called the anti-resonance frequency (f). The relationship between f, the stiffness coefficient, and the moment of inertia of the load is shown in Eq (3). According to Eq (3), the anti-resonant frequency is the natural oscillation frequency of the load and the spring, which is independent of the motor, but the system is difficult to operate under the anti-resonant frequency condition as all the energy input into the motor will flow quickly to the load. Although the motor is near to a calm state at the anti-resonant frequency, the load may be oscillating at a higher intensity at that time.

Fig 4

Frequency characteristic curve including the mechanical resonance frequency.

The peak frequency of the gain at the amplitude–frequency characteristic is called the resonant frequency (f). The expression of the resonant frequency is shown in Eq (4). At this resonant frequency, the motor and the load offer almost no obstacle to the movement. As the total inertia is small, the loop gain becomes large.

As the motor and the load are in free oscillation at the anti-resonance frequency, the focus is usually on the anti-resonance frequency rather than on the mechanical resonance frequency. The mechanical resonance frequency will lead to extreme instability of the motor and the load, imposing much greater difficulty on the controller design. In this study, the occurrence of the mechanical resonance frequency is suppressed by taking the following three measures:

Increasing the ratio of the motor’s rotational inertia to the load’s rotational inertia.

Increasing the ratio between the motor’s rotational inertia and the load’s rotational inertia is the most reliable method for improving the resonance problem. This is because a larger ratio of the motor’s rotational inertia to the load’s rotational inertia will result in less flexure of the system, and the prevalent resonance problems are caused by the flexure between the motor and the load. An increase in the ratio between the motor’s rotational inertia and the load’s rotational inertia can be achieved by either increasing the rotational inertia of the motor or decreasing the rotational inertia of the load. Increasing the motor’s rotational inertia will also increase the total inertia of the system. Assuming that the motor’s rotational inertia is increased by 25%, the available acceleration will be reduced by 25%; consequently, the torque would need to be increased by 25% to ensure the same acceleration as before. Increasing the motor’s rotational inertia will increase the cost of both the driver and the motor, which is contrary to the development goals of commercial satellites of high precision, low cost, low weight, and low power consumption. Therefore, reducing the rotational inertia of the load is the best way to increase the ratio of the motor’s rotational inertia to the load’s rotational inertia. As rotational inertia is a physical quantity that reflects the mass distribution of a rigid body, the rotational inertia of a load can be reduced by either reducing the mass of the load, or changing the dimensions of the load, or both. The magnitudes of the motor’s rotational inertia and the load’s rotational inertia of the SGCMG gimbal servo system designed in this study are shown in Table 1.

Table 1

Rotational inertia of the SGCMG gimbal servo system.

Parameter	JM	JL
Value	44340000 g·mm2	10058854.28 g·mm2

2. Increasing the stiffness of the system.

Increasing the stiffness of the SGCMG rotating mechanism is another important measure to prevent the occurrence of mechanical resonance. The SGCMG rotating mechanism is shown in Fig 5. It is designed to increase the stiffness by including the following features:

Fig 5

Schematic diagram of the SGCMG rotating mechanism.

The direct drive of the motor is used to avoid the introduction of additional flexible links.

A compact structure design is adopted to shorten the length of the transmission shaft.

Subject to the weight and size constraints, the diameter of the transmission shaft is increased as much as possible.

The rotating mechanism structures are designed using materials with high specific stiffness.

An axial preload is applied on the bearing to improve the supporting stiffness, rotation accuracy, and stability of the bearing unit.

Following the modal analysis of the entire machine construction of the SGCMG system, the frequency nephogram of each order is shown in Fig 6, and the natural frequency of each order is shown in Table 2. According to the modal analysis, the system possesses sufficient stiffness and avoids being detrimentally affected by the flywheel rotation frequency or the natural frequency of the gimbal structure; effectively, the design avoids the occurrence of resonance problems.

Fig 6

Frequency nephogram of each order.

Table 2

Natural frequency of each order.

Order	Frequency (Hz)
First-order	209.6
Second-order	226.21
Third-order	453.01
Fourth-order	461.13
Fifth-order	637.72
Sixth-order	662.03

3. Designing filters.

The resonance problem can be further attenuated in the SGCMG gimbal servo control system by designing filters, which are placed before the current loop and used to attenuate the gain variation caused by the flexibility between the motor and the load. The filters can be low-pass filters or hysteresis filters; adding filters can reduce the gain around the resonance frequency. The transfer function of the low-pass filter is shown in Eq (5). When the frequency of the low-pass filter is increased to the turning frequency ω, the attenuation becomes increasingly severe, and the low-pass filter causes a relatively large phase lag problem, which is equivalent to reducing the phase margin to some extent. The transfer function of the hysteresis filter is shown in Eq (6), where ω < ω, which has the maximum attenuation ω/ω at high frequencies. As the hysteresis filter has the advantage of a relatively small phase lag, a hysteresis filter is added before the current loop in this study.

3. Modeling the SGCMG servo system identification method

The SGCMG servo system identification process treats the SGCMG servo system as a black box [26]. A given input signal generates a specific output signal. The particular black-box model for a system is identified according to the characteristics of the input and output signals. According to the relevant identification algorithm, the frequency characteristic curve of the input and output signals is obtained, and then the transfer function of the black-box model of the system is identified. The system identification process is shown in Fig 7.

Fig 7

System identification flowchart.

First, an input of the corresponding excitation signal is made to the system. The speed response signal calculated by the circular grating of the SGCMG servo system is recorded. Secondly, the appropriate identification algorithm is used to obtain the spectrum characteristic response curve of the system. Then, the transfer function of the system is obtained by curve fitting.

The input excitation signal can take the form of either a white noise signal or a sinusoidal frequency sweep signal. Because the ideal white noise signal is difficult to generate, it will lead to poor repeatability of the identification experiment if the white noise signal is chosen as the excitation signal. In addition, because the frequency components contained in white noise cannot be controlled, there will be no excitation signal at certain frequencies, resulting in large errors in the calculation. Therefore, we select a sinusoidal frequency sweep signal as the input excitation signal in this study. The input signal frequency should be selected to ensure that the lowest frequency is at most 50% of the first corner frequency of the servo system.It is simultaneously necessary to ensure that the amplitude of the input current signal is reasonable. If the input current amplitude is too large, it will lead to an inappropriately high speed of the SGCMG servo system, which may cause permanent damage to the SGCMG. If the input current amplitude is too small, it would appear (because of the servo system) to overcome the friction of the shaft system, resulting in the entire system having no output response. As for the sampling frequency of the input and output signals, according to the Shannon sampling theorem, the sampling frequency of the signal should not be less than twice the highest resonant frequency of each link of the servo system.

The servo system frequency characteristics can be classified according to the type of input signal: deterministic (which can be described by an analytic expression), random, or pseudo-random (which can be described by an analytic expression while approaching a random nature). The signal processing method can be either a frequency response test, a Fourier analysis, or a correlation analysis. The frequency response test requires the evaluation of multiple frequency points to obtain the frequency response characteristics of the system; the accuracy of the result obtained with this test is relatively high, and the use of orthogonal correlation analysis is highly effective in this test. Fourier analysis can identify the frequency response of a linear system based on impulse and step responses; although this method is computationally less intensive, relatively simple to implement, and can be performed in a short time, it is only suitable for systems with a good signal-to-noise ratio. The correlation analysis method must be performed in the time domain, targeting linear systems that can be either continuous or discrete time signals, and it is equally applicable under poor signal-to-noise ratio conditions.

The input signal selected in this study is a deterministic signal, and the signal processing method is selected based on the frequency response test of the orthogonal correlation analysis method to obtain the amplitude–frequency and phase–frequency characteristic curves of the SGCMG servo system. We describe the following orthogonal correlation analysis for obtaining the frequency response processing method:

For a linear system, the autocorrelation function of the input signal is written as: where u(t) is the input response, T is the sampling period, and τ is the time delay.

The cross-correlation function can be written as follows: where y(t) is the output response.

Through a convolution integral, the two correlation functions (7) and (8) can be combined to obtain: where g(t’) is the impulse response. The frequency response of the system can be determined by the impulse response, and the Fourier transform of the impulse response can be performed to obtain the frequency response G(ω).

For sinusoidal input, the signal can be expressed as:

In this equation, u is the input signal amplitude, and ω is the input signal frequency.

According to Eqs (7) and (10), the autocorrelation function can be written as:

The expression of the output response y(t) can be given as: where φ(ω) represents the frequency response phase.

According to Eq (9), the cross-correlation function between the output response and the input signal is obtained as follows:

Eq (13) leads to the following:

The real and imaginary parts of the frequency response can be estimated from Eq (14).

When τ = 0, the real part of the frequency response can be expressed as: When , the imaginary part of the frequency response can be expressed as:

According to Eq (8), by multiplying and integrating the input signal and the output signal of the system, the cross-correlation function when τ = 0 can be calculated as: where n is the number of measurement cycles.

Similarly, the cross-correlation function when can be calculated as: where the phase shift of the velocity measurement signal is π/2 so that the sine signal becomes a cosine signal with the same frequency.

According to the following equation, the amplitude and phase relationships of the system’s frequency response can be determined as follows:

4. SGCMG servo system identification experiment

4.1 Principle

The system identification principle diagram is shown in Fig 8. The system identification experiment adopts the vector control strategy of i = 0. The computer provides an input signal to the SGCMG gimbal servo system. According to whether the current loop is a closed loop or not, there are two waysfor signal input:

Fig 8

System identification principle diagram.

If the current loop is open, a voltage signal is provided as an input to the system. The object of identification includes the motor part, the coupling part of the motor and the load, and the drive amplifier part.

If the current loop is closed, a current signal is provided as an input to the system. The object of identification includes the current loop controller, the motor part, the coupling part of the motor and the load, and the drive amplifier part.

As the current loop is the innermost loop of the SGCMG servo control system, the current loop is part of the open-loop model of the speed loop. Therefore, in this study, we choose the second scheme, which uses the closed loop of the current loop. Using the programmable characteristics of the controller, the current signal is input to the servo system at point I to drive the SGCMG gimbal servo system to rotate. Simultaneously, the frequency is increased from low to high. The system synchronously records the velocity response signal of the SGCMG servo system calculated by the circular grating at point O. The sampling can be stopped after the set sampling period is completed, and the data is stored in the computer.

4.2 Execution of servo system identification experiment

The experimental system includes the SGCMG system, a data acquisition circuit, a power supply, and a computer, as shown in Fig 9. The data acquisition circuit is used for transmitting the control signal and driving the motor rotation. The power supply provides a voltage of 28 V to the PMSM and a voltage of 5 V to the digital control chip. The computer is used to pass control instructions to the entire system and to write the control flow code. The software for data recording and processing is developed in-house. The current signal are provided as inputs to the SGCMG system, and the velocity signal, which is obtained through the circular grating differential, is transmitted by the controller area network (CAN) to the computer for recording.

Fig 9

System identification experiment.

The main flow of the frequency characteristic identification experiment for the SGCMG gimbal servo system is specified as follows:

Time-domain test: a current signal with appropriate frequency and amplitude is selected as the input signal of the system to drive the SGCMG gimbal to overcome the friction torque rotation, and the angular velocity response signal is collected and stored in real time.

Frequency-domain estimation: the orthogonal correlation analysis method is used to perform the calculation of the input and output time-domain signals collected in the system identification experiment. The frequency characteristic curves containing the amplitude–frequency and the phase–frequency characteristics are obtained using the previouslydescribed method.

Fitting transfer function: the transfer function model of the SGCMG gimbal servo system is obtained by the curve fitting method, in accordance with the frequency characteristic curve, and the frequency characteristic curve of the fitting transfer function is plotted and compared with that obtained by experiment.

An appropriate input signal needs to be designed as an excitation signal to ensure the experiment is appropriate and rigorous. Some scholars take a white noise [27] signal with an appropriate frequency band as the input excitation which can give a balanced excitation of each frequency point of the tested system and allows the dynamic frequency characteristics of the tested system to be obtained. However, a white noise signal with random characteristics does not exist in reality, and number of the frequency bands of white noise that can be realized physically are limited, so the repeatability of using white noise signal as input excitation signal for system identification is relatively poor. The sinusoidal frequency sweep method is a fast and effective measurement method of frequency characteristics. Compared with the white noise method, the biggest advantage is that the repeatability is better. In this paper, a sinusoidal frequency sweep signal is used as the input excitation signal of the system to test the frequency characteristics of SGCMG gimbal servo system. The sinusoidal frequency sweep signal is programmed in the controller according to Formulas (21)–(23).

Where I(t) is the sinusoidal frequency sweep signal, A is the amplitude, f and f are the starting point and end point of sinusoidal frequency sweep signal respectively, T is the period required for the servo system to complete the frequency sweep, and n is the polynomial order of the transfer function to identify the controlled object.

The parameters of the sinusoidal frequency sweep signal used in this paper are shown in Table 3 below.

Table 3

The parameters of the sinusoidal frequency sweep signal.

Parameter	A	f 0	f T	T	n
Value	0.08 A	1.2 HZ	100 HZ	30 s	3

5. Results and discussion

In order to verify the correctness of the experimental platform, a single frequency sinusoidal current signal is used as the input to the SGCMG gimbal servo system. The frequency and amplitude of the input signal, respectively, are 2 Hz and 0.08 A. The frequency and amplitude of the output velocity response signal, respectively, are 2 Hz and 0.4573 rad/s. The input and output curves are shown in Fig 10.

Fig 10

Single frequency input and output curves.

Because the above experiment only uses one frequency signal as the input, only one amplitude–frequency characteristic and one phase–frequency characteristic can be obtained. The frequency characteristic curve of the system obtained by repeatedly collecting a single frequency point in this way is not only time-consuming and labor-intensive, but the result obtained is also not accurate. Therefore, the input current is changed to a sinusoidal frequency sweep current signal. The frequency and amplitude of the input signal, respectively, are 1.2–100 Hz and 0.08 A. The input sine sweep signal and the output speed response signal have the same sampling frequency of 200 Hz. The resulting curves of the input sine sweep current signal, and the output velocity response, are shown in Fig 11.

Fig 11

Sinusoidal frequency sweep signal input and output curves.

Using the method of orthogonal correlation analyze to analyse the input sine sweep current signal and the output response signal, the frequency characteristic curve of the SGCMG servo system can be obtained, as shown in Fig 12. In addition to the orthogonal correlation analysis method used in this paper, the identification algorithm often used by other scholars is the Fourier transform [28-33]. In order to prove the superiority of the identification algorithm proposed in this paper, the frequency characteristic curve of the SGCMG gimbal servo system was obtained by the Fourier transform at the same time. The comparison of the frequency characteristic curve generated by the two methods is shown in Fig 13. Evidently, the noise on the frequency response characteristic curve calculated by the orthogonal correlation analysis method is smaller. As the computer is limited by measurement time, the transform can only be carried out in a limited range; this reduces the accuracy and produces the noise on the system output signal. The orthogonal correlation analysis method multiplies the input excitation signal and output response of the system to calculate the average value, so it can effectively suppress the interference harmonics and noise.

Fig 12

The frequency characteristic curve of the SGCMG servo system.

Fig 13

Comparison of frequency characteristic curves obtained by orthogonal correlation analysis and Fourier transform method.

Fig 12 indicates that the SGCMG gimbal servo system designed in this study can effectively prevent mechanical resonance by increasing the ratio of the motor rotational inertia to the load rotational inertia, increasing the stiffness of the system, and designing a filter. The system transfer function G(s) can be obtained by fitting the frequency characteristic curve of the SGCMG servo system, as shown in Eq (24). The result shown in Fig 13 displays the comparison between the frequency characteristics measured in the experiments and those calculated by the theoretical fitting. Fig 14 indicates that the amplitude–frequency and phase–frequency characteristics of the transfer function identified agree well with the measured curve. Fig 15 is the error curve between the measured frequency characteristic and the frequency characteristic of the identified transfer function. It is evident from this curve that the absolute errors in the amplitude–frequency and phase–frequency characteristics are less than 0.5 dB and 1°, respectively, in the common used frequency band.

Fig 14

Frequency characteristic curves of acquisition and identification.

Fig 15

Identification errors.

6. Conclusions

This paper mainly focuses on the study of mechanical resonance and transfer function using the system identification method. First, the SGCMG gimbal servo system was modelled as the equivalent of a two-mass block model, and the transfer function between the electromagnetic torque and the output angular velocity of the gimbal servo system motor was studied. Second, the typical mechanical resonance frequency in the general servo system was considered. By starting from the underlying cause of the mechanical resonance frequency, methods of preventing the mechanical resonance were proposed. Finally, an input sinusoidal current sweep signal was used to test the frequency characteristics of the system. The input sinusoidal current sweep signal and the output velocity response signal were analyzed and processed by the orthogonal correlation analysis method. The frequency characteristic curves of the system, including the amplitude–frequency and the phase–frequency characteristics, were obtained. According to the frequency characteristic curve, it can be observed that the SGCMG designed in this study has no mechanical resonance problem. The novelty of this research is in proposing the application of orthogonal correlation analysis to servo system identification and comparing the identification results with the widely used Fourier transform method. The results show that the frequency characteristic curve identified by the orthogonal correlation analysis method has less noise singal than that obtained by Fourier transform. In addition, this paper also proposes the following options to prevent mechanical resonance: increase the ratio of the motor rotational inertia to the load rotational inertia, increase the stiffness of the system, and design a filter; these steps can improve the stability of the servo system.

(RAR)

Click here for additional data file.

21 Jan 2022

A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.

If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: https://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols.

We look forward to receiving your revised manuscript.

Kind regards,

António M. Lopes, PhD

Academic Editor

PLOS ONE

Journal requirements:

When submitting your revision, we need you to address these additional requirements.

1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and

https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf”.

2. Thank you for stating the following financial disclosure:

“This work was supported in part by the Jilin Province Scientific and Technological Development Program, grant number 20180201111GX.”

Please state what role the funders took in the study.  If the funders had no role, please state: ""The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.""

If this statement is not correct you must amend it as needed.

Please include this amended Role of Funder statement in your cover letter; we will change the online submission form on your behalf.

3. In your Data Availability statement, you have not specified where the minimal data set underlying the results described in your manuscript can be found. PLOS defines a study's minimal data set as the underlying data used to reach the conclusions drawn in the manuscript and any additional data required to replicate the reported study findings in their entirety. All PLOS journals require that the minimal data set be made fully available. For more information about our data policy, please see http://journals.plos.org/plosone/s/data-availability.

Upon re-submitting your revised manuscript, please upload your study’s minimal underlying data set as either Supporting Information files or to a stable, public repository and include the relevant URLs, DOIs, or accession numbers within your revised cover letter. For a list of acceptable repositories, please see http://journals.plos.org/plosone/s/data-availability#loc-recommended-repositories. Any potentially identifying patient information must be fully anonymized.

Important: If there are ethical or legal restrictions to sharing your data publicly, please explain these restrictions in detail. Please see our guidelines for more information on what we consider unacceptable restrictions to publicly sharing data: http://journals.plos.org/plosone/s/data-availability#loc-unacceptable-data-access-restrictions. Note that it is not acceptable for the authors to be the sole named individuals responsible for ensuring data access.

We will update your Data Availability statement to reflect the information you provide in your cover letter.

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: No

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: This study focused on the establishment of a model for the SGCMG gimbal servo system, the suppression of mechanical resonance, and the identification of a transfer function. First, the SGCMG gimbal servo system was modelled as the equivalent of a two-mass block model, and the transfer function between the electromagnetic torque and the output angular velocity of the gimbal servo system motor was studied. Second, the typical mechanical resonance frequency in the general servo system was considered. I recommend the article to be accepted for publication in the PLOS ONE after providing following revision.

1. Further strengthen the motivation to write this article in the Abstract section;

2. Please describe in detail the identification and mechanical of problem in paper.

3. Further polish the English writing of the full article;

4. Rewrite the Abstract section, not to detail it, and to focus on the outlines and the task in the paper without going into details or writing letter

5. Please improve the quality of the drawn figures

6. Unify the writing format of all references base on the journal style.

7. More description on Examples should be given,

8. Please also check each reference, there are some typos.

9. The Conclusions section is not satisfactory; it should be further improved. Please emphasize the main novelty of this paper and the significance of the results in the conclusion.

10. Please check writing mathematical formulas in paper.

11. In general, the typeset equations should be regarded as parts of a sentence and treated accordingly with the appropriate grammatical convention and punctuation. More editing for writing is needed. At the end of all equations must be put ”COMMA” or ”POINT” according to the typing rules. Therefore, they need to pre-check all the equations.

12. It is advised to add the following refs related to the work

https://doi.org/10.1016/j.jare.2020.06.018

https://doi.org/10.1016/j.cnsns.2021.105755

https://doi.org/10.1016/j.jocs.2021.101394

Reviewer #2: 1) Considering the current state-of-the-art, what is the novelty of this work?

2) The relevance of the signal processing method (frequency analysis) would be better proven if you add a comparison with those obtained by Fourier and correlation analysis.

3) The experimental result is not clearly presented, with no links between the simulation results and their practical validation.

4) Line 359, “Fig 13 indicates that the fit of the transfer function is excellent”. It is not clear how the authors have concluded the relevance of the fit parameter.

5) The discussion part is poor, need more detail of comprising, you can use some of the related work and mention the references.

6) Line 379, “Therefore, the fitting effect of the transfer function is verified”. The fitting effect is not discussed in the paper. So, the fitting effect should be more detailed.

7) In the introduction part, lines 61, 62, 63, 70, 75, 76, and 79, the word “reference” should be rephrased.

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: Yes: LAHLOUH Ilyas

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.

29 Mar 2022

Response to Reviewer 1 Comments

Point 1: Further strengthen the motivation to write this article in the Abstract section;

Response 1: Thank you for pointing this out. We have added a sentence to further strengthen the motivation to write this article in the Abstract.

Point 2: Please describe in detail the identification and mechanical of problem in paper.

Response 2: Thank you for pointing this out. We have added a detailed description of the identification and mechanical of problem in lines 50-56 of the paper.

Point 3: Further polish the English writing of the full article;

Response 3: The manuscript has been edited to improve spelling, grammar and readability.

Point 4: Rewrite the Abstract section, not to detail it, and to focus on the outlines and the task in the paper without going into details or writing letter.

Response 4: Thank you for pointing this out. We have rewritten the abstract. The rewrite does not go into detail but focuses on the outlines and tasks to make the Abstract clearer.

Point 5: Please improve the quality of the drawn figures

Response 5: Thank you for pointing this out. We have improved the quality of the figures.

Point 6: Unify the writing format of all references base on the journal style

Response 6: Thank you for pointing this out. We have unified the writing format of all references base on the PLOS ONE journal style.

Point 7: More description on Examples should be given,

Response 7: We have added more description of the examples in lines 62 to72. In order to explain the principle and procedure of the experiment more clearly, we have added more detailed description in lines 298–302 and lines 332–349.

Point 8: Please also check each reference, there are some typos.

Response 8: Thank you for pointing this out. We have checked each reference and revised the typos.

Point 9: The Conclusions section is not satisfactory; it should be further improved. Please emphasize the main novelty of this paper and the significance of the results in the conclusion.

Response 9: Thank you for pointing this out. We have emphasized the main novelty of this paper and the significance of the results in the conclusion.

Point 10: Please check writing mathematical formulas in paper.

Response 10: Thank you for pointing this out. We have checked all the mathematical formulas and revised mistakes.

Point 11: In general, the typeset equations should be regarded as parts of a sentence and treated accordingly with the appropriate grammatical convention and punctuation. More editing for writing is needed. At the end of all equations must be put ”COMMA” or ”POINT” according to the typing rules. Therefore, they need to pre-check all the equations.

Response 11: Thank you for pointing this out. Following your comments, we have edited the full manuscript for this. We have checked each of the equations and added the missing comma or full stop as appropriate.

Point 12: It is advised to add the following refs related to the work

https://doi.org/10.1016/j.jare.2020.06.018

https://doi.org/10.1016/j.cnsns.2021.105755

https://doi.org/10.1016/j.jocs.2021.101394

Response 12: Thank you for this suggestion. We have added citations to the three references in the introduction.

Response to Reviewer 2 Comments

Point 1: Considering the current state-of-the-art, what is the novelty of this work?

Response 1: In this paper, orthogonal correlation analysis is applied to servo system identification for the first time. Other scholars previously only identified relevant parameters, making the identified model inaccurate, or use Fourier transform to identify the frequency characteristic curve of the system. The computer is limited by measurement time. Therefore, the Fourier transform can only be carried out in a limited range; this would reduce the accuracy and causes the noise to affect the system identified frequency characteristic curve. In order to verify the superiority of the orthogonal correlation analysis method proposed in this paper, the orthogonal correlation analysis method is compared with the Fourier analysis; it can be observed from the identification results that the noise on the frequency response characteristic curve identified by the orthogonal correlation analysis method is smaller. In addition, this paper also suggests the following to prevent mechanical resonance: increasing the ratio of the motor rotational inertia to the load rotational inertia, increasing the stiffness of the system, and designing a filter; these steps can improve the stability of the servo system. We have emphasized the main novelty of this paper and the significance of the results in the conclusion.

Point 2: The relevance of the signal processing method (frequency analysis) would be better proven if you add a comparison with those obtained by Fourier and correlation analysis.

Response 2: Thank you for pointing this out. We have added the identification result of Fourier analysis in lines 370–383, and compared with frequency analysis at the same time. The experimental results show that there is more noise in the frequency response characteristic curve identified by Fourier analysis. Therefore, the orthogonal correlation analysis method is used to identify the system.

Point 3: The experimental result is not clearly presented, with no links between the simulation results and their practical validation.

Response 3: We have remedied this by highlighting that the main purpose of this paper is to identify the frequency characteristic curve and transfer function of the system through experiment. In order to clearly explain the experimental result, we have added a detailed discussion of the experimental result in lines 391-398.

Point 4: Line 359, “Fig 13 indicates that the fit of the transfer function is excellent”. It is not clear how the authors have concluded the relevance of the fit parameter.

Response 4: Thank you for pointing this out. We have added the error curve between the measured frequency characteristic and the frequency characteristic of the identified transfer function (Fig 15). It can be seen from this that the absolute error of the amplitude-frequency and phase-frequency characteristics are less than 0.5 dB and 1 °, respectively, in the common frequency band.

Point 5: The discussion part is poor, need more detail of comprising, you can use some of the related work and mention the references.

Response 5: We have added more detailed discussion. In order to explain the principle and procedure of the experiment more clearly, we have added more detailed description in lines 298–302 and lines 332–349. In order to clearly explain the experimental result, we added detailed discussion of the experimental result in lines 370–383 and lines 391-398. We have also used the related work and mention the references in discussion.

Point 6: Line 379, “Therefore, the fitting effect of the transfer function is verified”. The fitting effect is not discussed in the paper. So, the fitting effect should be more detailed.

Response 6: Thank you for pointing this out. We have added more discussed about the fitting effect in lines 391–398.

Point 7: In the introduction part, lines 61, 62, 63, 70, 75, 76, and 79, the word “reference” should be rephrased.

Response 7: Thank you for pointing this out. We have checked the introduction and rephrased the word “reference” in the appropriate lines.

Submitted filename: Response to Reviewers.pdf

Click here for additional data file.

11 Apr 2022

System identification and mechanical resonance frequency suppression for servo control used in single gimbal control moment gyroscope

PONE-D-21-26743R1

Dear Dr. Yu,

We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.

Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.

An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

Kind regards,

António M. Lopes, PhD

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.

Reviewer #1: All comments have been addressed

Reviewer #2: All comments have been addressed

**********

2. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

3. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: No

Reviewer #2: Yes

**********

6. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: All comments were addressed. The revised paper can be accepted for publication in the PLOS ONE journal.

Reviewer #2: The new version is well prepared end organised and all the responses are very clear.

The novelty is now well detailed, but it would be better if you add it in the introduction part.

**********

7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

9 Aug 2022

PONE-D-21-26743R1

System identification and mechanical resonance frequency suppression for servo control used in single gimbal control moment gyroscope

Dear Dr. Yu:

I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

If we can help with anything else, please email us at plosone@plos.org.

Thank you for submitting your work to PLOS ONE and supporting open access.

Kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. António M. Lopes

Academic Editor

PLOS ONE