<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">IJCNS</journal-id><journal-title-group><journal-title>International Journal of Communications, Network and System Sciences</journal-title></journal-title-group><issn pub-type="epub">1913-3715</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ijcns.2017.105B011</article-id><article-id pub-id-type="publisher-id">IJCNS-76553</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Capacity Analysis and Optimization of Satellite MIMO System
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Haijin</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jianbo</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yuxin</surname><given-names>Cheng</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jianjun</surname><given-names>Wu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Institute of Advanced Communications, EECS, Peking University, Beijing, China</addr-line></aff><pub-date pub-type="epub"><day>26</day><month>05</month><year>2017</year></pub-date><volume>10</volume><issue>05</issue><fpage>116</fpage><lpage>126</lpage><history><date date-type="received"><day>March</day>	<month>15,</month>	<year>2017</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>May</month>	<year>23,</year>	</date><date date-type="accepted"><day>May</day>	<month>26,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   
   Multiple-Input Multiple-Output (MIMO) technology is widely applied in terrestrial wireless communication system, which greatly increases the system capacity. Satellite communication system has many advantages such as wide coverage and strong flexibility. Therefore, how to make a better use of MIMO technology in satellite communication system has become a research hotspot in recent years. The purpose of this paper is to analysis the relationship between satellite MIMO system capacity and parameters of terrestrial antenna such as angle and distance. The parameters of terrestrial antenna were derived and calculated to keep a higher capacity for satellite MIMO system. Numerical analysis of system capacity performance before and after optimization was obtained, which proved the correctness of the theory proposed in this paper. 
  
 
</p></abstract><kwd-group><kwd>MIMO</kwd><kwd> Satellite Communication</kwd><kwd> System Capacity</kwd><kwd> Capacity Analysis</kwd><kwd>  Capacity Optimization</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Multiple-Input Multiple-Output (MIMO) technology is widely applied in modern communication system. It provides users with a higher data rate, better service quality and greater network capacity [<xref ref-type="bibr" rid="scirp.76553-ref1">1</xref>]. Compared to terrestrial communication system, satellite communication system has many advantages such as wide coverage, flexible connection, and little influence by topography, etc. It is widely used in the fields of transoceanic television broadcasting and telephone, communication in remote areas, maritime and aeronautical communication, navigation positioning and military communication. Because of the huge potential of MIMO technology to increase the system capacity, the application of MIMO technology in satellite communication system becomes a hot topic. In recent years, scholars have carried on a large amount of researches about satellite MIMO communication.</p><p>In 2005, P. R. King provided a model of distributed satellite MIMO channel [<xref ref-type="bibr" rid="scirp.76553-ref2">2</xref>], and then in second year, he completed a similar analysis and research on the polarized satellite MIMO channel, which shows the benefits of using MIMO in satellite communication [<xref ref-type="bibr" rid="scirp.76553-ref3">3</xref>]. Sellathurai also proposed a simplified statistical model of polarized satellite MIMO channel [<xref ref-type="bibr" rid="scirp.76553-ref4">4</xref>]. Benjamin Ros proposed a system scheme for interactive satellite MIMO communication using the OFDM system. He compared the performance of satellite MIMO system with multi-antenna SISO system under DVB-SH standard, and the results indicate that MIMO technology is necessary for satellite communication [<xref ref-type="bibr" rid="scirp.76553-ref5">5</xref>].</p><p>Compared with the abundant research results of satellite MIMO, the investigations about optimization of satellite MIMO system capacity are limited. How to improve the system capacity of satellite MIMO system is worthy of research and analysis. In this paper, we provide theoretical calculation of the system capacity of GEO and LEO satellite MIMO system, and then complete the optimization of system capacity. Especially for LEO satellites, we carry out the numerical analysis, which proves the correctness of the theory proposed in this paper.</p></sec><sec id="s2"><title>2. System Model</title><p>Satellite system can be divided into GEO system and non-GEO system. Because the GEO system is relatively geostationary, MIMO system using GEO is more stable. The existing GEO systems are Inmarsat system and Thuraya system, but they both only have three satellites. The beam overlap area is not large enough since the number of satellites is too small. Thus, it is impossible to achieve MIMO system using multiple GEO satellites in most areas on the ground.</p><p>The orbit altitude of LEO satellite system is 700 to 1500 km, moving at a high speed relatively to the earth. Globalstar system and Iridium system are major existing LEO satellite communication systems until now, which respectively have 48 and 66 LEO satellites completing global coverage. The existing LEO satellites are abundant resources, and users can simultaneously connect multiple satellites in most areas of the Earth. Thus, there are enough LEO satellites to build MIMO system.</p><p>As shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>, we consider a complete closed-loop scenario. The whole satellite MIMO system contains three components as follows.</p><p>Terrestrial transmitting terminal: The transmitter is generally a terrestrial station. It is mainly responsible for controlling, scheduling, and connecting the ground network and the user terminal in the entire communication process.</p><p>Communication satellite: Generally, there are two or more satellites, playing a role as relaying and forwarding in the whole communication system. The signal transmitted by the ground station is delivered to the user terminal through transparent forwarding or decode-and-forwarding.</p><p>User terminal: It could be hand-held mobile terminals, or terminals on-board, ship-borne and air-borne moving at a high speed, or even another fixed ground station.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Satellite MIMO system model</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/76553x2.png"/></fig><p>General system includes two links: The uplink is from ground station to the satellite system and downlink is from satellite system to the user terminal. Maybe there also exist inter-satellite links among the satellites to assist communication.</p><p>In <xref ref-type="fig" rid="fig1">Figure 1</xref>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x3.png" xlink:type="simple"/></inline-formula>is the number of antenna of the ground transmitter, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x4.png" xlink:type="simple"/></inline-formula>is the number of antennas of the relay satellite and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x5.png" xlink:type="simple"/></inline-formula> is the number of antenna of the receiver on the ground. In general,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x6.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3"><title>3. Capacity Analysis and Optimization</title><sec id="s3_1"><title>3.1. Capacity Analysis</title><p>In satellite communication systems, the satellite channel is significantly different from the terrestrial channel, because the power of its direct path is much greater than that of the multipath. For example, according to the data from ITU-R M.1225 [<xref ref-type="bibr" rid="scirp.76553-ref6">6</xref>], the power of direct path is about 20 dB greater than the maximum multipath. In fact, because the signal energy in NLOS channel signal is too small, the multipath signal has little effect on the capacity of the system compared with the existence of direct path [<xref ref-type="bibr" rid="scirp.76553-ref7">7</xref>]. Therefore, in this paper we only consider the influence of direct path.</p><p>In consideration of the free space propagation of the direct path through satellite-ground channel, the channel function can be expressed as [<xref ref-type="bibr" rid="scirp.76553-ref8">8</xref>]</p><disp-formula id="scirp.76553-formula87"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x7.png"  xlink:type="simple"/></disp-formula><p>where r denotes the distance between the transmitter and the receiver, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x8.png" xlink:type="simple"/></inline-formula>denotes the signal amplitude attenuation when passing through the channel. Thus for 2*2 MIMO the capacity can be represented as</p><disp-formula id="scirp.76553-formula88"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x9.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x10.png" xlink:type="simple"/></inline-formula> denotes the system SNR removing the channel fading part. The channel matrix H is expressed as</p><disp-formula id="scirp.76553-formula89"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x11.png"  xlink:type="simple"/></disp-formula><p>We establish an unified coordinate system to analyse and calculate the difference between<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x12.png" xlink:type="simple"/></inline-formula>. It is shown in <xref ref-type="fig" rid="fig2">Figure 2</xref> and <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p><p>As shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>, we define the equatorial plane as XOY plane, the polar direction as Z-axis. Besides, OX axis intersects 0 degrees’ east longitude and OY axis intersects 90 degrees’ east longitude. The two red circles in <xref ref-type="fig" rid="fig2">Figure 2</xref> represent the location of ground terminal antennas. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x13.png" xlink:type="simple"/></inline-formula>represents the latitude of ground terminals and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x14.png" xlink:type="simple"/></inline-formula> represents the longitude of ground terminals (East longitude is positive and West longitude is negative). <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x15.png" xlink:type="simple"/></inline-formula>denotes the distance between two antennas on the ground, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x16.png" xlink:type="simple"/></inline-formula> denotes the angle between two antennas and east direction. Then the coordinates of the two antennas on the ground are respectively expressed as</p><disp-formula id="scirp.76553-formula90"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x17.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.76553-formula91"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x18.png"  xlink:type="simple"/></disp-formula><p>Because the coordinates of two satellites are not necessarily related, we analyze a single satellite coordinates respectively. Satellite coordinate is shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p><p>The dotted line circle in <xref ref-type="fig" rid="fig3">Figure 3</xref> is the orbit of satellite. A point is the ascending node (the intersection point of satellite orbit and the equatorial plane) of this orbit,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x19.png" xlink:type="simple"/></inline-formula>. S is the satellite position. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x20.png" xlink:type="simple"/></inline-formula>is the inclination angle of satellite orbit, which is a constant in one satellite system. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x21.png" xlink:type="simple"/></inline-formula>is the phase of satellite orbit ascending node in XOY plane. In the same satellite system, the difference of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x22.png" xlink:type="simple"/></inline-formula> between each orbit is a fixed value. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x23.png" xlink:type="simple"/></inline-formula>changes in the angular velocity of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x24.png" xlink:type="simple"/></inline-formula> with the earth’s rotation, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x25.png" xlink:type="simple"/></inline-formula> is the angular velocity of the earth’s rotation. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x26.png" xlink:type="simple"/></inline-formula>is the phase of satellite in its orbit. We define the phase of the ascending node A as 0˚, so</p><disp-formula id="scirp.76553-formula92"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x27.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x28.png" xlink:type="simple"/></inline-formula>is the angular velocity of the satellite and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x29.png" xlink:type="simple"/></inline-formula> is the initial phase of satellite in its orbit. Then the coordinates of the satellite can be expressed as</p><disp-formula id="scirp.76553-formula93"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x30.png"  xlink:type="simple"/></disp-formula><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Terrestrial terminal coordinate</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/76553x31.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Satellite coordinate</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/76553x32.png"/></fig><p>where</p><disp-formula id="scirp.76553-formula94"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x33.png"  xlink:type="simple"/></disp-formula><p>H is the orbit altitude of the satellite. Two satellites have different<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x34.png" xlink:type="simple"/></inline-formula>. If the two satellites are in different orbit, they also have different<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x35.png" xlink:type="simple"/></inline-formula>. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x36.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x37.png" xlink:type="simple"/></inline-formula> both change with time. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x38.png" xlink:type="simple"/></inline-formula>reflects the satellite motion in its orbit, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x39.png" xlink:type="simple"/></inline-formula> reflects the movement of the earth’s rotation.</p></sec><sec id="s3_2"><title>3.2. Parameter Optimization of GEO System</title><p><xref ref-type="fig" rid="fig4">Figure 4</xref> is a schematic diagram of the GEO satellite MIMO coordinate. It is necessary to discuss the optimization of distance and angle between antennas, so that we can always keep the system capacity in the optimal value.</p><p>First, put Formula (3) into Formula (2). Capacity can be expressed as</p><disp-formula id="scirp.76553-formula95"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x40.png"  xlink:type="simple"/></disp-formula><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> GEO satellite MIMO coordinate</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/76553x41.png"/></fig><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x42.png" xlink:type="simple"/></inline-formula> is only related to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x43.png" xlink:type="simple"/></inline-formula>, and it can be expressed as</p><disp-formula id="scirp.76553-formula96"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x44.png"  xlink:type="simple"/></disp-formula><p>In order to maximize C, the following expression should be guaranteed</p><disp-formula id="scirp.76553-formula97"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x45.png"  xlink:type="simple"/></disp-formula><p>In the GEO system, the satellite is fixed relatively to earth, and the rotation of the earth <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x46.png" xlink:type="simple"/></inline-formula> remains unchanged. Set<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x47.png" xlink:type="simple"/></inline-formula>, and keep the satellite orbit plane parallel to the equatorial plane, as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>. Therefore, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x48.png" xlink:type="simple"/></inline-formula>, put it into Formula (7), we obtain</p><disp-formula id="scirp.76553-formula98"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x49.png"  xlink:type="simple"/></disp-formula><p>Put Formula (12) into Formula (4), we obtain</p><disp-formula id="scirp.76553-formula99"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x50.png"  xlink:type="simple"/></disp-formula><p>Put Formula (12) into Formula (5), we obtain</p><disp-formula id="scirp.76553-formula100"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x51.png"  xlink:type="simple"/></disp-formula><p>Let</p><disp-formula id="scirp.76553-formula101"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x52.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.76553-formula102"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x53.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x54.png" xlink:type="simple"/></inline-formula>can be expressed as</p><disp-formula id="scirp.76553-formula103"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x55.png"  xlink:type="simple"/></disp-formula><p>Then by exploiting the small aqueous formula<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x56.png" xlink:type="simple"/></inline-formula>, we obtain</p><disp-formula id="scirp.76553-formula104"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x57.png"  xlink:type="simple"/></disp-formula><p>Put Formulas (14) and (18) into (11), it can be obtained that when place angle <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x58.png" xlink:type="simple"/></inline-formula> is fixed, the optimal value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x59.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.76553-formula105"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x60.png"  xlink:type="simple"/></disp-formula><p>It can be seen that the optimal value of the terrestrial antenna is periodic variation. In the same way, let</p><disp-formula id="scirp.76553-formula106"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x61.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.76553-formula107"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x62.png"  xlink:type="simple"/></disp-formula><p>Put Formula (20) and (21) into (19), It can be obtained that when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x63.png" xlink:type="simple"/></inline-formula> is fixed, the optimal value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x64.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.76553-formula108"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x65.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3_3"><title>3.3. Parameter Optimization of LEO System</title><p>In the LEO satellite system, there exists an angle between the satellite orbit plane and the equatorial plane, in addition, the angular velocity of satellites and the earth is different, so <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x66.png" xlink:type="simple"/></inline-formula> will change with time. Thus, it is necessary to discuss the optimization of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x67.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x68.png" xlink:type="simple"/></inline-formula> when time changes, so that we can always keep the system capacity in the optimal value.</p><p>First, we calculate the distance <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x69.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x70.png" xlink:type="simple"/></inline-formula> between two satellites and antenna No. 1 on the ground.</p><p>Let</p><disp-formula id="scirp.76553-formula109"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x71.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.76553-formula110"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x72.png"  xlink:type="simple"/></disp-formula><p>Make use of Formulas (4) and (7). After simplification, we obtain</p><disp-formula id="scirp.76553-formula111"><label>(25)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x73.png"  xlink:type="simple"/></disp-formula><p>Set the variable<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x74.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x75.png" xlink:type="simple"/></inline-formula>can be represented as</p><disp-formula id="scirp.76553-formula112"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x76.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.76553-formula113"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x77.png"  xlink:type="simple"/></disp-formula><p>Then by exploiting the small aqueous formula<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x78.png" xlink:type="simple"/></inline-formula>, we obtain</p><disp-formula id="scirp.76553-formula114"><label>(28)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x79.png"  xlink:type="simple"/></disp-formula><p>When place angle <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x80.png" xlink:type="simple"/></inline-formula> is fixed, the optimal value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x81.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.76553-formula115"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x82.png"  xlink:type="simple"/></disp-formula><p>In this formula, it is obvious that the optimum value of ground antennas changes periodically. Similarly, let</p><disp-formula id="scirp.76553-formula116"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x83.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.76553-formula117"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x84.png"  xlink:type="simple"/></disp-formula><p>When <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x85.png" xlink:type="simple"/></inline-formula> is fixed, it can be obtained that the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x86.png" xlink:type="simple"/></inline-formula> corresponding largest system capacity can be expressed as</p><disp-formula id="scirp.76553-formula118"><label>(32)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/76553x87.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s4"><title>4. Numerical Analysis Results</title><p>Through the analysis above, it can be seen that optimum channel capacity of GEO and LEO satellite MIMO system is related to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x88.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x89.png" xlink:type="simple"/></inline-formula>. Then numerical analysis is carried out to prove the correctness of the theory proposed in chapter III. Here we set LEO satellite system as an example of numerical calculation. For GEO satellite, results can be used in the same way.</p><p>In order to calculate LEO satellite MIMO system through numerical analysis, orbit parameters of Globalstar system and Iridium system are needed. We also need to determine the system SNR <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x90.png" xlink:type="simple"/></inline-formula> besides of the orbit parameters of the two systems. SNR depends on satellite transmission power<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x91.png" xlink:type="simple"/></inline-formula>, satellite antenna gain<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x91.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x92.png" xlink:type="simple"/></inline-formula>, terrestrial antenna gain <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x91.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x93.png" xlink:type="simple"/></inline-formula> and received noise <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x91.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x94.png" xlink:type="simple"/></inline-formula> of terrestrial antennas. These parameters are summarized in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>System capacity of satellite MIMO system is closely related to the arrangement of ground antennas. The system capacity can be maintained at a optimal value by using tracking antenna chrtrolling <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x95.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x96.png" xlink:type="simple"/></inline-formula>. <xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref> show the performance of the tracking antenna using Formula (25). It can be seen that</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Main parameters in numerical analysis</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Satellite transmit power</th><th align="center" valign="middle" >52 dBm</th></tr></thead><tr><td align="center" valign="middle" >Satellite antenna gain</td><td align="center" valign="middle" >48 dBi</td></tr><tr><td align="center" valign="middle" >Ground antenna gain</td><td align="center" valign="middle" >22 dBi</td></tr><tr><td align="center" valign="middle" >Received noise of terrestrial antenna</td><td align="center" valign="middle" >−90 dB</td></tr></tbody></table></table-wrap><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Tracking antenna performance (d = 0.2 m)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/76553x97.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Tracking antenna performance (d = 2 m)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/76553x98.png"/></fig><p>the system capacity of tracking antennas is significantly improved compared with the fixed antennas in both systems. The tracking antenna makes the system maintain a high and relatively fixed system capacity.</p><p>From <xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref>, we can see the different capacity performance when the antenna spacing <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/76553x99.png" xlink:type="simple"/></inline-formula> is 0.2 m and 2 m, when the antenna distance increases, we can find that the system capacity of the fixed antenna becomes more unstable with big fluctuation. It shows that tracking antenna is more significant for large-scale antenna system with larger antenna spacing.</p></sec><sec id="s5"><title>5. Conclusion</title><p>In this paper, we have analyzed the relationship between the capacity performance of satellite MIMO system and placement parameters of terrestrial antennas such the angle and distance. The numerical calculation verifies the possibility to optimize the parameters. Later, through formula derivation we obtain the time variation function of the optimal placement angle and distance of the ground antenna. Finally, numerical analysis is carried out according to the initial position of satellite in the practical system. The results demonstrate that the system capacity can be significantly improved when the antenna is traced using the theoretical analysis results.</p></sec><sec id="s6"><title>Acknowledgements</title><p>This paper is supported by the National Natural Science Foundation of China # 61371073.</p></sec><sec id="s7"><title>Cite this paper</title><p>Li, H.J., Li, J.B., Cheng, Y.X. and Wu, J.J. (2017) Capacity Analysis and Optimization of Satellite MIMO System. Int. J. Communications, Network and System Sciences, 10, 116-126. https://doi.org/10.4236/ijcns.2017.105B011</p></sec></body><back><ref-list><title>References</title><ref id="scirp.76553-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">G.Grillo, Gesbert, D., Shafi, M., Shiu, D.-S., et al. (2003) From Theory to Practice: An Overview of MIMO Space-Time Coded Wireless Systems. IEEE Journal on Selected Areas in Communications, 21, 281-302.  
https://doi.org/10.1109/JSAC.2003.809458</mixed-citation></ref><ref id="scirp.76553-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">King, P.R., Evans, B.G. and Stavrou, S. (2005) Physical-Statistical Model for the Land Mobile-Satellite Channel Applied to Satellite/HAP-MIMO. 11th European Proceedings of the Wireless Conference 2005-Next Generation Wireless and Mobile Communications and Services (European Wireless).</mixed-citation></ref><ref id="scirp.76553-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">King, P.R. and Stavrou, S. (2006) Capacity Improvement for a Land Mobile Single Satellite MIMO System. IEEE Antennas and Wireless Propagation Letters, 5. 
https://doi.org/10.1109/lawp.2006.872439</mixed-citation></ref><ref id="scirp.76553-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Sellathurai, M., Guinand, P. and Lodge, J. (2006) Space-Time Coding in Mobile Satellite Communications Using Dual-Polarized Channels. IEEE Transactions on Vehicular Technology, 55, 188-99. https://doi.org/10.1109/TVT.2005.861195</mixed-citation></ref><ref id="scirp.76553-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Ros, B., Lacoste, F., Scot, G., et al. (2009) Increasing Reliability of an Interactive Mobile Satellite Telecommunication System Using Diversity and MIMO Schemes. 2009 International Conference on Proceedings of the Wireless Communications &amp; Signal Processing. https://doi.org/10.1109/WCSP.2009.5371395</mixed-citation></ref><ref id="scirp.76553-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Recommendation, I. (1997) Guidelines for Evaluation of Radio Transmission Technologies for IMT-2000. International Telecommunication Union.</mixed-citation></ref><ref id="scirp.76553-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Knopp, A., Schwarz, R.T., Hofmann, C.A., et al. (2007) Measurements on the Impact of Sparse Multipath Components on the LOS MIMO Channel Capacity. 4th International Symposium on Proceedings of the Wireless Communication Systems 
https://doi.org/10.1109/iswcs.2007.4392301</mixed-citation></ref><ref id="scirp.76553-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Oh, C.-I., Choi, S.-H., Chang, D.-I., et al. (2006) Analysis of the Rain Fading Channel and the System Applying MIMO. International Symposium on Proceedings of the Communications and Information Tech-nologies.  
https://doi.org/10.1109/iscit.2006.340001</mixed-citation></ref></ref-list></back></article>