IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 11, NO. 1, JANUARY 2003
3
Direct Versus Indirect Line of Sight (LOS) Stabilization Peter J. Kennedy, Member, IEEE, and Rhonda L. Kennedy, Member, IEEE
Abstract—Two methods are analyzed for inertially stabilizing the pointing vector defining the line of sight (LOS) of a two-axis gimbaled laser tracker. Mounting the angular rate and acceleration sensors directly on the LOS axes is often used for precision pointing applications. This configuration impacts gimbal size, and the sensors must be capable of withstanding high angular slew rates. With the other stabilization method, sensors are mounted on the gimbal base, which alleviates some issues with the direct approach but may be less efficient, since disturbances are not measured in the LOS coordinate frame. This paper investigates the impact of LOS disturbances and sensor noise on the performance of each stabilization control loop configuration. It provides a detailed analysis of the mechanisms by which disturbances are coupled to the LOS track vector for each approach, and describes the advantages and disadvantages of each. It concludes with a performance comparison based upon simulated sensor noise and three sets of platform disturbance inputs ranging from mild to harsh disturbance environments. Index Terms—Angular rate sensor, gimbal, rate gyro, stabilization.
I. INTRODUCTION
T
HE REQUIREMENT to point a laser from a dynamic platform continues to surface in both defense and communication applications. The laser and track sensor (i.e., camera, detector array, etc.) are mounted on the inner axis of a multiaxis mechanical gimbal; or below the gimbal, coupled to the line of sight (LOS) via a stabilized pointing mirror. Pointing control is implemented via two servo loops, the outer track or pointing loop and an inner stabilization or rate loop. The track sensor detects the laser returns from the target location. The track processor uses this information to generate rate commands that direct the gimbal bore-sight toward the target LOS. The stabilization loop isolates the laser and sensor from platform motion and disturbances that would otherwise perturb the aim-point. The track loop must have sufficient bandwidth to track the LOS kinematics. The stabilization loop bandwidth must be high enough to reject the platform disturbance spectrum. A typical configuration is illustrated in Fig. 1. Disturbances that affect the pointing vector arise from platform angular motion or maneuvers, and external loads such as wind and air-stream induced torque. Platform linear motion and vibration generate disturbance torques due to mass Manuscript received November 30, 2000; revised March 15, 2002. Manuscript received in final form October 1, 2002. Recommended by Associate Editor D. W. Repperger. The authors are with David H. Pollock Consultants, Inc., Sensor Technology Specialists, Westwood, NJ 07675 USA (e-mail: dpollock@dhpconsultants. com). Digital Object Identifier 10.1109/TCST.2002.806443
imbalance and gimbal geometry. Conventional stabilization techniques employ rate gyros, rate integrating gyros, or rate sensors to sense rate disturbances about the LOS. They are often mounted on the inner gimbal axis, as close as possible to the gimbal bore-sight (i.e., in Fig. 1), or in some application at the gimbal base. The mounting configuration can have a significant impact on pointing performance. This paper focuses on two methods of implementing the sta), bilization servo loop design. A two-axis gimbal ( over as shown in Fig. 1, is considered. Direct LOS stabilization, with the angular rate sensors mounted on the LOS axes, is normally recommended for precision pointing applications [1], [2]. The technique is simple to implement, requiring only two rate sensors and the associated rate loop compensation. Gimbal size is impacted, since a larger payload volume is required to mount the sensors on the inner gimbal. For high-performance aircraft applications, the penalty in drag and turbulence induced by a large gimbal penetrating the air-stream may be severe. In addition, the sensors must be capable of functioning, or at least withstanding, high angular rates generated during slew. With indirect LOS stabilization, currently being pursued for several applications, the rate sensors are mounted on the gimbal base. This can alleviate some of the problems associated with the direct stabilization approach. The disturbances, however, are not measured in the LOS coordinate frame, which can degrade performance. The base mounted sensors must accurately measure the rates at the gimbal base. A processor transforms them to an equivalent disturbance about the LOS. A rate loop feedforward compensator applies a control signal proportional to the estimated rates, attempting to cancel the disturbances. The relationship between the actual LOS disturbances and those measured at the base of the gimbal is dependent upon gimbal geometry, structural rigidity, resolver accuracy, tachometer and/or encoder accuracy, and the processor sampling rate (if the control signal is processor generated). This method is more complex than the direct approach.
II. GIMBAL KINEMATIC RELATIONSHIPS Three coordinate frames are used to define gimbal motion, starting with the gimbal base frame ( ), followed by the gimbal outer frame ( ), and finally the gimbal inner frame ( ). LOS elevation control is exercised via the inner gimbal -axis, while azimuth control is applied via the outer gimbal -axis. It is important to note that stabilization is required about the inner gimbal -axis, called cross elevation. Therefore, the cross-elevation axis is actually controlled indirectly. The coordinate trans-
1063-6536/03$17.00 © 2003 IEEE
4
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 11, NO. 1, JANUARY 2003
Fig. 1. Two-axis tracker configuration.
formation between the base and outer coordinate frame is (note a “bar” denotes a matrix, and a “hat” denotes a vector)
The relationship between the angular acceleration of the inner and outer gimbal coordinate frames is obtained by differentiating (2) or (3) where
—relative azimuth angle between the outer and base coordinate frames. The base coordinate frame rate vector by outer coordinate frame rate vector
is related to the
(1) where
III. GIMBAL DYNAMICS MODEL The gimbal dynamics model can be derived from the torque relationships about the inner and outer gimbal body axes based on rigid body dynamics. The expression for the inner gimbal angular dynamics is given by (“ ” denotes the cross product operation) (4)
The coordinate transformation between the outer and inner gimbal coordinate frame is given by
—relative elevation angle between inner and outer gimbal coordinate frames.
where sum of the kinematic torques about the inner gimbal, , is
with
The outer coordinate frame rate vector is related to the inner , by coordinate frame rate vector, (2)
Friction and cable restraint torques. Inner gimbal inertia matrix. Reaction torques exerted by inner gimbal on outer gimbal. Mass imbalance torques about each gimbal axis. Elevation stabilization control torque.
KENNEDY AND KENNEDY: DIRECT VERSUS INDIRECT LINE OF SIGHT (LOS) STABILIZATION
The friction and cable restraint torques generate both linear and nonlinear disturbances, [3]. To simplify the analysis, it is assumed that the gimbal rotation axes are aligned with the principal axes of inertia so that the inertia matrix is diagonal and the inner gimbal dynamics equation can then be expanded as
5
As with the inner gimbal, alignment is assumed between the principal inertia axes and the axes of rotation, so that the outer inertia matrix is also diagonal. To express the dynamics in terms of the inner body -axis (cross elevation) one can solve (3) for and substitute the resulting expression and (10) into (9) to obtain
(5) (11) (6) (7)
The cross elevation axis dynamics are generated by the third element of vector (11). Expanding the equation for this component results in
where viscous friction coefficient; cable restraint coefficient; ; ; nonlinear friction torques; nonlinear cable restraint torques. In (5)–(7), and those equations that follow, the explicit dependency on time, , has been suppressed in the expanded expressions where no confusion will arise. The reaction torques on the outer gimbal are determined from (5) and (7). The elevation control axis is the inner gimbal -axis, whose dynamics are generated by (6). The inputs to this equation are the base rates and acand . Subcelerations, while the controlled variables are stituting (1) into (2) and then into (6), the elevation axis dynamics can then be expressed in terms of the base inputs as
(12) denotes the third element of the vector. Also the last where term on the left-hand side of (11) is zero. From (3), the angular acceleration about the inner -axis can be obtained. Substituting into (12), expanding the cross product terms and substituting the kinematic torques leads to
(13) where
(8) In a similar manner, the rigid body torque dynamics for the outer gimbal body are derived from (9) where sum of the kinematic torques about the outer gimbal is
outer gimbal viscous friction coefficient; outer gimbal cable restraint coefficient; ; ; ; nonlinear friction torque; nonlinear cable restraint torque. The cross elevation axis dynamics are described by (13). As with the elevation axis, this equation needs to be expressed in , , the base disturbances terms of the controlled variables , and the elevation and azimuth angles and . The modified expression is derived in Appendix A as
with outer gimbal reaction torques on base; outer gimbal mass unbalance torques; azimuth axis stabilization control torque; friction and cable restraint torque;
(10)
(14) The impact of the elevation angle on the cross elevation dynamics is evident from (A2)–(A4). Control of this axis diminishes as the elevation angle increases. A secant gain function can compensate to some degree, however, in practice the gain gets very high for angles beyond 70 . The inability of the two axis system to track through the “cone of occlusion” [4] is commonly termed the NADIR problem.
6
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 11, NO. 1, JANUARY 2003
Fig. 2.
Control configuration for direct stabilization.
Fig. 3.
Control configuration for indirect stabilization.
IV. STABILIZATION LOOP SERVO CONTROL The elevation axis control torque and the cross elevation can be generated by various algorithms. Recontrol torque gardless of the method, the primary objective is to suppress the (or its derivative). This can be effect of any base motion, achieved by direct feedback (direct LOS stabilization) and/or disturbance reconstruction and cancellation (indirect LOS stabilization). Significant differences exist between the two approaches, as will be discussed. The control algorithms presented are common practical implementations, derived from classical control methods and bandwidth considerations. With direct LOS stabilization, a rate sensor measures the rate and about the LOS. These are used to gendisturbances , which attempts to erate a feedback compensating torque null the measured disturbance. The typical control configuration is shown functionally in Fig. 2 for either axis. With indirect stabilization, rate sensors mounted on the base of the gimbal
measure platform motion and are used to generate a compensating torque , which attempts disturbance cancellation by matching the measured base rates, transformed into the LOS coordinates, to the LOS disturbance rates. In this method, disturbances induced by relative motion are compensated for by a , as tachometer feedback loop, which generates a torque in the direct method. The typical control configuration is shown functionally in Fig. 3 for either axis. A. LOS Stabilization of the Elevation Axis The elevation axis dynamics were described by (8). With the direct approach a rate sensor mounted on the LOS axes mea. To shape the loop response a proportional plus intesures gral (PI) control feedback compensator is often used, or (15)
KENNEDY AND KENNEDY: DIRECT VERSUS INDIRECT LINE OF SIGHT (LOS) STABILIZATION
. (measured elevation LOS rate). Elevation rate sensor measurement error. The sensor measurement error includes sensor noise. The PI compensator can also be designed to operate on the rate error, which is the difference between the commanded rate from the track loop and the measured rate. In this case the compensator is in the feedforward path, however the results obtained are essentially the same. Substituting the control expression into (8), one obtains
(16) where (inner axis related disturbances); (base disturbances); (sensor noise). The compensator coefficients are often determined as
where elevation rate loop bandwidth; elevation rate loop lead frequency [ to ]. , (16) can then Using differential operator notation be expressed as
7
It can be seen that disturbance rejection is derived from the freand the quency dependent gain (or attenuation factor) . In a point design, and are chosen such that inertia the magnitude of the denominator is sufficient to attenuate a measured torque or acceleration disturbance at a specific frequency(s). The base disturbances are coupled via friction and a cable restraint, but are attenuated by the loop gain. Normally, the desired rate loop bandwidth is such that and This essentially results in a high gain feed back approach. It is somewhat robust to parameter variations in the sense that exact , , and is not required as long as the knowledge of above inequality is satisfied. Sensor noise, however, is only marginally affected by the gain, so that it will still have a significant impact. Sensor noise is coupled via the compensator and will be . Since attenuated (or amplified) roughly as disturbance attenuation is inversely proportional to the servo loop bandwidth, it is limited by stability and noise. If the required bandwidth is too high, an alternative algorithm is needed. Identifying the dominant disturbance sources, such as nonlinear friction torque [3], [5], is necessary. In a severe platform motion and environment, base sensors can be used to measure and construct a feedback control torque to compensate for the (often termed Strong Stabilizabase disturbance input, tion [4]). The basic servo algorithm can be modified to include , (or ), and the nonlinear terms that cancel . This requires the additional measurement of the term in . outer axis roll rate With the indirect stabilization method, a tachometer is used to measure the relative elevation angular rate . A tri-axis rate sensor is mounted at the gimbal base to measure . The control torques for this configuration are (19a) and (19b)
(17) where The elevation axis response is coupled to the cross elevation axis rate and the outer axis roll rate by the fourth term on the left side of (17). In addition, the base rate inputs are also coupled to the azimuth axis response by the relative azimuth angle . Equation (17) can be rewritten as
(18) where
and
measured relative elevation angular rate; elevation tachometer noise; measured body , rates; ;
(measured resolver angle); ; ; azimuth resolver noise. It is also likely that the command for the compensation , will be computer generated and contain samtorque, pling errors associated with the base rate sensor and resolver measurements. For a fast enough sampling rate, this error comThis ponent will be proportional to the sampling interval, analysis, ignores this source of error and treats the command
8
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 11, NO. 1, JANUARY 2003
Again, a simple proportional plus integral control compensator is used, or
torque as an analog signal. Expanding the expression for and assuming the resolver error is small results in
where
(22) where
The combined inner gimbal compensation torques are then obtained as
measured cross-elevation rate. cross-elevation rate sensor measurement error. As with the elevation axis, the sensor measurement error includes noise. Substituting (22) into (14), and restricting the analysis to small elevation angles, then
(20) Substituting (19) into (8), and using the same algorithm to choose the servo loop compensator coefficients as described in Section IV-A, one has
(23) where .
(21)
sensor noise. The PI coefficients are chosen using the same algorithm as and are again the loop bandwidth for the elevation axis ( and lead frequency), or and
Comparison of (18) with (21), indicates that the loop gain for both approaches will be equivalent. The major difference is in the noise terms. Whereas the direct approach contains only the elevation rate sensor noise, the indirect approach is perturbed by both tachometer noise as well as base rate sensor noise. The difference in the noise terms arises from a fundamental difference in the control technique. The direct approach suppresses the disturbances via gain attenuation. The disturbance is never totally applied to the rate loop since it is always working to null it. Its magnitude as seen by the loop, is dependent on the loop bandwidth. The indirect approach requires a cancellation of the so that it must actually be measured. To estimate perrate formance of the two approaches, one can compare the root mean of the direct LOS method square (rss) rate sensor error with an estimate of the indirect LOS method noise as given by the root sum square (rss) of the rms noise of the tachometer, resolver, and base rate sensors, or
It is important to note that a noise contribution from the azimuth resolver, dependent on the outer gimbal axis roll rate, now impacts the elevation axis performance. Indirect method performance could be improved by mounting a rate sensor on the outer directly. gimbal to measure
The cross-elevation response is then approximated as
(24) where
and . With these approximations, rate loop performance depends upon loop gain and axis inertia, as with the elevation axis. Again, coupling of this axis with the elevation axis is seen to be dependent upon the outer axis roll rate. The indirect stabilization control configuration again includes a tachometer to measure the relative azimuth angular rate, , and the tri-axis rate sensor, specified previously, mounted on the gimbal base. The control torques for this case are given by (25a)
B. LOS Stabilization of the Cross-Elevation Axis The cross elevation axis dynamics were described by (14). is measured With the direct approach the cross elevation rate directly and fed back to suppress disturbances about this axis.
(25b)
KENNEDY AND KENNEDY: DIRECT VERSUS INDIRECT LINE OF SIGHT (LOS) STABILIZATION
where
where
measured relative elevation axis angle; ; . Again, issues regarding the base torque being computer generated and having an associated sampling error will be ignored. can be expanded (again asThe expression for the rate suming the resolver errors are small) as
(28) With the small elevation angle assumption, the response can be approximated as
where
For small elevation angles, mated as can be expressed as
9
(29)
in (25b) can be approxiand the base control torque
and are similar to those defined The transfer functions for the elevation axis. The primary difference between the direct and indirect approaches is again the noise term. The rate , can be compared with an upper bound sensor noise, of the rss of the rms noise associated with the indirect approach, where the or is approximated as shown in the equation at rms of the bottom of the page. Again the dependence of the rms noise on the resolver angle noise contribution and the outer axis roll and pitch rates can be observed. As with the elevation axis, in the presence of bandwidth limitations, enhanced performance may be possible by and as well as canceling the terms associated with direct compensation of friction and cable restraint torques. C. General Stabilization Servo Loop Requirements
(26) The combined outer gimbal azimuth compensation torque can finally be generated as
The preceding stabilization loop analysis focused on specific control algorithms based on known design techniques. More insight into stabilization loop design requirements may be obtained by formulating the problem in a general structure. The models describing the gimbal dynamics for the elevation and cross elevation axes, given by (8) (EL) and (14) (Cross EL), can be expressed in a differential vector form as
(30)
(27) where Substituting (27) into (14), one obtains
.
where
10
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 11, NO. 1, JANUARY 2003
for the indirect LOS stabilization case. As suggested earlier to improve performance, a strong stabilization approach [4] might be utilized in which case the feedback and base torques could be modified as ; where is an estimate of the nonlinear feedback vector the estimated base disturbance input. With these conand trol inputs, the response becomes
and
where
and
The direct LOS stabilization feedback design produces a control torque given by (31) is given in the equation at the bottom of the page. where Substituting (31) into (30), one obtains (32)
The difference vectors are the residual errors resulting from the cancellation of the respective terms. Implementation of the base disturbance cancellation torque requires knowledge of the base rates, gimbal inertia, linear friction and cable restraint constants, and gimbal angles as well as a model of the nonlinear friction torque [5], [6]. Another stability measure for the feedback design can be obtained by examining the linear part of the response to determine matrix eigenvalues have negconditions under which the ative real parts. A first-order estimate of the feedback gain requirements can be obtained by treating the time varying terms as constants. The characteristic equation is then given of by
where
or
Equation (32) represents the gimbal response to a disturbance input vector for the direct LOS stabilization approach. The servo loop gains, which attenuate the input disturbance, are now in. A similar expression can be derived herent in the matrix
Using the Routh criterion to establish stability conditions, a quick estimate can be had by assuming that the bandwidths and lead frequencies of each loop are equal, or
With this assumption, it can be shown that the following stability condition results:
KENNEDY AND KENNEDY: DIRECT VERSUS INDIRECT LINE OF SIGHT (LOS) STABILIZATION
This lower bound on loop bandwidth indicates it must exceed a value dependent upon the inertia terms and the maximum ex(see Appected roll rate, assuming is small and the term pendix A) is bounded. Conditions under which the solution to the homogenous part of (32) approaches zero as time tends to infinity can also be derived directly from the state-space repcan be expressed as: resentation. The matrix where contains only the constant elements contains and eigenvalues with negative real parts and the two time varying roll elements (element 2,4 and element 4,2). Equation (32) can then be written as (33)
11
TABLE I
vector norm is exponentially decaying. To satisfy the condition [ for [ [ and all , it is necessary that such that to select
This vector equation has the solution or for (34) is the state transition matrix. The following where positive bounds are assumed:
The bound is a function of the initial state, due to the nonlinear feedback term. Although the practical significance of this result may be limited, it does provide some insight into the lower bound on the loop bandwidth and gain required for a reasonably stable response. V. SIMULATION
The nonlinear term is a product of the state variables. For small, it will satisfy the following condition [ is a max] imum value of the inertia terms associated with
The gimbal model derived in previous sections and defined by (8) and (14), was simulated in order to provide a quantitative comparison of the two stabilization methods. The gimbal parameters assumed were the following. Inner Axis Inertia: lb-in-s
lb-in-s
lb-in-s
Inner and Outer Axis Friction: (viscous) lb-in/rad/s (Coulomb) lb-in/rad/s
since
inner and outer axis cable restraint: lb-in/rad
Therefore, for and the following condition . This bound on the state will be met: vector norm is dependent on the initial state vector. Taking the norm of both sides of (34), one has
outer axis inertia: lb-in-s
lb-in-s
lb-in-s
A 30-Hz servo bandwidth was obtained on the inner (elevation) axis and an approximate 20 Hz bandwidth on the outer (cross elevation) axis using a PI compensator given by
(35) where
assuming that for
It is shown in Appendix B, that this bound can be expressed . If the roll rate is as: , then the state bounded such that
Three cases were simulated, representative of different platform types, with the following sinusoidal disturbance rates [i.e., ] as listed in Table I. The rate sensor noise was assumed to be normally distributed with a zero mean and a standard deviation of 0.005 rad/s. The resolver noise was also assumed zero mean normal with a standard deviation of 0.001 rad. Sampling of the resolver angles was not modeled so that this source of error is not included.
12
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 11, NO. 1, JANUARY 2003
Fig. 4. Direct LOS stabilization Case One.
Fig. 5. Indirect LOS stabilization Case One.
The simulation results for Case One, representing a severe disturbance environment (with a significant platform roll rate), are shown in Figs. 4–6. Fig. 4 is a plot of the LOS angular and rate disturbance about the inner axis for the direct approach. Fig. 5 is the same plot, but for the indirect approach. Comparing the rate disturbances for each approach, it is evident that the density of the variations observed with the indirect approach is higher than for the direct approach. The disturbance rate distribution is plotted in Fig. 6 to quantify this observation. It can be seen that for the indirect approach, the distribution has a much broader spectrum than for the direct approach. This supports the analysis that indicated more error sources were associated with the indirect approach. Similar results were obtained for the outer axis. Figs. 4 and 5 also show the LOS angular devi-
ation. With the indirect approach, the drift in the deviation due to the low frequency rate bias resulting from the composite rate disturbances, such as the integrated roll/resolver error disturbance product [see (28)], is much higher. The inner axis disturbance rate distribution plot for Case Two is shown in Fig. 7. The average amplitude of the rate disturbances for this case is lower than for Case One, however, the same trends can be observed. The disturbance distribution resulting from the indirect control method is again significantly greater than with the direct method. Again similar results were obtained for the outer axis. The final case represents very low disturbance conditions with a correspondingly low roll rate. The rate distribution plot for the inner axis is shown in Fig. 8 and again the trend is similar to the previous two cases.
KENNEDY AND KENNEDY: DIRECT VERSUS INDIRECT LINE OF SIGHT (LOS) STABILIZATION
Fig. 6.
Y -axis rate distribution Case One.
Fig. 7.
Y -axis rate distribution Case Two. VI. SUMMARY
A two-axis gimbal model was derived to evaluate the performance of two methods for stabilizing the LOS pointing vector. With the direct approach, the rate sensors are mounted on the LOS axes and directly measure rate disturbances about the elevation and cross-elevation axes. With the indirect approach, the rate sensors are mounted at the gimbal base, and compensation requires that these rates be transformed into equivalent disturbances about the LOS. In this approach, relative rate disturbances such as friction, and imbalance are compensated for by tachometer feedback about the torque motor. The analysis showed that more error sources were associated with the indirect approach. The simulation results verified this and indicated that the performance of the indirect approach would
13
not be as efficient as with the direct approach. In the indirect approach, several of the noise inputs are coupled to the platform roll disturbance. For both methods, it was shown that a lower bound on stabilization loop bandwidth was proportional to the maximum expected roll rate disturbance, coupled in via gimbal geometry. The analysis of the indirect method did not account for sampling rate errors or those associated with the gimbal structural rigidity. Essentially, all signals were assumed to be analog. Before choosing a particular method, one must tradeoff the advantage gained with the indirect approach from the payload volume reduction versus the penalty paid in LOS disturbance attenuation performance. Although, indirect LOS stabilization performance may be satisfactory, this approach will not attain the same performance level as the direct approach, given an equal design effort.
14
Fig. 8.
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 11, NO. 1, JANUARY 2003
Y -axis rate distribution Case Three. APPENDIX A
Equation (13), given in Section III, can be modified by exand in terms of the controlled pressing the variables variables and base disturbances. The following expressions for and are used to generate this representation:
Substituting each of these expressions into (13) results in the following representation for the cross elevation dynamics:
(A1) Using these equations, the first term on the right-hand side of (13) can be expressed as
where
the second term as
(A2) the third term as
(A3) and (A4)
and finally, the friction and cable restraint term as
APPENDIX B This derivation follows that presented using the Gronwell– Bellman Lemma [7, Th. 16]. The vector norm inequality was given by (35), or
(B1)
KENNEDY AND KENNEDY: DIRECT VERSUS INDIRECT LINE OF SIGHT (LOS) STABILIZATION
Multiplying by
and dividing by , one obtains
(B2) Dividing the left side of this expression by the right side and results in multiplying by
The left hand side can be expressed as the derivative of the natural log of the quantity in the denominator, or
Integrating and raising each side to the power results in
From (B2), one then has
15
REFERENCES [1] A. K. Rue, “Stabilization of precision electro-optical pointing and tracking systems,” IEEE Trans. Aerosp. Electron. Syst., vol. AES-5, pp. 805–819, Sept. 1969. [2] , “Precision stabilization systems,” IEEE Trans. Aerosp. Electron. Syst., vol. AES-10, pp. 34–42, Jan. 1974. [3] B. Friedland and D. Haessig, Jr., “On modeling and simulation of friction,” ASME J. Dyn. Syst., Meas., Contr., vol. 113, pp. 354–362, 1991. [4] W. J. Bigley and F. Schupan, “Wide-band base motion isolation control for a mobile platform,” in Proc. 1987 American Control Conference, Minneapolis, MN, June 10, 1987. [5] C. D. Walrath, “Adaptive bearing friction compensation based on recent knowledge of dynamic friction,” Automatica, vol. 20, pp. 717–727, 1984. [6] B. Li, D. Hullender, and M. DeRenzo, “Nonlinear induced disturbance rejection in inertial stabilization systems,” IEEE Trans. Contr. Syst. Technol., vol. 6, pp. 421–427, May 1998. [7] N. Balabanian and T. Bickart, Electrical Network Theory. New York: Wiley, 1969, pp. 727–767.
Peter J. Kennedy (M’76) received the B.S. degree in mechanical engineering from Lafayette College, Easton, PA, in 1971, the M.S. degree in industrial engineering from Lehigh University, Bethlehem, PA, in 1973, and the M.S. degree in electrical engineering from Fairleigh Dickenson University, Teaneck, NJ, in 1978. He has worked in the area of pointing, tracking, and stabilization for the U.S. Army Electronic Warfare Laboratory, Bell Laboratories, SciTec, Inc., Northrop Grumman Corporation, and is presently employed with David H. Pollock Consultants, Inc., Westwood, NJ. Mr. Kennedy is a licensed professional engineer and a member of the National Society of Professional Engineers, and the Instrumentation, Systems, and Automation Society (ISA).
or
ACKNOWLEDGMENT The authors would like to thank the reviewers whose critical analysis and suggestions brought insight into the final editing of this paper.
Rhonda L. Kennedy (M’02) received the B.S. degree in mathematics from Phillips University, Enid, OK, in 1979 and the M.S. degree in operations research from Texas Tech University, Lubbock, TX, in 1981. She was with Bell Laboratories from 1981 to 1983. Currently she is working in the area of statistical analysis, modeling and interface design for simulation of EO/IR sensor systems, and real-time simulations of end-to-end countermeasure systems. She is presently employed with David H. Pollock Consultants, Inc., Westwood, NJ.