<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Space Technol.</journal-id>
<journal-title>Frontiers in Space Technologies</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Space Technol.</abbrev-journal-title>
<issn pub-type="epub">2673-5075</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">755520</article-id>
<article-id pub-id-type="doi">10.3389/frspt.2021.755520</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Space Technologies</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Multi-Satellite Cooperative Beamforming ALOHA for LEO Satellite&#x20;IoT Networks</article-title>
<alt-title alt-title-type="left-running-head">Liu et&#x20;al.</alt-title>
<alt-title alt-title-type="right-running-head">Cooperative Beamforming ALOHA</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Rui</given-names>
</name>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Hong</surname>
<given-names>Tao</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1299646/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ding</surname>
<given-names>Xiaojin</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Yunfeng</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Gengxin</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/978340/overview"/>
</contrib>
</contrib-group>
<aff>Nanjing University of Posts and Telecommunications, <addr-line>Nanjing</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1386340/overview">Chunguo Li</ext-link>, Southeast University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1354916/overview">Qiong Wu</ext-link>, Jiangnan University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1439144/overview">Gang Wang</ext-link>, Southeast University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1029050/overview">Kanglian Zhao</ext-link>, Nanjing University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Tao Hong, <email>hongt@njupt.edu.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Aerial and Space Networks, a section of the journal Frontiers in Space Technologies</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>12</day>
<month>10</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>2</volume>
<elocation-id>755520</elocation-id>
<history>
<date date-type="received">
<day>09</day>
<month>08</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>20</day>
<month>09</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2021 Liu, Hong, Ding, Wang and Zhang.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Liu, Hong, Ding, Wang and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these&#x20;terms.</p>
</license>
</permissions>
<abstract>
<p>In this paper, we proposed a cooperative beamforming ALOHA (CBA) scheme based on linearly constrained minimum variance criterion for low Earth orbit satellite (LEO) IoT networks to solve the problem of &#x2018;deadlock&#x2019; in multi-satellite scenario. In multi-satellite overlapping coverage areas, packets can be received by multiple satellite receivers by sending them only once, which forms the concept of spatial diversity. The cooperative beamforming collision resolution technique combined with successive interference cancellation scheme is design to efficiently resolve packet collision by iteration way at the gateway station. The performance of cooperative beamforming ALOHA scheme is evaluated via mathematical analysis and simulations. Simulation results show that the proposed CBA scheme can effectively solve the problem of &#x2018;deadlock&#x2019; and improve the performances of random access compared with benchmark problems.</p>
</abstract>
<kwd-group>
<kwd>internet of things</kwd>
<kwd>satellite IOT</kwd>
<kwd>random access</kwd>
<kwd>cooperative beamforming</kwd>
<kwd>collision resolution</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>Internet of Things (IoT) is one of the main research directions in the future 5G wireless network applications such as industrial assets monitoring, smart city, logistics and environmental data sensing (<xref ref-type="bibr" rid="B4">Stefano et&#x20;al., 2018</xref>). Nowadays, the information transmission between IoT terminals and the gateway station mainly depends on terrestrial wireless network, Wi-Fi and other local area networks, which need to construct numerous ground base stations to support massive IoT terminals. However, there will exist a conflict between service capability and network QoS requirement when IoT is applied in mountainous area, ocean and desert, because ground base stations are not suitable to construct and maintain in these remote areas (<xref ref-type="bibr" rid="B3">Wei-Che et&#x20;al., 2019</xref>). Due to the nature of wide coverage, satellite communication systems have become an important supplement network to meet the requirement of the global coverage for 5G IoT applications (<xref ref-type="bibr" rid="B7">Mauro et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B18">Zhicheng et&#x20;al., 2017</xref>; <xref ref-type="bibr" rid="B23">Tingting et&#x20;al., 2019</xref>). Therefore, low Earth orbit (LEO) satellite based IoT networks has become a hot research topic in IoT&#x20;area.</p>
<p>Compared with traditional wireless broadband access, the IoT service has the characteristics of massive connections, strong burstiness, and short packets. Therefore, RA is an effective access scheme from point view of channel utilization compared with fixed allocation and on-demand allocation multiple access methods. The RA scheme is to share channel resources in a competitive way. Multi-terminals occupying the same time-frequency resource will cause packets collision, which results in increasing the system packet loss ratio (PLR) and decreasing throughput of system. The channel utilization rate of the pure ALOHA (P-ALOHA) protocol is only 18.4% (<xref ref-type="bibr" rid="B1">Abramson, 1970</xref>). On this basis, the time slot ALOHA (S-ALOHA) protocol is proposed to improve the throughput for satellite networks by dividing the transmission time into several time slots (<xref ref-type="bibr" rid="B19">Roberts, 1975</xref>). The terminal can only send packets at the beginning of each time slot. In this way, the channel utilization rate improves to 36.7% by reducing the collision window of packets. However, S-ALOHA protocol requires network synchronization, which increases the complexity of the system. In (<xref ref-type="bibr" rid="B2">Enrico et&#x20;al., 2007</xref>), E. Casini proposed a contention resolution diversity slotted ALOHA (CRDSA) to resolve contention of the packet collision, which combines multiple copies in the time or frequency domain with SIC technology. The collision packets have a chance to be recovered by this scheme. To improve the performance of CRDSA, R. De Gaudenzi presented a CRDSA&#x2b;&#x2b; scheme by transmitting more than two replicas of packet in (<xref ref-type="bibr" rid="B6">Riccardo and Oscar, 2007</xref>). This proposed RA scheme reveals that the increased replicas allow better resolve collision performance by SIC processing, although the more replicas can increase the probability of collision. G. Liva optimized the CRDSA scheme in (<xref ref-type="bibr" rid="B15">Gianluigi, 2010</xref>), called as irregular repetition slotted ALOHA (IRSA). In IRSA, the terminal generates a different number of copies according to the optimized allocation rules, which makes the throughput higher than the CRDSA protocol. However, terminals need to transmit signal copies in the time or frequency domain for these CRDSA based schemes, which not only consume the more power of terminals but also trigger &#x2018;avalanche effect&#x2019; in high load condition of IoT networks. These time-frequency diversity improved RA schemes at the terminals fall far short of the initial design goal of low power consumption and massive access.</p>
<p>Inspired by transmitting signal copies in the time or frequency domain, massive multiple input multiple output (MIMO) is an effective spatial diversity measure to improve the system capacity in terrestrial wireless networks (<xref ref-type="bibr" rid="B21">Mingjin et&#x20;al., 2019</xref>). However, compared with rich-scattering channel environment in terrestrial wireless network, light-of-sight (LOS) channel plays an important role with characteristic of non-rich-scattering channel in satellite communication channel. Therefore, the technique of massive MIMO applied in satellite communication system is difficult to achieve ideal diversity performance. To improve the system capacity in LOS channel, X.Q. Gao&#x2019;s research team proposed a beam division multiple access (BDMA) transmission technology to reuse time-frequency resources in orthogonal beam domain for massive MIMO communications (<xref ref-type="bibr" rid="B20">Chen et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B24">Xin et&#x20;al., 2017</xref>). On this basis, <ext-link ext-link-type="uri" xlink:href="https://elkssl0060f01354d8aa7bea8c889b638382aflib.webvpn.njupt.edu.cn:4443/author/37086290858">R.D. Jia</ext-link> combine BDMA with power domain non-orthogonal multiple access (NOMA) to increase the number of concurrent access terminals in the same time-frequency resource for the terrestrial IoT networks in (<xref ref-type="bibr" rid="B12">Rundong et&#x20;al., 2019</xref>). However, the spatial resolving ability of antenna array in LEO satellite is much less than the same aperture of antenna array in terrestrial wireless networks because of the long distance between LEO satellite and IoT terminals, such as 150&#xa0;km covering radius for antenna array with meter aperture in L band. J Ding proposed a novel machine learning based framework solution to address the preamble collision problem by joint processing distributed massive MIMO signal in terrestrial wireless networks in (<xref ref-type="bibr" rid="B9">Jie et&#x20;al., 2021</xref>). This joint processing algorithm clustered the neighboring access points (APs) of collided RA terminals and organize each AP cluster to decode the received data individually. Motivated by this distributed cooperative scheme in terrestrial wireless networks, based on the precondition of synchronous mode, B. Zhao proposed a novel multi-satellite cooperative random access (MSC-RA) scheme combine with CRDSA, CRDSA&#x2b;&#x2b;, and IRSA, respectively, to improve the RA performances in (<xref ref-type="bibr" rid="B26">Bo et&#x20;al., 2019</xref>), which utilized multi-satellite overlapping coverage condition for the IoT terminals, such as OneWeb, Starlink and Telesat LEO satellite constellation development plan (<xref ref-type="bibr" rid="B17">Nils et&#x20;al., 2021</xref>; <xref ref-type="bibr" rid="B10">Tong and Venkata, 2021</xref>). P. Li proposed an Asynchronous Cooperative ALOHA (ACA) scheme to form a spatial replica of the packet by multiple receiving satellites in (<xref ref-type="bibr" rid="B13">Pengxu et&#x20;al., 2017</xref>). For this distributed spatial diversity scheme, the non-collision packet in a certain satellite is used to resolve contention in other satellites by SIC technology, which improves the system channel utilization rate. The techniques in ACA schemes that utilize multiple distributed nodes to receive signal separately and then process signal centrally can be categorized as coordinated multipoint (CoMP) transmission techniques for uplink. For example, Chen Zhang in the (<xref ref-type="bibr" rid="B25">Chen et&#x20;al., 2021</xref>) used 3-dimensionally placed unmanned aerial vehicles as receiving nodes to process uplink data information under complex geographical conditions. Liang Liu used UAVs as CoMP nodes for interference cancellation in (<xref ref-type="bibr" rid="B14">Liang et&#x20;al., 2019</xref>) to enhance the sum-rate capacity of the cellular uplink network. Compared with the time-frequency replicas schemes, the terminals do not need to send packet replicas, which not only saves the power of terminals but also alleviates network load at the expense of LEO multi-satellite spatial domain resources. However, the author considers that at least one of the receive packets by multiple satellites must not collide with other packets in a certain satellite. In this way, the SIC technology can be used to iteratively recover more collision packets. For this SIC cooperative scheme among multi-satellite, a &#x2018;deadlock&#x2019; phenomenon will occur when the receive packet collide in all satellites, which will degrade the advantages of multi-satellite spatial diversity.</p>
<p>To solve this &#x2018;deadlock&#x2019; phenomenon of multi-satellite spatial diversity scheme, we propose a cooperative beamforming ALOHA (CBA) scheme based on linearly constrained minimum variance (LCMV) criterion for LEO satellite IoT networks. This scheme combines multi-satellite receive diversity of RA packets with cooperative beamforming technique to resolve contention of the packet collision at the gateway station. The power of the collision packet will be suppressed in the spatial domain by cooperative beamforming technique compared with SIC cooperative scheme for ACA scheme in (<xref ref-type="bibr" rid="B13">Pengxu et&#x20;al., 2017</xref>), cluster cooperative scheme in (<xref ref-type="bibr" rid="B9">Jie et&#x20;al., 2021</xref>) and MSC MIMO detection in (<xref ref-type="bibr" rid="B26">Bo et&#x20;al., 2019</xref>). The proposed CBA technique combined with SIC scheme is also designed to efficiently resolve packet collision for multi-satellite receive signal by iteration way. In this way, the multiple access interference (MAI) in multi-satellite overlapping coverage areas is changed into the available spatial resource for LEO satellite IoT networks. Furthermore, we derived the theoretical performance of proposed CBA scheme. Simulation results show that the proposed CBA scheme can effectively solve the problem of &#x2018;deadlock&#x2019; and improve the performances of RA compared with benchmark problems.</p>
<p>The rest of this paper is organized as follows:<list list-type="simple">
<list-item>
<p>&#x2022; In <italic>System Model</italic>, a RA system scenario with multi-satellite coverage is established.</p>
</list-item>
<list-item>
<p>&#x2022; In <italic>Multi-Satellite Cooperative Beamforming ALOHA Scheme</italic>, the process of the CBA scheme is proposed and analyzed; the performance expression of the entire system is derived.</p>
</list-item>
<list-item>
<p>&#x2022; In <italic>Simulation</italic>, the simulation parameters, and results of the CBA scheme are given, and the performance comparison with the CRDSA and ACA scheme under the same conditions is carried&#x20;out.</p>
</list-item>
<list-item>
<p>&#x2022; In <italic>Conclusion</italic>, the conclusion of the paper is presented.</p>
</list-item>
</list>
</p>
</sec>
<sec id="s2">
<title>System Model</title>
<p>
<xref ref-type="fig" rid="F1">Figure&#x20;1</xref> shows a typical multi-satellite coverage scenario in LEO satellite IoT networks with terminals, LEO multi-satellite and gateway station. The area bounded by the black line denotes the multi-satellite overlapping coverage areas. Thus, the packets received by multi-satellite are forwarded to the gateway station transparently with independent backhaul channels which are assumed to be error free. Packets are performed at the gateway station. Without loss of generality, we supposed that the IoT network load follows Poisson distribution with parameter <inline-formula id="inf1">
<mml:math id="m1">
<mml:mtext>&#x3bb;</mml:mtext>
</mml:math>
</inline-formula> in one packet transmission duration and each transmitted packet length is fixed with the same duration <inline-formula id="inf2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. Thus, the normalized load <inline-formula id="inf3">
<mml:math id="m3">
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi mathvariant="normal">lo</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">g</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> (bit/symbol) denotes the amount of data received by the gateway station per unit time, where <inline-formula id="inf4">
<mml:math id="m4">
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:math>
</inline-formula> is the signal modulation order and <inline-formula id="inf5">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the total number of cooperative satellites. The normalized throughput <inline-formula id="inf6">
<mml:math id="m6">
<mml:mrow>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">PLR</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> represents the amount of data can be correctly demodulated by the gateway station per unit time, where <inline-formula id="inf7">
<mml:math id="m7">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PLR</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the PLR under a specific network load. Due to the difference distance between each terminal and multi-satellite, a specific packet is transmitted once without any copy and it can be received by multi-satellite with different propagation delays in the receiving window <inline-formula id="inf8">
<mml:math id="m8">
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula id="inf9">
<mml:math id="m9">
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:math>
</inline-formula> denotes the minimum propagation delay and <inline-formula id="inf10">
<mml:math id="m10">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the maximum difference of propagation delays. For the typical LEO satellite constellation, such as OneWeb, Starlink and Telesat, <inline-formula id="inf11">
<mml:math id="m11">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is about 10&#x2013;100&#xa0;ms, while the typical packet length of IoT networks <inline-formula id="inf12">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is approximately equals to 1&#xa0;ms. Therefore, we consider <inline-formula id="inf13">
<mml:math id="m13">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> in the following sections.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Multi-satellite coverage scenario in LEO satellite IoT networks.</p>
</caption>
<graphic xlink:href="frspt-02-755520-g001.tif"/>
</fig>
</sec>
<sec id="s3">
<title>Multi-Satellite Cooperative Beamforming ALOHA Scheme</title>
<p>To describe proposed CBA scheme in detail, the packet and receiving window status is represented in <italic>Packet and receiving window structure</italic>, the processing flow of proposed Multi-satellite CBA scheme is introduced in <italic>Packet processing flow of proposed CBA scheme</italic>, and the throughput and PLR performance of proposed CBA scheme is analyzed in <italic>Performance Analysis</italic>.</p>
<sec id="s3-1">
<title>Packet and Receiving Window Structure</title>
<p>
<xref ref-type="fig" rid="F2">Figure&#x20;2A</xref> shows a collision-free packet structure, which contains two parts: one is preamble and the other is data. The preamble part is used to estimate the channel state information (CSI) of the packet, whose design method can reference in (<xref ref-type="bibr" rid="B5">Riccardo et&#x20;al., 2014</xref>) and (<xref ref-type="bibr" rid="B16">Andrea et&#x20;al., 2015</xref>). Due to the fact that a specific packet can be received by multi-satellite, we consider that these receive packets as the spatial copies are processed by sliding window at the ground gateway station. The process mode of sliding window for each LEO receive signal is schematically shown in <xref ref-type="fig" rid="F2">Figure&#x20;2B</xref>, where <inline-formula id="inf14">
<mml:math id="m14">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is sliding step in time order. Generally, the window length <inline-formula id="inf15">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> should be three times the length of the physical layer packets and the sliding step <inline-formula id="inf16">
<mml:math id="m16">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is taken as half of the window length.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>
<bold>(A)</bold> Single packet structure. <bold>(B)</bold>. Structure of sliding window and receiving window.</p>
</caption>
<graphic xlink:href="frspt-02-755520-g002.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>Packet Processing Flow of Proposed CBA Scheme</title>
<p>RA scheme does not perform resource allocation in the time-frequency domain, which leads to packet collision at the receiver. At low load area <inline-formula id="inf17">
<mml:math id="m17">
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, packet collision has little impact on the system performance due to the fallback transmission mechanism (FTM). However, when the load spikes in a short period of time <inline-formula id="inf18">
<mml:math id="m18">
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, a large number of packet collision are generated at the receiver causing the FTM performance and the normalized throughput <inline-formula id="inf19">
<mml:math id="m19">
<mml:mi mathvariant="bold-italic">T</mml:mi>
</mml:math>
</inline-formula> degradation. For the multi-satellite spatial diversity scheme, packet collision happens in all the receiving satellites under the medium and high load areas, which lead to the SIC cooperative scheme failure, called as &#x2018;deadlock&#x2019;. Therefore, we proposed a CBA scheme to address the problem of low throughput performance under medium and high load areas in multi-satellite spatial diversity scheme. Packet collision can be resolved by cooperative beamforming technique when packet collision status satisfied certain requirement. The principle of this requirement is that the preamble can be identified for the desired demodulation packet in different packet collision statuses, and the other collision packets are regarded as interference signals. The gateway station utilizes the cooperative beamforming algorithm by multi-satellite receive signal to constructively enhance the desired demodulation packet, simultaneously, suppress the other collision packets. In this way, the signal to interference plus noise power ratio (SINR) of the desired demodulation packet is greater than the demodulation threshold to achieve the demodulation function of collision packet. The proposed CBA scheme is divided into four steps as follows.</p>
<sec id="s3-2-1">
<title>Packet Status Detection</title>
<p>According to the collision situation of preamble, the packet status is divided into two categories: one is preamble collision, the other is packet date part collision while packet preamble detection. If a packet preamble cannot be detected, it is considered as a collision packet and its signal is temporarily reserved for subsequent SIC operations. Conversely, if the packet preamble can be detected, the gateway station uses sliding correlation detection to locate the copies in other receive windows, and its copy&#x2019;s location are used for subsequent copies SIC operation. Then, the correctness of demodulation packet can be determined by the cyclic redundancy check (CRC). <xref ref-type="fig" rid="F3">Figure&#x20;3</xref> shows the packet status detection&#x20;flow.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Packet status detection flow&#x20;chart.</p>
</caption>
<graphic xlink:href="frspt-02-755520-g003.tif"/>
</fig>
</sec>
<sec id="s3-2-2">
<title>Packet Collision Type Classification</title>
<p>According to the CRC result in step (1), the packet collision type can be classified into three cases: 1) The CRC results of this packet and its spatial copies are correct in all receiving windows. This means that this packet can be correctly demodulated. Therefore, the gateway station continues to process the next packet signal in time order and return to step (1); 2) The CRC result is correct at least in one receiving window. For this case, the gateway station utilizes the demodulation information to reconstruct the packet with the corresponding CSI and go to step (4); 3) The CRC results of this packet and its spatial copies are incorrect in all receiving windows, the cooperative beamforming algorithm is performed in step&#x20;(3).</p>
</sec>
<sec id="s3-2-3">
<title>Cooperative Beamforming Algorithm Based on LCMV Criterion</title>
<p>According to the reference (<xref ref-type="bibr" rid="B22">Zhiqiang et&#x20;al., 2019</xref>), the gateway station makes use of multiple CSI to estimate the angle of arrival (AoA) of the desired demodulation packet. Therefore, a cooperative beamforming algorithm based on the LCMV criterion can be used to suppress the interference signal power in the undesired directions. The optimization problem is formulated as follows:<disp-formula id="e1">
<mml:math id="m20">
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="normal">min</mml:mi>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">w</mml:mi>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>R</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">w</mml:mi>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where <inline-formula id="inf20">
<mml:math id="m21">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold">t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>,<inline-formula id="inf21">
<mml:math id="m22">
<mml:mrow>
<mml:mi mathvariant="bold-italic"> x</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> denotes the receive signal, <inline-formula id="inf22">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the steering vector <inline-formula id="inf23">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the AoA angle. This optimization problem can be solved by the Lagrange multiplier method:<disp-formula id="e2">
<mml:math id="m25">
<mml:mrow>
<mml:mi mathvariant="bold-italic">L</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3bb;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mi mathvariant="bold-italic">w</mml:mi>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi mathvariant="bold-italic">w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3bb;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">w</mml:mi>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>The optimal weight vector can be written as:<disp-formula id="e3">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">w</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">opt</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mo mathvariant="bold">&#x5e;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mo mathvariant="bold">&#x5e;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>In summary, the cooperative beamforming method can enhance the desired signal <inline-formula id="inf24">
<mml:math id="m27">
<mml:mrow>
<mml:mi mathvariant="bold-italic">SINR</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>value. If <inline-formula id="inf25">
<mml:math id="m28">
<mml:mrow>
<mml:mi mathvariant="bold-italic">SINR</mml:mi>
<mml:mo>&#x2265;</mml:mo>
<mml:mi mathvariant="bold-italic">SIN</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the desired packet can be demodulated. If <inline-formula id="inf26">
<mml:math id="m29">
<mml:mrow>
<mml:mi mathvariant="bold-italic">SINR</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mi mathvariant="bold-italic">SIN</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the desired packet cannot be demodulated. The packet bit information is retained and return to step 1. <xref ref-type="fig" rid="F4">Figure&#x20;4</xref> shows a schematic diagram of applying beamforming technique to resolve the collision problem.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Schematic diagram of collision resolution using beamforming.</p>
</caption>
<graphic xlink:href="frspt-02-755520-g004.tif"/>
</fig>
</sec>
<sec id="s3-2-4">
<title>Successive Interference Cancellation</title>
<p>If a packet is successfully demodulated in step (2) or step (3), SIC technology can be used to resolve collision packets in certain receiving windows. The principle is that the reconstructed packet is identical to the original packet, and its colliding spatial copies can be eliminated by using the reconstructed packet, thus restoring the packets colliding with its spatial copies to a collision-free state. When this collision is eliminated, a new collision-free packet is generated, which means that the new preamble can be detected. Then the gateway station will restart the execution of step (1) for a new iterative&#x20;loop.</p>
<p>
<xref ref-type="fig" rid="F5">Figure&#x20;5</xref> shows a schematic diagram of the SIC. For packet 1, although it collides with packet 2 at receiving window B, the packet 1 signal at receiving window A can be used for cooperative reconstruction to eliminate the packet 1 in receiving window B. When packet 1 in receiving window B is eliminated, packet 2 is restored to a collision-free state. Then in the next iteration, packet 2 of receiving window B can be used to eliminate the collision packet 2 signal in receiving window A. And so on, after several iteration loops, packets 1, 2, 3 and 4 can be successfully demodulated. <xref ref-type="fig" rid="F6">Figure&#x20;6</xref> shows the flow chart of the proposed CBA algorithm.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Schematic diagram of the SIC.</p>
</caption>
<graphic xlink:href="frspt-02-755520-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>The flow chart of the CBA algorithm.</p>
</caption>
<graphic xlink:href="frspt-02-755520-g006.tif"/>
</fig>
<p>Compared to traditional ALOHA schemes, improved ALOHA schemes are using complex algorithms in exchange for higher system throughput. Since both use a three-level cyclic structure for processing packets, the time complexity of both the ACA scheme algorithm and the CBA scheme algorithm is approximately <inline-formula id="inf27">
<mml:math id="m30">
<mml:mrow>
<mml:mi mathvariant="bold-italic">O</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</sec>
</sec>
<sec id="s3-3">
<title>Performance Analysis</title>
<p>In this subsection, the probability of packet arrival in a receiving window during the time-period <inline-formula id="inf28">
<mml:math id="m31">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> follows a Poisson traffic process. If the length of a receiving window is <inline-formula id="inf29">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3bb;</mml:mi>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula id="inf30">
<mml:math id="m33">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> denotes the maximum delay difference of packets arriving at different satellites. The probability of <inline-formula id="inf31">
<mml:math id="m34">
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:math>
</inline-formula> packets arriving in a window is written as:<disp-formula id="e4">
<mml:math id="m35">
<mml:mrow>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msup>
<mml:mi mathvariant="normal">exp</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>!</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>Thus, the probability of collision packets number <inline-formula id="inf32">
<mml:math id="m36">
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:math>
</inline-formula> in the arriving packets number <inline-formula id="inf33">
<mml:math id="m37">
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:math>
</inline-formula> follows the binomial distribution as:<disp-formula id="e5">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">x</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo> &#x2217;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic"> p</mml:mi>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:msubsup>
<mml:mo> &#x2217; </mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">x</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where <inline-formula id="inf34">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the probability of collision in different parts of the packet. It can be divided into two cases: one is preamble collision <inline-formula id="inf35">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>; the other is data part collision <inline-formula id="inf36">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. We assume that the length of preamble accounts for <inline-formula id="inf37">
<mml:math id="m42">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mi>&#x3b5;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3b5;</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> of the total packet length <inline-formula id="inf38">
<mml:math id="m43">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. Take <inline-formula id="inf39">
<mml:math id="m44">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> as an example, <inline-formula id="inf40">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be written as according to the full probability formula:<disp-formula id="e6_A">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mi mathvariant="bold">&#x3b5;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
<label>(6-A)</label>
</disp-formula>
<disp-formula id="e6_B">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo mathvariant="bold">/</mml:mo>
<mml:mi mathvariant="bold">&#x3b5;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6-B)</label>
</disp-formula>
<disp-formula id="e6_C">
<mml:math id="m48">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6-C)</label>
</disp-formula>
<disp-formula id="e6_D">
<mml:math id="m49">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3b5;</mml:mi>
<mml:mo mathvariant="bold">&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold">&#x3b5;</mml:mi>
<mml:msup>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(6-D)</label>
</disp-formula>
</p>
<p>Therefore, we obtain the probability of packets collision and different parts collision in the receiving window as shown in equation (6) and (7), respectively. Meanwhile, we also need to consider the influence of multiple packets on the desired modulation packet in the receiving window. In a time-slotted synchronous RA system, packets collision means that these packets transmitted in the same time slot, while in an asynchronous RA system, packets collision has a different degree of collision. In general, the degree of collision is a time-varying with interference signal, which is a function of the number of interference packets. To facilitate modeling this time-varying interference, we focus on the average interference on the desired demodulation packet. From the reference (<xref ref-type="bibr" rid="B11">Hall, 1927</xref>), it can be derived according to the average collision ratio of a Poisson process which approximately follows the Irwin-Hall distribution. The probability density function (PDF) of the average collision ratio is written as:<disp-formula id="e7_A">
<mml:math id="m50">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi>&#x3c7;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="bold">1</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold">2</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>!</mml:mo>
</mml:mrow>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mi>a</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi mathvariant="italic">n</mml:mi>
</mml:msup>
</mml:mrow>
</mml:mstyle>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3c7;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">n</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mi mathvariant="bold-italic">sign</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7-A)</label>
</disp-formula>
<disp-formula id="e7_B">
<mml:math id="m51">
<mml:mrow>
<mml:mi mathvariant="bold">sign</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold">1</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:mi mathvariant="bold">0</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">0</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mi mathvariant="bold">0</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7-B)</label>
</disp-formula>where <inline-formula id="inf52">
<mml:math id="m63">
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
</mml:math>
</inline-formula> is the collision ratio satisfied with <inline-formula id="inf53">
<mml:math id="m64">
<mml:mrow>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> (<inline-formula id="inf54">
<mml:math id="m65">
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">0</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> means no collision, <inline-formula id="inf55">
<mml:math id="m66">
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> means complete collision, and <inline-formula id="inf56">
<mml:math id="m67">
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
</mml:math>
</inline-formula> in fractional form means partial collision).</p>
<p>Based on the analysis in <italic>Packet and receiving window structure</italic>, the PLR can be divided into following three cases: (a) the preamble of desired demodulation packets collide in all receiving windows; (b) the weight of multi-collision packets is greater than the number of receiving satellites; (c) the SINR of output signal of cooperative beamformer by iteration way does not meet the requirement of demodulation threshold.</p>
<sec id="s3-3-1">
<title>Packet Loss Rate for Case (a)</title>
<p>Due to the preamble of the desired packets collide in all receiving windows, the bit information corresponding to these packets cannot be identified and located in step 1) of proposed CBA scheme. Thus, the PLR of a desired demodulation packet multi-collision with <inline-formula id="inf57">
<mml:math id="m68">
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:math>
</inline-formula> packets can be expressed as:<disp-formula id="e8">
<mml:math id="m69">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2245;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mi mathvariant="bold">&#x221e;</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mi mathvariant="bold">&#x393;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">SINR</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mi mathvariant="bold-italic">d&#x3c7;</mml:mi>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>where <inline-formula id="inf58">
<mml:math id="m70">
<mml:mrow>
<mml:mi mathvariant="bold">&#x393;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">SINR</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is packet error rate function in (<xref ref-type="bibr" rid="B8">Oscar and Riccardo, 2014</xref>) and<inline-formula id="inf59">
<mml:math id="m71">
<mml:mrow>
<mml:mi mathvariant="bold-italic"> SINR</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c7;P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">N</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>. The probability <inline-formula id="inf60">
<mml:math id="m72">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> of a desired demodulation packet multi-collision with <inline-formula id="inf61">
<mml:math id="m73">
<mml:mi>l</mml:mi>
</mml:math>
</inline-formula> packets can be written as:<disp-formula id="e9">
<mml:math id="m74">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="italic">k</mml:mi>
<mml:mo mathvariant="bold">&#x3d;</mml:mo>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:mrow>
<mml:mi>&#x221e;</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2217;</mml:mo>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>where<inline-formula id="inf62">
<mml:math id="m75">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic"> p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2217;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:msubsup>
<mml:mo>&#x2217;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> denotes the probability of collision packets number <inline-formula id="inf63">
<mml:math id="m76">
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:math>
</inline-formula> in the arriving packets number <inline-formula id="inf64">
<mml:math id="m77">
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:math>
</inline-formula>, and <inline-formula id="inf65">
<mml:math id="m78">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the probability of preamble collision for a single packet, which can be obtained from <bold>Eq.7</bold>. Therefore, in the presence of <inline-formula id="inf66">
<mml:math id="m79">
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:math>
</inline-formula> collision packets, the PLR <inline-formula id="inf67">
<mml:math id="m80">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> can be approximately expressed by the cumulative sum of the desired demodulation packet PLR <inline-formula id="inf68">
<mml:math id="m81">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and the probability <inline-formula id="inf69">
<mml:math id="m82">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> as follows:<disp-formula id="e10">
<mml:math id="m83">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">l</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi>&#x221e;</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev </mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x2217;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
</sec>
<sec id="s3-3-2">
<title>Packet loss rate for case (b).</title>
<p>In this case, the desired demodulation packets collide, while we can detect the non-collision preamble in all receiving windows. The gateway station can perform cooperative beamforming to suppress collision packets. However, according to the antenna array theory, <inline-formula id="inf70">
<mml:math id="m84">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> receiving satellites can only generate <inline-formula id="inf71">
<mml:math id="m85">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> nulls to resolve <inline-formula id="inf72">
<mml:math id="m86">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> weight of multi-collision packets. According to the analysis of the case (a), the probability that <inline-formula id="inf73">
<mml:math id="m87">
<mml:mi>n</mml:mi>
</mml:math>
</inline-formula> packets out of <inline-formula id="inf74">
<mml:math id="m88">
<mml:mi>k</mml:mi>
</mml:math>
</inline-formula> packets collide with a particular packet in a receiving window, but their preamble does not collide, can be expressed as the formula: <disp-formula id="e11">
<mml:math id="m89">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="italic">n</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2217;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:msubsup>
<mml:mo>&#x2217;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>where <inline-formula id="inf75">
<mml:math id="m90">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the probability of data part collision for a single packet, which can be obtained from <bold>Eq.7</bold>. In a receiving window, the PLR of the desired demodulation packet in the presence of multi-collision with <inline-formula id="inf76">
<mml:math id="m91">
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:math>
</inline-formula> packets can be written as:<disp-formula id="e12">
<mml:math id="m92">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2245;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mi>&#x221e;</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mi mathvariant="bold">&#x393;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">SINR</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi>&#x3c7;</mml:mi>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>Then, the probability that the weight of multi-collision packets is greater than the number of receiving windows can be formulated as:<disp-formula id="e13">
<mml:math id="m93">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi>&#x221e;</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>The PLR for case (b) can be written as:<disp-formula id="e14">
<mml:math id="m94">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi>&#x221e;</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
</p>
</sec>
<sec id="s3-3-3">
<title>Packet loss rate for case (c).</title>
<p>Due to the fact that the packets performed by beamforming operation may also fail to reach the demodulation threshold by iteration SIC operation, which will result in packet loss. In this case, the PLR is related to with the number of iterations <inline-formula id="inf77">
<mml:math id="m95">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of cooperative beamforming. According to the reference (<xref ref-type="bibr" rid="B5">Riccardo et&#x20;al., 2014</xref>), the PLR of a desired demodulation packet with <inline-formula id="inf78">
<mml:math id="m96">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> iterations can be expressed as:<disp-formula id="e15">
<mml:math id="m97">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msup>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:mrow>
<mml:mi>&#x221e;</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mi mathvariant="bold-italic">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>where <inline-formula id="inf79">
<mml:math id="m98">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the PLR that there is still <inline-formula id="inf80">
<mml:math id="m99">
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:math>
</inline-formula> packets collision with <inline-formula id="inf81">
<mml:math id="m100">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> iterations. However, in the process of iterative SIC technology, it is possible that some copies of the desired demodulation packet can be eliminated. The collision probability is variable with time, which is related to the number of iterations <inline-formula id="inf82">
<mml:math id="m101">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and the number of copies. Thus <inline-formula id="inf83">
<mml:math id="m102">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> can be written as:<disp-formula id="e16">
<mml:math id="m103">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">r</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold">k</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>where <inline-formula id="inf84">
<mml:math id="m104">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the binomial distribution with the number of experiments <inline-formula id="inf85">
<mml:math id="m105">
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:math>
</inline-formula>, and <inline-formula id="inf86">
<mml:math id="m106">
<mml:mrow>
<mml:mi mathvariant="bold-italic">q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msup>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo mathvariant="bold">&#x2212;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo mathvariant="bold">&#x2212;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the binomial distribution success probability. The value of <inline-formula id="inf87">
<mml:math id="m107">
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:math>
</inline-formula> is the PLR with the previous round of SIC operation, and the initial value <inline-formula id="inf88">
<mml:math id="m108">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msup>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msup>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> when <inline-formula id="inf89">
<mml:math id="m109">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. The <inline-formula id="inf90">
<mml:math id="m110">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> denotes the probability that there is still <inline-formula id="inf91">
<mml:math id="m111">
<mml:mi mathvariant="bold-italic">r</mml:mi>
</mml:math>
</inline-formula> packets cannot be correctly demodulated with <inline-formula id="inf92">
<mml:math id="m112">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> iterations, which can be approximately expressed as:<disp-formula id="e17">
<mml:math id="m113">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2245;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mi>&#x221e;</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mi mathvariant="bold">&#x393;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">SINR</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3c7;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">r</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mi mathvariant="bold-italic">d&#x3c7;</mml:mi>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
</p>
<p>Therefore, the total PLR of proposed CBA scheme can be written as:<disp-formula id="e18">
<mml:math id="m114">
<mml:mrow>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">CBA</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">iter</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
</p>
<p>According to the relationship between total PLR and throughput, the throughput of proposed CBA scheme can be expressed as follows:<disp-formula id="e19">
<mml:math id="m115">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">CBA</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>&#x2217;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">1</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">PL</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">CBA</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>
</p>
</sec>
</sec>
</sec>
<sec id="s4">
<title>Simulation</title>
<p>In this section, we analyze the throughput and PLR performance of proposed CBA scheme compared with traditional CRDSA scheme in (<xref ref-type="bibr" rid="B2">Enrico et&#x20;al., 2007</xref>) and ACA algorithm in (<xref ref-type="bibr" rid="B13">Pengxu et&#x20;al., 2017</xref>). The physical and MAC layer is jointly performed in the MATLAB simulation platform and the Monte Carlo method is adopted to verify the theoretical performance. These simulation conditions are similar to the benchmark problems. The simulation conditions are similar to the benchmark problems and the simulation parameters are given in <xref ref-type="table" rid="T1">Table&#x20;1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Simulation parameters.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">
<italic>variable</italic>
</th>
<th align="center">
<italic>significance</italic>
</th>
<th align="center">
<italic>value</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">
<inline-formula id="inf98">
<mml:math id="m121">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">rev</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">Number of satellites</td>
<td align="center">[2,3]</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf99">
<mml:math id="m122">
<mml:mi mathvariant="bold-italic">&#x3bb;</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="center">Traffic arrival parameter</td>
<td align="center">[0:0.1:2]</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf100">
<mml:math id="m123">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">Packet length</td>
<td align="center">100 bits</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf101">
<mml:math id="m124">
<mml:mi mathvariant="bold-italic">r</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="center">FEC Code rate</td>
<td align="center">0.5</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf102">
<mml:math id="m125">
<mml:mi mathvariant="bold-italic">b</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="center">Modulation order</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf103">
<mml:math id="m126">
<mml:mi mathvariant="bold">&#x3b5;</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="center">Preamble to packet ratio</td>
<td align="center">1/10</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf104">
<mml:math id="m127">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mrow>
<mml:mi>iter</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">Maximum iteration number</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf105">
<mml:math id="m128">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">&#x3c4;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">Maximum propagation time delay difference</td>
<td align="center">
<inline-formula id="inf106">
<mml:math id="m129">
<mml:mrow>
<mml:mn mathvariant="bold">50</mml:mn>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>
<bold>Figure&#x20;7</bold> shows the normalized throughput versus normalized load for the proposed CBA scheme compared with the CRDSA scheme and ACA scheme. As a baseline, we consider that the number of receiving satellites for the proposed CBA scheme and ACA scheme is the same as the number of replicas for the CRDSA scheme to generate the packet copies in spatial and time domain, respectively. From <xref ref-type="fig" rid="F7">Figure&#x20;7A</xref>, we can find that throughput curves of different schemes linearly increase to a peak value in the low load area and then steeply degrade in the high load area. For proposed CBA scheme and ACA scheme, the throughput performances increase with the number of receiving satellite. The maximum normalized throughput reaches to <inline-formula id="inf107">
<mml:math id="m130">
<mml:mrow>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi mathvariant="bold">5</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> (bit/symbol) and <inline-formula id="inf108">
<mml:math id="m131">
<mml:mrow>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi mathvariant="bold">9</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> (bit/symbol) using two and three receiving satellites, respectively, for ACA scheme. Meanwhile, the maximum normalized throughput approximately reaches to <inline-formula id="inf109">
<mml:math id="m132">
<mml:mrow>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi mathvariant="bold">6</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>(bit/symbol) and <inline-formula id="inf110">
<mml:math id="m133">
<mml:mrow>
<mml:mi mathvariant="bold">1</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> (bit/symbol) using two and three receiving satellites, respectively, for the proposed CBA scheme. However, for CRDSA scheme, the throughput performance seriously degrades with the packet copies increase in the high load area, because packet copies generated in the time/frequency domain increase the network load, which will result in &#x2018;avalanche effect&#x2019; in high load area. The maximum normalized throughput approximately reaches to <inline-formula id="inf111">
<mml:math id="m134">
<mml:mrow>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi mathvariant="bold">55</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> (bit/symbol) and <inline-formula id="inf112">
<mml:math id="m135">
<mml:mrow>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi mathvariant="bold">7</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> (bit/symbol) using two and three packet copies, respectively, for the CRDSA scheme. For a fixed normalized load <inline-formula id="inf113">
<mml:math id="m136">
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, the throughput performance of proposed CBA scheme increases <inline-formula id="inf114">
<mml:math id="m137">
<mml:mrow>
<mml:mi mathvariant="bold">10</mml:mi>
<mml:mo>%</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf115">
<mml:math id="m138">
<mml:mrow>
<mml:mi mathvariant="bold">80</mml:mi>
<mml:mo>%</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> performance compared with ACA scheme and CRDSA scheme, respectively, when the number of packet copies equals to <inline-formula id="inf116">
<mml:math id="m139">
<mml:mi mathvariant="bold">3</mml:mi>
</mml:math>
</inline-formula>. From <xref ref-type="fig" rid="F7">Figure&#x20;7B</xref>, the simulation results show that the theoretical Throughput curves are well matched with those simulated by Monte Carlo method of both two and three receiving satellites.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>
<bold>(A)</bold> Throughput performance of proposed CBA scheme compared with CRDSA and ACA scheme. <bold>(B)</bold>. Theoretical Throughput and Monte Carlo simulation Throughput performance of proposed CBA scheme.</p>
</caption>
<graphic xlink:href="frspt-02-755520-g007.tif"/>
</fig>
<p>
<bold>Figure&#x20;8</bold> shows the PLR performance versus normalized load for the proposed CBA scheme compared with the CRDSA scheme and ACA scheme. From <xref ref-type="fig" rid="F8">Figure&#x20;8A</xref>, we can find that the PLR of CRDSA scheme is lower than proposed CBA scheme in low load area, because CRDSA scheme adopts time slot synchronization. However, the PLR of proposed CBA scheme is lower than CRDSA scheme in the medium and high load area. For a fixed normalized load <inline-formula id="inf117">
<mml:math id="m140">
<mml:mrow>
<mml:mi mathvariant="bold-italic">G</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, PLR value for the proposed CBA scheme reaches to <inline-formula id="inf118">
<mml:math id="m141">
<mml:mrow>
<mml:mi mathvariant="bold">1</mml:mi>
<mml:msup>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold">2</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> when the number of copies equals to <inline-formula id="inf119">
<mml:math id="m142">
<mml:mi mathvariant="bold">3</mml:mi>
</mml:math>
</inline-formula>, while PLR value for CRDSA and ACA scheme approximately degrades to <inline-formula id="inf120">
<mml:math id="m143">
<mml:mn>1</mml:mn>
</mml:math>
</inline-formula> and <inline-formula id="inf121">
<mml:math id="m144">
<mml:mrow>
<mml:mi mathvariant="bold">1</mml:mi>
<mml:msup>
<mml:mi mathvariant="bold">0</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold">1</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula>, respectively. From <xref ref-type="fig" rid="F8">Figure&#x20;8B</xref>, the simulation results show that the theoretical PLR curves are well matched with those simulated by Monte Carlo method of both two and three receiving satellites.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>
<bold>(A)</bold> PLR performance of proposed CBA scheme compared with CRDSA and ACA scheme. <bold>(B)</bold> Theoretical PLR and Monte Carlo simulation PLR performance of proposed CBA scheme.</p>
</caption>
<graphic xlink:href="frspt-02-755520-g008.tif"/>
</fig>
</sec>
<sec sec-type="conclusion" id="s5">
<title>Conclusion</title>
<p>To increase the throughput performance of RA system, we presented a CBA scheme to solve the problem of &#x2018;deadlock&#x2019; in ACA scheme under the multi-satellite common-view scenario. The cooperative beamforming collision resolution technique combined with SIC scheme was designed to efficiently resolve packet collision by iteration way at the gateway. Simulation results show that the proposed CBA scheme can effectively increase the RA performance compared with traditional ACA scheme. Further works on cooperative beamforming algorithm for packet collision resolution need to be discussed such as the robust design and collaboration receiver selection strategy of beamforming algorithm.</p>
</sec>
</body>
<back>
<sec id="s6">
<title>Data Availability Statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>LR and HT led the main contribution and writing of the manuscript. DX, WY, and ZG supervised the findings of this work. All authors contributed to the manuscript and approved the submitted version.</p>
</sec>
<sec id="s8">
<title>Funding</title>
<p>This work presented was partially supported by the National Science Foundation of China (No. 91738201&#x3001;62171234 and 61971440), the Jiangsu Province Basic Research Project (No. BK20192002), the China Postdoctoral Science Foundation (No. 2018M632347), and the Natural Science Research of Higher Education Institutions of Jiangsu Province (No. 18KJB510030).</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abramson</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>1970</year>). <article-title>The ALOHA System-Another Alternative for Computer Communications</article-title>. <source>Afips Conf. Proc.</source> <pub-id pub-id-type="doi">10.21236/ad0707853</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Casini</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>De Gaudenzi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Herrero</surname>
<given-names>O. R.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Contention Resolution Diversity Slotted ALOHA (CRDSA): An Enhanced Random Access Schemefor Satellite Access Packet Networks</article-title>. <source>IEEE Trans. Wireless Commun.</source> <volume>6</volume> (<issue>4</issue>), <fpage>1408</fpage>&#x2013;<lpage>1419</lpage>. <pub-id pub-id-type="doi">10.1109/TWC.2007.348337</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chien</surname>
<given-names>W.-C.</given-names>
</name>
<name>
<surname>Lai</surname>
<given-names>C.-F.</given-names>
</name>
<name>
<surname>Hossain</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>Muhammad</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Heterogeneous Space and Terrestrial Integrated Networks for IoT: Architecture and Challenges</article-title>. <source>IEEE Netw.</source> <volume>33</volume> (<issue>1</issue>), <fpage>15</fpage>&#x2013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.1109/MNET.2018.1800182</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cioni</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>De Gaudenzi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Del Rio Herrero</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Girault</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>On the Satellite Role in the Era of 5G Massive Machine Type Communications</article-title>. <source>IEEE Netw.</source> <volume>32</volume> (<issue>5</issue>), <fpage>54</fpage>&#x2013;<lpage>61</lpage>. <pub-id pub-id-type="doi">10.1109/MNET.2018.1800024</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>De Gaudenzi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>del Rio Herrero</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Acar</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Garrido Barrabes</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Asynchronous Contention Resolution Diversity ALOHA: Making CRDSA Truly Asynchronous</article-title>. <source>IEEE Trans. Wireless Commun.</source> <volume>13</volume> (<issue>11</issue>), <fpage>6193</fpage>&#x2013;<lpage>6206</lpage>. <pub-id pub-id-type="doi">10.1109/TWC.2014.2334620</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>De Gaudenzi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>del Rio Herrero</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Advances in Random Access Protocols for Satellite Networks</article-title>. <source>Int. Workshop Satellite Space Commun.</source>, <fpage>331</fpage>&#x2013;<lpage>336</lpage>. <pub-id pub-id-type="doi">10.1109/IWSSC.2009.5286351</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>De Sanctis</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Cianca</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Araniti</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Bisio</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Prasad</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Satellite Communications Supporting Internet of Remote Things</article-title>. <source>IEEE Internet Things J.</source> <volume>3</volume> (<issue>1</issue>), <fpage>113</fpage>&#x2013;<lpage>123</lpage>. <pub-id pub-id-type="doi">10.1109/JIOT.2015.2487046</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>del Rio Herrero</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>De Gaudenzi</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Generalized Analytical Framework for the Performance Assessment of Slotted Random Access Protocols</article-title>. <source>IEEE Trans. Wireless Commun.</source> <volume>13</volume> (<issue>2</issue>), <fpage>809</fpage>&#x2013;<lpage>821</lpage>. <pub-id pub-id-type="doi">10.1109/TWC.2013.121813.130435</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ding</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Qu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Machine Learning Enabled Preamble Collision Resolution in Distributed Massive MIMO</article-title>. <source>IEEE Trans. Commun.</source> <volume>69</volume> (<issue>4</issue>), <fpage>2317</fpage>&#x2013;<lpage>2330</lpage>. <pub-id pub-id-type="doi">10.1109/TCOMM.2021.3051202</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Duan</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Dinavahi</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Starlink Space Network-Enhanced Cyber-Physical Power System</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>12</volume> (<issue>4</issue>), <fpage>3673</fpage>&#x2013;<lpage>3675</lpage>. <pub-id pub-id-type="doi">10.1109/TSG.2021.3068046</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hall</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>1927</year>). <article-title>The Distribution of Means for Samples of Size N Drawn from a Population in Which the Variate Takes Values between 0 and 1, All Such Values Being Equally Probable</article-title>. <source>Biometrika</source> <volume>19</volume> (<issue>3</issue>), <fpage>240</fpage>&#x2013;<lpage>245</lpage>. <pub-id pub-id-type="doi">10.2307/2331961</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jia</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Qi</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Massive Beam-Division Multiple Access for B5G Cellular Internet of Things</article-title>. <source>IEEE Internet Things J.</source> <volume>7</volume> (<issue>3</issue>), <fpage>2386</fpage>&#x2013;<lpage>2396</lpage>. <pub-id pub-id-type="doi">10.1109/JIOT.2019.2958129</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Asynchronous Cooperative Aloha for Multi-Receiver Satellite Communication Networks</article-title>. <source>IEEE Commun. Lett.</source> <volume>21</volume> (<issue>6</issue>), <fpage>1321</fpage>&#x2013;<lpage>1324</lpage>. <pub-id pub-id-type="doi">10.1109/LCOMM.2017.2672985</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Multi-Beam UAV Communication in Cellular Uplink: Cooperative Interference Cancellation and Sum-Rate Maximization</article-title>. <source>IEEE Trans. Wireless Commun.</source> <volume>18</volume>, <fpage>4679</fpage>&#x2013;<lpage>4691</lpage>. <pub-id pub-id-type="doi">10.1109/TWC.2019.2926981</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liva</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Graph-Based Analysis and Optimization of Contention Resolution Diversity Slotted ALOHA</article-title>. <source>IEEE Trans. Commun.</source> <volume>59</volume> (<issue>2</issue>), <fpage>477</fpage>&#x2013;<lpage>487</lpage>. <pub-id pub-id-type="doi">10.1109/TCOMM.2010.120710.100054</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Munari</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Clazzer</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Liva</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Multi-receiver Aloha Systems - a Survey and New Results</article-title>. <source>IEEE Int. Conf. Commun. Workshop</source>, <fpage>2108</fpage>&#x2013;<lpage>2114</lpage>. <pub-id pub-id-type="doi">10.1109/ICCW.2015.7247493</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pachler</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>del Portillo</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Crawley</surname>
<given-names>E. F.</given-names>
</name>
<name>
<surname>Cameron</surname>
<given-names>B. G.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>An Updated Comparison of Four Low Earth Orbit Satellite Constellation Systems to Provide&#x20;Global Broadband</article-title>. <source>ICC Workshops</source>, <fpage>1</fpage>&#x2013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1109/ICCWorkshops50388.2021.9473799</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>LEO Satellite Constellation for Internet of Things</article-title>. <source>IEEE Access</source> <volume>5</volume>, <fpage>18391</fpage>&#x2013;<lpage>18401</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2017.2735988</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roberts</surname>
<given-names>L. G.</given-names>
</name>
</person-group> (<year>1975</year>). <article-title>ALOHA Packet System with and without Slots and Capture</article-title>. <source>SIGCOMM Comput. Commun. Rev.</source> <volume>5</volume> (<issue>2</issue>), <fpage>28</fpage>&#x2013;<lpage>42</lpage>. <pub-id pub-id-type="doi">10.1145/1024916.1024920</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Matthaiou</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Beam Division Multiple Access Transmission for Massive MIMO Communications</article-title>. <source>IEEE Trans. Commun.</source> <volume>63</volume> (<issue>6</issue>), <fpage>2170</fpage>&#x2013;<lpage>2184</lpage>. <pub-id pub-id-type="doi">10.1109/TCOMM.2015.2425882</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>An Overview of Enhanced Massive MIMO with Array Signal Processing Techniques</article-title>. <source>IEEE J.&#x20;Sel. Top. Signal. Process.</source> <volume>13</volume> (<issue>5</issue>), <fpage>886</fpage>&#x2013;<lpage>901</lpage>. <pub-id pub-id-type="doi">10.1109/JSTSP.2019.2934931</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ng</surname>
<given-names>D. W. K.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>NOMA for Hybrid mmWave Communication Systems with Beamwidth Control</article-title>. <source>IEEE J.&#x20;Sel. Top. Signal. Process.</source> <volume>13</volume> (<issue>3</issue>), <fpage>567</fpage>&#x2013;<lpage>583</lpage>. <pub-id pub-id-type="doi">10.1109/JSTSP.2019.2901593</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xia</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>You</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Satellite Machine-type Communication for Maritime Internet of Things: An Interference Perspective</article-title>. <source>IEEE Access</source> <volume>7</volume>, <fpage>76404</fpage>&#x2013;<lpage>76415</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2019.2919253</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>You</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Beam-Domain Channel Estimation for FDD Massive MIMO Systems with Optimal Thresholds</article-title>. <source>IEEE Trans. Wireless Commun.</source> <volume>16</volume> (<issue>7</issue>), <fpage>4669</fpage>&#x2013;<lpage>4682</lpage>. <pub-id pub-id-type="doi">10.1109/TWC.2017.2701371</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>X.-G.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>3D Deployment of Multiple UAV-Mounted Base Stations for UAV Communications</article-title>. <source>IEEE Trans. Commun.</source> <volume>69</volume> (<issue>4</issue>), <fpage>2473</fpage>&#x2013;<lpage>2488</lpage>. <pub-id pub-id-type="doi">10.1109/TCOMM.2021.3049387</pub-id> </citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Multisatellite Cooperative Random Access Scheme in Low Earth Orbit Satellite Networks</article-title>. <source>IEEE Syst. J.</source> <volume>13</volume> (<issue>3</issue>), <fpage>2617</fpage>&#x2013;<lpage>2628</lpage>. <pub-id pub-id-type="doi">10.1109/JSYST.2018.2870416</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>