Tattooing and body piercing guidance Toolkit - Public Health England
Full description
ocean engineering
The Software Configuration Index (SCI) functions as a master list for the configuration of items under configuration control for the Guidance and Control Software (GCS) project. The software Life C...
ocean engineeringFull description
Data acquisition, processing and interpretation in the field of Geophysics.Full description
Modem Navigation, Guidance, and Control Processing
b y Ching-Fang L i n
Modeling, Design, Analysis, Simulation, and Evaluation ( M D A S E ) Modern Navigation, Guidance, and Control Processing Advanced Control Systems Design Integrated, Adaptive, and Intelligent Navigation, Guidance, and Control Systems Design Digital Navigation, Guidance, and Control Systems Design
Modem Navigation, Guidance, and Control Processing
Ching-Fang Lin American Gh'C Corporation
gig
Prentice Hall, Englewood Cliffs, New Jersey 07632
L i b r a r y of C o n g r e s s Catalo~lng-ln-Publlcatlon D a t a
Lln. Chlng-Fang. M o d e r n n a v l g a t l o n , guldance. and c o n t r o l processing / by C h l n g -Fang Lln. p. CI. ( S e r l e s In a d v a n c e d n a v l g a t l o n , guldance. a n d c o n t ~ o l . a n d t h e l r a p p l l c a t l o n s : b k . 2) I n c l u d e s b l b l l o g r a p h t c a l r e f e r e n c e s and index. I S B N 0-13-596230-7 1. F l l g h t c o n t r o l . 2. G u i d e d ntsslles--Control s y s t e m s . I. T l t l e . 11. S e r l e s . TL589.4.L55 1991 90-43009 629.1--0C20 CIP
--
Editoriallproduction supervision and interior design: Brendan M . Stewart Cover design: Bruce Kenselaar Manufacturing buyers: Kelly Behr and Susan Brunke Acquisitions editors: Bernard Goodwin and Michael Hays
O 1991 by Prentice-Hall, Inc. A Simpn & Schuster Company Englewood Cliffs. New Jersey 07632
Thii book can be made available to businesses and organizations at a special discount when ordered in large quantities. For more information, please eontad: PrenticeHall, Inc., Special Sales and Markets, College Division, Englewood Cliffs, NJ 07632 All rights reserved. N o part of this book may be reproduced, in any form or by any means, without permission in writing from the publisher Printed in the United States of America 1 0 9 8 7 6 5 4 3 2 1
ISBN 0-13-59b230-7 Prenticc-Hall International (UK) Limited, Lotidon Prentice-Hall of Australia Pty. Limited, Sydney Prentice-Hall Canada Inc., Toronto Prenticc-Hall Hispanoamericana. S.A.. Mexico Prentice-Hall of India Private Limited, ,Vru, Drlhi Prentice-Ha11 ofJapan. Inc.. Tokyo Simon & Schuster Asia Pte. Ltd.. Singapore Editora Prenticc-Hall do Brasil. Ltda.. Rio de Jatreiro
To my family for their love, understanding, and support throtrghout.
2 Modeling-Design-Analysis-SimulationEvaluation (MDASE) of NGC
Processing 2.1 LinearlNonlinear Intercept NGC System, 13 2.1.1 LinearlNonlinear Intercept N G C Processing, 14
Contents
2.1.2 2.1.3 2.1.4
Modeling and Simulation, 18 Guidance System Classification, 21 Flight Control System (FCS) and FCS Sensing, 22 2.2 Target Signal Processing, 23 2.2.1 Targeting, 30 2.2.2 Kinematic/Relative Geometry, 34 2.2.3 Targeting Sensor Dynamics, 37 2.3 N G C System Design and Analysis, 45 2.3.1 Guidance Filtering/Processing, 45 2.3.2 N G C Stability and Pevfovmance Analysis, 52 2.4 Target Tracking State Modeling, 68 2.4.1 Target Noise and Target Maneuver Modeling, 68 2.4.2 Two-Dimensional Target Tracking State Modeling, 72 2.4.3 Three-Dimensional lntercept State Modeling, 74
3 Modem Multivariable Control Analysis Singular Value Analysis, 78 Sensitivity and Complementary Sensitivity Functions, 82 Design Requirements, 84 Structured Singular Value, 86 General Robustness Analysis, 89 Robustness of Real Perturbations, 90 3.6.1 State-Space Model for Additive Uncertainty, 91 3.6.2 Slate-Spare Sing~rlarI/allres, 92 3.6.3 Root Loctrs, 92 3.6.4 A Monte Carlo Comnputation, 94 3.6.5 A Monte Carlo Analysis.for SecondOrder Systetns, 95 3.6.6 Stability Margin Comttp~rtatiot~, 98 Monte Carlo Analysis, 100 Covariance Analysis, 102 Adjoint Method, 106 3.9.1 Adjoint P/~ilosophy,106 3.9.2 Applications, 108
Contents 3.10 Statistical Lzncarization, 11 1 and Tools for Statistical 3.10.1 l'c~clrtriq~rc~s L~trcartzation, 116 3.10.2 Statistical Litlearization with Adjoint hlrtlrod, 1 18 3.10.3 Applicatiot~s, 120 3.11 Qualitativc Comparison, 121 3.12 Other l'erformancc Analysts, 124
5 Inertial Navigation 5.1 Introduction, 176 5.2 Common Requirements for Inertial Navigators, 181 5.3 Navigation Computation and Error Modeling, 184 5.3.1 Coordinate Systems, 184 5.3.2 Position and Velocity Generation, 186 5.3.3 Data Processing, 188 5.3.4 IA'S Error Analysis and Modeling, 188 5.4 Gimballed Inertial Navigation System (INS), 192 5.4.1 Ginzbal Mechanism, 192 '5.4.2 Iitertial Setisors on the Stable PlatJorrn, 193 5.4.3 Plarfornt Aligrlriient Modes, 193 5.4.4 nhvigatioit Mode, 194 5.4.5 Systein Ititernal arid External Integaces, 195 5.4.6 Navigatioii Mechanizatiorz and Error Model, 195 5.5 Strapdown Inertial Navigation System (INS), 199 5.5.1 Generalized Strapdowti Mechai~izatiori, 199 r, 5.5.2 Strapdouwi Navixatiori ~ o r n ~ u t e 201 5.5.3 Strapdown Coniputer Requireineiits, 201 5.5.4 Strapdowrl Error Model and Attalysis, 204 5.5.5 Operation Flow Diagram, 208 Comparison and Analysis of Gimbal Versus 5.6 Strapdown, 209 5.6.1 Fcatitrc Coiiiparisoii, 209 5.6.2 h'ai~igatiori Errors Cori~parisoil,209 5.7 External Navigation Aids, 213 5.7.1 Aided liierti~l,Vavkatioir Meckatiizatior?, 213 5.7.2 Global Positioi~iiigSysterri ( G P S ) , 215 5.7.3 Tactical A i r hlavigatiorl ( T A C A N ) , 277 5.7.4 L o r ~ gRange Navigation ( L O R A N ) , 220
176
Contents
5.7.5
Terrain Contour Matching (TEHCO.M), 220 5.7.6 Doppler Radar, 225 5.7.7 Star Trackers, 225 5.7.8 Kaln~anFiltering, 226 5.7.9 Kaltnon Filtering Performance, 228 5.8 Integrated Inertial Navigation System (IINS). 229 5.8.1 I t ~ t e ~ ~ r a rSensittcq/Flight ed Control Reference Systrtn (ISFCRS), 231 Sensory Subsystetn 5.8.2 it~te~qrated ( I S S ) , 231 5.8.3 Itltecqrared Inertial Sensing Assembly ( I I S A ) , 232 5.8.4 Helicopter Integrated Inertial h'avigation Systern (HIINS), 241 5.8.5 integrated Missile Guidance Systerns, 245
Defense and Offense Systems, 392 6.7.1 Performance Paratneters, 392 6.7.2 Low-Altitude Air Defense Systems, 397 6.8 Future Guidance Processing, 404 6.8.1 Areas for Teclztiological Advartces in Signal Processitig, 405 6.8.2 Future Signal Processing for Missile Guidatlce, 407
State Estimation, 414 Target Tracking Filter Suiizmary, 416 The Wiener Filter, 420 Kalmatl Filter, 420 The Sitiiplified Kalnlari Filter, 420 Alpha-Beta-Canima Filter, 421 Modified Maxitn~rm-Likelihood Filter, 424 7.1.7 Two-Point Extrapolator, 425 7.1.8 Coriiparisori of Targct Trackirrg Filters, 425 7.2 Practical Navigation and Guidance Filtcr Design, 427 7.2.1 Guidance Tracking Filtcr, 427 7.2.2 Navigation and Guidance Filteringfor Position Estiriiatc, 430 7.2.3 Navigatiorr arid Cuidatice Filteritlgfor Position and Velocity Estitilate, 431 7.2.4 Navigation atrd Cuidarrcc Filterit~gfor Position, I/clociry, and Accrleratiot~ Estiriiatc, 434 7.2.5 Advatrccd Guidance Filter, 440 7.3 Radar Tracking, 442 7.4 Spacccraft Attitudc Estimation, 444 7.5 Advanccd Navigation System Design, 445 7.5.1 Global Positioning Systetir (GPS) Acc~rracyItnpro~~etttcrit, 445 7.q.2 Integrated CPSIIR'S, 457
Contents
8 Advanced Guidance System Design 8.1
8.2
8.3 8.4
8.5
8.6
Advanced Guidancc Laws, 468 8.1.1 O~itirttnlCrridatrce Law Survey, 469 8.1.2 Atrnlyticnl Solrrtiotr of Optirnal Filters and Optitnal G~ridatrceLaw, 474 Complenicntary/Kalman Filtered Proportional Navigation: Biased PNG and Complcmcntary I'NG, 481 8.2.1 Corttliitted Seeker-Czridatrce Filtering it1 a Cottrplettret~taryFilter, 483 8.2.2 Biased PAYC ( B P N G ) Algorithm, 486 8.2.3 Cort~plettret~tary PNG (CPNG) Alyoritl~trr, 492 8.2.4 Terrtr inal G~ridatlceSystem Analysis, 501 Other Terminal Guidance Laws, 510 Radome Error Calibration and Compensation, 517 8.4.1 Radorne Error Cotnpensntion, 51 7 8.4.2 Gtridowce Perforrnat~ceAnalysis with In-Fl!qht Radotne Error Calibratiotr, 520 8.4.3 Des[gtt Equations for P N C with Parasitic Feedback, 528 Command Versus Semiactive Homing Guidance System Design and Analysis, 544 8.5.1 Modeling, 544 8.5.2 Miss-Distance Analysis for Command Guidance, 550 8.5.3 Analysis of Optimal Command Ccridance Versus Optimal Semiactive Homirlg C~ridance,557 Analytical S o l ~ ~ t i oonf Optimal Trajectory Shaping for Combined Midcourse and Terminal Guidance, 562 8.6.1 Optitnal Trajectory Sllaping Guidance, 562 8.6.2 Problem Formulation, 563 8.6.3 Attalytic Optimal Gtridance Law, 566 8.6.4 Real-Time lmplemetztation and Pevformance, 578 8.6.5 Discussiotls, 583
Series Foreword The role played by modern navigation, guidance, and control (NGC) in the development and advancement of such areas as commercial and military aviation, to name only two, has continually expanded since its earliest inception in the 1950's. As this field began to grow and take on added importance, many books were written dealing mainly with the theoretical aspects of NGC, but most of these were confined to the earlier years of N G C development. Currently, although N G C system applications continue to take on an ever-increasing importance, the availability of reference books, especially textbooks suitable for graduate level and advanced undergraduate students as well as those who practice in the field, has not kept pace. This series emphasizing NGC systems and their applications is long overdue; in fact, it has been 30 years since such a series dealing with N G C systems has been written and available to the academic and ~rofessionalcommunities. Moreover. it is the first ever such series to thoroughly discuss the advanced control system design (modern multivariable control analysis; robust control; estimation; adaptive control; nonlinear control; intelligent control; etc.). It comes at a time when concern over issues such as the status of education and the decreasing number of trained, qualified professionals in this country is a t an all-time high. It is against such a background that the present series was conceived to assess state-of-the-art systems and control theories, and engineering applications of advanced N G C systems. Another purpose of the series is to develop future research agenda and at the same time encourage A
xuiii
Series Foreword
discussions in those areas that do not always find the systems and control community in complete agreement. The series provides a comprehensive coverage of the latest N G C technology as follows. The first book begins by introducing the various applications of N G C systems, after which it provides a thorough, fundamental treatment of what is considered the five most important stages in N G C system development: modeling, design, analysis, simulation, and evaluation (MDASE). The second book in the series takes up the subject of advanced estimation and guidance systems design, as well as N G C processing. The third book is concerned with the subject of advanced control system design, with particular emphasis placed on the topic of flight control system (FCS) design. The topics that constitute the fourth book include integrated, adaptive, and intelligent N G C systems design, while the fifth book is devoted completely to digital N G C systems design. Although most of the material in these five books is self-contained, there is a natural progression in the series as a whole toward more advanced topics. For example, much of the material in the second book actually serves as a prelude to the third, fourth, and fifth books. These books are the result of several years of experience gained on the part of the authorleditor both as a professor at the university level and as a practitioner in the field. It is believed that this series will provide invaluable insight and instruction to students, mainly at the graduate level but also to advanced urldcrgraduate students, as well as to those engineers who work directly or indirectly in the field of N G C system design and applications. In addition, this series is intended to provide both engineers and managers with the advanced N G C knowledge and concepts necessary to make correct decisions concerning the best N G C system design in a particular situation.
Preface It is very likely that few people who labor in any scientific discipline are unaware of the contributions of advanced navigation, guidance, and control (NGC) theory to aerospace-related programs. It is, however, equally unlikely that many are aware of the dramatic impact of this field on such diverse areas as medicine, industrial manufacturing, energy management, and chemical engineering. While its broad range of applications would at first appear to indicate that advanced N G C theory is enjoying an immense popularity in scientific and academic settings in general, this unfortunately turns out not to be the case. It is felt by many experts that many of those in the aerospace field in particular either are content to rest on the laurels surrounding the success of NGC theory developed in the 1950's and 1960's or have become so conservative in their design philosophies as to be unduly apprehensive about using advanced N G C theory. The latter appears to be especially true in the area of aviation. The author is quick to point out, however, that the NGC field itself is somewhat responsible for many of the misperceptions on the part of those w h o are not convinced of the usefulness of advanced N G C theory. More than a mere shadow of doubt has been cast on the usefulness of this theory as a result of its having taken
xx
Preface
a much too mathematically-oriented turn almost immediately after the theory was first applied in the solution of practical problems. It is the author's opinion that, while N G C theorv is built around a rather beautiful framework of mathematics. its primary emphasis'must nonetheless always bk placed on solving engineering problems of great practical importance. N G C technology has always been the focal point of aerospace engineering and automation research and development. A combination of theoretical concepts, the rapid evolution of computer and microelectronics technology, and the continued refinement of sensor and actuator technology has contributed to its advances. This book examines the role of modern N G C processing in the design of advanced N G C systems. This volume places major emphasis on the practical applications of advanced NGC systems, treating the subject more from an engineering than a mathematical perspective. Nevertheless, theoretical and mathematical concepts are introduced and adequately developed to make the book a self-sufficient source of instruction for readers. The intent of this book is to enable readers to achieve a level of competence that will permit their participation in the practical applications of modeling, design, analysis, simulation, and evaluation (MDASE) to advanced N G C systems. The book presents basic as well as advanced algorithms. A wide range of examples culled from various applications are provided to meet the needs of the different levels and types of readers, extending from issues requiring only a rudimentary knowledge to those involving avant-garde research. Morever, problems in the text span from those that concern only N G C to those that are interdisciplinary, and ultimately to those that encompass the entire systems and control field. An outline of the topics presented in this book is given in Fig. 1. Following the Introduction (Chap. I ) , the text is organized according to five principal topics: MDASE of N G C Processing (Chap. 2), Modern Multivariablc Control Analysis (Chap. 3), Design Algorithms for Advanced NGC Systems Design (Chap. 4). Fundamentals of N G C Processing (Chaps. 5 and 6), and Advanced N G C Systems Design (Chaps. 7 and 8). Each of these is in turn divided into a number of subtopics that discuss the relevant theories, algorithms, and computing tools related to their applications in advanced NGC systems. Referring to the outline organizing the five principal topics covered in this book, only those subtopics that are connected by solid lines are treated in this book. Those subtopics that are connected by broken lines are treated specifically in the books referenced under them. The numerous examples that are included in this book are supplemented by a liberal use of rcferenccs, thus making it easier for the reader to get access to a tremendous body of literature it] this field. As noted by one rcvicwer, there are n o comparable books currently on the market that prcscnt in a usable format such a complete collection of practical tools that are applicable to real world problems. Moreover, the organization of thc material coupled with comprchensive examples make this book well suited to self-teaching, bridging the gap betwccn thc theoretical and the practical.
Preface
Introduction
Modeling-Design-Analysis-Simulation-Evaluation (MDASE) of NGC Processing
-----
I
Hodern Multivariable Control Analysis
1
Performance
-----
rn Advanced
Modern Filtering and Estimation Techniques (Chap. 4)
Design for Advanced NGC Systems Design
Adaptive NGC I Systems: Detection I and Identification (Book 4)
Multivariable Control Design Techniques
1
Knowledge-Based and Neural-Network Approaches to Systems Control (Book 4)
- - I- - - - 1 --Fundamentals of NGC Processing
Guidance Processing
1
Flight Control Processing (Bogk 3)
--7
Navigation and Guidance Filtering Design (Chap. 7)
+
Advanced Guidance Systems Design
Design
+
Advanced Flight Control System Design (Book 3)
I
r
7
- - -C - - - 1
Integrated, Adaptive, and Intelligent NGC Systems Design (Book 4)
Digital NGC Systems Design (Book 5)
Figure 1. Modern Navigation, Guidance, and Control (NGC) Processing
Acknowledgments The number of people who have lent assistance in one way or another to this book is too great to allow me to thank each one individually. 1 am indebted to each person who has contributed to making this project a success, but I feel that I must acknowledge a few individuals by name for their support and effort during the countless hours spent on this project. My first acknowledgement goes out to my staff at American GNC Corporation (AGNC) who have assisted me as follows: Chun Yang and Jerry Juang have provided important technical input, while Jim Wright has spent many hours editing the manuscript. Special thanks are due to Kylie Hsu and Janet Young for their involvement in designing, planning, coordinating, and preparing the entire manuscript, including word processing and numerous original pieces of computer and manual artwork. I am also grateful to William R. Yueh for his technical advice. Several of my colleagues and friends have provided invaluable technical assistance. I specifically want to acknowledge William H. Gilbert of Martin Marietta Electronic Systems who, through his ardent support and encouragement of this project, made this impossible task appear less so at times. I am deeply touched by his unselfish giving of his time in the overall guidance and input to this project. I am
xxiv
Acknowledgments
am also indebted to J. Stanley Ausman of Litton Guidance and Control Systems Division who has greatly improved the quality of Chapter 5. Dr. Ausman is an authority on the subject of this chapter and has co-authored a book on Inertial Guidance (Wiley, 1962). I would also like to thank him for granting me permission to adopt his lecture notes in Tactical Aircraft Weapon Delivery Systems. I also want to express my gratitude to Gary Hewer of Naval Weapons Center, Hsi-Han Yeh of Wright-Patterson Flight Dynamics Lab., and Keqin Gu of Southern Illinois University for their painstaking effort in reviewing Chapter 3. T. Sen Lee of MIT Lincoln Lab., is also appreciated for his effort in reviewing Chapters 4 and 5. A very special note of thanks goes out to Bernard Goodwin, Vice-President, Professional and Technical Reference Publishing of Prentice Hall, who has taken on the task of publishing this book and has continued to support enthusiastically throughout this project. I would also like to extend my appreciation to Michael Hays, Executive Editor and Assistant Vice-President, Professional and Technical Reference Publishing of Prentice Hall, for his constructive advice and his special interest in this project. Also deserving special mention is Brendan Stewart, Production Editor of Prentice Hall, for his excellent job throughout the production process. I am equally obliged to other Prentice Hall personnel who have participated behind the scenes in the different stages of this book project. Finally, I would like to acknowledge those individuals and organizations that have permitted me to reprintladapt portions of their outstanding work in order to make this book a well-rounded source of information. They include (in alphabetical order): AACC, AGARD, AIAA, Arthur Gelb of TASC, Frederick W. Hardy of Hughes Missile Systems Group, Robert J. Heaston of GACIAC, Kaz Hiroshige formerly of General Dynamics Convair Division, IEEE, Johns Hopkins APL Tcchnical Digest, Robert J. Kelly of Allied-Signal Aerospace Co., Litton. James A. McLean of the U.S. Army Missile Command, Pergamon I'rcss, Rockwell International Collins Avionics Division, and SCS International.
Modern Navigation, Guidance, and Control Processing
Introduction This chapter provides an overview of advanced navigation, guidance, and control (NGC) design. The text is oriented to the applied rather than the theoretical aspects of the subject matter. Although advanced techniques are discussed, the contents are presented in such a manner as to also provide a simple and interesting picture of the central issues underlying both classical and advanced control theory and the practice of the NGC modeling, design, analysis, simulation, and evaluation (MDASE) process.
1.1 OVERVIEW The objective of this book is to provide both engineers and managers with the advanced NGC knowledge and concepts necessary to make sensible decisions concerning the best NGC system design in a particular situation. It will hopefully serve as a useful source of information for NGC systems designers by providing them with ideas for the solutions of current problems and future designs, as well as information about the problems encountered with microprocessor-based systems, microelectronics, aerodynamics, structures, propulsion, sensing, actuation, target acquisition, and weapon systems. In an attempt to achieve such an objective, this book demonstrates practical tools that are realistically applicable in the work area.
2
Introduction
Chap. 1
The development and applications of present-day systems and control theory were spurred on by the challenge of unsolved aerospace problems, especially by the series of events that has occurred since the late 1950s. Therefore, it is beneficial to review the development of systems and control theory since that time. The jolting success of the Soviet Union's satellite technology during that time inspired the United States to excel in aerospace technology, thus giving birth to an entire new generation of support for the field of systems and control. This in turn led to success with regard to solving urgent aerospace problems. The emergence of the Apollo program in the 1960s restored confidence in systems and control research and provided opportunities for conceptualized systems and control theories to be transformed into actual practical NGC system designs. Among the more well-established concepts that found great applicability to solving real engineering/control problems in the 1960s are the recursive minimum variance estimator, often referred to as the Kalman-Bucy filter, the LQG technique, the time domain concepts related to the fields of linear algebra and probability, and the Kelley-Bryson variational optimization procedure. Research into this last area proved particularly fruitful, as it resulted in the successful design of optimal trajectories for several space missions. In particular, the Apollo program and space shuttle flights that followed benefitted greatly from the application of optimal control theory. The contributions of Kalnlan filtering to the Apollo program represent yet another nlilestone for the emerging control theory of the late 1950s and early 1960s. The Kalman filter was implemented using a square-root algorithm in the measurement of stars taken by the astronauts with the aid of a sextant. Both optimal control theory and Kalman filtering provide tangible examples of the cruc~alrole played by control science in important programs such as the Apollo project. As a result of its early success in solving mainly aerospace-related problems, NGC theory soon found application in such diverse areas as mcdicine, industrial manufacturing, and energy management. For example, two French researchers werc able to formulate the problem of treating certain cerebral edemas or malignant brain tumors by the simultaneous administration of vasopressin and cortisone into a nonlinear rnultivariable control problem. The treatment resulting from their work is currently used in the neurosurgery clinic of the Hbpital de la Pitie in Paris [Fleming, 19881. Although examples of this type demonstrate the profound impact of control theory in general, and are not altogether rare, it remains true that the majority of progress in N G C research and development continues to find its greatest application in aerospace. The preceding historical account is provided to illustrate the emergence of advanced NGC systems and their applications. The emergence of the various N G C methodologies contributed significantly to progress achieved in the development of state-of-the-art systems and control theories in the 1960s. Unfortunately, it was not long before the development of the NGC system became too mathematically oriented in spite of the pressing need to solve many remaining practical problems. It would appear reasonable that the obligation of the NGC engineer should be first to commit to solving practical NGC problems [Ho, 19871. In the course of carrying
Sec. 1.1
Overview
3
out such an obligation, thc cngincer will oftcn discover relevant theorics that surface spontaneously. While theoretical discovcries of this type often allow thc engineer to conduct furthcr systcm research, it remains impcrativc for professional engineers to commit themselves to solving what are considered to bc rcal-world problems. Very oftcn, practicing engineers encounter difficulty in just understanding the theorics presented at conferences and in journals, not to mention applying them. What is at issue here is not a debate of the merits of theory versus those of applied technology. Rathcr, thc issue is onc of determining an effcctivc way to build a practical N G C systcm; understanding an csotcric N G C theory is of sccondary importance. While there is currently an abundance of N G C problems waiting to be solved, along with a whole assortment of tools capable of attacking these problems, diminished funding in this area hinders further development of N G C technology. In order to recapture the strong financial support for NGC research of the 1960s, the N G C engineering community must first be seriously committed to attacking those immediate ~racticalNGC ~ r o b l e m sthat remain outstanding. " In the absence of such a first step, any practical system and control theory is not likely to be developed, as evidenced by the slow progress of N G C technology in aerospace. T o highlight this fact, N G C tcchnologv used in the Apollo project was later directly applicd to the space shuttle with hardly any new development. As a result of the diminishing research effort in the NGC discipline since the Apollo era, scarcely any systems or control theories have proven valuable enough for practical use in N G C systems. The extent to which system and control theories have facilitated the clarification of various N G C problems and issues cannot be overemphasized. However, the N G C engineering field still faces several practical problems, one of which involves the question of identifying the basic reason for the gap between theory and practice and recommending a means to close this gap for the sake of advancement in the N G C field. In light of this, the impetus behind this book is an attempt to bridge the gap between theoretical arguments and the practical needs of the N G C community by providing detailed discussions and practical examples of MDASE in this area. A more detailed treatment of the MDASE cycle is given in Book 1 of the series. Because experimental and theoretical aspects of NGC system research both represent integral parts of systems and control science, it is essential to consider N G C theories and physical systems together. If an N G C system of a flight vehicle or robotic system is carefully modeled, it usually serves as a strong foundation for the design process which, in turn, leads to analysis and simulation for validation. The final step of the procedure involves evaluation. Thus, the MDASE cycle is an important process in yielding practical results. Through demonstrations of the various techniques and examples in this book, the author hopes to shed light on and unravel complicated systems and control theories, and to translate them into practice. Before they can be brought into the operational stage, advanced N G C algorithms must pass through a developmental stage involving several years of studies and experiments. Through organizing sessions, panel discussions [Lin, 1983-1985; Lin et al., 1986; Lin and Speyer, 19851, and workshops [Lin, 1987; Lin and Franklin,
4
Introduction
Chap. 1
19901 on aerospace vehicle N G C systems for the Amcrican Control Conference and the IEEE Conference on Decision and Control over a period of several years, the author has been able to keep abreast of the current trends in N G C technology and to understand its practical needs in different situations. Thus, the author is able to provide a diverse range of advanced NGC techniques in the book which the reader can then apply to various areas of NGC with competence, developing the best intuitive design. Since all flight vehicles and robotic systems share the same N G C techniques, the techniques presented in the book are applicable to these complex dynamic systems. Applications presented throughout the book and techniques for these applications are summarized in Book 1 of the series. In addition, appropriate applications are embedded throughout the text to enhance the reader's understanding of the techniques presented.
Importance of advanced NGC concepts and their impact.
Advanced N G C systems, although designed so that machines might interact effectively with the elements of nature, rely on human elements for their continued progress and success. The beginning of this chapter highlights the role played by control theory in the success of the Apollo project. The successful development of modern fly-by-wire aircraft such as the F-16 fighter jet can also be credited largely to control theory. And while modern manufacturing is becoming increasingly dependcnt on highly accurate process and machine control, it is an unfortunate circumstance that the continuous contributions that are made in N G C technology and the resulting achievements go largely unappreciated by the broader scientific and engil~ccring communities. The way it1 which N G C technology is isolated from the othcr technological disciplines in governmental, industrial, and academic settings is vcry apparent [Speyer, 19871. Consequently, the acrospacc industry has not fully rcalizcd the potctltial of N G C tech11010gy it1 acrospace systems, and progrcss in this area especially continues to be delayed. This situation will continue unabated so long as those involved fail to grasp both the importance of making progress in the area of advanced N G C technology and the relation of this tcchnology to vital functions in certain acrospace projects. For example, in certain commercial aircraft projects, advocates of new methodologies for enhancing safety and ride quality performances are for the most part ignored. While this can be partially attributed to thosc who prefer conventional dcsign approaches such as the root locus or othcr classical techniques, some of the blame must be laid on inadequately trained pcrsonncl, many of whom lack a complcte understanding of thc basic principles of modcrn and advanced control tcchniqucs. At~othcrreason for this is that thc designers of these commercial transports dcpcnd entirely on the aerodynamic dcsign to improve vehicle stability and efficiency, and totally ignorc the merits of advanccd NGC systems design. Those who are usually quick to criticize techt~ologicaladvance have always been skeptical about the complex and unreliabic nature of new technology, doubtirlg the wisdom of investing so much trust in the capabilities of ncw tcchnology. However, competition and survival dcpend on technological advatlcc. It xvould bc a grave
Sec. 1.1
Overview
5
error to halt technological advance just becausc some inconvcnicnces arise at the beginning of each new era. Despite the fact that any criticisms can generally be found to be true in the short term, the long-term positive contributions that advanced technology has to offer should be emphasized. After a new era of technology has set in, new systems perform more complex tasks and become less expensive than their archaic counterparts. The reliability factor, that is, the mean time between failures, also improves several times over that of older technology of similar complexity. The cost does not increase as much as the cr~ticsargue, as evidenced in the recent shift between the technologies of the 1970s and 1980s. Spccifically, in the weapon systems field, the gains have measured up to the cost. Therefore, critics can discard the notion that new technological concepts are not important [Deitchman, 19871. In the realm of aerospace, unlike other technologies such as aerodynamics, structures, and propulsion, advanced NGC technology has yet to distinguish itself as an essential discipline. There are two main factors that account for this predicament. First, N G C technology in aircraft development has been traditionally accorded a secondary role. This is because many of those who direct large aerospace programs are unacquainted with the important role of advanced NGC technology and its impact on aerospace. Second, no hardware or other tangible entities are produced by N G C technology; whereas, shapes are formed by aerodynamics, delivered by structures, and driven by propulsion [Speyer, 19871. Therefore, NGC technology has barely been utilized in any of the aerospace vehicle design processes with the exception of the Apollo project. According to the author's knowledge, only a few newly developed European commercial airplanes employ advanced control laws. Faulty instrumentation and inadequate computational capabilities reinforced the negative view of NGC technology in the past, but even with the rapid improvements existing in the present this view still persists. T o reap the benefits of advanced technology, timely investments in the pursuit thereof are crucial. Most people cite only the recent progress in military technology. However, new technology must be applied to the commercial arena as well if economic dominance is to be possible. For example, airlines tend to purchase ultramodern airplanes from companies employing more advanced technical sophistication. From this viewpoint, Airbus, the European company, is a commercial success story. It is actively involved in investigating N G C issues and incorporating new N G C technology into its products as well as its human resources. Thus, research and development (R&D) of advanced N G C concepts cannot be stressed enough. It is important that advanced NGC design concepts enter early in the design phase. The earlier engineering errors are detected, the less severe future costs and schedule delays will be. An example of this type of situation is brought into focus in the now famous article entitled, "Probing Boeing's Crossed Connections" by Karen Fitzgerald in the May 1989 issue of IEEE Spectrum, in which the author writes: "By the end of Uanuary], the [FAA] had received 12 more reports of crossed wires, not only on 757s, but on 737s, and not only on cargo fire-extinguishers, but on engine fire extinguishers, and one case of crossed wires from engine temperature sensors. . . . As of March 17, the FAA had received 66 more reports of bungled ~
~
Introduction
6
Chap. 1
wiring and plumbing in fire-protection systems on all four classes U737, J747, J757, J7671 of aircraft."' If it is possible to make such errors in those systems which are by n o means as complicated as an NGC system, then errors in more complicated systems will certainly result unless advanced system concepts coupled with the proper quality control techniques are utilized. As shown in Figure 1-la, the overall cost resulting from engineering error increases dramatically as time progresses. Therefore, errors detected early can drastically reduce the overall cost. Although engineering R&D design and analysis constitutes only approximately 15 percent of the entire project cost, as high as 90 percent of the overall cost of the project may rely on R&D results, as illustrated in Figure 1-lb. Therefore, inadequacy or errors in the R&D design and analysis due to ina)
Cost
Detection of )Engineering Error 1 2 (Normalized Year) Cost Resulting From Engineering Error Increases Dramatically As Time Progresses Therefore Errors Detected Early Can Drastically Cut Down The Overall Cost
b,
Cost
OPS: Operating Cost MAN: Manpower cost R&D: Research 8: Development Cost
Life Cycle CDst = OPS t MAN t R&D Figure 1-1 I
Cost lnvolvemcnt
Karen Fitzgerald. "Probing Bocing's Crosscd Connections." IEEE Spccfrurn. May 1989. p 33.
experience or other conservative factors can often lead to a drastic increase in the overall time and cost. Quoting from the same article, the author writes, "A Boeing engineer who asked that he not be identified said a too ambitious schedule for the new 747-400 aircraft has caused wiring errors so extensive that a prototype had to be completely rewired last year, a $1 million job. . . . In a technology that can tolerate few errors, the crossed connections raise the specter of undetected errors in other parts of aircraft. Each new II&D assignment presents additional challenges which will tax the skill and ingenuity of ~ ~ C - N G design C engineer and will only be solved through hard work, experimentation, and tests. Hence, each new assignment must be analyzed on its own merits, and past practices and methods must not be allowed to stifle new ideas and concepts. Both the past and the present must never forget the common denominator for all designs, which is the human element. Unfortunately, in advanced technology, undue conservatism often hinders progress. Ironically, the objective of any technology should be to go beyond the current state o f design concepts. The NGC field, although relatively new, is currently in an excellent state. I t is responsible ā¬or substantial contributions to engineering, science, and economics as well as the standard of living in the United States and other countries. One of the main contributions of this text is a new perspective of the N G C system whereby the NGC system's interaction with other disciplines will establish the basis for truly innovative problem formulations and methodologies. Methodologies used throughout the book are practical and have been employed in the designing a n d testing of NGC systems. In the past, classical designs were generally used for NGC of unaugmented dynamic systems. Thus, some designs covered in the book are inevitably related to the classical approach. However, modern state-of-the-art N G C system design for highly augmented dynamic systems, which is rapidly gaining popularity, plays the principal role in the design methodologies used in the book. Integration of N G C designs with various engineering automation and signal processing systems are the ultimate goals of this advanced technology. The intention, therefore, of this book and Books 3-5 of the series is to meet this challenge.
"'
Advanced NGC systems. It is apparent that the usefulness of advanced NGC systems has infiltrated the modern world. In the modern aerospace field, NGC designers must be thoroughly equipped to handle N G C systems problems and able to draw on a rather sophisticated knowledge base encompassing many diverse fields; in particular, they must draw on their in-depth knowledge of modern dynamics and MDASE techniques (see Book 1 of the series). They must also understand the finer interactions between systems and components in their nonlinear operational range. Thus, it is natural for them to appreciate the fact that an NGC system should be designed as an integrated system. Additionally, they must consider the wide range of applications of modern digital, analog, and hybrid computers. Very often,
* Ibid.. p 34-35.
8
Introduction
Chap. 1
more sophisticated analytical and computational tools are required to properly model N G C systems. The advantage of new technologies can be utilized to the fullest extent through the advances of new mathematics, analysis, and computation. For example, N G C system scientists and engineers are now more actively involved with the designing and manufacturing of microprocessor/microelectronic chips and computers in addition to interfacing with t h m d y n a m i c s , structure, propulsion, sensing, actuation, target acquisition, and weapon systems.
1.2 OUTLINWSCOPE As shown in the Preface, in Figure 1 labeled Modern Navigation, Guidance, and Control (NGC) Processing, the five primary topics that make up this book are: (1) MDASE of N G C processing (Chapter 2); (2) modern multivariable control analysis (Chapter 3); (3) design algorithms for advanced N G C system design (Chapter 4); (4) fundamentals of N G C processing (Chapters 5 and 6); and (5) advanced N G C systems design (Chapters 7 and 8). O f the subtopics that appear alongside each of these five main categories in the figure, only those that are connected by solid lines are treated in this book. Those subtopics that are connected by broken lines are treated specifically in the books referenced after them.
MDASE of NGC processing. For the purpose of continuity, the first of the five main topics, as shown in Figure 1, deals with MDASE of N G C processing in Chapter 2, highlighting some ofthe key features of Book 1. In particular, Chaprer 2 summarizes the major stages in the evolution of the design and development process. Modem multivariable control analysis. Modern multivariable control analysis, which is the second of the five primary topics as shown in Figure 1, is mainly concerned with the analyses of robustness, performance, dynamics, and stability. Chapter 3 deals specifically with robustness and performance analysis. Robustness analysis is particularly important for examining multivariable control systems design. Performance analysis, which is taken up in the later part of Chapter 3, includes the various statistical methods used in the analysis and synthesis of any modern NGC system. The subject of stability and dynamics analysis, which is threaded throughout this book, is one of the principal topics in Book 1 of the series. The performance predicted by advanced control theory is rarcly achieved when designing and developing an N G C system. A detailed analysis of N G C systems must therefore include a treatment of the actual hardware equipment, citing those characteristics that tend to limit system performance. The most difficult aspects of analyzing N G C loops (at least, in missile applications) include, but are not lirnitcd to, estimating their stability, determining their accuracy, and finding the trajcctory of the pursuer vehicle along with the normal and lateral accclcrations necessary to achieve that trajcctory for varying types of targct motion. Thc theoretical and cx-
perimental techniques for treating these problems are presented in Chapters 2 and 3 in this book, and in Book 1 of the series. The methods that together constitute the third of the five main topics listed in Figure 1, design algorithms for advanced N G C systems design, arc listed with this category in the same figure. The first of these, modern filtering and estimation techniques, is covered in Chapter 4, while advanced multivariablc control systems design techniques are the subject of Book 3. The last two topics in this category, adaptive N G C systems (detection and identification) and knowledge-based and neural network approaches to system control are given special treatmcnt in Book 4 of the series. I t is worthwhile at this point to examine the hierarchical and interactive relationships among the different design techniques mentioned previously that make up the design algorithms for advanced NGC systems. These techniques provide a methodology for selecting successful N G C systems and integrating them according to the mission requirements of a particular vehicle. The primary focus of this particular treatmcnt of design algorithms is to describe present-day components, systems, and synthesis techniques from the system-integration point of view. The goal is to demonstrate how to analyze and select an N G C system to meet a set of performance requirements when the vehicle maneuver and dynamic environment are specified. The various techniques which make up the design algorithms are arranged in hierarchical order in Figure 1-2a. From the figure, it can be seen that modern
Design algorithms for advanced NGC systems.
Modern Filering and Estimation
+
Advanced Multivariable
Control Design Techniques (Book3)
(Chap. 4)
+
Adaptive NGC System Design Twhniques
CBook 4)
v v Knowledge-Eased and Neural-Network Awroacha to Svstems ConUol
Figure 1-2 Design Algorithms for Advanced NGC Systems (a) Hierarchical Relationship (b) Advanced NGC Systems
10
Introduction
Chap. 1
filtering and estimation techniques represent the starting point in the use of design algorithms. Results from the application of these techniques are then used for both advanced multivariable control and adaptive N G C system designs; however, adaptive N G C techniques are also used to iterate or improve upon the multivariable control design. The results from the application of these last two techniques serve as input to the design of a knowledge-based and neural network intelligent system. The interactive relationships between components and techniques are depicted in Figure 1-2b. This figure illustrates how the design algorithms are made more intelligent by incorporating a knowledge-based and neural network approach to N G C system design. Ultimately, an N G C system is designed to accurately control the outputs of a system whose dynamics contain significant uncertainties. This involves the following fundamental processes and their associated algorithms: (1) modeling ofthe system based on physical laws (see Book 1 ofthe series); (2) systern identification based on experimental data (presented in Book 4 of the series); (3) signal processing of the output by filtering, prediction, state estimation, and detection (modern filtering and state estimation are presented in Chapter 4, while adaptive N G C systems detection is presented in Book 4 of the series, and digital processing is presented in Book 5 of the series); and (4) synthesizing the control input and applying it to the system (advanced control theory is presented in Book 3).
NGC processing. Referring again to Figure 1 of the Preface, the fourth of the five main subjects which comprise this book deals with N G C processing. As indicated in the figure, N G C processing includes as subtopics navigation. guidance, and flight control processing. The first of these, navigation processing, is covered in Chapter 5 . Guidance processing is treated in Chapter 6, while flight control processing is presented in Book 3 of the series. In the course of the last three decades, N G C systems have evolved from a state in which they existed only in the imagination to their current state in which they are very much a part of reality. Both military and commercial aviation, as well as other fields including automation and manufacturing, owe a good deal of their success and technical advancement to the parallel advancement in the area of N G C systems. It can readily be said that the development of N G C systems and techniques, which in their ,--infancy often led to unpredictable accuracy and reliability performance, has been most successful, especially when examined in light of the often extreme environmental conditions to which the components of these systems are typically subjected. Figure 1-3 gives a block diagram representation of all the signal processing elements that are needed to perform the functions listed in the previous paragraphs. Each block in the diagram additionally lists the corresponding chapter number where the individual elements are treated in more detail. It can be seen from the figure that N G C proccssing involves several processors, each of which must function individually as well as in unison with others. Advanced NGC systems design. This constitutes the last of the five main topics covered in this book. Looking once again at Figure 1 of the Preface, it
Feedback Sensing Information
I
I
I Tracking Sensing Information
+
+
I
Target Tracking (Chap. 7) I
Guidance Algorithms (Chap. 6)
---
Advanced Navigation Guidance System 4 System (Chap. 5) (Chap. 8)
can be seen that this topic is concerned with the following subjects: navigation and guidance filter design (Chapter 7); advanced guidance systems design (Chapter 8); and advanced flight control system design (Book 3). These subjects prepare the reader for the actual design of integrated, adaptive, and intelligent N G C systems design, presented in Book 4 of the series, and of digital NGC systems design, presented in Book 5 of the series.
Modeling-Design-
Analysis-SimulationEvaluation (MDASE) of NGC Processing It is brought to the reader's attention at this point that a very broad picture of the modeling-design-analysis-simulation-evaluation (MDASE) cycle of the NGC system has already been put forth previously in Book 1 of the series. In order to allow for smooth discussions of the modern NGC processing tasks and for highlighting the way in which NGC processing gets involved in the MDASE cycle, advanced NGC concepts must first be examined. The examples chosen in the following discussions serve to introduce various parameters and functions of NGC processing with respect to particular applications.
2.1 LINEARINONLINEAR INTERCEtT NGC SYSTEM The conventional approach taken by design engineers to the avionics has been to break up the avionics into three distinct and independent systems: the navigation, guidance (midcourse andlor terminal), and flight control systems.
Navigation system. The navigation system functions to provide position, velocity, and attitude of the vehicle with respect to a reference coordinate frame. Using high-accuracy gyros and accelerometers, it is conventionally configured as an inertial system in either a gimballed or strapdown mode.
14
MDASE of NGC Processing
Chap. 2
Guidance. From the perspective of a control system, guidance is a matter of finding the appropriate compensation network to place in series with the plant in order to accomplish an intercept. In order for the pursuer to impact a maneuverable target with little miss distance, guidance uses the principles of feedback control. The purpose of the guidance system is to determine appropriate pursuer flight path dynamics such that some pursuer objective might be achieved efficiently. The guidance system decides the best trajectory (physical action) for the pursuer based on its knowledge of the pursucr's capability, target capability, and desired objectives. In many applications, the guidance system is designed so that it makes use of an inertially stabilized tracker (for example, seeker) that directly measures the angular rates between the pursuer and its target in a fixed coordinate frame. The function of the guidance computer is to mathematically integrate the separate functions of navigation and the flight control system (FCS).
FCS. The function of the FCS is to control the pursucr in pitch, yaw, and roll motion. The FCS executes the guidance commands and stabilizes the pursuer in flight. The FCS, upon receiving commands from the guidance law, then issues its own commands to the appropriate aerodynamic andlor thrust controls of the pursuer so that the guidance command can be properly executed. It is usually configured as a system equipped with low-accuracy inertial components (gyros and accelerometers). Small tactical missiles and most large transport vehicles often have open-loop control instead of the more complex FCS control. 2.1.1 LinearlNonlinear Intercept NGC Processing
Guided weapons or missilcs are normally guided from shortly after launch until target intcrccption. The NGC systcrn supplics stccring commands to aerodynamic control surfaces or corrccting clemcnts of the thrust vector subsystem to maneuver the weapon to its targct and to make it possible for the weapon to intercept moving targets. The guidance process, which is ruled by signal proccssing algorithms implcmented in the N G C system, a prcset flight program, or both, is essentially a feedback control systcrn where the pursuer-target engagcmcnt is considered part of the N G C loop. Thc N G C loop consists of the guidance systcnl together with dynamic controls. Elc~ncntsof this loop include an information subsystem, control clcmcnts, an opcrator, and pursucr dynamics. The components for the N G C systcrn arc shown in Figurc 2-1 in which the ovcrall control of thc pursucr is dividcd into two or lnorc loops. The main control loop in thc diagram is thc guidancc loop, which is thc outcr loop that controls translational dcgrccs of frccdom, whilc the inner controlIFCS loop controls pursucr attitude. ,The guidancc loop contains guidancc scnsors for sensing pursucr motion, targct motion, or rclativc motion of thc targct with rcspcct to the pursucr. This information is uscd in thc guidance computer or the corrccting nctworks, togcthcr with information conccrning thc intended flight profilc, to gcncratc guidance (latcral acceleration) commands for the FCS. Thc FCS and actuation in turn direct control surface dcflcctions to altcr thc pursucr's trajcctory. The body ratcs and accclcrations arc fed back to thc incrtial sensors to closc the FCS loop.
16
MDASE of NGC Processing
Chap. 2
The guidance and control laws used in current tactical missiles are based largely on classical control design techniques. These control laws took birth in the 1950s and have evolved into fairly standard design procedures. Proportional feedback is generally used to correct missile course in the guidance loop, which is commonly referred to as proportional navigation guidance (PNG), and is quite successful against nonmaneuvering targets. The controller for a homing missile, in general, is a closedloop system known as an autopilot, which is a minor loop inside the main guidance loop. In addition to the control surface and servon~echanism,the autopilot consists of mainly the acceleron~etersand/or (rate) gyros to provide additional feedback into the missile servos for missile motion modifications. Advanced sensors may measure other variables. N o explicit state estimators are used and the signals are filtered to reject high-frequency noise. Broadly speaking, autopilots control either the motion in the pitch and yaw planes (lateral autopilots), or the motion around the missile axis (roll autopilots). In general, the roll: pitch, and yaw channels are uncoupled and are typically controlled independently of each other. All commands are amplitude or torque constrained to ensure autopilot and n~issilcstability. Classical controllers have two major advantages, simplicity in design and simplicity in implementation, but they also have several problems. A pursuer is designed to complete the basic homing loop, which requires a sensor (a seeker in the case of a guided missile) to track the targct, a noise filter to reduce the effect of noise, a guidance law to generate thc dcsircd guidance (acccleration) commands to home on the target. an FCS to rcccibc thc dcsircd accelcration commandfrom the guidance system and to generate the required acceleration capability to insure interception of the maneuvering targct, and, finally, a good understanding of the physics of the homing engagement itsclf. Figure 7-2 depicts a miss-distance analysis model conlposed of elemenrs rcprcscnting the intercept kinematics plus clclncnts affecting the intcrccpt guidancc dynamics (that is. pursuer dynamics and the NGC system), and illustrates the intcrdcpcndcncc bct~veensystcnl elements. Inasmuch as thc total system perfornlancc is affcctcd by the individual characteristics of evcry elcmcnt, the system cnginecr is grcatly conccrned with this system interaction. The kinematic part inside the lowcr left dottcd box generates the LOS angle in terms of pursucr and target motion. Based on the LOS anglc, the N G C system generates a latcral acceleration for the pursuer, whilc trying to bring the projectcd miss to zero. In the model, a rclativc position ) I , is added to a glint noisc y,, thc result of which is multiplied by the inverse of thc range R to produce a true LOS anglc corruptcd by glint noise. u, in inertial coordinates. A mcasured LOS angle u is thcn obtaincd by adding ro (5. a radomc crror coupling u,and the following angular noisc tcrms: rangc-depcndcnt (thermal) noisc o,,,, rangc-indcpcdcnt (fading) noisc a,, and cluttcr noisc o,.It thcn scrvcs as input to thc scckcr dynamics to producc a mcasured LOS ratc 6 that is itsclf corruptcd by bias u l , and g sensitivity. The output from thc scckcr u is the input to thc guidancc filter. Many signal processing tcchniqucs arc uscd to discrin~inatcthctargct from its background, other targcts, and decoys. The guidancc law functions to gcncratc midcoursc and/or
MDASE of NGC Processing
18
Chap. 2
terminal steering commands u, based on the guidance filter output. Before being sent to the autopilot, u, is limited. In general (except for coordinated turn vehicle), the methods used to generate the guidance command signals u, for guiding a pursuer in each of two mutually perpendicular planes are identical. Consequently, the guidance algorithm needs to be determined forwnly one of the planes. In the treatment presented here, this plane will be the missile s horizontal plane. In Figure 2-2, u, is normally taken to be the lateral acceleration command. As long as saturation effects are ignored, the model shown in Figure 2-2 is a linear time-varying system driven by stochastic inputs. The inclusion of acceleration saturation effects causes the model to become nonlinear.
2.1.2 Modeling and Simulation The digital simulation technique in which the equations of motion are represented by' a set of first-order nonlinear differential equations is widely recognized. The modeling vector form is given as (see Figure 2-3): x(t) = f(x(t), u(t), t) with initial conditions x(tO) = xg
(2- 1a)
where the function f(x, u ) includes all of the model equations and NGC algorithms. Equation (2-la) corresponds to Figure 2-3, in the absence of any noise inputs. Depending on the specific application, the initial conditions are given, and the terminal time tf is to be determined. For intercept flight simulation, rf is determined by either the point of closest approach of the two vehicles or the time at which the range rate passes through zero for the first time. Numerical integration techniques are then used to integrate these equations with respect to time so that the time history of the state from an initial condition x(t0) to a terminal state x(tf) can be obtained. The output of the simulation is typically the time history of the state, from which a performance criterion can be formulated as a function of the terminal state. The flight vehicle to be designed determines the performance criteria in the simulation. In the case of a guided missile, good performance is represented by short miss distance or short range between the two vehicles at t,. The overall simulation block diagram, shown in Figure 2-4, defines all of the state equations. Referring again to
Nonlinear Function
t
-
f(x(t).uo),t) Figure 2-3
Nonlinear Dynamic Systems
MDASE of N G C Processing
20
Chap. 2
Figure 2-3, the modeling vector form, when a noise input vector is considered as shown in the figure, can be written k(t) = f(x(t), 4 t h t )
+
w(t)
(2-1 b)
where w(t) usually, although not necessarily, .denotes a white-noise input vector. A more general case is considered in Figure 2-4, in which the model includes certain error sources resulting from imperfect measurements or unmodeled dynamics. Incorporating these error sources in Equation (2-la) gives k = f(x,
U)
+ b + n ~ ( x t), + m(t)
(2- lc)
where b is the constant bias term; n l (x, t) is a state-dependent random vector; and nz(t) is a state-independent random vector. Equation (2-lc) could also apply to Figure 2-3, in which case w ( t ) is simply the sum of 6, nl(x, t), and n2(t). If M I or n2 represents a colored-noise process, a shaping-filter technique can be used to model either or both so that the new state vector can be augmented to Equation (2-lc). The resulting augmented state equation is driven by the white-noise characteristics only. A Monte Carlo technique, covariance analysis method, or an adjoint method can be used to properly evaluate the selected performance criteria. These techniques are described in Chapter 3. In closed-loop control system analysis, it must be decided what level of detail will be used to represent each element in the guidance loop. If it is desired to study high-frequency system instabilities, then one must employ a representation that is accurate at those frequencies. The elements that may be required to construct a sufficiently accurate representation in this case include modeling the dynamics of the nonrigid system, knowledge of the full nonlinear aerodynamic characteristics, and detailed modeling of the seeker track loop, guidance signal processing, autopilot, and noise sources. Conversely, a complicated model of this type would not be required in the case of a preliminary parametric study of homing trends since the main effccts of interest are found at low frequencies. Rathcr, a simplified representation of the more complex system capable of incorporating the available trim aerodynamic data and the low-frequency approximations to the different subsystcms that make up the guidance kinematic loop would be more appropriate. After being validated, this simplified representation can be used to study the relative performance of various missile configurations as well as the result of varying system paramctcrs such as time constants and limits [Reichert, 19811. In problems involving dynamic modeling and controls, a state model is usually specified in addition to thc nieasuremcnt equation, relating the state and measurement vcctors as shown in Table 2-1. Bcfore implctnenting thc preceding N G C simulation, it is important to understand the roles of thc pursucr system, the targets, and thc cr~virot~ment so that appropriate modcls in Figure 2-4 can be devclopcd for a particular pursuer-targct engagcmcnt scenario. Essential guidance and control softwarc modulcs arc as follows: intcrfacc, flight path planning, strapdown device, navigation, guidance, stabilization and control, fin actuation, target tracking, and self-test. Their relationships are dcscribcd in Figure 2-5.
LinearlNonlinear Intercept NGC System
Sec. 2.1
TABLE 2-1
-
21
STATE MODELING
D t f ~ n xg ( m l ) state vcclor. y
-
(pxl) system output vcclor. z = (1x1) measured output vcetor.
u = (mxl) mnlml vmor, I, i(qil) cxogenram inpu v m o r .?(I)w d a i n s all thc variables that can be measured and arc fed back to thc controller. Thc enuies of y(t) arc the variables lo be conmllcd. up) is h c output of ihs mnuollsr. lc(l] wnsisu of thc reference inputs, the disturbance inpuu (e.g.. process noise w(t)), and the rcnror M i s s ~(1).Note that y(t) and z(t) may have some wmponenu in common. The wnlinuous-lime nonlinear system is dcscribcd by
i t
=x
I
I t
I
Z ( I ) s h(x(1). I) t
,
X
"(I),
I
=x
Q ~ I= 0. I ~ 4q1)*TI
T
I = at)b(t-r)
T Qv(t)] = 0. Qv(l) "(1) ] = R(t) b(1-z)
The continuous-time nonlinear system with d*crete.time nonlinear mcasuremcnls is repnscnled by: i(1) = f(X(1). Nt), 11 + *I) T
(la) ., (lb)
(2=4
q k ) = h[x(tkh kl + v(k); qv(kl1 = 0. Qv(k) v(k) I = Q
Gb) w h e n the parameter k is considered to be a timc t and the measurements an sampled at discrete timc intcwak. The discrete form of the state equation (la) is as follows: k =k l1 k l Qqk)l = 0, q W k ) wT(k)l = or (2c) h applying many MDASE uchniqucs. the system quations (1)need to be linearized as:
where u, w, and v arc defined in Eqs. (1). Eq.(3) is equivalently expressed as i(t) = &(I) + Bqf) t Br i(1) and y = Cx + h ~,
The state equation in Egs. (4) is generally a nonhomogenwm differential quation with forcing term ~ .t f. (where ~ ) f(l) is a command vcctor anUor an anticipated forcing function mulling from a system nonlinearity andior an estimated dinurbance. In the output quation of Eq. (4), h consisls of the measurement nonlinearilia and uncenainties. The discrete form of the system quations (3) is: x(k) = Wk.1) x(k-1) + W(k-1) qk.1) + wk-1) qk-1)
2.1.3 Guidance System Classification
There are several solutions that can be applied to the problem of guiding a missile to an eventual near miss or collision with a target using only the observed motion of the missile-target line of sight (LOS). In one of these, the missile is constantly flying . - directly toward the current target location. This method, termed simple pursuit, suffers from the fact that the missile must generally undergo rather severe maneuvers as intercept is approached. A second solution concerns a constant-bearing trajectory wherein the missile leads the target, much like the way in which a projectile is traditionally fired at a moving object. A collision invariably occurs as long as the missile flies a course that keeps the relative missile-target velocity aligned with the LOS. As shown in Figure 2-2, the missile trajectory becomes straight line when the missile speed, target speed, and course are constant. PNG is a steering law that is designed to produce a constant-bearing course even in the case of a maneuvering target. T o accomplish this, one makes the rate of change of the missile heading directly proportional to the LOS rate. While this technique presupposes that the missile is able to respond instantaneously to changes in the LOS direction,
MDASE of NGC Processing
22
Chap. 2
-
StrapdownAFtU Guidance Navigation o initialization o Midcourse & Terminal o State Estimation o Euler o Position Transformation _* 4 + o Power On Mode o Power Off Mode o Velocity o Bias Compensation o Vertical Launch o Acceleration o Coordinate E m r s o Separation Mode o Estimate Transformation
Autopilot o Rate Stabilization o Acceleration Control o Aerodynamic Coupling Stability
*
t I
Seeker
I
t Relative
Interface o Aim Point o Initialization 0 Status
4-
Self Test o Test Sequence 4-W
-
1
Right Path Planning o Stored Mission Profile o Change Profile
Actuation o Fin Rate Limitation o Fin Angle Limitation
RF Environment & Clutter ECM: Target Kinematics
Figure 2-5
Sensors o Rate Gyro o Accelerometer
,
Dynamics o Equations of Motion o Engagement o Environment
there is in practice some delay in returning to a constant-bearing trajectory subsequent to a target maneuver as a result of finite missile responsiveness and filtering introduced to diminish the effects of noise. Figure 2-6 defines the conventional PNG configuration, in which the measured LOS rate u is filtered through a noise filter to produce a smoothed LOS rate estimate u, and u, is made proportional to u A similar idea is used in modern guidance approaches based on optimization techniques. If measurement noise is to be suppressed, then consideration must be given to a low-pass noise filter. PNG systems generally employ a first-order lowpass filter, which has a transfer function as shown in Figure 2-6. The missile guidante computer, using the preceding relationships, computes steering commands. Besides noise filters, com~ensationsare introduced to offset such factors as variations in the missile's velocity and radome effects. Also, in order to maintain aerodynamic stability along with structural integrity, steering command limits are imposed [Witte and McDonald, 19811.
2.1.4 Flight Control System (FCS) and FCS Sensing That part o f the system comprising the FCS also includes the inertial reference unit (IRU)to measure A,,,. A,,, is the acceleration of the pursuer in response to the guidance command, which is applied to the FCS and results in operation of the pursuer's control surfaces. The control surfaces can be caused by aerodynamic cross coupling, dynamics nonlinearity, and structural vibration to deflect at different orientations from the desired guidance input, thus contributing to miss distance. Au-
Target Signal Processing
Set. 2.2
vmf
Mi& velocity vector
23
--
MiSc without
, , R
Missile maneuver
~ V target ~ velocity vector
The basic proportional navlgation equation
Seeker Guidance filter
Am
a
V,,,? = AV,b
/
Autopilot
-*-
/ 1\
Gs(S)
GF(S)
GA(S)
Missile ~uidance.res~onse / \ maneuver Navigation clos/nng Line sight ratio velocity = Velocity I line of sight
,
rate
Range
g(S)
* k = Ik, Noise Ntcr bandwiith
V, Closing vebcily A Effcaive navigalioo mtw
s+t
Figure 2-6
Proportional Navigation (Courtesy o j K . Hiroshige, 1984)
topilot, actuator, structural aeroelasticity, sensors, and nonlinear airframe dynamics are part of the actual mode of an FCS. Details of the FCS are presented in Book 3 of the series. The useful models for NGC design are summarized in Figure 2-7, Table 2-2. and Table 2-3. 2.2 TARGET SIGNAL PROCESSING
A mathematical model is used to embody pursueritarget kinematics in the form of parameters of the LOS as shown in Figure 2-8, which shows what type of processing is required. A complete simulated engagement of pursuer target tracking begins with the launching of the missile, is followed by a midcourse and a terminal phase, and ends in an intercept and the scoring of missile-to-target miss distance. In active homing, after launch the missile seeker transmits radio-frequency (RF) energy to the target. The target converts this to a synthesized return and radiates it back to the seeker radar. Multiple target returns, which include different ranges, velocities, amplitudes, angular noise, and environments, can be generated simultaneously with or without ECM. The returns are tracked by the seeker. As shown in the relative geometry box in Figures 2-2 and 2-8, the homing loop is closed by calculating the missile-to-target geometry as well as using the RF array to update the relative angular
. MDASE of NGC Processing
Chap. 2
N [NORMAL FORCE1 X AXIS
(PITCH TOTAL MOMENT) REFERENCE Fx
(YAW TOTAL MOMENT1
(AXIAL FORCE) V
Z AXIS (A)
XCG XACC
-
center of gravity lacation arcelermeter location
(B)
6
-
-
c~ TOTAL CN 0 NORMAL FORCE COEFFICIENT
MAX (TRIM1
0
MOMENT COEFFICIENT
-
MACH =CONSTANT ALTITUDE = CONSTANT
k@ya MACH CONSTANT ALTITUDE =CONSTANT
.30°
OMAX ITRIM)
& = . 3 0 ° = 6MAX
6 = 00 b = -100
6 = .20°
TOTAL MOMENT COEFFICIENT CM VERSUS ANGLE OF ATTACK n (D)
Figure 2-7 (a) Horning Missile Axis System and Forces and Moments (b) Forces Due t o Angle of Attack a and Control Fin Uetlection S in a Tail-Controlled Missile (c) Normal Force Coefficient CN VS. Angle of Attack a (d) Total Moment Coefficicnr CM 1,s. Angle of Attack a (Froin [n'erlincs, 1984(a)] u,irk prnnisrio~tfi.otnAACC)
position. As shown in the endgame computation and miss-distance computation boxes, immediatcly before interception, the characteristics of the endgame including the miss distance are determined by extrapolation. Finally, the LOS rates in the pitch and yaw axes are computed and sent to the missile seeker as pitch and yawhead radar errors to fully close the guidance loop. E C M J E C C M modeling involves digital computer simulations of an operating radarlsensing system including the effects from an enen~v'sE C M . This allows the sensine " unit to be assessed under field conditions as to whether it would perform satisfactorily so as to support the tactical control system nlission. A sensing detection model is required to obtain signal-to-noise ratio data as a function of target range. The model must be able to allow these data to be obtained in both clear and noise-jamming environn~cnts. Jamming can be produced by a self-protecting or stand-offjammer operating against the radar signal in either the main lobe or the side lobe. Enemy threats may be taking placc in ascending, level, or descending flight profiles. Multiparh and terrain (terrain-masking and clutter) effects may also need to be considered.
TABLE 2-2
LINeARIZED AIRPRAME RESPONSE
Hjrmal Auxlemlion: A,
= v,?,
? = Aa + I36
"
I1ilch Anxllrr - Aoxlcrdlion: 0 = -(h-1%
6
(
6
s2
-= -
'liansfcr l'maions:
-=
G
Dl
~- I i~ -C HI^
- lib; u = 0 - y - 1*i2)
+ (A+#>&+ <: + All
-
(8)
(1)
kl
w m lag
(2)
c')
-=
-0= .
(5)
6
v
m
- I SkI(~+~lls+~12s2) (4) (6) s2 t (A+l>)s + C + N) k ( I + )&I' (W. - HC + 1%) - 3 (6) s2 + (A+l))s + (: + N) "1 (")
-
Ik,r 1;nil mnlml or c;tn:~rrlnirfr;mw with negligihlc liII from ctmlq,l s u t t ~ i.c.. ~ , 1% 0. knrdynamic lift: A,,, = Vm?
A = -7ff
- -at
-
Vm Atz
CN, aftcr boas1 (thrwl = 0 )
n = -i.,= -a, cNa C
= A Vm (0-Y);
= -Ma = -a,CM,
O =(I
+ s/A)Y
(8)
-
-
D = -Mi, = -ao (
6 = fin &fleaion or thnnt A, =achieved normal ameleralion a1 missile Iznfer o f gravity prrprndicllar lo missile body A,,," = auzleration measured al the amelemmeler 6 rate P = Vm ( H ( z - a l k, = acceleration gain umstant = k 2 = (J((!+/U)), 2 ( 8 1 ~ k 3 = L I P m = MaI(l&u)d) Ul = 2 W w d = (A+I))/(C+AD) if M a < ( ) &- h m p i n g ralio of airframe o natual frequrncy of airframe = Ma > 0 and 2Sdwd i s small d B., = 11o:" = lllC+r\o, . Ta = 4 missile l m i n g rate time mnslanl = -P/(BC-AE) = 111% lIA(if B 0 ) u b = low f r e w n c y lead
4
--
.
& .
-
-
-
rA Courtesy
of [ H i r o s h i g e , 19861 and from [Neslines, 1')84(a)l w i t h pernlissioll f r o t n AACC
-
AIRFRAME IS T H E HISULT O F TRADE-OFFS TO REDUCE 1/A &INCREASE MANEUVERABILITY WHILE MEhTlNG PHYSICALCONSTRAINTS BO kn &n!s%h W c w Wnpw-5253.
I
TABLE 2-3
HIOH ~ R E Q ~ E~ N ~~ ~1 l From - s [Neslines. 19851 with permission from AACC
a(~) Aaua10r: -4(5)
-
1
25,
s2
"A
"1
I+-st-
4 ~ s ) 1 Acaleromaer: -= 2t*m %(S) lt-st"xc
-
-
Rate Gyro: ~ G ( S )-
(1)
0 )
(3)
s2
Suuclural Filter:
1
2h
1+-s+"G
02 1t- s2
6 (s)
"sf
&= US)
1 +2 hsf + "sf
w& ~~~
(2)
s2
~
Nominal Values for a Typical Missile: UG = UW) &% LU, = 220 radkec,
= 0.65.
fA=
(4)
s2
"fF
0.65,
The mode shape cnvclopes a n shown below.
.
MODE
5 .-
2
MODE
FREOUENCY
5
GENERALtYD MASS t l b / ~/in) 0.0517
SECOND MODE
0.6. 0.4 i 0.2
4.6
0.0373
Gym Location
)X
-
LOCATION
bl Normalized Bodv Bendinn Sham
a) Fint and E m n d Mode Shapes The detailed flcxiblc body transfer Function from b
LO body
rae is
(av% &
- 24
1
( I v l s2
-
~ ~ I Y Y &
-. &B
6
-6 <-
Xffi
,1 2Ei l-+-s+l) w1' Wi
Nominal values typically used in thc above equalion are k damping of i t h v i b r a t i o n a l mode 9 t h e v i b r a t i o n a l frequency of t h e i t h .ode n o r n a l i z e d .lop a t r a t s g y r o a t a t i o n 10.00503 rad/in., -0.01986 rad/in.) % n o r u l i a e d airvane deflection 10.315 i n . / i n . , -0.1485 i n . / i n . ) normalized a i r v a n e s l o p e 1-0.02 r a d / i n . , 0.03308 r a d / i n . ) lt. e f f e c t i v e mass 40.0514 l b - s e c 2 / i n . , -0.03131 lb-sec / i n . )
Defining
%mi Klsi
XBL
I n
Xm
%
% I,
-
-*
m i s s i l e molaent of i n e r t i a 11.75 (lb-ft-sec2) l o c a t i o n of missile ffi (5.885 f t ) l o c a t i o n of missile hinge l i n e (11.3108 f t l mass of one f i n (0.0489 slupm) d i s t a n c e between ffi of on8 f i n and hinge l i n e 10.084 I t ) m e n t of i n e r t i a of m e f i n (0.00683 lb-ft-s.c2)
g a i n of s t r u c t u r a l p a t h of the I t h moda
qu*dIatlC a s r o Of the i t h .Dde a t r u c t u r s l path
the general transfer funaion from b lo body rate for each mode is given by
404%
- 24 41"
Sec. 2.2
Target Signal Processing ECM Target Dynamicl&
+ Sewor bra
HITMU
1
Relative ~ c o r n c n y Compurarion
-Target Stale
4 Sellsor Models -el Filtering &A Estimation
Esbres
t Target Starc Pmdicnon
T
rare Esnmator lnrduaoo
Target Tracking Processing
Figure 2-8
The most common measurements of target tracking systems consist of the range R, range rate R, azimuth u, or elevation u, angles (see Figure 2-9). The three coordinates to specify the position of the target are:
1. The x , y , z components of the tracking sensor-target range R 2. The azimuth u, elevation u,, and range R. 3. Direction cosines of the vector R and its length. A tracking sensor or seeker mounted on a homing missile or an aircraft is a device to detectltrack the relative position vector of a target with respect to the missile or the aircraft. The function of the tracking sensor is to receive external signals that are sent to the pursuer vehicle for the purpose of directing its course. Because of its role in providing an essential communication link to the external environment, the sensor becomes the key element of the entire vehicle system. In fact, a vehicle is given a particular identity by the choice ofsensor. The source ofenergy responsible for the signal that is picked up by the sensor may be energy reflected from the target, energy emitted from the target, or energy from the target's own environment. Assuming at this point that the sensor is capable of recognizing hot targets, or targets that emit infrared energy, then the sensor will be selected on the basis of a certain frequency or frequency band emitted by the target and will be optimized to respond to the target signature in question. Given the relation of the sensor to the successful operation of a vehicle, it becomes apparent that the aforementioned optimization to a great extent affects the design of all smart weapons. Although there are other important components besides the sensor'in the vehiclelmissile, the ability to neutralize the sensor reduces the vehiclelmissile to something less than a smart weapon [Heaston and Smoots, 19831. One of the requirements of a sensorltracker is to acquire and track all modes of target day or night. T o do this, an air-to-air tracker uses a laser-auadrant tracker and a centroid track. whereas an air-to-ground tracker uses a scene track, target track, and track adjust. A sensorltracker must accurately update target position relative to IP, OAP, and VRP,
Tracking sensor/guidance comparisons. L
MDASE of NGC PrOCesSlng
I
ELEVAT 1ON VIEW
Chap. 2
PLAN VIEW
(A)
ALTITUDE (BAR0 OR RADAR ALTIMETER) DIRECT
ANGLE RATE ( E - J TPACKER)
Figure 2-9 (a) Visual Line-of-Sight (LOS) Angle Measurements Determine T w o Coordinates of Aircraft-to-Target Position (b) Determination of Slant Range (c) External Grid Reference Targeting (Courfesy 0jJ.S.Ausman, 1986)
permitting target identification day or night. Low-altitude ranging sensors for NGC algorithms include redundant range measurement and accurate range measurement during terminal maneuver. Sensors used in the capacity of seeker, depending on their operating wavelengths, are typically characterized, depending on the sensor's operating wavelength,
SIX. 2.2
Target Signal Proceulng
29
MULTllATERATlON OR TDOA FROM 3 OR MORE GROUND-BASED, AIRBORNE, OR SATELLITE STATIONS GROUND-BASED RADAR
wave (MMW), laser radar, television (TV)ivideo, microwave, antiradiation homing (ARH), acoustic sensors, and multiple sensors. While . ~. of good .quality . . . . .and high resolution, an optical sensor nonetheless can suffer degraded-performance due to weather conditions or time-of-day effects. An IR sensor is designed for use at day or night, in conditions of rain and/or smoke, and is capable of hot-spot detection. I t is also well suited for use against armor. and high-value and ship targets, and has reasonably good resistance to ECM. However, this all comes at a high relative cost. Like the SAR, the MMW sensor is designed for all-weather operation with the capability of cloud penetration, and is particularly well suited to bad weather conditions. MMW sensors can be used against armor and high-value targets, but are not designed tor use against ship targets. They have a moderate resistance to ECM, and are not unduly expensive to build. The laser radar, which utilizes three-dimensional distance information and surface reflectance, is desigaedfouse.atday or night against armor and high-value targets. While not designed for use against ship targets, it has a high resistance to ECM and comes at a relatively low cost. Acoustic sensors are designed to function in underwater-erivironments. A T V sensor cannot function at night and is designed primarily tor use against armor targets. It .,
,
MDASE of NGC Processing
30
Chap. 2
is not very useful against high-value or ship --- targets. -.Although TV sensors have a low resistance to E C M , - & ~are ~ relitiVdy inexpensive. While a forward-looking IR (FLIR) imaging sensor provides a different mode of target detection and recognition, a laser provides a different mode of target range and track. A video tracker microprocessor FLlR consists of sensors and ancillary electronics as well as video processing. Microwave sensors can function in all types of weather conditions and are designed for use against high-value and ship targets, but not against armor targets. They have low resistance to ECM and are relatively expensive to build. ARH sensors are designed for all types of weather conditions and are only moderately expensive to build [Hardy, 19861. Finally, multiple sensors act to process multisensor information. The interceptor/missile uses sensor data that are processed by the seeker, which also acts to ensure that the sensor is receiving as much target information as possible. Related to this last function, the seeker actually orients the sensor so that it can survey, acquire, and then lock-on and track the target. The important components making up the seeker include [Heaston and Smoots, 19831: energy-gathering system, that is, antenna for radio frequency or lens/mirror for IR, visible, and ultraviolet spectral regions, including radome o r IRdomel~vindow stable platform and its associated control system (body-fixed seekers employ inertial sensors for an equivalent function) sensor to convert the received energy into a more usable signal for processing signal processing system used to produce a signal to point the antennalmirror at the target and/or to provide signals to ultimately control the missile's trajectory Rather than focusing too much attention on any particular tracking sensor or seeker, it is more advantageous to consider their fundamental properties. When, however, the problem is that of conducting a precise system analysis, then a specific dynamic model is assumed for a particular seeker, and this is integrated with that o f the pursuer. Three homing configurations are assumed: passive, active, and semiactive. In the passive case, only the LOS angle is available through an 1R seeker. In the active homing case, both range and closing-velocity radar information are available through measurements, but in the semiactive mode, these measurements are bistatic.
2.2.1 Targeting The presentation on targeting in this section follows an excellent lecture by [Ausman, 19861. Position sensor data must be available for targeting so as to be able to locate the pursuer with respect to the target. One or more coordinates of pursuer-to-target relative position is measured with a targeting sensor. This device may also be used
Target Signal Processing
Sec. 2.2
31
to locate the pursuer with respect to a pseudo-target o r offset aim point, the position of which with respect to the real target is known beforehand. Various types of targeting sensors with examples of each are listed in Table 2-4. A particularly common typc of targeting sensor is one that measures the azimuth angle and depression angle of the pursuer-to-target LOS. Examples of this typc of sensor system include gunsights and headup displays (HUD) in combination with an attitude reference. As shown in Figure 2-9a, these sensor systems measure just two angular position coordinatcs. A slant-range sensor such as an air-to-ground ranging radar or a laser range tinder can be used to directly measure the third coordinate of pursuer-totarget position. This is illustrated in Figure 2-9b, which also shows two additional more indirect techniques for determining the third (slant-range) coordinate. Pursuer altitude h can be measured using a barometric altimeter. Subtracting from this the target's altitude hT, assuming knowledge of the latter, results in the quantity h - h ~ A. radar altimeter, however, can measure h - hT directly in a target overfly. If it is not possible to realize a preparatory target overfly, the pursuer must accept thc difference in terrain elevation between the target and the spot below the aircraft at release or pickle as an additional measurement error. Slant range can also be obtained indirectly from LOS angle ratc measurements. One example of an angle ratc measuring device is the Norden bomb sight used in World War 11, which relied on manual tracking during a long, straight, and level bombing approach. By using EO trackers, modern angle rate sensors avoid the long, straight and level runin. The last entry in Table 2-4 corresponds to the external reference type of targeting sensor, which is illustrated in Figure 2-9c. The typical application of this sensor type is a blind bombing of known coordinates. As will be shown in Chapter 5, examples of external references that can locate the pursuer with respect to the target, assuming the same reference system is used to specify target coordinates, include ground-based radar, LORAN, TACAN, DME, and GPS.
TABLE 2-4
TARGETING SENSORS
Generic Tyve Line-of-Sight ( L a ) Angles Azimuth o and Elevation o,
Range R Altitude h LOS Rate b, bc Range Rate k = -Vc External Reference Courtesy 0fJ.S. Ausman, 1986
Exam~les unsight, T kpping ~ a d a r , ' ~ Sonar. ~. Optical, Video Tracker Radar, Laser, Sonar, Optical Instruments, FCS, Passive Ranging Processor, Onboard Active Sensor Baro-Altimeter, Radar-Altimeter EO Tracker Doppler Radar Ground Radar Director LORAN, TACAN, DME, GPS
MDASE of NGC Processing
Chap. 2
Line-of-sight (U)S)sensors.
Common to all LOS sensors is that some type of image of the target area is displayed to the weapon system operator. The operator then selects the target by positioning a marker or pair of crossed reticles over the target and depresses a target designating button. The computer responds by recording the azimuth and depression angles of the target marker (pipper). What the operator sees may be the real visual scene, or it may be an electronic reproduction using TV, IR, or radar sensor subsystems. The resolution and accuracy ofthe shorter wavelength systems (optical) are generally better than that of other systems. Accuracy may be improved by using scene magnification, which can also aid in target identification.
Range sensors. A range sensor, as shown in Table 2-4, provides t w o means of measuring slant range from the pursuer to the target. For weapon delivery, it is not necessary to have extremely accurate range measurement since the bomb miss distance for typical dive deliveries is on the order of one fourth of the range error responsible for the miss. That is, a range measurement error of 100 feet causes a miss distance of only 25 feet, typically. What may at first appear to be an inconsistency is actually a consequence of the fact that the range error causes correlated errors in ground range and in altitude that partially compensate each other. As an example, a very large slant-range measurement causes the computer to estimate the aircraft to be further from the target (short impact) but at a higher altitude (long impact). It should be noted that the range R, the quality of information concerning which would be very important if it were being directly utilized in the solution of the NGC problem, is not an easily measured quantity. For missile applications, the explicit range information may be obtained from either the launching platform or launching shiplground defense site, as when the missile and target are illuminated by the launching apparatus, or from the target sensors listed in Table 2-4. It still may happen, however, that range information is not readily available, or that the price in terms of space andlor dollars associated with the tracking and guidance hardware to acquire reliable range information is prohibitively great. Such can be the case with smaller missiles. Thus, it is advantageous to derive guidance algorithms that do not require range information. Altitude sensors.
Inasmuch as it is a component of the pursuer's flight instrument package, the barometric altimeter provides a convenient targeting sensor. However, by reason of its rather doubtful accuracy, it usually serves as a backup to be used in the event that the slant range sensor ceases to function correctly. Assuming that local barometric pressure is accounted for, the largcst altimetry errors are typically those arising from uncompensated stat~cpressure defect, which is a bias type error that is on the order of a few pcrccnt of the dynamic pressure, and to uncertainty in lapse rate, which is a scale factor typc error that is on the order of a few percent of the pursuer's altitude above mean sea level. Besides altimetry errors, any unccrtainty in the target altitude will result in an error in the computed value of the pursucr's altitude above the target. While this unccrtainty is removed
Sec. 2.2
Target Signal Processing
33
on a bombing range, it can be significant in actual operational usage. Consequently, one expects a certain amount of degradation in accuracy between bombing range results and operational results with this type of bombing system. It becomes unnecessary to know target altitude with a radar altimeter since this measures altitude above terrain directly. Given that the pursuer is not usually situated directly above the target at the instant of release, the flatness of the tcrrain between the release point and the target to a large extent determines the bombing accuracy associated with a radar altimeter. For this reason, it is advantagcous to make a dry pass over the target for a relative altitude update before proceeding to the attack. This, however, is only possible if the threat environment permits.
LOS angular rate sensors. In an angular rate targeting scheme, the relative accuracy of the slant range computation is directly proportional to the relative accuracy of the angular rate measurement. A 1 percent error in measuring angular rate therefore results in a 1 percent error in slant range. For a typical depression angular rate of 20 milliradians per sec, this accuracy translates to 0.2 milliradians per sec, a level not easily achieved with current state-of-the-art trackers. The great majority of angular rate sensing systems that have been built so far are of the T V contrast tracker type. In addition to requiring three-axis stabilization, these sensors have a response time that is scan rate limited. Owing to this last limitation, these devices can measure angular rates only so accurately. The T V contrast tracker is also characterized by the requirement of a single, high-contrast point to track.
Range rate (closing velocity). Determination of the closing velocity V , is done with the aid of a Doppler tracking device. Doppler frequencies, which are one of the principal outputs of the radar sensor, are proportional to V,. Continuouswave-illumination homing systems can exploit the Doppler frequency of the target return to generate a good estimate of V,.
External references. Each of these targeting sensors is a navigation aid, which locates the pursuer in an electronically generated grid or coordinate system. The rest of the targeting task centers around locating the target in this same electronic grid or set o f coordinates. One way in which to implement this type of targeting is to use radar ground control to direct the pursuer's flight, which requires that the pursuer be within LOS of the ground radar. The pursuer, in order to maintain this LOS, must fly at relatively high altitude, particularly when it is very far from the radar. Moreover, accuracy degrades in proportion to the distance from the radar. Other LOS positioning systems include TACAN and DME, which generally rely on at least two ground stations and use triangulation and/or multilateration to establish the pursuer's position. T o a large extent, accuracy is determined by the relative geometry between the pursuer and the several ground stations, also known as geometric dilution of position (GDOP). G D O P has a greater influence on the altitude coordinate than on the horizontal coordinates, unless the pursuer is situated directly overhead a ground station. This is not likely to occur, however, during a
MDASE of NGC Processing
34
Chap. 2
bombing run. The LORAN system is capable of obtaining coverage beyond the LOS and at low altitudes because it uses ground-wave transmission oflow-frequency radio waves. However, since it is only a two-dimensional positioning system, LORAN must rely on another method such as a baro-altimeter to obtain vertical targeting information. When GPS is operational, it will be able to position lowflying pursuers in all three dimensions since its reference stations are satellites. With the GPS system, each satellite can maintain LOS to almost half the earth. 2.2.2 Kinematic/Relative Geometry
Because only kinematic equations are considered in the kinematic method, the pursuer may be treated as a geometric point as illustrated in Figure 2-10. The corresponding guidance is therefore ideal, and the character of the pursuer trajectory can be determined approximately. Moreover, the required accelerations may be estimated. The kinematics of a pursueritarget scenario are contained in the change in time of the vector from the pursuer to the target R. In other practical applications, the kinematics are characterized by two pairs of kinematic equations, one of which determines the relative motion of the centers of mass of the pursuer and target in the horizontal plane, and the other of which determines this motion in the vertical plane. The geometry of the homing guidance in the horizontal plane is described by the following relations:
where, for convenience of discussion in this section, the subscript true in u,,, is omitted in this section from this point on. The nomenclature is displayed in Figure 2-10 and the measurable variables are 2 I12
R = (x? + y,)
,
o = tan-' y , / x ,
(2-3)
The closing velocity V , = - R of the pursuer and target is determined by Equation (2-2a), while the LOS rate u is related to the parameters of motion of the pursuer and target by Equation (2-2b). A set of analogous equations is used for the vertical plane. Differentiation of Equation (2-2b) and substitution of Equation (2-2a) yield R6
+ 2db
=
-
VT
sin(u
- yT) +
V T cos(u
- ~T)+T
(2-4)
+ i', sin(u ,, - y) - V,,, COS(U - y)+ Because kinematic Equations (2-2) and (2-4) are nonlinear differential equations, they render it a difficult process to study and analyze guidance loops. The analysis can be facilitated by linearizing the kinematic equations relative to the reference trajectory of the pursuer, using approximation methods. With initial values given by YT, and yo, theangles y~ and y are perturbed during the engagement by y, and y, where y, and y,,, are very small compared to yr,, and yo, respectively. Consequently, the instantaneous angles of the velocity vectors are y r = y ~ + , y, and
MDASE of NGC Processing
36
Chap. 2
+ yn,. A
l y simde -case, which is of course linear as well, occurs when sin a zz u , cos u = 1, s i n y,, cos y, = 1, sin y, = y,,, and cos y,, = 1. The L O S angle in Figure 2-10 can be approximated as:
y = yo
where t, = tf - t = time to go
and
t f =- flight time
For analysis purposes, assuming the intercept geometry is simplified to a head-on engagement with constant closing velocity, the range is determined from R
= Ro -
VCtZZ V,(tf - t) = I/,!,
(2- 5b)
where R, = V C t fThe relative velocity and acceleration are j/r
=
YT y,
- jrn = VT Y I cos Y T ~- V,n ~m cos =
yT-
YO
y,,, = AT" - A m y = V T j T - V,,,?
(2-SC) (2-Sd)
where
AT^
=
y ~ =. actual target lateral acceleration = V T jT cos y ~ , = ! VT jT
A,,
=
,
=
actual pursuer acceleration = V,,,j cos yo
--
V,, j
(2-5el
Equations (2-5) are linear equations representing the dynamics of the engagement. The signal that results from subtracting A,,,,, from the simulated AT, is integrated twice to generate y,. With this information, Equation (2-Sa) can then be used to obtain the true L O S angle a,,,,. Assuming VT and V,,, are constant, and hence V, is constant, differentiating Equation (2-Sa) leads to the following formula for the L O S rate u: Rewriting the equation with the help of the geometry yields yr = Vr&u - Vru =
~ J T
jm = VT y, cos y ~ "- V,,, y,,, cos yo
(2-7)
or in block diagram form shown in Figure 2-lla. Differentiating Equation (2-7) yields
y,
= V,t,U
- 21/,u = AT,. - A,., = VT jr cos y n , - V,,,j cos yo
(2-8a)
or, in a more general form, y , = RU
-
~ V , U= AT" - A","
(2-8b)
By applying the assumption made so far, Equation (2-4) can be reduced to Equation (2-8b), while Equation (2-7) can be derived from Equation (2-2). The component of the system for y , to I?, based on Equation (2-8b),can be modeled with an in-
tegrator contained in the positive feedback loop, as shown in Figure 2-11b. The kinematic element, however, becomes unstable on account of the presence of positive feedback. That is, the output signal u grows continuously for a constant input y,. A negative feedback (A,,,) loop is consequently added to the model in order to counter the effects of the positive feedback. The guidance law, FCS, and the pursuer dynamics and propulsion are all contained in this loop. From Equation (2-8), AT, is given by
AT, = (Rs
- 2V,)u + A ,,,, =
V<(t,?s - 2)u
+ A,,
(2-9)
2.2.3 Targeting Sensor Dynamics For the active RF seeker and the passive IR seeker, two gimbals (that is, azimuth gimbal and elevation gimbal) are required. In order to describe the generalized seeker tracking dynamics, the RF seeker with an antenna dish is used as an example. Figure 2-10 shows the essential seeker angles and geometrical quantities. As shown in Figure 2-12, due to the presence of lags in the seeker tracking loop and radome refraction effects, the seeker is not pointing directly at the actual target [Murray,* 19841. Figure 2-13a illustrates the tracking loop that describes the dish dynamics in the absence of radome effects and having the seeker bandwidth w, = 117,. In the case of small track loop time constant T , , the seeker dish rate command u; and the LOS rate are approximately equal. Care must be taken to ensure that T, is adequately small to be capable of tracking highly maneuvering targets, but not so small as to induce excessive noise transmission or a stability problem. The stabilization represents a closed-loop system encompassing the dish rate feedback loop.
MDASE of NQC Processing
-
Chap. 2
RADOh ERROR
Figure 2-12 Radome Error Distorts Boresight Error (From [Munay, 19841 with permission . t o m AACC)
Radome coupling and compensation loops. The increased miss distance associated with radar-guided air defense missiles can be blamed on radome refraction error. Miss-distance predictions for nonlinear radome error are historically less than that predicted for a linear radome error with constant slope [Yost et al., 19801. As shown in Figure 2-12, the effect of radome error, a,,is to distort boresight error measurement E. a,,which is the amount by which the direction of the apparent Total LOS Angle Noise v,
Gyro Acceleration
(B) Figure 2-13 Effect
rsuer laation
Sensitivity
Seeker Block Diagram (a) Without Kadomc Effect (b) With Radomc
Sec. 2.2
Target Signal Processing
39
target and that of the true target differ, is a nonlinear function of the gimbal angle q,. Differentiating u, with respect to u, yields the radome slope r, which is not a constant but varies inp predictably with angle-off and from dome to dome. u, itself depends on both pitch and yaw gimbal angle [Murray et al., 19761. Here, the effect of u, on miss distance can be shown quite well using only a single-plane guidance model [Murray, 19841. The radome error a, arises because RF energy passing through the radorne to the seeker antenna is refracted by the dome material. The magnitude of refraction during homing depends on many factors including radome shape, fiticness ratio, thickness, material, tcnipcrature, operating frequency, and polarization of the target 'echo signal. . Other sources of error such as standing waves inside the dome, reflection from gimbals, receiver phasing, and processing errors also contribute to the radome error. All of these factors differ with each flight and therefore cannot be specifically accounted for in the computer algorithm. As a result, the radome error magnitude can neither be precisely measured nor predicted. This is not to say that it cannot be modeled. The most serious aspect is that the magnitude of refraction error is also a function of where the target energy impinges the radome. This causes in-plane and cross-plane radomc errors to become look-angle dependent with the slopes also depending on the roll attitude of the in~erceptor.Specification of radome characteristics involves the radome slope r, which . ~. is ~shown to be-the main parameter in the homing loop. The most optimistic situation is o i i in which the designer would be able to specify the manufacturing tolerances and the limits on the allowable variations of r. As shown in Figure 2-2, if the nonlinear kinematic guidance loop is to include the radome error coupling loop, then the missile body pitch angle 0 must be generated by Equation (7) of Table 2-2. In addicion to being used to generate 0 in Figure 2-2, A,,,l is also applied to both the guidance filtering and guidance algorithms. In this study of targeting sensor dynamics, a PNG system is assumed where the missile's antenna provides a measure of the target bearing u as follows. Looking at Figure 2-10, there is an angle a d called seeker dish made by the seeker antenna's center line with the inertial reference. The geometry tracking error without radome error is ~ , , r f e ~ ~= ulrue- ud (see Figure 2-13a). The radar converts this error into a voltage, which is then filtered. This is done in order to achieve the proper driving signal that will precess the seeker antenna in the direction that acts to reduce ~,,rfe,, to zero. The seeker dish angle and the LOS angle are equal in the case the seeker tracks the target precisely. However, the actual angle tracking error E = ~ ~ ~ +f u, v, (see Figure 2-13b, v, = measurement noise) is the angle between the apparent target LOS and the seeker center line measured by the gimballed seeker dish. The seeker center line is positioned by the executive computer system and periodically updated based on the reconstructed LOS angle 6. The pitch angle 0, measured by the IRU, has added to it the actual gimbal angle u2measured by the gimballed pickoff. This provides UJ (that is, ud = 0 + a,), which is then added to E to obtain the measured LOS angle u as follows: ~~
~
+
u = e
+ ud = brmr+ u, +
= LOS angle measurement
(2-10)
i . ~
MDASE of NGC Processing
40
Chap. 2
where the apparent LOS angle u as seen by the homing eye is equal to the true LOS angle plus the aberration angle u,, as shown in Figure 2-10. In order to simplify the analysis here, it is assumed that u, is a linear function of the gimbal angle u, such that
Making a correction to include radome effects in Equation (2-lo), the tracking error measurement (see Figure 2-10) then becomes - ~ -
r , e - e - ~ g + ~ r + V , = ~ ~ - u , + ~ , + ~ ,
Equation (2-11) is substituted into Equation (2-12) to yield
The angular relationship of Equation (2-13) and the dish dynamics from Figure 213a are employed to develop a seeker block diagram with radome effects shown in Figure 2-13b. Since the radome-induced error here depends not only on the LOS angle but also on 0 as well, 0 is used in the radome compensation algorithms for stabilizing the (negative) radome error slope. T o enlarge the stability region for the positive error slope and for cross-plane coupling slopes, the radome compensation algorithm is inserted at the output rather than the input side of the guidance processing algorithms [Yueh, 1983(a); Yueh and Lin, 1984, 1985(a)]. Integrating 0 to get 0, the latter quantity is then multiplied by the radome error r and used in the radome coupling block of Figure 2-2 to compute the LOS angle measurement. Combining Equations (2-10) and (2-ll), another simplified equation that can be used to calculate the LOS angle measurement is
where a,,,,,is given by Equation (2-5a). The look angle is defined as Assuming radome aberration error u, is explicitly dependent on the look angle U L , then the time derivative of the measured LOS [Equation (2-lo)] to the first-order linearization approximation can be written as
where the radome local slope is The importance of the seeker as a missile subsystem becomes evident when one considers what it must be able to do. The seeker must be able to track a target, remain spatially stabilized or go into a search mode if target track is lost, and measure
Sec. 2.2
Target Signal Processing
41
a guidance signal. Most seekers have a stabilization loop (see Figure 2-13b) to isolate body rotations and track loop to maintain the seeker axis along the missile-target direction. These loops, in conjunction with the appropriate signal processing, permit the measurement of the guidance signal. The knowledge of the seeker dynamics together with information about the radome characteristics and noise sources would be required to develop a detailed model of this measurement process. In order for the tracking and stabilization functions to work properly, the seeker slew-rate capability for a BTT missile must be compatible with thc missile's banking motion, being ablc to deal with potctltially large look angles. It is important also to assess the impact of banking maneuvers on track loop requirements. A detailed seeker analysis must include modeling of the inertial coupling between the two seeker channels and the effects of friction. When analyzing missile banking maneuvers, one must also consider the effects of polarization. In the case of a semiactive, linearly polarized seeker with linearly polarized radiation, the seeker gain becomes zero when the polarizations are 90 deg out of phase. This situation can arise, for example, when the missile undergoes a 90-dcg bank. Also, as the relative polarization of the seeker and incoming radiation change, there can be serious problems having to do with changes in radome aberration error. That is, if the radome is designed for a particular polarization, a change in this angle may produce larger aberration errors. An active seeker may also be subject to polarization effects when used against a large target such as a ship, for example. If the target's reflectivity is a function of the signal's polarization, the aim point may be altered as the missile banks. Should the guidance system mistakenly interpret this rapid change in aim point as a spike in the LOS rate, it may cause the missile to miss the target [Riedel, 19801
Target tracker (seeker) modeling. The LOS rate u is corrupted by bias ub and gyro sensitivity as shown in Figure 2-2. It is convenient at this point to consider the tracking loop time constant T,, the stabilization loop gain K,, and the gyro sensitivity gain K?,as shown in Figure 2-13b. In general, if the tracking loop time constant T, is too large, the seeker will be unable to track highly maneuvering targets. If, on the other hand, this constant is too small, there will be problems associated with excessive noise transmission and, hence, a stability problem. T h e tracking loop time constant must therefore be chosen somewhere between these two limiting values. he stabilization loop gain K, which is the loop crossover frequency should be large (that is, perfect seeker dynamics) so as to speed up the guidance system, with a practical upper bound placed on it so as to stabilize the seeker well. Because it will have almost no influence on the seeker track loop dynamics if K, is large, the seeker dynamics can be approximated by
+
Without radome effect (see Fig. 2-13a) u'dlu = W,S/(S w,) With radome effect (see Fig. 2-13b)
u'd
= [(I
+ r)u,,,
- re] w,/(s
+ w,)
(2- 16)
that is, a,,,, is passed through a target tracker (for example, seeker), composed of a differentiator followed by a first-order lag, to produce a seeker dish rate command which is approximately the LOS rate I?.
MDASE of NGC Processing
42
Chap. 2
Tracking error guidance system seeker model. The transfer function of the seeker as it appears in Equation (2-16) is more representative of a tracking error guidance system shown in Figure 2-13b. Figure 2-14a illustrates a typical tracking error guidance seeker loop in which an error signal is produced by the angle error detector and is proportional to the misalignment of measured target angular tracking error E. This error signal'is then input to a proportional plus integral amplifier whose transfer function is A rate gyro feedback must be included in the design in order to improve the stability margin in the Nyquist sense and to help reject outside disturbances [Garnell, 19801. The output of gimbal dynamics plus pitch rate 8 is the antenna (dish) rate ed in body coordinates, and, since the feedback is defined as dish angle and not rate, the
Total LOS
Angle Noise
v-
Anele Error
Power
Amalifier
Rate Gym
1
-,---~ Sensitivity
(B) Figure 2-14 (a) Typical Tracking Error Guidance Seeker Loop (b) LOS Rate Guidance System Seeker Model
Sec. 2.2
Target Slgnal Rocesslng
43
loop is closed by the implied integration of u d . An evaluation is now made of the quantitics T I and K I K2K3 which fall under the designer's control. The generalized tracking guidancc system seeker model is illustrated in Figure 2-13b in which the tracking loop time constant T, is proportional to the torque loop gain TI of the seeker. The PNG guidance system transfer function is approximated by u
=
V [ (
+ T , s ) ( ~ + TJ)]
(3-17b)
LOS rate guidance system seeker model. An alternative to the tracking error guidance system sccker model is the LOS rate guidancc system seeker model which measures LOS rate u. The target is tracked by the seeker antenna assembly, which is stabilized with respect to missile body motion. The measured LOS rate is given by the sum of the time derivative of the angular tracking error E and the antenna rotation rate u d throughout Equation (2-10). That is, Thus, to obtain a measurement of u, the signal processor provides i ( s u,) and seeker-mounted rate gyros provide u d . The measured LOS rate u is actually very close to the geometric u,,,,. There are roll-stabilized radar seekers that generate the LOS rate u for each axis using Equation (2-183). While it can be shown that the measurement of u becomes independent of the seeker track loop dynamics under this operation, it is also true that this method may lead to noise amplification owing to the process of differentiation. Moreover, the differentiation is electronic and is performed on a signal measured in seeker coordinates. The process of generating u using Equation (2-18a) together with all seeker low-frequency signal processing must be examined in the context of being suitable when the missile and seeker cease to be roll stabilized, but can bank at a rate of a few hundred degrees per second. A LOS rate guidance system seeker model is presented in Figure 2-14b. The PNG system transfer function is approximated by
Parasitic feedbacks.
As illustrated in Figures 2-2 and 2-13b, bodyrotation rates always have the effect of corrupting LOS rate measurements. Imperfect sensor stabilization or unmatched signal processing gains are only two factors that may act to produce undesirable perturbations on the LOS rate. A parasitic loop may develop in which body-rate induced perturbations on the guidance signal cause a perturbation in the commanded and achieved accelerations and, hence, additional body rate. The result of this situation is that the system can become unstable and may suffer degraded performance [Riedel, 19801. Both the body rate and body acceleration are parasitic feedbacks, owing to the necessity of an aerodynamic missile to pitch to an angle of attack to be able to maneuver. Because of the radome refraction effects, the autopilot and seeker dynamics are coupled through the missile bodyrate feedback signal. This creates an attitude feedback loop in which the missile responds to a target LOS change by maneuvering. The resulting rotation of the
44
MDASE of NGC Processing
Chap. 2
missile body causes an apparent additional change in the LOS angle, which closes the loop. This feedback increases as the altitude increases. When the values o f the radome slope I. are large, they represent the worst effect of an undesirable feedback path through the gimbal angle as shown in Equation (2-13). Beside the radome effects, the stabilization loop is fed with a false angular rate due to missile maneuvers. The gyro acceleration sensitivity drift produces a direct unwanted feedback path with a gain K, from the achieved acceleration to the measured antenna rate in space. If the gyro mechanism is not perfectly balanced about the output axis, any linear acceleration normal to that axis to which the gyro is subjected will cause an erroneous rate signal to be generated. Such a signal will result in a tracking error that in turn commands an additional lateral acceleration to close the loop. This effect was first brought to light during a flight test in which the missile flew a helical path to an excessive miss distance. After conducting an analysis, improved mass balance procedures were put into effect. In addition, the gyro was mounted in a direction that minimized the miss-distance impact [Fossier, 19841. The preceding two parasitic feedbacks become important at high closing velocities occurring in short-range or high-altitude flights. Thus, the design of the seeker must take into account the seeker's own performance characteristics and provide sufficient stability margins to handle unwanted feedback paths. As will be shown later in Examples 3-3 and 2-4, the effect of radome error feedback on the control system is determined by the sign of the radome error. The effect of a degenerative error (positive), which is one that tends to reduce the input LOS rate, is to slow the system response and to reduce the autopilot's stability margin. A regenerative (negative) error may bring about guidance instability at very low frequency. Although both can result in increased miss distance, a regenerative error is of greater concern inasmuch as low-frequency oscillations are characterized by large displacements of thc missile. In a situation somewhat similar to what occurs with the radome error, a feedback loop is formed by in~perfectionsin the antenna stabilization system (sec Figure 2-13b). This causes a change in the antenna direction as the missile body attitude, or missile acceleration, changes. In early systems. the stabilization loop consisted of rate gyros mounted on thc back of the antenna, with their outputs electronically integrated to drive a hydraulically actuated gimbal system. Owing to the servo's finitc frcqucncy response, the antenna is unable to kccp perfectly stationary for body motions a t frcqucncies of primary concern, and the resulting motion gcncratcs a borcsight error cquivalcnt to that caused by positivc radon~eslopc. One aspcct of the solution to this problem involved improving thc hydraulic valve rcsponsc so that the head stabilization loop could be closed to a high enough frcqucncy to rcducc this effect to acceptable icvcls. Another aspect of thc solution was clcctronic in nature and centered around the closed-out frequcncy of the stabilization loop. It was realized that this frequency mattcrcd only to the extent that it affects thc gain in the stabilization loop at frequencics of interest, typically 1 to 3 Hz, and that it could bc increased by changing the elcctronic integrator to provide a - 2 slopc instcad of a - 1 slopc, keeping in mind loop stability considerations. The resemblance to radome slope was brought into grcatcr focus when it was discovcred that purposely limiting the DC gain of
Sec. 2.3
NGC System Deslgn and Analysk
45
the stabilization loop, as opposed to allowing perfect integration of the rate gyro, created the same cffect as a positive radome slope. Thus, one could compensate for the negative radome slope at the very low frequencies that produced large miss distances. A rather serious accuracy problem was alleviated by the implementation of this simple compensation scheme. Chapter 8.4.3 presents the design cquations for PNG with thc previously-mentioned parasitic feedback. Thc cffect of body bending is yet another type of parasitic feedback. In contrast to bcing a separate parasitic loop, this cffect is simply a high-frequency autopilot instability in which body bcnding is detected by the autopilot as a motion ofthe missile. In early systems, the autopilot rate gyros were mounted as closely as possible to the nonrotational point for the first bending mode. This in conjunction with the electronic filtering employed made it possible to avoid autopilot instability [Fossier, 19841. Details of this structural vibration control are presented in Book 3. 2.3 NGC SYSTEM DESIGN AND ANALYSIS The N G C analyst works with all members of the design team. He or she must be able to work from simple to very complex analyses and simulations, and to provide preliminary and firm requirements and designs.
2.3.1 Guidance FilteringlProcessing
Simplified pursuer dynamics. A functional block diagram of the PNG kinematic loop is shown in Figure 2-13a with the guidance computer defined previously in Figure 2-6. Seeker dynamics and corruption of the LOS rates due to body motion coupling through the radome are modeled in Figures 2-lja(iii) and 2-15a(iv). Guidance Jdnematic loop. Consider the block diagram corresponding to Figure 2-15b which is simplified from Figure 2-2. The target position measurement is (2- 19) "' = Yr + Yg where y , corresponds to the target maneuver, which is generally considered as a random variable, and y, is the target glint/scintillation noise (generated by the input equivalent radar noise), which is a property of target radar return signal and arises from sources that physically enter the system. Noise inputs. The two target noise inputs in Figure 2-15b are denoted by y , and U ~ (angle N noise of the tracker). These noises, which corrupt the observation, induce miss distance, and impact pursuer acceleration requirements, are assumed to be uncorrelated white-noise processes. In Figures 2-15b and 2-2, u ~ ~ ( t ) is equal to u,,,(t) + u f ( t ) + uc(t).Represented as a linear displacement at the target, y , has a PSD of and, like miss distance, is divided by rhe range to become an angular glint noise u,,i,,, with PSD
+,,
MDASE of NGC Processing
Chap. 2
a) Simplified P u r s ~ Dynamics r TNE U S Angle
i)
I I qne I I
Guidance & N D ~ Gs ('1 Total LOS Angle N o i s
Guidance Conlpuvr G , ($1 Air+c
I *kc,
Dirh
FCS & Pmpukion
I I
Ratc Command
b;l-a
I
Pursuer mnamicr
Simplified Punuer D y o a k (See F i e a)
Figure 2-15 (a) Sirnplificd I'ursuer Ilynarnics (h) Trajectory Ilynarnics Model (c) Equivalent Trajectory Dynamics Model
to be input to the scckcr. It is pointed out here that, evcn under the assumption o f stationary scintillation noisc at the target, this becomes nonstationary at the seeker by virtue of division by R ( t ) . Target scintillation is commonly assumed to be aptly represented as stationary band-limited noise, thc amplitude and bandwidth being functions of the physical dimensions and motion of thc target. Inasmuch as band-
NGC System DesIgn and Analysis
Sec. 2.3
47
limited noise and filtered white noise possess the same statistical properties, the simulation can be performed by introducing a simple lag filter with ~ C L T= l/wg as a correlation time constant between w, and y , as shown in Figure 2-2. The PSD of y , is +,\.
= 27CLTrCLT, ~ C L T= glint noise variance
(2-2Ob)
The filter was excluded in Figure 2-15b as a matter of simplicity because it falls outside thc homing loop and has no effect on the system adjoint [Peterson, 19611. The net return signal from a target with many radar reflectors can be modeled as a movement of the apparent radar target position. It is typically modeled as a Gaussian random variable only with zero mean and variance rcLT. that is,
where (r~gET)is 112 to 118 the target's dimension perpendicular to the LOS. That IS, rg;T = WJa, where a = 5 is the rule o f thumb and can vary between 2 a n d 8 according to data, and W , is the wing span of the target [Alpert, 19881. Rangeindependent noise (KIN) uf is the noise inherent in the receiver. It is independent of target return power and is caused by many sources, only some of which include signal processing effects, cross polarization of returns, quantization effects, and gimbal servo drive inaccuracies. uf is modeled as of
- iV(0, rf) with equivalent white RIN PSD +f
= 2?frf
(2-20d)
where ~f is the RIN correlation time constant. Thermal noise a,,,is the receiver noise, which is a function of the strength of the target return, and hence the signalto-noise ratio. u,,, is modeled as
- N(0, r,,)
u,,,
with equivalent white receiver noise PSD +,,, = 27,,r,,,
(2-20e)
where T,,, is the correlation time constant of u,,,and r,, is the range-dependent noise variance governed by (Semiactive)
~,,,sA
= r,,,~(RIRO)',
(Active) r , , ~ = r,,~ ( R I R o ) ~ (2-200
where Ro is referred to as the target-to-pursuer reference range. The clutter noise u, is modeled as u,
- N(0, r,) with equivalent white clutter noise PSD +, = 2~,r,
(2-20g)
where 7, is the correlation time constant of a,. Glint noise, receiver noise, and RIN are the important contributions to tracking accuracy when radar is used. Decreasing the range to go increases the glint noise but decreases the receiver noise. The variance ~ T N and the spectral input + T N for the angle noise arNare determined by the characteristics of the seeker and are of the following forms, respectively:
MDASE of NGC Processing
48
Chap. 2
The total LOS angle measurement noise is then Vm
=
UTN
+
uglinr
=
urn
+ uf
f uc
+
Ug*llittr
(2-22)
The angle noise variance and spectral input are of the following forms: v, = v,,,
+ rf + r, + rCLTIRZ
In addition, the tracker noise, or.., can be added to target position noise, y,, so that the measurement equation is modified from Equation (2-19) as: z(t) = y T
+n
=
target position measurement
where n = target position noise = o , R = y, of n are of the following forms, respectively: r,, = +tt=
YTN
R2
+ R U T N , and the variance and PSD '
+ YCLT = [rm + r, + v , ] R 2 +
~ T , v R+~ 4 s =
[@r,z
(2-24)
+ 4f
f
+c]R2
+
ICLT
(2-25)
Thus, Figure 2-15b has been modified to become equivalent to Figure 2-15c.
Nonnoise inputs. The initial condition at the start of homing (launch, transient, and so on) is the nonnoise inputs to the system. The initial seeker error and the initial turning rate of the missile are both assumed to be constant random variables. In certain instances, either of these two errors may be zero. For example, the seeker initial error input becomes zero if the seeker is lockcd on and settled out at the start of homing. Similarly, the initial turning rate input becomes zero if the missile is not maneuvering when homing begins. Other quantities assumed to be constant random variables include initial heading error and detcrn~inistictarget maneuvering inputs [Peterson, 19611. NGC system design requirements and design considerations. Consider the possibility of choosing a guidance and control law such that the system within the dotted box of Figure 2-1% looks like a guidance filter. Then, pursuer lateral position output y,,, is controlled by the guidance computer and FCS to minT y , ) T ( y T - Y , ) ] , that is, minimum miss. The guidance and control imize E [ ( ~ system must get the missile close enough to kill the target when the warhead explodes. Therefore, an appropriate measure of guidance performance is miss distance. which is defined as the minimum distance between the missile and target. Besides minimizing miss, other goals that the NGC system must attempt to accomplish include maximizing the following: the range of operating conditions, the allowable environmental disturbances, the reliability, and the effectiveness. At the same time, the N G C system must be designed so that both cost and sensitivity to target tactics are minimized. When viewed from the standpoint of meeting these objectives, the guidance and control problem becomcs a matter of designing a compensation filter that optimizes system performance with respect to accuracy, stability, speed of response, and external or environmental disturbances. A high relative stability is
Sec. 2.3
NGC System Deslgn and Analysls
49
measured in terms of good gain and phase margins. Concerning speed o f response, the system must be suitably fast so as to wipe out the nonnoise inputs with small guidance time constant T,,. As for environmental disturbances, the filter must be able to filter whatever noise inputs arc present to acceptable levels [Axelband, 19861.
Advanced multivariable control system design. T o begin with, consider Figure 2-16 which describes a standard procedure for designing a robust, integrated control law (see Books 3, 4, and 5 of the series). After a nonlincar simulation modcl is developed and the requirements for robustness and performance are determined, the results are used to conduct an open-loop analysis. As shown in Book 1 of the series, the purpose of the open-loop analysis is to determine the significant dynamic characteristics and interactions of the airframe, propulsion, and control surfaces. It consists of nonlinear and linear simulation, as well as stability, covariance, observability, controllability, and frequency response analyses. The results of these analyses are used to define the design model. Once the design model is defined, control law synthesis can be performed to develop control law algorithms. Control law synthesis consists of: (1) feedforward controller design to enhance command response; (2) feedback controller design to satisfy stability and robustness requirements; and (3) controller order reduction to minimize implementation complexity. These results are then used to conduct closed-loop analysis which, as Figure 2-16 shows, also considers directly robustness and performance requirements. The closed-loop analysis, which makes use of the control law combined with the openloop simulation model, consists of nonliner and linear simulation, as well as stability, covariance, and robustness analyses. This is done to ensure that the control law
.
P
Model helupmnent ' Wind tunnel test L CFD. CSM * Physical laws SyPtem identification Model validillim
Design Specifieatim Military speciliatiow * Flyingqualities Stability requirements 4-M i u distance 4Performance
'
I c
Open L m p Analysh Design model Model reduction * System limitation
C
,
Conlml Syslem D e i g n CACSD Advanced guidanceand control laws 'Design knowledge h w
Simulation and Performance Evaluation l i m e and frequency responses I.inearand nonlinear Simulatim
Implcmentatia, Fast prototype Envirol~mmlled
Figure 2-16 Robust Integrated Control Law Design Procedure
50
MDASE of NGC Processing
Chap. 2
design satisfies control functional requirements. If these requirements are not adequately met, it may happen that further open-loop analysis is required. However, it may be that the algorithms developed in the control law synthesis stage simply need to be refined, in which case the design model does not have to be altered. Some of the methods available for advanced multivariable control law design are shown in Book 3 of the series.
Classical versus modem guidance and control. As shown in Figs. 2-6 and 2-17a, low-pass filters attenuate the noise component of the sensor signal, given the frequency characteristics of target signal and noise. Also, as can be seen from Figure 2-17a, separation principle allows one to design the optimal estimator and the advanced guidance law independently. Optimal estimators such as Kalman filters ideally separate the target signal from the noise by using information about target (and pursuer) dynamics and noise covariances. The objective of an optimal estimator is to provide accurate estimates of all states and model parameters required for the advanced guidance law without significantly increasing the sensor requirements (and therefore cost). The advanced guidance law provides the optimal steering commands assuming a noise-free environment. Figure 2-17b is a functional comparison of how the current modern approach differs from the classical approach to missile guidance and control. Classical approaches have become firmly entrenched in guided missile designs because they worked well in the rather benign environment of past target engagements and were easily implemented with analog circuitry. Because of such precedence, missile designers have often tried to satisfy the increased performance requirements of modern missiles by increasing the complexity of associated hardware such as seekers, gyros, accelerometers, airframes, and engines. Improving performance nearly always involves a trade-off between more sophisticated hardware or more sophisticated software. While such approaches in many cases have improved performance, the increased hardware sophistication almost always results in increased costs. In fact, the resulting cost has often been so high that the systems either were never developed for operational use or were purchased in small quantities [Gonzalez, 1979(a), (b)]. In spite of many classical terminal guidance and control laws which are still being used for tactical missiles, each characterized by varying degrees of performance, complexity, and seeker-sensor requirements, the increascd accuracy requirements and morc dynamic tactics of modern warfare render contemporary guidance and control laws unsatisfactory in many applications. Whilc many modern guidance and control system des~gntechniques have yiclded better solutions than those obta~nedfrom a classical approach, these remain nevertheless overly complex and are not well undcrstood by industry as a whole. As a result, advanced guidance and control sysfetn design, as extensively discussed in this book, is an ideal approach which combines the most attractive features of both classical and modern guidance and control system design techniques. This advanced approach, coupled with the advent of new theoretical methods and low-costlhigh-speed microprocessing techniques, makes it possible to provide better effectiveness, reliability, tremendous increases in missile performance, and greater leeway in other subsystem performance for the same or less power, weight, and cost.
See. 2.3
NOC System Design and Analysls
c
LINEONIOW1 U T E
Xl'W
C
FILTERED LINEOF. LlOn M T Z
Low C A N FILTER
C
CRWAV LAW
KCELERATIW COMMAND
-
FILTERED LlNEOFalOHT
MTE
OTHER SENSORS
LINEOFJIOHT RATā¬
-
/L=.
L1
XltW
nus rson~orr~~
EIIIMATU
u
EITIMATOR
bDVU(CED OUlOANCE
U X E6ETlER LERATIW COYUAND
U W
TARGET
C L W K A L LIIROhCH NOISE
OUIDANCE LAW ~mo-wAv)
SEEKER
AUTOPILOT
AIRFRAME
MISSILE KINEMATICS
TARGET KINEMATICS
ESTIMATOR
MISSILE KINEMATICS
Figure 2-17 Classical vs. Modern Guidance and Control (a) Guidance Filtering and Processing (b) Functional Comparison ( ( a ) j o r n [Conzalez, 1979(b)] withpennissionjom AGARD ( b j j o m [Conzalez, 1979(a)], CC 1979 IEEE)
52
MDASE of
NGC Processing
Chap. 2
2.3.2 NGC Stability and Performance Analysis The guidance analysis model shown in all parts of Figure 2-15 is composed of elements representing the homing kinematics plus other systems that affect the homing guidance dynamics. The root-locus method and frequency response analysis can be applied to guidance system analysis. Performance analysis techniques are used to evaluate the NGC performance. The guidance system transfer function from a,,,, to A,,, is shown in Figure 2-15a. The stability analysis pertains only to the guidance system transfer function, whereas the miss-distance analysis will apply to the entire homing loop of Figure 2-15b. The homing loop shown in Figure 2-15b can be used to obtain a clear understanding of the requirements of an analysis model, and to obtain quantitative comparisons of miss distance as a function of the significant guidance system paran~eters.Such information is then used to evaluate the effects of N G C parameters on missile performance. Miss distance is brought about by the dynamics of the guidance system transfer function, the simplest representation for which is shown to be a fifth-order binomial [Nesline and Zarchan, 1984(a)]. The desir.ed/actual effective navigation ratio A, which varies due to changes in the flight condition, is one major determinant of homing system performance. Another major determinant of homing performance is the guidance system tinte constant T,, which is the sum of system constants. The ejfertivegrridance systenl titne constant, T;, is defined as the guidance system time constant with radome effect. Examples 2-3 through 2-5 later will discuss radome error and its severity on the guidance and control responsiveness of homing missiles. As will be shown in Example 2-4, the overall system stability needs to be studied before computing the time constant of the FCS because the timc constant is the major factor that affects the miss distance. The trend of the effect of cach feedback loop on the linear, time-\.arying gain systcm is clarified by the preceding discussion using root locus. The adjoint simulation is required because all the transient behavior in the miss-distancc response cannot bc predicted by only the simple root structure. The most efficient method of parametric study can be realized through a combination of the root-locus and the adjoint simulation techniques. The former is initially used to idcntify the region of parameters for system stability, while the lattcr is used to check thc conclusions dcrived from the root-locus study, concentrating on the optimum choice of N G C gains. -
Example 2-1: Guidance kinematic loop stability analysis. A vcry good guidancc stability analysis has already bccn put forth prcviously in Book 1 of thc scrics. Thc opcn-loop phase and gain frcqucncy rcsponsc charactcristics arc gcneratcd so that the stability charactcristics of thc guidance kincrnatic loop can be examincd. Once specific values arc assigned to cach paramctcr, the contribution to this responsc by the linear clcmcnts in this loop arc easily evaluatcd. The opcn-loop transfer function in Figurc 2-l5b with the simplified pursucr dynamics defined in Figurc 2- 15a(i) is
Substituting jo for s in Equation (2-26) for sinusoidal input yields for magnitudc
Sec. 2.3
NGC System Deslgn and Analysls
53
and phase
where GA(w) is the autopilot gain and L$A(o) is the autopilot phase response. Although the linear mode of the autopilot can easily be solvcd, it is difficult to analytically obtain the contribution to this response by the nonlinear autopilot. Thus, the nonlinear autopilot phase and gain characteristics are obtained through computer simulation.
Example 2-2: Comparison of statistical digital simulation methods. This example begins with the simplified model of an intercept homing loop shown in Figures 2-15a(ii) and 2-15b. The FCS and the airframe and propulsion dynamics are approximately represented by a single tag (7, = llw,) transfer function CA(S) = A,,lu, = 1/(7,s + 1) (although a second-order representation seems more reasonable) to include their effects in the guidance gains. The guidance system trans~ 5 T) , = T, + fer function without radome effect is A,,lu = AV, s/(l + ~ ~ where 7, + T,. Included in this simplified model are: (1) effective navigation ratio A = 3; (2) closing velocity V, = 3,000 ftlsec; (3) pursuer velocity V,, = 2,300 ftlsec; (4) 7, = l / w y = 0.1 sec; (5) seeker tracking time constant 7, = l/w, = 0.05 sec; (6) noise filter bandwidth k = l / i , = 10 radlsec; and (7) flight time t f = 3 sec [Zarchan, 1979(a) & 19881. The target maneuver model is assumed to be a Poisson jinking maneuver shown later in Figure 2-21e where B = RMS target acceleration = 161 ft/sec2, and A, = the target maneuver bandwidth = 0.2 sec-'. The PSD of the white glint noise u, and the white range-independent fading noise uy are denoted = 4 (ft2/Hz) and +f = 1 x 10-16 (rad2/Hz), respectively. This example is by used in Chapter 3 to demonstrate the utility of each of the performance analysis techniques via digital simulation.
+,
Example 2-3: Radome error-induced miss-distance predictions. Following Example 2-2 and Figures 2-15a(iii) and 2-15b, included in this simplified model are: (1) effective navigation ratio A = 3.5; (2) closing velocity V, = 4,500 ftlsec; (3) pursuer velocity V, = 2,500 ft/sec; (4) T, = l / o , = 0.45 sec; (5) seeker tracking time constant T, = l/w, = 0.1 sec; (6) T, = l l k = 0 sec infinity bandwidth; (7) aerodynamic turning rate constant T , = 2 sec; (8) flight time tf = 5 sec. The model developed by Murray [I9841 and presented here begins with the fact that system disturbances that cause miss distance result also in gimbal angle motion. The effective slopes r are reduced by this motion, while the effective noise q is increased. This example, largely fdllowing the excellent article by Murray [1984], looks at a linear radome error model that spans the gap between nonlinear and linear radome error miss-distance predictions. Included in this model are an effective slope and an effective noise variance. Gimbal angle motion, which is a function of guidance system disturbances such as target maneuver and range-independent noise (RIN), determines the effective slope and effective noise. In the treatment presented here,
MDASE of NGC Processing
54
Chap. 2
analytical expressions for effective slope and effective noise are derived for a sinusoidal radome error. The guidance system transfer function with radome effect is: 7;
'2 4 rn U
= effecti4e guidance system constant
A' Vcs l + Tis ' A! =
=
=1 + ~ A V / V , ) ' ( T T , ) I I (+~ A V A = effective navigation ratio 1 + rAV,/V,
V (2-28)
+ +
where T~ = T, T, T~ Figure 2-15a(iii) illustrates how the radome slope changes the missile system dynamics. The radome slope r appears twice in A,,,/u. A' is altered in terms of r effect by a factor of 1 over the quantity 1 + rAV,/V,,. For positive radome slope, A' is reduced and the response time of the system is increased, the effect of which is to increase the miss distance due to target maneuver. T h e opposite occurs for negative radome slope, the result of which is to increase the effect of A' and thereby increase the miss distance due to RIN [Fossier and Hall, 19671. System damping is reduced in the case of a negative radome error when higher-order terms are concerned. Figure 2-18a shows a typical miss-distance sensitivity to r. The figure shows that there exists a negative threshold value of r, below which miss rapidly escalates. This value depends on the flight condition through a/-$ and on the control system design through the value of T,. Given the variation of a/-$ with altitude, the value of T, must be increased appropriately so as to maintain acceptable performance [Fossier, 19841. The radome slope r also enters in the denominator of Arn,lu, where the time constant (s coefficient) T; is modified by a factor of r T,A V,/ V,,,. The effect of a positive r is to increase 7i, thereby increasing the miss distance due to target maneuver. A negative r., if it is sufficiently large, will cause T: to go negative, resulting in an unstable guidance system. The cffect of r on 7; is of greater concern at high altitude where T , is large. Given the parameter values, the missile system in Figures 2-lja(iii) and 2-15b gocs unstable for r = -0.045. This second-order missile system remains stable for all values of positivc r, a fact whlch is not significant for the following reason. By adding to the model a typical third-order autopilot as will be given later in Examples 2-4 and 2-5, together with a first-order noise filter, r = 0.095 is required for missile system instability. However, fcr positive slope values much less than 0.095, the miss distance due to target maneuver becomes unduly high.
Miss due to radome slope. Target maneuver and RIN are thc two disturbanccs evaluated. As Figure 2-18b shows, the miss distance increases with raFigure 2-18 (a) Representative Effect of Radome Boresight Error Slope on Miss Distance (b) Miss Increases with Radonie Slope (c) Miss Is Less with l'eriodic Radoinc (d) Target Maneuver Miss Is Less with Periodic lladome (e) Gimbal Angle Motion Deterniines Effective Slope and Effective Noise (f) Target Maneuver and Noise Cause Gimbal Angle Motion (g) Effective Slopes and Effective Noise Are Derived by Mininiizing Mean Square Error ((aJ,fmitr /Fossirr, 1984/, 8 1984 AIAA (h)-(g) jam IMuwoy, 19841 urith permissiotr j o m A A C C )
.-
~
89
CONSTANT SLWE I?
8. 4.
8 2. 5 0.
.2. -4.
3
2
1
TARGET MANEUVERTIME TOGO i
3 9 TARGET MANEUVER 7. mred RANGE INDEPENDENT NOISE
4
rl
(D)
TOTAL
-----.O.W
.0.02
0
0.02
0.04
0.06
0.08
RADOME SLOPE
q = r,
SINUSOIDAL RADOME ERROR:
-
TARGET MANEUVER:
- ...-
- . .-
RAWME
0, !0.08.o.m
- 0 . ' ~ .o.b2
o
0.b2
0.04
MAXIMUM RADOME SLOPE
0.05
I 0.08
RANGE INDEPENDENT NOISE:
(C)
t
Effective Slope for Target Maneuver Gimbal Angle Due to Target Manewer
Gimbal Angle Due to
Range Independent Noise
Oggf
Effec~iveSlope for Range Independent Noise
sin (+(up
Effective Bias b
*))
d
'
MDASE of NGC Processing
56
Chap. 2
dome slope. This miss due to target maneuver increases with positive slope, whereas the miss due to RIN increases with negative slope. For this example, the initiation time of the 3-g step target maneuver is uniformly distributed with 25 samples over the 5-sec terminal flight time. The amplitude standard deviation of the RIN is 2 mrad within a bandwidth of 1.8 radfsec (7 = 0.55 sec). In Figure 2-18b, the RMS value for zero radome slope is 1.0 because all RMS miss distances in this example are normalized by the same factor.
Miss i s less with periodic radome. Actual radome is periodic with gimbal angle [Murray, 19841. As shown in Figure 2-18c, periodic radome error results in smaller miss distance than occurs with constant maximum radome slope. The figure shows that miss is significantly reduced for a sinusoidal radome error with a 30-deg period. Specifically, the miss distance due to target maneuver is less for a periodic radome error. Miss distance is plotted against target maneuver time in Figure 2-18d for three types of radome error. The first corresponds to r = 0, the second to r = - 0.045, and the third to a 30-deg period sinusoidal radome with a maximum slope of 1 0.045 1. At the initiation time of the target maneuver, the gimbal angle is such that the radome slope has the maximum negative value of -0.045. While the RMS target-induced miss is 0.7 for r = 0, it is only slightly higher with a value of 1.1 for the case of the 30-deg period sinusoidal radome with a 10.045 1 maximum slope. These ~ a l u e scan be compared with the RMS targetinduced miss for the case of constant -0.045 radorne slope, which is much higher at 3.4. Actually, this adjoint continues to increase with increasing flight time owing to the fact that the missile system is marginally stable for r = -0.043. Effective slope and effective noise. As demonstrated in Figure 2-18e, the effective slope and effective noise of a radonle error are determined by gimbal angle motion. In this figure, the effective slope r due to gimbal angle motion is less than the maximum slope. The difference between the nonlinear radome error, u,, and the effective slope approximation, r u,, gives the effective noise, q. Gimbal angle motion is brought about as a result of missile system disturbances such as target maneuver and RIN. Gimbal angle motion, particularly that component which owes itself to target maneuver, can be a substantial fraction of a radome error. This example looks at a sinusoidal radome error, target maneuver, and RIN disturbances as shown in Figure 2-18f. The sinusoidal radome is characterizcd by amplitude, period, and phase of ro, P, and @, respectively. The gimbal angle is incremented as a rcsult of target maneuver by the value u,l over the flight time. A uniform distribution with a (1/3)"2u,,/2 standard deviation is assumed. Figure 2-18g shows the two effective slopes, one for each of the two nlissile system disturbances. Minimizing the mean-square error, q, yields r~ = ur(aE)* u ~ ~ l ( u , ~ ~ l l 2r, ) ,= u,(ux) *
02
= u,Z
-
~,~f14.
r+u:I/12
- r34,
(1-29)
b =i
Ssc. 2.3
NGC System Design and Analysls
57
Evaluating Equations (2-29) for the sinusoidal radome given in Figure 2-18f results in thc following normalized effective slopcs and effective noisc variance:
-r-T ~.WAX
-0s(
- rf- -
)
(y)
[[sin
COS
~.\IAX
where r ~ f . 4= ~2 ~ r , ~ and /P
(J,ro
sin
2 ~ 4 ~ (Y) sin
P
Because it is not significant for long flight times, the bias is neglected. Equations (2-30b and 2-30c) are plotted in Figure 2-19a for zero phase a. As the target maneuver and RINs increase. the effective sloves decrease, whereas the effective noise increases. The effective slopes are much less than the maximum slope, M MAX, as the gimbal angle increment, r,~, due to target maneuver approaches and exceeds the radome error period, P. As this occurs, the effective noise approaches 0.7 of the radome error amplitudes. The greater gimbal angle oscillation, + ' I 2 , due to RIN accelerates this effect of reduced effective slopes and increased effective noise.
Linear radome miss prediction. Equation (2-31) represents the linear radome miss-distance prediction, which is the RSS of the components due to system disturbances and the component due to effective noise RangeRMS Miss
SlNUSOlDAL RADOME DENOTES EFFECTIVE SLOPE V A L U E GIMBAL ANGLE STANDARD D E V I A T I O N ldegl I
-
EFFECTIVE
SLOPE, rf
TRAJECTORY DYNAMIC MODEL WITH 2 mrad RANGE INDEPENDENT NOISE
TARGET MANEUVER EFFECTIVE
1---
7.5 deg .0.06 -0.04
.0.02
0
0.02
0.04
R A D O M E SLOPE
(B) 3 0 d q PERIODSINUSOIDAL R A D W E 0 --180d.p INDEPENDENT EFFECTIVE SLOPE. rf
EFFECTIVE NOISE STANDARD DEVIATION 1rnr.d)
:;;,,
.0.0(5
.
o
.
m
h
,
0
0
5
10
15
0
GIMBAL ANGLE STANDARD DEVIATION Id41
5
10
15
GIMBAL ANGLE STANDARD DEVIATION 1-1
(C)
CONSTANT SLOPE RAWME
30dW PERIODSINUSOIDAL RADOME 2 rnrad RANGE INDEPENDENT NOISE
'J
EFFECTIVE I SLOPE AND \
rmr MISS
4-
3-
(A)
SlNUSOlDAL RADDME
0
1
.o.a . 0 . b
.o.b2
o
0.b2
0.k
o.'a 0.h
MAXIMUM RADOME SLOPE
Figure 2-19 (a) Effective Slopes Decrease and Effective Noise Increases with Target Maneuver and Range Independent Noise (b) Effective Slopes are Reduced with Range Independent Noise (c) Effective Slope and Noise are Defined by Gimbal Angle Standard Deviation (d) Effective Slope and Effective Noise Predict Range Independent Noise Induced Miss (e) Target Maneuver Causes Gimbal Angle Motion (1) Effective Slope Due to Noise Is Reduced with
Sec. 2.3
NGC System Design and Analyski 3 9 TARGET MANEUVER
SIMULATION
ANALYTICAL
1
Y
%
SIMULATION
1
TARGETMANEUVERAT
5wTOW
-40 0
TIME TO GO lul
TARGET MANEUVER AT 3ucTOGO
O
1 5
TIME TOGO (ucl
(E) 30 deg PERIOD SlNUSOlDAL RADOME 3 9 TARGET MANEUVER . 2 mrad RANGE INDEPENDENT NOISE
30 deg PFRIOD SlNUSOlDAL RADOME DENOTES EFFECTIVE SLOPE VALUE GIMBAL ANGLE STANDARD DEVIATION (deg)
I
-CONSTANT 1 RADOME
SLOPE
SlNUSOlDAL RADOME
TRAJECTORY
INDEPENDENT SLOPE A N 0 NOISE RADOME I I - 0 . E -0.@4
I
I
-0.02
I
0
0.02
1
I
0.04
0.06 0
1
rf = .0.004
RAOOME SLOPE
0.02 0.04 0.06 0.08 MAXIMUM RADOME SLOPE
(H) 30 dog PERIOD SlNUSOlDAL RADOME (0.045( MAXIMUM SLOPE
Target Maneuver and Range Independent Noise (g) Effective Slope Due to Target Maneuver Is Reduced and Effective Noise Increased with Target Maneuver and Range Independent Noise (h) Effective Slopes and Effective Noise Predict Miss (i) Effective Slopes and Effective Noise Predict Target Maneuver Miss (From [Murray, 19841 with permissionfrom AACC)
5
60
MDASE of NGC Processing
Chap. 2
Range-independent npise (RIN) miss.
The result of RIN is gimbal angle motion, and the oscillation amplitude is greater for a negative radome slope. here fore, in addition to the miss distance being greater with negative slope ( ~ i ~ i r e 2-18b), the gimbal angle standard deviation is larger as well. The result is a lower effective slope and thus a reduction in miss with a periodic radome from that with a constant slope radome as demonstrated later. Figure 2-19b shows the gimbal angle standard deviation. The figure also shows the effective radome slope versus gimbal angle standard deviation for a 1-0.045 1 maximum slope radome. These effective slopes for periods of 7.5, 30, and 120 deg are obtained from Figure 2-19a. The value of the effective slope due to RIN is given by the intersection of these two sets of curves. For @ = 180 deg, corresponding to a -0.045 maximum slope (Figure 2-18f), the effective slope, rf, is -0.028 for a 30-deg period sinusoidal radome as shown in Figure 2-19c. The resultant effective noise standard deviation is 1.0 mrad. The RIN induced miss with a sinusoidal radome is predicted by the effective slope and effective noise, as shown in Figure 2-19d. When Equation (2-31) and the effective slope and effective noise are used to compute miss, this correlates with the miss for the sinusoidal radome. The reason for which the miss computed this way agrees so much better than the unrealistically large miss with a constant slope radome is not difficult to understand. In addition to increasing RIN miss, negative radome slope also increases the gimbal angle standard deviation. The result of this is to lower the magnitude of the effective negative slope, which in turn reduces the miss. Consider a linear radome miss vrediction model for RIN-induced miss as follows: a 30deg period sinusoidal radome, a -0.045 maximum slope at 0-deg gimbal angle, and a 2-mrad RIN. Based on this model, the RMS RIN-induced miss is 1.5, while the RMS miss component due to effective noise is 0.3. The RSS for this model is therefore 1.5. This can be compared to a value of 1.2 for the RSS of the sinusoidal radome error model. In this examvle the RIN comvonent for a -0.028 effective slope is dominant. The effective noise component is relatively small.
Noise and target maneuver miss. Gimbal angle motion is induced as a result of target maneuver as shown in Figure 2-19e. Here, both simulation and analytical results are displayed for target maneuvers initiated at 3 and 5 sec to go. The gimbal angle increment, a,,, due to target maneuver can be expressed analyticall y as
An arialytical expression for a missile acceleration, A,,, due to target maneuver. AT,, for a zero time lag missile system can be found in Fossier and Hall [I9671 and Jerger [1960]. Equation (2-32) is obtained for computing 8, given A,,,, and assumirlg perfect seeker stabilization and therefore that q, = 8. This equation is not highly
Sec. 2.3
NGC System Design and Analysls
61
sensitive to radome slope, r. There is a variation of 12 deg, from 22 deg to 34 deg, in u , ~for , r varying from -0.045 to +0.045 with a value of 27 deg for r equal to 0.0. The results that follow therefore neglect r in Equation (2-32). That is, A, = A. As can be seen in Figure 2-19f, the effective slope due to RIN is dramatically reduced with target maneuver and RIN. This reduction is primarily a consequence of the 3-g target maneuver, which results in a 27-deg gimbal angle increment, u,~,. This caused the effective slope, r,, to decrease from -0.028 (no maneuver) to -0.004 for the 30-deg period sinusoidal radome. As shown in Figure 2-19g. there is also a big falloff in the effective slope due to target maneuver as target maneuver is increased. Thc 27-deg gimbal angle increment, q,,, owing to a 3-g target ma, fall off from a value of 0.038 to 0.014 for neuver causes the effective slope, r ~ . to the 30-dcg period sinusoidal radome. The effective slopes and effective noise predict the miss with a sinusoidal radome as shown in Figure 2-19h. Using Equation (2-31) together with the effective slopes and effective noise to compute miss gives values that correlate well with the miss with the sinusoidal radome. These miss distances are the RMS of four sets of 25 flights each, with the sets having the initial gimbal anglc at the radome phase angle, 0,of 6.0, 13.5, 21.0, and 28.5 deg. Thus, these miss distances are the RMS over RIN, uniformly distributed target maneuver, and gimbal angle. Consider a linear radome miss-distance prediction model as follows: a 30-deg period sinusoidal radome; a 1 0.043 / maximum slope; a 3-g target maneuver; a 2-mrad RIN; and an RMS miss of four radome phases. Using this linear radome miss-distance prediction model, the RMS miss-distance component due to target maneuver is 0.8, while that due to RIN is 0.8 and that due to effective noise is 0.8. With these values, the RSS is 1.4, which can be compared with a value of 1.5 for the sinusoidal radome. In this example, the effective noise component is equal to the target maneuver and RIN components. The effective slopes and effective noise predict the miss with increasing target maneuver as shown in Figure 2-19i.
Empirical radome slope calculation. It is not difficult to find in the literature empirical relationships for radome slope [Donatelli and Fleeman, 1982; Giragosian, 1979; Travers, 1982; Youngren, 19611. That found in Travers [I9821 is examined here as an example to illustrate certain key points. where p is the radome fineness ratio (radome length divided by radome diameter), E is the dielectric constant, f~ is the deviation from design frequency, A is the signal wavelength, K , is the performance measure, and Dd is the dish diameter. The dielectric constant is a function of both temperature and radome material. Walton [I9701 gives the dielectric constant for different radome materials at various temperatures. Generally, radome slope error decreases with increasing seeker dish diameter, with decreasing fineness ratio, and with decreasing signal wavelength. The thickness of the radome is usually constructed to be one half of the wavelength to avoid multiple reflections. The performance measure k, reflects how sophisticated the manufacturing process is, higher numbers indicating an advanced radome grind-
62
MDASE of NGC Processing
Chap. 2
ing process. While a hemisphere ( p = 0.5) possesses ideal electromagnetic properties, it is possible to achieve lower drag with shapes having a higher fineness ratio. Therefore, a compromise must be struck between these opposing factors. Radome slope is related to the total radome swing by r = + r ~ / 2which is treated as a maximum expected value for 90 percent of the slopes for a given radome [Nesline and Zarchan, 1984(b)].
Higher-order system analysis. chan [1984(b)] consider
In Figure 2-15a(iii), Nesline and Zar-
The guidance system transfer functions are
The effective guidance time constant, T;, can be approximated as the coefficient of the s term in the denominator of ~ ~ u a t i o(2-35) n anh is the same as that in Equation (2-28). Here, as in the case of the first-order system, the effective navigation ratio, A', is modified by radome slope according to Equation (2-28). It is easily seen that rand T, have an effect on the guidance time constant and effective navigation ratio the same as occurs in the first-order system of this example. It is possible to achieve larger stability regions if one can keep small any or all of the following: the navigation ratio, the ratio of the turning rate to guidance time constant, the ratio of closing to missile velocity, and the radome slope [Nesline and Zarchan, 1984(b); Peterson, 19611. Because for a given radome, increasing A acts to destabilize the system, the designer must strike a balance between having A to be sufficiently small so as to satisfy stability requirements and at the same time sufficiently large so as to ensure effective homing. For a given higher Ta, radome effects were responsible for stability problems as well as an increase in miss distance. Higher T, is found at higher altitude and lower missile velocity. While T, is usually increased at higher altitudes to improve stability, this brings about a corresponding increase in miss distance. Consequently, the designer is faced with attempting to make 7 ,as small as possible to reduce the miss distancc without sacrificing too much stability. O f course thosc components of miss distancc caused by saturation and radomc/parasitic effects will not be affected by lowering 7,. The minimum achievable guidance time constant, (~~),i,,,given T, and the expected r, can be determined on the basis of a useful stability criterion originally developed by Nesline and Zarchan [1984(b)] as - 0.79 < r h,(VJ V,,,)(Ta/rg) < 2.07. The left-hand side of the preceding inequality gives the lower limit of the allowable negative radome slope, while the right-hand side establishes the positive upper limit of radome slope. If the guidance time constant is increased, there is a slight degradation in miss accompanied by a decrease in the system scnsitivity to radome swing. If the guidance time constant is decreased below
Sec. 2.3
NGC System Design and Analysis
63
the minimum predicted by the stability analysis, it becomcs impossible to attain satisfactory performancc over the expected radomc swing range. When thc resultant performance estimatc is not compatible with the allowable warhcad size, eithcr r or T, must be reduccd. Reducing thc former may involvc blunting the nose, using compensation, or somc other techniquc, while reducing the latter can be done by increasing the lifting surface area, reducing the maximum operational altitudc, increasing missile velocity, reducing missilc weight, and so on [Ncslinc and Zarchan, 1984(b); l'etcrson. 19611.
Example 2-4: Homing missile guidance and control analysis. Trimming thc missile at high altitudc may require that the fin be deflected through large angles. If the fin angle is limitcd, the missile's acceleration capability will also be limited. Determining the fin dcflection required for deterministic inputs can be aided through the use o f zero-lag guidance system results. Neglecting guidance dynamics in Equation (4) o f Table 2-2, fin angle can be expressed according to Nesline and Zarchan [1984(a)] as Achieving in a rapid fashion the large control fin angles necessary to trim the missile at high altitudc requires that the control fin move at high rates. The result of an inability on the part of the control fin actuator to attain these rates can be instability. Combining Equation (2-35) and Equation (4) of Table 2-2, the fin rate response due to LOS rate for the fifth-order binominal guidance system becomes
As was true in the case of fin angle, zero-lag guidance system results are very useful for determining fin rate response due to deterministic inputs. When guidance system dynamics are neglected, Equation (2-37) becomes [Nesline and Zarchan, 1984(a)]
Following Examples 2-2 and 2-3, if the structural filter, actuator, accelerometer, and gyro dynamics are relatively fast, the FCS can be represented by the ideal thirdorder transfer function described by
where T, 5, and o are the FCS desired time constant, desired damping, and desired natural frequency, respectively. The values for the zeros are determined by the aerodynamic configuration (the way in which these values are determined i~ pre-
MDASE of
64
NGC Processing
Chap. 2
sented in Book 1 of the series). The FCS gain is represented by K.This third-order model of the autopilot dynamics is useful in generating a nominal system design. The model supposes only two homing loop disturbances with white glint noise with PSD +, in addition to a uniformly distributed sinusoidal target maneuver shown later in Figure 2-21d. These are two important factors that contribute to high miss distance in a radar-guided PNG system. An adjoint model similar to that to be presented in Chapter 3 for Example 2-2 is constructed based on the homing loop shown in Figures 2-15a(iii) and 2-ljb, which follows the adjoint rules presented in Chapter 3. From such an adjoint model, an adjoint simulation is created using representative values for error sources and system parameters, the results of which are shown in Figure 2-20. Neglecting parasitic effects, the guidance system transfer function is given by
TOTAL rmr
TARGET MANEUVER CONTRIBUTION
0
0.2
0.4
0.6
0.8
NOISE CONTRIBUTION
1.0
FLIGHT CONTROL SYSTEM T l M E CONSTANT k c )
(B)
o
i
2
6
E F F E C T I V E NAVIGATION RATIO, A (A)
NORMALIZED RADOME SLOPE-
NORMALIZED RADOME SLOPE
(D)
Figure 2-20 (a) Effective Navigation Ratio Influences System Performance (b) Flight Control System Time Constant Influences System Performance (c) Decreasing Flight Control System Time Constant Makes System Performance More Sensitive to Radome Slope (d) Decreasing Flight Control System Damping Makes System Performance More Sensitive to Radomc Slope (From /Neslines, 1984(a)] and /Nesline G Zarchan, 19821 utith permissions from AACC and ACARD)
NaC System Design and Analysls
Sec. 2.3
65
-
where A,, = KA is the desiredlactual effective navigation ratio. Since the gain K is varied around unity, the actual and the desired navigation ratios are almost identical. 140, which varies due to changes in the flight condition, is one major determinant of homing system performance. Typical performance results displayed in Figure 220a show that if A,, is too high, miss distance increases due to noise and, if it is too low, miss distance increases due to target maneuver. Therefore, A. has an optimal value based on the error sources and system parameters, and in no case should exceed certain bounds. This imposes limitations on the allowable variation in K,as a variafrom tion in gain is equivalently a variation in A,]. For example, a variation in unity along with the corresponding variation in A,, makes it impossible to attain optimal homing performance. Neglecting parasitic effects, T , ~can be approximated as the coefficient of the first-order term in the denominator of Equation (2-40a) as
Therefore, increasing the FCS time constant T increases T,,, while decreasing T de~ and Zarchan, 19821. The system performance is also substantially creases T , [Nesline affected by T.
Autopilot response characteristics necessary to achieve homing. The speed of response of an FCS can be measured to first-order accuracy by 7. The maximum system time constant permissible to achieve successful homing depends primarily on the maneuverability of the targets the missile is expected to engage and the reauired miss distance necessarv to intercept these targets. Seven-tenths of " a second would be an acceptable maximum for a moderately difficult homing engagement. Figure 2-20b shows the effect of T on miss distance [Neslines, 1984(a)]. As shown in this figure, in which it is assumed that the time constants T,, T,, and 2510 are all known, if 7 is decreased, the miss due to target maneuver will also decrease, but that due to noise will increase. Increasing 7 produces the opposite result. For the set of inputs considered it appears from Figure 2-20b that T can be made arbitrarily small to optimize the miss distance, but this cannot be done because of radome refraction effects. When the parasitic radome loop is considered, the guidance transfer function becomes A*,,
where 7; and A' are the same as those in Equation (2-28) except that T, is given by Equation (2-40b). The effect of radome slope error on system stability is itself strongly dependent on T,. At high altitudes, where T, is large, the guidance transfer function can become unstable due to either excessively positive or negative radome slopes. Thus, small T yields a large guidance system sensitivity to radome slopes as shown in Figure 2-20c [Nesline and Zarchan, 19821. Therefore, the optimal value
MDASE of NGC Processing
66
Chap. 2
of 7 is bounded from above by target maneuver (especially at low altitude) considerations, and from below by radome refraction effects (especially at high altitude). The minimum amount of damping required in order to maintain performance, given a reasonable radome. should be determined usingu a detailed complete missile model: however, a damping ratio of 0.3 is a reasonable estimate of the minimum value that a radar-guided missile could tolerate. Figure 2-20d shows the effect of reduced damping and radome slope on miss distance [Neslines, 1984(a)]. The figure also shows typical high-altitude performance for a FCS with high and low damping. As shown in this figure, the FCS damping becomes more crucial when radome effects are considered. A system with low damping is far more sensitive to radome effects than one with higher damping. This means that if the unaugmented airframe is to be used in a radar homing application, a FCS must improve the low damping.
Example 2-5: Radome error compensation. 2-4 and Figures 2-15a(iv) and 2-15b, assuming 7, and assuming
7y
-
=
Following Example Ilk = 0 sec infinity bandwidth,
effective FCS time constant
;---
27A
(2-42a)
Combining the guidance equation (2-16) and autopilot Equation (2-42a) then gives
A,, -
1
+
A V
Substituting for
6 yields the guidance system transfer function
~ where 7; and A' are the same as those in Equation (2-28) except that 7, = 2 7 + 7,. With AV,/V,,, = 6, r = 5.02, and T , = 10, A V,rT,I V,,, can be & 1.2 sec and hence the guidance system becomes unstable (that is, 7,;<0 ). The radome crror slope can be the dominant factor in the guidance response time [Hiroshige, 19861. A guidance and control system is designed to compensate for the stability problcm caused by the radome error. To alleviate the severe radome stability problem for a high-altitude engagement scenario, one of the commonly used engineering "fixes" is to introduce pitch rate compensation G,0 with gain G , feedback to either before or after the guidance filtering network as shown in Figure 2-2. It can somewhat symmetrize the positive and negative slope stability region at the cost of slowing down the autopilot response [Yueh and Lin, 1984, 1985(a)]. Adding the pitch ratc compensation and the guidance filtering (1/(1 7.5)) to the guidance cquatiorl
+
NOC System Design and Analysis
Sec. 2.3
67
(2-42c). The guidance system transfer function bccomcs
-A,,,, --
(I
urn,<. s(l
+ r)AV
+ T ~ ~ ) [ ~ / G , ,+( SG,AV,(I ) + Tas)lV,,,llGs(s) + r A V 4 + ?,s)/Vm
The transfcr function for thc FCS, G.d(s) can be modeled as a third-order transfcr function as described in Equation (2-3'9) or (3-42a). If thc nonlinear elcmcnts are not included in the analysis, Equation (2-43a) based on the linear FCS model is easily solved. Howevcr, if it bccomcs necessary to employ the nonlinear FCS model, it is difficult to obtain an analytical solution using the classical frequency-response approach. Assuming the input to the saturated nonlinear elements is a sinusoid, then the nonlinearity can be represented by a variable gain N obtained from the describing function of the nonlinear element to be presented in Chap. 3. The transfer function ) N(ao + als + azr2)l(bo + bls + bzs2 for the FCS would then be given by G . 4 ( ~ = + b3s3) where the parameters are specified according to a particular condition. If G,(s) = G , = constant, C s ( s ) = 1/(1 7,s) have been assumed, then from Equation (2-43a).
where the effective autopilot time constant 7, -- ~ T Athe , guidance filter constant T,, and the radome compensation gain G , can be selected to compensate for the stability problem caused by the radome error. In comparing Equations (2-28) and (2-43), one can see that the radome compensation G,0 has reduced A' and increased 7;, both of which tend to stabilize the system of a given radome but at the expense of decreasing the guidance response time. The root-locus technique is applicable only to a linear, time-invariant system. The root locus for a linear, time-varying system of Equation (2-43a) can be obtained with the roots at different t,. The evaluation can be performed on the transfer function given in Equation (2-43a) for T,, T,, and G,, which vary with respect to t,.
Example 2-6: Gyro sensitivity gain design. In this section, a method of approximating aeroelastic effects on the dynamic behavior of a rate gyro is developed. More detailed derivations of a flexible body transfer function from control surface deflections to body rotation rates have been presented in Book 1 of the series. A gyro sensitivity gain expression of body-mounted rate gyro coming from body acceleration due to the aeroelastic vibration is considered here. The following assumptions are made in this derivation: (1) the first mode of vibration considered here is excited by body accelerations, (2) only longitudinal and lateral vibrations are considered, and (3) superposition of the aeroelastic angular rate at any given point
MDASE of NGC Processing
68
Chap. 2
on the rigid-body angular rates is valid. The transfer function from the generalized acceleration f to the normalized body-bending curve can be written as
where 5 1 is a structural damping rate, w l is the natural frequency, and Z is the reference displacement to which the mode shape 1 shown in Table 2-3 is normalized. The location of the body-mounted rate gyro on the normalized body bending mode shape 1 shown in Table 2-3 is arbitrary, and it is desired to derive an expression for the aeroelastic body angular rate at this station. T o this end, the slope of thc normalized bending curve at the rate gyro station is determined first. Then, the time rate of change of this slope, due to acceleration, is approximately equal to the aeroelastic body angular rate at the rate gyro station. The product of the normalized displacement Z N and the reference displacement Z located at the rate gyro body station give the instantaneous reflection Z D at this station. That is, Z D = ZAjZ. The slope of the normalized bending curve h at the rate gyro station is h = ZN/XD and the actual slope at the same station is K = Z D / X D= hZ. The rate of change of the actual instantaneous slope is K = h i . For small-amplitude oscillations. the slope of the tangent to the curve at the gyro station can be approximated by the angle between this tangent and the reference x-axis. In this case, the aeroelastic angular body rates 8 , at the rate gyro station is
8,
=
K
=
/ z ~
(2-43)
Thus, the transfer function from the generalized accelerations f to the aeroelastic angular body ratcs a t the gyro station is found by substituting Equation (2-45) into Equation (2-44) giving
8,/f
=
s / ~ / ( +s ~251wls
+~
f )
(2-46)
Because the first-ordcr bending frcqucncy is typically very Iargc, thc gyro sensitivity gain to the step accelcration can be expressed to an approximation as
K, = hlw:
(2-47)
For missile applications, the gyro sensitivity gain is a very .small value since the normalized slope h is a small value and the bending frequency w l is a large value of first bcnding modc.
2.4 TARGET TRACKING STATE MODELING
2.4.1 Target N o i s e and Target Maneuver Modeling Targct noise has been presented in Section 2.3. I , and is applied hcrc to gencratc the target mancuver signal-to-noise ratio. Figure 2-21 shows thc cquivalcnt signal noisc diagram, where the signal, representing a target maneuver, is in turn represcntcd by a shaping filter. As was done in Section 2.3.1, it is possiblc to convert thc total
Target Tracklng State Modeling
Sec. 2.4
Target . Position Noise Target Maneuver Model T ~position ~ ~n ~ , Target Position Spectral Density 7 Signal yT Measurement z Shaping Filter H(s)
-
a) De~ermineManeuver: b(l) b)
Brownian Target Position: u,@ --L + J
C)
Zero-Mean White Gaussian Targel Acceleration:
,,
d) Uniformly Distributed Targel Maneuver: us e) Poisson Jinking Maneuver: u, St),,.
Note: u, = white noisc with PSD 9, Figure 2-21
Target Maneuver Model
seeker angular noise to a positlon noise by multiplying the angle noise by the range to go. In general, it is possible to generate a number of random processes by passing a white noise w through a suitable filter. If it is assumed that the target noise n in Figures 2-1% and 2-21 is a first-order Gauss-Markov process with correlation time constant T,, and white noise input w , then the corresponding equivalent white noise PSD of the position noise is where r , and 4, are given by Equation (2-25). Figure 2-21a shows a deterministic model, which is the simplest case. In most applications, target position or velocity or acceleration is assumed to be filtered white noise (see Figures 2-21b through 2-21d). Uniformly distributed target maneuver. For a target acceleration of step form with r a g o m starting time varying uniformly over the flight time, the autocorrelation function of this model is identical to that of the integrated white noise (see Figure 2-21d). The uncertainty in target position y + is described by a third-order integrator with a constant PSD 4, white-noise input us as:
where n T is the magnitude of the step maneuver, and tr is the flight time over which the maneuver is uniformly distributed. Two integrators are used in tandem in order to convert the acceleration to the target position y T , as shown in Figure 2-21d.
MDASE of NGC Processing
70
Chap. 2
With regard to the uniformly distributed target maneuver, it can be seen by examining Figure 2-21d that the noise spectral density level is constant, but the signal spectral density passes through three integrators. As a consequence, the target signal is frequency dependent while the target noise is not at a particular range to go. Moose et al. [I9791 give a short exposition on the evolution of the relationship of target acceleration modeling to stochastic processes. This evolutionary development is graphically depicted in Figur62-22 [Cloutier et al., 19881. Initially, as shown in Figures 2-22a and 2-21c, target acceleration was simply accounted for with white noise (depicted by a correlated process whose spectrum is flat over the system bandwidth). Consequently, tracking could be maintained only for those maneuvers remaining within the envelope of the white noise, and even then performance suffered due to the incorrect assumption of uncorrelated acceleration.
Correlated acceleration process. A correlated acceleration process, Figure 2-22b, can be achieved using Gauss-Markov models [Fitts, 1974; Pearson
(D)
Figure 2-22 Historical Development of Maneuvering Target Acceleration Models (a) Zero-Mean White Gaussian Process (b) Zero-Mean Correlated Gaussian Process (c) White Gaussian Process with llandomly Switching Mean (d) Correlated Gaussian Process with Randomly Switching Mean (Frotn /Clourirr et a!., 19881 with prmisrion porn AACC)
Sec. 2.4
Target 'Itacking State Modeling
71
and Stcar, 1974; Singer, 1970; Vcrgez and Liefer. 19841. The well-known Singer model [Singcr, 1970) has been used with various Kalman-type filters to achieve excellent tracking characteristics over a broad class of large scale maneuvers. Its primary drawbacks arc the need to specify apriori both the acceleration time constant and power of the driving noise and the inability of the model-based filter to track target motion resulting from abrupt changes (jumps) in the acceleration process. The velocity a t which the target is moving is constant, but thc target's lateral acceleration is described as a Poisson jinking maneuver. Such a maneuver, which is one whosc magnitude B (RMS target acceleration) is constant but whosc sign is randomly switching, can be modeled as white noise through a single-polc filter. This is because the acceleration autocorrelation function (RA.,.)I(7) = ~ ~ e - ~ .' f' ' = average frequency) corresponding to this maneuver process is identical to the Poisson process. Thus, the target maneuver models are as shown in Figure 2-21e
AT, = -2fAr,
+ 2B d? w ,
= -A,AT,
+ h,u,
(2-50)
where u, is a zero-mean Gaussian noise, that is,
-
~v,
with PSD
=
-
(2-51)
w , is a white noise with unit PSL), and A, = 2f = the target maneuver bandwidth. Gholson and Moose [I9771 modeled rapid, major changes in target motion by a semi-Markov process. The mean of their process, Figure 2-22c, randomly switched to a finite number of states according to a Markov transition probability matrix, with the time duration in each state itself being a random variable. Successf~~l tracking with this method necessitates a large number of preselected states, and the erroneous assumption of uncorrelated acceleration is still retained. Moose et al. [I9791 extended this work by employing the Singer model to represent correlated acceleration within each of the states, Figure 2-22d. Larimore and Lebow [I9871 developed a similar model based on a parameterized Gauss-Markov process with the parameters changing according to a point process [Cloutier et al., 19881. Lin and haf froth [1983(a)] concluded that oneof the most versatile representations of target acceleration is that of the sum of a continuous-time Gaussian process and a finite basisljump process. The continuous-time process can be obtained by means of a classical shaping filter, while the jump process which is characterized by its basis, jump size distribution, and interjump time distribution, can be obtained by driving a shaping filter with a white sequence.
Other target acceleration modeling.
As stated in Cloutier et al. [1988, 19891, another aspect of target acceleration modeling is a consideration of the aerodynamic characteristics of the target. A winged aircraft accelerates mostly orthogonal to its velocity vector, in particular, orthogonal to the plane of its wings. The magnitude of the acceleration also has asymmetric bounds (positive or negative g) which are set by pilot and/or aircraft limitations. Modeling based solely on target point-mass motion fails to take many of these characteristics into account. Kendrick
MDASE of
72
NGC Processing
Chap. 2
et al. [I9811 modeled normal load acceleration with a random variable whose probability density function (PDF) was asymmetrically distributed between small negative g and large positive g. Nonnormal accelerations were modeled as first-order Gauss-Markov. It was assumed that the orientation of the three-dimensional target could be obtained by means of an EO seeker in conjunction with pattern recognition. Bullock and Sangsuk-Iam [I9841 developed a nonlinear Cartesian model of target dynamics by using polar coordinates to model planar, circular turns under the assumption of constant velocity magnitude. Similarly, Hull et al. [I9831 used polar coordinates to model both acceleration magnitude and angle as random processes. This led to a Cartesian model characterized by state-dependent noise. Finally, the parameterized Gauss-Markov model developed by Larimore and Lebow [I9871 took into account aerodynamic parameters such as bank angle, lift force, and thrust minus drag. A particular maneuver has a given set of parameters, while any abrupt change in the parameters corresponds to the initiation of a new maneuver. Merging a point process, or jump process, with a continuous-time Gaussian process, either additively, as recommended in Lin and Shafroth [1983(a)], or through parametric embedding as performed in Larimore and Lebow [1987], is an excellent way to model the acceleration of a highly maneuverable target. However, of the models of this type developed to date, only one [Larimore and Lebow, 19871 has given any consideration to the aerodynamic characteristics of the target. The natural next step in this technology area is to merge the concepts of those models based on aircraft flying characteristics with those symbolized by Figure 2-22d.
2.4.2 Two-Dimensional Target Tracking State Modeling Consider the two-dimensional target tracking problem illustrated in Figure 2-10. For simplicity, let velocity VT be in the x-direction while V,, = 0 and y = 0 for a fixed tracking station. From Equations (2-2),
R =
v7-cos u,
6
=
- (VTIR) sin a
(2-52)
Differentiating these equations yields the nonlinear state equations
R =
VT
cos u - l,/6 sin g ,ii = - [ ( v ~ R- v ~ R ) / R sin ~ ] u - (VT/R)6 cos cr (2-53)
Following the treatment of Lewis [1986], assume that the target is not accelerating, that is, VT = 0. Then, assuming also very small changes in R and u during each 1 to 10-scc radar scan, one has that R ;= 0 and & = 0. Hencc, random disturbances O R ( [ ) and ow(() are introduced to account for the small changes in R and b during the scan time T, giving
Sec. 2.4
Target nacklng State Modellng
73
The processes w R ( t )and w,(t) are known as maneuver noise. Because the Kalman filter is robust in the presence of process noise, R and & can be tracked. The addition of process noise is done to offset unmodeled dynamics introduced as a result of excluding nonlinear terms. The state equations then become linear as
The measurement model, given that the radar measures range R and azimuth a is defined in Equation (2-3) and is given by
Suppose independent disturbance accelerations are uniformly distributed between t n T,then w R ( t )is a radial acceleration with variance n$/3. Since w,(t) is an angular acceleration, w,(t)R is distributed according to the maneuver noise PDF and w,(t) has variance n$/(3R2). The covariances of process noise w ( t ) and the measurement noise v ( t ) are given respectively by
where :u and a : denore the variances of the range and azimuth measurements, and a :
where