Research Article  Open Access
Xinsheng Wang, Danwei Wang, Senqiang Zhu, Eng Kee Poh, "Fractional Describing Function Analysis of PWPF Modulator", Mathematical Problems in Engineering, vol. 2013, Article ID 287040, 5 pages, 2013. https://doi.org/10.1155/2013/287040
Fractional Describing Function Analysis of PWPF Modulator
Abstract
Pulsewidth pulsefrequency (PWPF) modulators are widely used in spacecraft thruster control. Their dynamic characteristic is still lack of effective analysis tools. This paper presents a fractional describing function method to describe the frequency characteristics of PWPF. A frequencydependent gain and phase shift are clearly described by fractionalorder expression, and the fractionalorder behaviors depict the nonlinear properties of PWPF modulators. This fractional describing function method can also be applied to other kinds of modulators.
1. Introduction
Spacecrafts commonly deploy thrusters as actuators for attitude control [1]. As thrusters work on on/off switching, there is a conversion from continuous desired torque control command to an on/off signal for spacecraft thrusters. There are two common approaches for thruster activation. The simplest converter is a bangbang controller which is, however, vulnerable to noise. Several improved controllers have been proposed, such as bangbang controller with dead zone or timeoptimal bangbang controller [2]. Another kind of the converter is pulse modulators, which produce a pulse command sequence to the thruster valves by adjusting pulse width and/or pulse frequency according to the level of input [3]. Among the known pulse modulators, the pulsewidth pulsefrequency (PWPF) modulators are the most common and enjoy advantages over bangbang control systems [4]. But the inherent nonlinearity of PWPF has inhibited the dynamic analysis of attitude control system.
Describing function method is a wellknown analysis tool for a kind of nonlinear system with certain structure. If the output signal of a nonlinear device can be approximated by the fundamental harmonics, the fundamental harmonics can be used to define the frequency characteristics of the nonlinear element and is called a describing function [5]. The describing function method is primarily used to analyze the stability and limit cycle of nonlinear control systems and is inherently approximate. In the last decades, many new criteria were established to improve the accuracy. These new types of describing function methods include areamatching method, rootmeansquare method, and corrected RMS method [3].
Recently, fractional calculus has been increasingly applied to mechanical systems, electricity, and bioengineering [6, 7]. Fractional calculus studies derivatives and integrals of fractional order. It is shown that fractionalorder system works more accurately than integerorder system [8, 9]. With the development of computational tools, more and more fractionalorder models and controllers are studied [10, 11]. Noticeably, the describing function method is also considered in the framework of nonlinear fractionalorder systems [12, 13]. In this paper, fractional describing function will be introduced to describe the dynamics of PWPF actuators.
2. PulseWidth PulseFrequency Modulator
Pulse modulators are commonly used in thruster control of fuel valves. There are various kinds of pulse modulators, such as pulsewidth modulator, pulsefrequency modulator, pseudorate modulator, and pulsewidth pulsefrequency (PWPF) modulator [14, 15]. PWPF modulator is preferred for its operation has almost linear input/output relationship [16, 17]. A PWPF modulator mainly comprises two components: a firstorder lag filter and a Schmitt trigger inside a feedback loop, as shown in Figure 1. A Schmitt trigger is an onoff relay with a dead zone and hysteresis. It differs from bangbang controller in that there are two thresholds: one onvalue and one offvalue . These values define a hysteresis as . The output of the Schmitt trigger is compared with the reference signal, and the error is fed to the firstorder filter whose output is the input of the Schmitt trigger. PWPF modulator operates in a quasilinear mode by modulating the width of the output pulses and the distance between them simultaneously. And it can produce pulses in two directions: positive and negative pulses.
With a constant input , the PWPF modulator drives the thruster valve with an onoff pulse sequence having a nearly linear duty cycle. The time interval in which the modulator has a nonzero output is denoted , and the time interval with a zero output is denoted . Static characteristics variables of PWPF are collected in Table 1. These variables are considered in the modulator design. Paper [18] shows the relationship between the static characteristics of PWPF modulator and selection of its parameters.

The static analysis of PWPF modulator shows that it operates near linear to the constant input over a large range between the deadband and saturation level . The modulator’s operations are independent from the spacecraft’s parameters and allow easy parameter tuning, especially when there are different requirements through different phases of operation. In addition, this modulator has the superiority in fuel consumption and pointing accuracy in the presence of vibrations. However, attitude control systems usually operate on dynamic mode, and PWPF modulators will introduce phase lag to the attitude control systems, which can cause instability. Dynamic analysis is necessary for the attitude control system design. But as to this nonlinear device, effective tools are limited. This work will propose the use of fractionalorder describing function method and develops some useful techniques.
3. Derivation of Fractional Describing Function
Describing function method is an approximation method for analyzing nonlinear dynamics, because only the first harmonic of the output of a nonlinear element is considered. But here describing function method is beneficial for PWPF modulator analysis for the following reasons. Firstly, the firstorder filter which is in series with the nonlinear element serves as a lowpass filter in PWPF modulator. Secondly, the attitude control system is a highorder system, where high harmonics are attenuated substantially. In this section, the PWPF modulator is considered as a single device and its fractional describing function is developed.
3.1. Describing Function of Hysteresis
A nonlinear element to a sinusoidal input in general does not generate a sinusoidal output, but the output is periodic. If the nonlinearity is symmetrical with respect to the variation around zero, the output signal can be decomposed to the Fourier series: where and the frequency of the fundamental harmonic is the same frequency as the input with amplitude and phase shift :
With the assumption that the amplitude of fundamental harmonic is much larger than the amplitude of other harmonics, the describing function is defined as the complex ratio of the fundamental harmonic component of the output and the input, that is, where is the describing function of the nonlinear element, is the sinusoidal input amplitude, is the frequency of input sinusoid, is the amplitude of the fundamental harmonic, is the describing function gain, and is the describing function phase.
Imposing this harmonic balance principle, can be used as frequency characteristic of nonlinear element to analyze the dynamics of a linear closedloop system. The Schmitt trigger in PWPF modulator exhibits hysteresis nonlinearity, whose describing function of hysteresis can be deduced as shown above directly:
Remark 1. The describing function defined above is based on fundamental harmonic equivalence. The influence of high order harmonics can be considered as system uncertainty.
3.2. Frequency Characteristics of PWPF
The PWPF modulator is a unit in the attitude control system. Considering, the frequency characteristics of PWPF as a whole is convenient for analysis and design of the attitude control system. In the structure of Figure 1, the frequency characteristic of the firstorder filter is
Then the frequency characteristic of the PWPF closedloop system is
If the describing function is independent of frequency , it is plotted with varying nonlinear input . But when it is a function of both amplitude and frequency , certain values are selected to view the plot of frequency characteristics. Consequently, this kind of plot lacks frequency dependent information. For control system dynamic analysis, information about how the amplitude and phase of a nonlinear element change with frequency is important because the nonlinear part usually adds phase lag to control system, which would destroy closeloop stability of the control system.
Now, a fractional describing function is introduced to give a direct relationship between the nonlinear characteristic and frequency. For illustration purpose, set the parameters of PWPF as follows: , , , and . Figure 2 shows the Nichols plot of the closedloop describing function . It can be seen that with the amplitude changing from 0.5 to 10 the magnitude of decreases, and the phase lag increases. To show the characteristic with frequency, several curves corresponding with different frequencies are plotted. Though it is clear that higher frequency induces larger phase lag, but the curves are separated and cannot show continuous information about frequency.
To reveal the relationship between the describing function and the frequency of interest, the real part and imaginary part of are studied, respectively. For attitude control, the interested frequency is relatively low. Figures 3 and 4 show the loglog plots of the real part and imaginary part of versus the exciting frequency , respectively, where the nonlinear input is from 0.5 to 1. Approximately, the curves can be considered as straight lines. Then a fractionalorder behavior is investigated.
The fractionalorder behavior is described by power functions, so it is called fractional describing function. The new function can be written as
Figure 3 shows that the real part of can be considered as a constant, that is, and . Figure 4 shows that the imaginary part of is fractional with and . These values will vary with system parameters, but the fractional property remains the same. It can be seen that the range of parameters does not vary much. We can choose one set of values to simulate the frequency characteristic of the PWPF modulator. Here, we set the values . With the help of available fractional calculator software [19, 20], frequency characteristic diagrams can be plotted. Figure 5 shows the Bode diagram of . It predicts a phase lag of about 10 degrees to 25 degrees between frequency 1 rad/s to 5 rad/s, which relates well with the analysis result in [3]. This information is crucial to the stability analysis of attitude control system.
(a)
(b)
Remark 2. The deduced fractional describing function is useful for system stability analysis for it is considered around the crossover frequency.
4. Conclusions
PWPF modulator is a nonlinear actuator in spacecraft attitude control system. The nonlinear dynamic behavior of PWPF modulator is investigated by the fractional describing function in this paper. The nonlinear element of PWPF is a Schmitt trigger, and its frequency characteristic can be described by describing function, and the fractional behaviors are caused by nonlinear element in PWPF. The frequency characteristic of the actuator is frequency dependent. The loglog plots of the real part and imaginary part of a modulator clearly reveal the fractionalorder behavior. The imaginary component is described by fractionalorder power function over a certain frequency range. With fractional calculus, these frequencydependent gain and phase information can be plotted in Bode diagram and used for control system design. Furthermore, the fractional describing function method should be an effective tool for other kinds of modulators.
Acknowledgments
This work is funded in part by the China Scholarship Council (CSC) and in part under Project Agreement no. DSOCL10004 with DSO National Laboratories, Singapore.
References
 M. J. Sidi, Spacecraft Dynamics and Control: A Practical Engineering Approach, Cambridge Aerospace Series, Cambridge university press, New York, NY, USA, 1997.
 B. Wie and C. T. Plescia, “Attitude stabilization of flexible spacecraft during stationkeeping maneuvers,” Journal of Guidance, Control, and Dynamics, vol. 7, no. 4, pp. 430–436, 1984. View at: Google Scholar
 T. C. Anthony, B. Wie, and S. Carroll, “Pulsemodulated control synthesis for a flexible spacecraft,” Journal of Guidance, Control, and Dynamics, vol. 13, no. 6, pp. 1014–1022, 1990. View at: Google Scholar
 B. N. Agrawal, “Spacecraft vibration suppression using smart structure,” in Proceedings of the 4th International Congress on Sound and Vibration, pp. 563–570, St. Petersburg, Russia, 1996. View at: Google Scholar
 H. Khalil, Nonlinear System, Prentice Hall, Upper Saddle River, NJ, USA, 2002.
 R. L. Bagley and R. A. Calico, “Fractional order state equations for the control of viscoelastically damped structures,” Journal of Guidance, Control, and Dynamics, vol. 14, no. 2, pp. 304–311, 1991. View at: Google Scholar
 C. M. Ionescu, J. A. T. Machado, and R. De Keyser, “Modeling of the lung impedance using a fractionalorder ladder network with constant phase elements,” IEEE Transactions on Biomedical Circuits and Systems, vol. 5, no. 1, pp. 83–89, 2011. View at: Publisher Site  Google Scholar
 I. Podlubny, “Fractionalorder systems and $P{I}^{\lambda}{D}^{\mu}$controllers,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 44, no. 1, pp. 208–214, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 A. Oustaloup, J. Sabatier, and P. Lanusse, “From fractal robustness to the CRONE control,” Fractional Calculus & Applied Analysis, vol. 2, no. 1, pp. 1–30, 1999. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 C. A. Monje, Y. Chen, B. M. Vinagre, D. Xue, and V. Feliu, FractionalOrder Systems and Controls, Springer, New York, NY, USA, 2010.
 Y. Luo and Y. Chen, “Fractional order [proportional derivative] controller for a class of fractional order systems,” Automatica, vol. 45, no. 10, pp. 2446–2450, 2009. View at: Publisher Site  Google Scholar
 M. S. Tavazoei and M. Haeri, “Describing function based methods for predicting chaos in a class of fractional order differential equations,” Nonlinear Dynamics, vol. 57, no. 3, pp. 363–373, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. B. Duarte and J. T. MacHado, “Fractional describing function of systems with Coulomb friction,” Nonlinear Dynamics, vol. 56, no. 4, pp. 381–387, 2009. View at: Publisher Site  Google Scholar
 B. Wie, Space Vehicle Dynamics and Control, AIAA Educational Series, Washington, DC, USA, 1998.
 Q. Hu and G. Ma, “Variable structure control and active vibration suppression of flexible spacecraft during attitude maneuver,” Aerospace Science and Technology, vol. 9, no. 4, pp. 307–317, 2005. View at: Publisher Site  Google Scholar
 G. Song, N. V. Buck, and B. N. Agrawal, “Spacecraft vibration reduction using pulsewidth pulsefrequency modulated input shaper,” Journal of Guidance, Control, and Dynamics, vol. 22, no. 3, pp. 433–440, 1999. View at: Google Scholar
 K. H. Kienitz and J. Bals, “Pulse modulation for attitude control with thrusters subject to switching restrictions,” Aerospace Science and Technology, vol. 9, no. 7, pp. 635–640, 2005. View at: Publisher Site  Google Scholar
 T. D. Krovel, Optimal tuning of PWPF modulator for attitude control [M.S. thesis], Norwegian University of Science and Technology, 2005.
 Y. Q. Chen, I. Petráš, and D. Xue, “Fractional order control—a tutorial,” in Proceedings of the American Control Conference (ACC '09), pp. 1397–1411, St. Louis, Mo, USA, June 2009. View at: Publisher Site  Google Scholar
 I. Petráš, “Fractionalorder memristorbased Chua's circuit,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 57, no. 12, pp. 975–979, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Xinsheng Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.