Research on formation of microsatellite communication with genetic algorithm.
Compared to using large, heavy, single-mission satellites, resource integrations of several smaller and smarter satellites have significant benefits in engineering applications like distributed aperture remote sensing due to its advantage in performance, cost, and so forth . To make use of recent advances in MEMS, robust microsatellites can be built . In order to integrate the resources of autonomous and formation-flying groups of microsatellites effectively, the satellites must have good ability to communicate with each other. Autonomy represents minimal dependence on ground stations for communication purpose, so ISLs can be used to allow each satellite to share their individual information and use their integrated resources in order to achieve much more complex application goals. A novel way to perform space missions is to utilize the concept of using satellite clusters, which can cooperate with each other in order to achieve the function of single large-scale satellite. Every microsatellite can communicate with any other one to share processing, communications, payloads, and some other mission functions. Hence, the required function is to spread across all the satellites of the cluster, of which the aggregation can form a virtual satellite.
Many investigations of satellites flying in formation have been presented. A brief review includes ESA missions Proba 3, Darwin, SMART-3, TechSat 21, SSTL SNAP-1, and NMP ST5 [5-8]. The satellite formation consists of a group of satellites performing a uniform mission, which is usually difficult to be achieved by a single larger-scale satellite. Wang et al.  proposed an optimal virtual center selection for formation flying maintenance. Zeng et al.  proposed a method to evaluate the safety of formation flying satellites.
Formation stability is a key technology for microsatellite formation communication system. For the problems of formation design, Becerra et al.  investigated the possibility that can obtain periodic or quasiperiodic natural relative motion; Schaub and Alfriend  used J2 theory to design the same order formation, but found that their design produced inconvenient constraint. Kasdin and Koleman  used the epicyclic orbital elements theory to investigate bounded, periodic orbits in the presence of various perturbations. Inalhan et al.  found the analytical expression of the initial conditions which can result in periodic motion based on the classical Tschauner-Hempel equations. Sabol et al.  discussed four kinds of formation designs based on Hill equation and analyzed their stability under various perturbations. The space circular formation and nadir circular formation designed with the passage of time could produce deviation, which cannot satisfy the mission requirements. In this paper, an improved optimizing design method based on genetic algorithm is proposed for space circular orbit formation in order to satisfy the performance requirements of intersatellite communication.
Firstly, the satellite links are designed, the equations of intersatellite link are presented, and a simulation is provided to confirm the correction of the equations. Secondly, the derivation process of an improved genetic algorithm is presented, and a simulation is provided to confirm the feasibility of the proposed algorithm. Finally, a formation microsatellite intersatellite communication system is proposed under the design formation orbit based on the improved genetic algorithm, and we discuss the results of simulation.
2. Links Equations
The basic task for intersatellite communication is transferring information among satellites. So the design of the satellite links must satisfy the communication quality demands. The guide of links equation is given by
[E.sub.b]/[N.sub.0] = [P.sub.t] [L.sub.1] [G.sub.t] [L.sub.s] [L.sub.a] [G.sub.r]/k[T.sub.s]R. (1)
If (1) takes dB as unit, the calculation of the links can be convenient. By this way, if the designers calculate the links parameter, it is only needed to do some addition and reduction operations:
[E.sub.b]/[N.sub.0] = [P.sub.t] + [L.sub.1] + [G.sub.t] + [L.sub.s] + [L.sub.a] + [G.sub.r] + 228.6 - 10 lg [T.sub.s] - 10 lgR, [N.sub.0] (2)
where: [E.sub.b]/[N.sub.0], [L.sub.i], [G.sub.t], [L.sub.s], [L.sub.a], and [G.sub.r] take dB as unit; [P.sub.t] takes dBW as unit; [T.sub.s] takes K as unit; R takes bit/s as unit. BPSK modulation style is chosen here. And the space environment of microsatellite formation-flying is assured. At the beginning, the change of power along with different code rate is simulated. Then, an appropriate transmitted power can be chosen for microsatellite transmitter.
The intersatellite RF links take S frequency. The transmitted frequency is 2 GHz, the error code rate request of links is 10E-5. The full direction antenna is adopted. The communication distance can vary from 100 m to 100 km. The code rate must achieve 1 Mbit/s within 10 km and 10 kbit/s within the range of 10 km to 100 km, respectively. The feedback loss is 3 dB. The loss led by ionosphere is 0.7 dB. The remaining quantity of links is 5 dB. The simulation is given according to the links budget. The simulation results of the change of the transmitted power relative to distance are shown in Figure 1. It can be seen that, when the code rate is 10 kbps and the communication distance varies from 100 m to 100 km, the range of the power change is 250 mw, and that, when the code rate is 1 Mbps, the range of the power is 30 W. If the communication distance is within 10 km, the 250 mw transmitted power is enough. If 30% efficiency of the power output is adopted, the power of transmitter is less than 1 W. So the burden of satellites is not increased.
3. Design Method
The change of intersatellite distance for microsatellites formation-flying can result in the change of space attenuation and satellite channel. Thus, the design of a steady formation order under various absorb impetus can reduce the complexity and the difficulty of intersatellite communication technology. The design of formation-flying order using Hill equation can appear as a biggish excursion under [J.sub.2] absorb impetus. In order to solve this problem, an optimizing design method of space circular formation order based on an improved generic algorithm is proposed in this paper. Genetic algorithm is initially built from problems which could represent the solution set of the beginning process of a population. This process will result in a population-like natural evolution as kid generation population is more adapted to the environment than the previous generations. So after decoding of the previous generations, the optimal individuals can be treated as approximate to the optimal solution . The method can keep formation order steady for a long time under various absorbs impetus. According to the task requirements, the restrict condition of the relative position and velocity can be gained under the Hill coordinate series when the formation order is designed. The optimizing objective function is the apiece average orbit root number of principal and subordinate satellites, with its relative excursion low enough under [J.sub.2] absorb impetus. The relative position and velocity should be transformed to relative orbit root number during the optimizing process. The average orbit factors of satellite under [J.sub.2] absorb impetus, which occur in long-term excursion, include perigee breadth angle, rise point ofintersection equator longitude and aclinic close point angle. The size of excursion item is relative to semilong axes, eccentricity rate, and orbit obliquity:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)
Therefore, the optimizing objective function can be given as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (5)
If a formation includes many satellites, the optimizing objective function is given by
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (6)
In this paper, the formation includes three satellites which consist of one host satellite and two subordinate satellites. The presented simulation is a generic algorithm based optimizing design method for space circular formation order. And the radius of space circular formation order is 1 km. During the optimizing process, the relative position, velocity, and orbit root numbers are:
[[theta].sub.0] = -3.300l040e + 002, [[??].sup.l.sub.0] = [[??].sup.2.sub.0] = -0.8328342. (7)
Table 1 presents the orbit elements of formation flying satellites including several satellites. The results are shown in Figures 2 and 3 (i.e., in Figures 2 and 3:(a) therelativeposition in v-direction; (b) the relative position in y-direction; (c) the relative position in z-direction; and (d) the relative distance of intersatellite), where the former figure represents the first subordinate satellite, and the latter one represents the second subordinate satellite. The simulation time is 14 days, and the number of iterations is 120,960 (10 s per step). The simulation includes [J.sub.2], [J.sub.3], [J.sub.4], and atmosphere resistance absorb impetus. It can be seen from Figures 2 and 3 that the relative position and relative distance excursion of the two satellites on the three axes are very small. These simulation results demonstrate that the proposed optimizing design method can design a long-time steady formation order under various absorb impetus.
Additionally, a comparison among the proposed generic algorithm and other optimizing methods has been taken in the presented work. Similar as what was found by Shi et al. , the presented method based on generic algorithm can obtain accuracy optimizing results with acceptable time requirement compared to the method based on another optimizing algorithm.
4. Simulation Results
4.1. The Parameters of Direct Spread Spectrum Communication System. Supposing that a formation including three microsatellites are flying on a same orbit and performing three dimension solid mapping for terra. Three microsatellites can get three solid images and process the error estimations of solid mapping. In such applications, the available load of each satellite is a camera. Mapping elevation data can be integrated with three satellites in a different perspective to get images of the same area. Compared with a single time-sharing satellite which maps the same area, a formation with three satellites has an advantage that the same area can be measured at the same time. So, the present microsatellites formation can measure not only the elevation of static objects, but also the elevation of dynamic objects. These features can provide a convenience to gain a long measurement baseline and a help to improve the accuracy of the survey and mapping. After the satellite gains the image data, its data transmission also has its own characteristics.
Since every satellite has its own data-transmission system, the mapping data which solid mapping needs can transmit at the same time or in sequence according to the demands. It can decrease pressure on the data transmission and power consumption. Because the formation is fixed, free space loss can be calculated by it. Formation of the microsatellite system using satellite links between the S-band takes center frequency 2 GHz. Each satellite can link with any other satellites of formation. The distance between every two adjacent satellites is 1 km. The satellite antenna is an omnidirectional antenna. The intersatellite link parameters and budgets are shown in Tables 2 and 3.
4.2. Formation of Microsatellite Intersatellite Communication System Model. The model of intersatellite communication system for formation microsatellites model is shown in Figure 4, where data source signal uses 1 MHz frequency random series, and spreading codes use 10 MHz frequency gold code.
The (512, 256) LDPC code is investigated under Gauss channel based on the platform. The code rate is 0.5, and the maximum number of iterations is 50. The bit error rates can be gotten under a different value of [E.sub.b]/[N.sub.0] through changing the [E.sub.b]/[N.sub.0] value ofAWGN channel. The simulation results are shown in Figure 5.
It can be seen from Figure 5 that, if the LDPC code is applied in the system, the bit error rates of the system can be improved obviously. This result can satisfy the requirements of intersatellite communication of a formation including three satellites. It can be seen from Figure 6 that the input wave matches well with the output wave. So, the result can also satisfy the requirements of the communication. Although the usage rate of available frequency is reduced under the transmission, the structure of transmitter can be predigested greatly, and the precision requirement of the fake code generator reduces greatly too.
It can be concluded from the simulation results of the inter-satellite direct sequence spread spectrum communication system that LDPC code can improve the communication performance; the correct capability of (512,256) LDPC code is better than (2,1,7) convolution code; and when the distance is 1 km and the data rate is 1 Mbps, the input wave matches well with the output wave. It can demonstrate that the proposed design system can satisfy the communication requirements of formations including many microsatellites.
The present work not only proposes an optimizing design method of formation order based on an improved generic algorithm, but also provides an intersatellite direct sequence spread spectrum communication system under the designed formation order with encouraging simulation results. So the presented work can provide a feasible theory foundation of the intersatellite communication for formation-flying microsatellites. For engineers, the satellite formation order and intersatellite communication system can be designed effectively and accurately with the convenient method proposed in this paper, which has significant potential benefits in aerospace engineering applications.
Nomenclature a: Semilong axes BPSK: Binary phase shift keying c: Host satellite CDMA: Code division multiple access d: Slave satellite e: Eccentricity rate [E.sub.0]: Received energy per unit ESA: European Space Agency [G.sub.t]: Gain of transmitting antenna ISLs: Intersatellite links [J.sub.2].[J.sub.3].[J.sub.4]: Absorb impetus [K.sub.[OMEGA]], Adding right coefficients [K.sub.[omega]], [K.sub.M]: [L.sub.1] Feed line loss [L.sub.a]: Transmission path loss [L.sub.s]: Propagation loss LDPC: Low density parity check code  LEO: Low earth orbit  M: Aclinic close point angle MEMS: Micro electromechanical systems n: [square root of (u/[n.sup.3])] [N.sub.0]: Noise spectrum density P: a[square root of (1 - [e.sup.2])] [P.sub.t]: Transmitter power R: Data rate [R.sub.e]: A constant equal to 6378.136 km RF: Radio frequency [T.sub.s]: System noise temperature [OEMGA][??]: Rise point of intersection equator longitude and its time derivatives [oemga]: Perigee breadth angle [delta]M, [delta]w, [delta][OMEGA]: The relative orbit root number based on aclinic close point angle, orbit obliquitous, and rise point of intersection equator longitude, respectively.
This work was supported by the National High Technology Research and Development Program of China (no. 2012AA120601), National Natural Science Foundation of China (no. 11202044, no. 11072044), and the Fundamental Research Funds for the Central Universities.
 W. Zhang, Z. M. Deng, and J. S. Li, "Dynamics analysis of satellite cluster non-coplanar orbit transfer based on geometric method," Aircraft Engineering and Aerospace Technology, vol. 84, no. 4, pp. 252-259, 2012.
 N. Sengil and F. O. Edis, "Atomistic simulation of micro-scale adiabatic piston problem," Aircraft Engineering and Aerospace Technology, vol. 81, no. 6, pp. 499-507, 2009.
 A. Ben Amar, M. Faucher, B. Grimbert et al., "Bias dependence of gallium nitride micro-electro-mechanical systems actuation using a two-dimensional electron gas," Applied Physics Express, vol. 5, no. 6, Article ID 067201, 2012.
 J. Shi, L. Miao, M. Ni, and J. Shen, "Optimal robust fault-detection filter for micro-electro-mechanical system-based inertial navigation system/global positioning system," IET Control Theory and Applications, vol. 6, no. 2, pp. 254-260, 2012.
 A. Das and R. Cobb, "TechSat 21: Space missions using collaborating constellation of satellites," in Proceedings of the 12th AIAA/USUAnnual Conference on Small Satellites, vol. 1, pp. 1-5, Logan, Utah, 1998.
 L. Kim, "University nanosatellite distributed satellite capabilities to support TechSat 21," in Proceedings of the 13th AIAA/USU Conference on Small Satellites, vol. 8, pp. 1-9, Logan, Utah, 1999.
 J. J. C. M. Bik, P. N. A. M. Visser, and O. Jennrich, "LISA satellite formation control," Advances in Space Research, vol. 40, no. 1, pp. 25-34, 2007.
 M. Sabatini and R. Bevilacqua, "Special inclinations allowing minimal drift orbits for formation flying satellites," Journal of Guidance, Control, and Dynamics, vol. 31, no. 1, pp. 94-101, 2008.
 J. H. Wang, X. B. Cao, and J. X. Zhang, "Optimal virtual center selection for formation flying maintenance," Aircraft Engineering and Aerospace Technology, vol. 84, no. 4, pp. 260-270, 2012.
 G. Q. Zeng, M. Hu, and J. L. Song, "Collision monitoring and optimal collision avoidance manoeuvre for formation flying satellites," Aircraft Engineering and Aerospace Technology, vol. 84, no. 6, pp. 413-422, 2012.
 V. M. Becerra, J. D. Biggs, S. J. Nasuto, V. F. Ruiz, W. Holderbaum, and D. Izzo, "Using Newton's method to search for quasi-periodic relative satellite motion based on nonlinear hamiltonian models," in Proceedings of the 7th International Conference On Dynamics and Control of Systems and Structures in Space (DCSSS '06), pp. 193-202, Cranfield University Press, London, UK, 2006.
 H. Schaub and K. T. Alfriend, J2 Invariant Relative Orbits For Spacecraft Formations, Nasa Conference Publication, NASA, 1999.
 N. J. Kasdin and E. Koleman, "Bounded, periodic relative motion using canonical epicyclic orbital elements," Advances in the Astronautical Sciences, vol. 120, pp. 1381-1398, 2005.
 G. Inalhan, M. Tillerson, and J. P. How, "Relative dynamics and control of spacecraft formations in eccentric orbits," Journal of Guidance, Control, and Dynamics, vol. 25, no. 1, pp. 48-59, 2002.
 C. Sabol, R. Burns, and C. A. McLaughlin, "Satellite formation flying design and evolution," Journal ofSpacecraft and Rockets, vol. 38, no. 2, pp. 270-278, 2001.
 Q. Zhu, F. Ye, and G.-Z. Feng, "Performance analysis of channel estimation for LDPC-coded OFDM system in multipath fading channel," Acta Electronica Sinica, vol. 23, no. 6, pp. 801-805, 2006.
 T. Shuai, Q. Y. Ren, H. J. Liu, X. W. Liang, and G. Q. Yang, "Hybrid parallel PN code acquisition technique for LEO satellite spread spectrum communication," Journal on Communications, vol. 28, pp. 79-84, 2007.
Guoqiang Wu, (1,2) Yuguang Bai, (1,2) and Zhaowei Sun (3)
(1) State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China
(2) School of Aeronautics and Astronautics, Dalian University of Technology, Dalian 116024, China
(3) Research Center of Satellite Technology, Harbin Institute of Technology, Harbin 150006, China
Correspondence should be addressed to Yuguang Bai; email@example.com
Received 29 April 2013; Accepted 24 July 2013
Academic Editors: D. Missirlis and A. F. B. A. Prado
TABLE 1: Orbit elements of formation flying satellites including several satellites. Orbital Host Difference of Difference of element satellite orbit element 1 orbit element 2 a (Km) 7555 -3.1816E -7 3.2870E -7 e 0.0001 -4.5999E -5 6.0762E -5 i (deg) 48 1.76130E -6 -1.7613E -6 [omega] (deg) 30 3.7782E +1 -1.1876E +1 [OMEGA] (deg) 30 -8.8383E -3 8.8383E -3 M (deg) 0 -3.7776E +1 1.1870E +1 TABLE 2: The parameters of intersatellite link. Multiple methods CDMA Baseband filtering Raised cosine Modulation BPSK Data rate (Mb/s) 1 Center frequency (GHz) 2 Satellite launch power (mW) 250 Satellite [G.sub.t] (dBW) 0 Satellite [G.sub.r] (dBW) 0 Satellite receiving system 300 noise temperature (K) Coding and decoding manner LDPC TABLE 3: The budgets of intersatellite link. Power (mW) 250 Satellite [G.sub.t] (dBW) 0 Satellite [G.sub.r] (dBW) 0 [L.sub.a] (dBW) 0.7 [L.sub.1] (dBW) 3 Noise temperature (K) 300 Link [E.sub.b]/[N.sub.0] (dB) 9.6 Margin (dB) 5 Channel model Gauss
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Wu, Guoqiang; Bai, Yuguang; Sun, Zhaowei|
|Publication:||The Scientific World Journal|
|Date:||Jan 1, 2013|
|Previous Article:||Structural health monitoring of a reinforced concrete building during the severe typhoon Vicente in 2012.|
|Next Article:||Nonlinear analysis of electrocardiography signals for atrial fibrillation.|