#### Abstract

We evaluate a new approach for multisatellite state modelling: the Master-Slave approach. By this concept slave satellites are modelled according to an existing master, whereas the correlation between multiple slaves is omitted. Master-Slave is therefore a generic name for a state modelling concept, for which different realisations are possible. As a possible realisation we present the Conditional Assembling Method. For modelling of only two satellites (one master and one slave), the Conditional Assembling Method enables an accurate resimulation of the correlation coefficient between the satellites and the probabilities of single and combined states. Based on this condition, the performance of Master-Slave for three, four, and five satellites is evaluated in terms of state probability modelling. Therefore, the correlation coefficients and the all bad-state probabilities with Master-Slave are compared with the measurements for different elevation angles and azimuth angle separations of the multisatellite system. Master-Slave has a high modelling error in case of small azimuth separation between the slave satellites (except that one slave has a small azimuth separation to the master). Furthermore, a master satellite with a high elevation provides a lower probability error compared to a master with low elevation.

#### 1. Introduction

Satellites play an important role in today’s commercial broadcasting systems. In cooperation with terrestrial repeaters they can ensure uninterrupted service of multimedia content (e.g., audio and video streaming) to stationary, portable, and mobile receivers. However, in case of mobile reception fading regularly disrupts the signal transmission due to shadowing or blocking objects between satellite and receiver. To mitigate these fading effects, diversity techniques such as angle diversity (multiple satellites) and time diversity (interleaving) are attractive. For link-level studies of the land mobile satellite (LMS) channel, statistical channel models are frequently used that are able to generate time series of the received fading signal [1–3].

Statistical LMS channel models describe several fading processes of the received signal: slow variations of the signal are caused by obstacles between the satellite and the receiver, which induce varying shadowing conditions of the direct signal component. Fast signal variations are caused by multipath effects due to static or moving scatterers in the vicinity of the mobile terminal. For short time periods these two components (*slow and fast variations*) are usually modelled by a stationary stochastic process, for example, as a Loo-distributed fading signal [1]. For longer time periods the received signal cannot be assumed as stationary. Therefore, statistical LMS channel models describe different receiving states to assess the large dynamic range of the received signal. The states (very slow variations) correspond to slowly varying environmental conditions (e.g., line of trees, buildings, zones with line-of-sight condition) in the transmission path. Figure 1 shows a state-of-the-art LMS narrowband model for single-satellite reception according to Prieto-Cerdeira et al. [4]. It describes two states: a “good” state (corresponding to line-of-sight/light shadowing) and a “bad” state (corresponding to heavy shadowing/blockage). Within the states this model assumes a stationary Loo-distributed fading signal. It includes a slow fading component (lognormal fading) corresponding to varying shadowing conditions of the direct signal and a fast fading component due to multipath effects. To enable a realistic modelling over the full dynamic range of the received signals, a random set of Loo parameters (, , and ) is generated after each state transition.

Focus of recent activities is the extension of this model for dual-satellite and multisatellite reception. It is realised by introducing correlation between states and between statistical parameters for fast and slow variations. However, the dominating part of the overall signal correlation stems from the state correlation.

Focus of this paper is state sequence modelling for multisatellite systems, assuming two states per satellite: “good” state and “bad” state. For this purpose, a new method for state sequence generation is introduced: the Master-Slave approach. Within Master-Slave it is assumed that each slave satellite depends only on a master satellite, whereas the correlation between different slaves is not described.

The paper is organised as follows: Section 2 gives an overview on state modelling approaches for single- and dual-satellite reception. In Section 3 the concept Master-Slave and a first implementation are introduced. Afterwards, in Section 4 a validation of Master-Slave for multisatellite reception in terms of state probability modelling is presented. Finally in Section 5 the conclusions are drawn.

#### 2. State Modelling for Single-Satellite and Dual-Satellite Reception

In this section available state modelling approaches for single- and dual-satellite reception are summarised.

##### 2.1. Channel State Models for Single-Satellite Systems

###### 2.1.1. First-Order Markov Model

A Markov model is a special random process for generating discrete samples corresponding to channel states of a predefined sample length. For a first-order Markov model, each state depends only on the previous state. The conditional probabilities of state given the state are described by state transition probabilities . Therefore, the only parameter of the Markov chain is the state transition probability matrix (STPM) with being the number of states. The main characteristic of a first-order Markov chain is that it enables an exact modelling of the state probability and the average state duration. The state duration probability density function (SDPDF) of the first-order Markov chain follows an exponential distribution, which is found as not realistic for the LMS channel [5, 6]. The first-order Markov model is used in early LMS models [2, 3].

###### 2.1.2. Semi-Markov Model

To improve the state duration modelling for the LMS channel, semi-Markov chains were proposed in [5]. In contrast to the first-order Markov model, the state transitions do not occur at concrete time intervals. In fact, the time interval of the model staying in state depends directly on its SDPDF. As with the first-order Markov model, the state transitions are described with the state transition probability , but with . Assuming a single-satellite model of only two states, the state transition probability is . The semi-Markov model allows some options to describe the SDPDF of each state, which defines the number of required parameters for the semi-Markov model. An overview is given in [7].

###### 2.1.3. Dynamic Markov Model

A further method to improve the state duration modelling are dynamic Markov chains introduced in [6]. For dynamic Markov chains the state transition probability depends on the current state duration with . Therefore, the two-dimensional STPM is extended to a three-dimensional state transition probability tensor (STPT) , where corresponds to the maximum state length. The dynamic Markov model enables an exact reproduction of the state probabilities and also an exact remodelling of the measured SDPDF. As a disadvantage, a high number of parameters are required to describe the STPT. Model approximations to reduce the number of parameters are presented in [6].

##### 2.2. Channel State Models for Dual-Satellite Systems

###### 2.2.1. Straightforward Method: Extension to a Multistate Model

The first-order Markov model, semi-Markov model, and the dynamic Markov model can be easily adapted for dual- or multisatellite modelling. This can be achieved by combining the single-satellite states “good” and “bad” from two satellites into joint states: “good good,” “good bad,” “bad good,” and “bad bad”. In case of the first-order Markov model, the STPM becomes a STPM. For the dynamic Markov model a STPT is required for the state series simulation. In case of the semi-Markov approach a STPM and four separate state duration statistics are required for the dual-satellite modelling.

An obvious problem with this straightforward approach is the exponential growth of the number of parameters with the number of satellites. A multisatellite model has combined states with being the number of states per satellite and being the number of satellites.

###### 2.2.2. Lutz Model

In [8] a state model for two correlated satellites based on first-order Markov chains was developed. This algorithm is based on first-order Markov chains and generates a joint STPM ( elements) from two independent single-satellite STPMs (each with elements), only by the use of one correlation coefficient. Using the joint STPM, a joint sequence of four combined states can be generated. The high flexibility of this algorithm becomes clear, since it requires only single-satellite parameter sets that are easy to derive from measurements and are already available in the literature for different elevation angles and a high number of environments. Databases for correlation coefficients are available as well for different environments, elevation angles, and angular separations of the azimuth and elevation angles. In contrast to this Lutz model, the above mentioned straightforward methods need complete datasets for any combination of elevation angles, azimuth angle separations, and environments to achieve the same variability. The Lutz model accurately resimulates the state probabilities of each single satellite as well as of the combined states. A limitation is that there is no flexibility in describing the state duration distribution. Furthermore, the current approach in [8] is limited to the usage of two satellites with two states per satellite.

#### 3. State Modelling for Multisatellite Reception: The Master-Slave Approach

In general, all state models mentioned in the previous section are extendable to multiple satellites by assuming a multistate model with combined states. However, the number of required parameters for these models grows exponentially with the number of satellites. In practise, already the description of a satellite system with two satellites is challenging [7]. To avoid excessive complexity of a state model parametrisation with more than two satellites, the Master-Slave approach was proposed in [9]. Within the Master-Slave approach it is assumed that each slave satellite depends only on a master satellite, whereas the correlation between the slaves is not described.

Figure 2 demonstrates the Master-Slave approach with four satellites. Within Master-Slave only three satellite links in terms of correlation coefficients are described therefore. A full description of a four-satellite system would require the description of six correlations. In case of satellites the Master-Slave approach describes satellite links, whereas the full description requires links. As a result, Master-Slave reduces the complexity of a -satellite model by .

It is not surprising that the missing correlation branches between the slaves could be a reduction of relevant information in case of multisatellite modelling. Therefore, a performance analysis of Master-Slave for multisatellite modelling (with at least three satellites) is found in Section 4. A first realisation of Master-Slave is described in Section 3.1.

##### 3.1. First Realisation of State Modelling with Master-Slave: Conditional Assembling Method

The Master-Slave approach in its simplest form can be reduced to a two-satellite modelling problem, that is, the modelling of the master and one slave. The challenge within Master-Slave is the generation of a slave state sequence based on an existing master state sequence by providing, for example, a certain state probability and a correlation coefficient between the satellites. For this purpose, the conditional state sequences of the slave in case of a constant master state have to be analysed.

According to this condition, various realisations of Master-Slave are possible. Thus, the term Master-Slave describes a concept—not a concrete implementation.

In the following we present one possible realisation of Master-Slave that we define as *Conditional Assembling Method* (cf. Figure 3). For parametrisation, the principle of the Conditional Assembling Method is first to concatenate all parts of the slave state sequence for which the master is in “good” state. Afterwards, this state sequence is parametrised following an arbitrary Markov model (e.g., first-order Markov, semi-Markov). Same procedure is done for the slave state sequence in case master is in “bad” state.

The simulation is done in reversed order. First, the master state sequence is modelled. Second, two conditional slave state sequences with respect to the master state are generated independently (i.e., the conditional slave sequence for master is “good,” and the conditional slave sequence for master is “bad”). For master and conditional slave states an arbitrary Markov model can be used. The final state sequence of the slave is composed of parts from both conditional sequences according to the current master state. To provide a sample-by-sample generation of the slave state sequence, the current conditional slave state generator is kept in “standby” mode when a state transition of the master is obtained.

In case that the chosen Markov model accurately describes the state probabilities of the master as well as the conditional state probabilities of the slave in case of a constant master state, then the state probabilities of the slave, the joint state probabilities, and the correlation coefficient between the satellites are perfectly remodelled. Furthermore, the state durations of the master sequence can be modelled with arbitrary accuracy (by, e.g., taking a complex model such as dynamic Markov).

However, the Conditional Assembling Method has some limitations in state duration modelling of the combined states (“good good,” “good bad,” etc.) and of the slave states (cf. Section 4.4).

##### 3.2. Exemplary Evaluation of Master-Slave

To validate the Master-Slave concept for multisatellite reception, Figure 4 shows an exemplary constellation of three satellites including correlation coefficients as obtained during SDARS measurements from the project MiLADY [9, 10]. Since only the correlation between the master and slaves is described, the simulation results depend highly on the definition of the master satellite. For this purpose, in Figure 4 two mappings (Mapping 1 and Mapping 2) are defined by assuming different positions of master and slave satellites. Based on these two definitions, Figure 5 shows the resimulation results for the Master-Slave approach in terms of correlation coefficients, state probabilities, and mean state durations.

**(a)**

**(b)**

In case of Master-Slave the slave satellites are modelled independently. In theory, the correlation coefficient between two independent state sequences is zero. In fact, within Master-Slave the correlation coefficient between the slaves depends on the individual correlation coefficients between master and the slaves (cf. in the appendix). It always holds .

For comparison, the results of a first-order Markov model (using a STPM for 8 joint states) are presented. It is obtained that the first-order Markov model enables an accurate description of correlation coefficients between all satellites, of state probabilities, and of average state durations for single satellite states and combined states. The Master-Slave realisation with Conditional Assembling Method in this example is based on the first-order Markov model as well (i.e., both the master state sequence and the conditional state sequences of the slaves are modelled with first-order Markov chains). The result is that for Mapping 1 and Mapping 2 the state probabilities of the individual satellites, the correlation coefficients between master and slaves, and the mean state durations of the master (which has a different position in the examples) are modelled accurately.

Mapping 1 was defined such that the correlation between the slaves is low (with ). The higher correlation coefficients between master and slaves ( and ) provide some correlation of the slaves as well. In Mapping 1 is remodelled accurately. As a consequence, also the joint state probabilities (“good good good”,…,“bad bad bad”) are remodelled accurately.

In Mapping 2 the correlation between the slaves is high (with ). Due to low correlation coefficients between the master and the slaves ( and ), *≈* 0 and deviates strongly from the measurements. As a result, the joint state probabilities are not accurately remodelled. For the application of system planning it should be noted that in case of describing an insufficient correlation a higher diversity gain for the multisatellite constellation will be predicted.

Hence, Mapping 2 suggests a lower probability of the combined state “bad bad bad” than obtained from measurements. Such a modelling error should be avoided. A solution for Master-Slave is an appropriate definition of master and slave satellites. Therefore, in Section 4 the performance of Master-Slave for different positions of master and slaves is evaluated.

#### 4. Validation of State Probability Modelling with Master-Slave

This section will focus on the modelling accuracy when using the Master-Slave approach incorporating three satellites. Quality criterion is the deviation of the joint “bad” state probability (), depending on the azimuth and elevation constellation. At the end of this section, we will derive constraints for appropriate satellite constellations and how to map them into the Master-Slave context.

##### 4.1. Analysis of a Three-Satellite System

Basis of the analysis is a three-satellite constellation according to Figure 6. In the sequel, the master satellite is always referred to as index 1, while indices 2 and 3 refer to the slaves. The constellation (and the mapping) is described by the elevation angles , , and and the azimuth angle separations , , and . Following the Master-Slave approach, only the correlations between master and slave are described ( and ), while the correlation between the slaves () is neglected. In the following, the consequences of the neglected correlation are investigated.

###### 4.1.1. Influence of the Azimuth Angle Separation

Figure 7 shows the state-correlation coefficient between two satellites with dependency on their azimuth angle separation . The elevation angles of both satellites are ; the driving direction can be assumed as uniformly distributed between 0 and . These data were derived from analysing dual-satellite constellations of GNSS measurements in urban environments during the project MiLADY and are presented in detail in [7]. It can be observed that low azimuth angle separations lead to a high correlation, while for the satellites are almost uncorrelated. A slight increase of the correlation is obtained for the antipodal constellation (). Due to the uniformly distributed driving direction the graph is symmetrical to and .

Based on these results for the dual-satellite case we now analyse the consequences of neglecting the slave correlation in the Master-Slave approach. Following the configuration according to Figure 6, we compare after resimulation using the Master-Slave approach versus the measured results. Figure 8 shows both cases with dependency on each of the azimuth angle separations , , and . The error is given by . Clearly, this error is large when the two slaves come close to each other, that is, —except when one slave is close to the master.

Besides correctly remodelling the state correlation, for designing mobile satellite communication systems a very important parameter is the probability of all satellites being in “bad” state, that is, . It characterises the probability of service unavailability for the entire system and requires special attention, therefore. According to Figure 9, in the measured scenario if all three satellites are colocated, and if two of the three satellites are colocated. Considering the probability error (the dashed red line in Figure 9) as a quality measure of the modelling approach, we can extend the analysis to further parameters. Therefore, in the next section the elevation-angle dependency is investigated.

For the sake of completeness, we show the consequences when neglecting not only the correlation between the slaves, but neglecting a further correlation as well as neglecting all correlations in the three-satellite system: as can be observed from Figure 10, for azimuth angle separations of *≈*30° the correlation cannot be neglected. Furthermore, the correlation should be considered for antipodal constellations, as can be seen for .

**(a)**

**(b)**

**(c)**

###### 4.1.2. Influence of the Elevation Angles

Based on the validation of Master-Slave with dependency on the azimuth angle separation, we now incorporate the elevation angles as further variables. The Figures 11 and 12 show the probability error with Master-Slave with dependency on the azimuth separations and for different combinations of elevation angles between master, slave satellite 2, and slave satellite 3. Figure 11 shows results with a constant master elevation , Figure 12 for , respectively. The figures include also the results from Figure 10 for the special case . Since the slave satellite 2, and slave satellite 3 can be exchanged, only the lower triangle is shown. From Figures 11 and 12 appropriate and inappropriate constellations (and mappings) for the three-satellite Master-Slave approach can be determined. The following is obtained.(i)Assuming constant elevation angles (represented by one subfigure), the error is maximal for a low azimuth separation between the slaves, while some azimuth separation to the master exists. This error can be avoided by redefinition of the system such that the low azimuth separation is between the master and one of the slaves. (ii)If one of the slaves has a higher elevation than the master and the second slave (e.g., ), then the error for low slave separations is reduced with respect to . The reason is that the correlation coefficient and therefore the correlation error between slaves are lower due to the elevation angle separation. (iii)The error with Master-Slave is high when both slaves have a higher elevation than the master. (iv)Comparing Figures 11 and 12, it is obtained that a master with a high elevation has great benefits. For all positions of the slaves the error for elevation is much lower than for . (v)Comparing the constellations and it is seen that the higher elevation angles provide a lower probability error . It should be noted that the correlation error in case of some elevation angles is similar. This result is therefore based on the satellite state probabilities, which depends on the current elevation angle. An analysis of the relation between single satellite state probabilities and the probability error with Master-Slave is found in Section 4.2. (vi)A great performance difference is seen between the constellations and . Although the satellites have the same position on the hemisphere, for the second mapping nearly no limitations of azimuth angle separations are obtained. An examination on the influence of different elevation angles of master and slave satellites is done in Section 4.2.

##### 4.2. Estimation of Worst-Case Probability Error for Master-Slave with Three, Four, and Five Satellites

From the results in the previous section it is known that a high probability error occurs when master and slave satellites are weakly correlated, but the correlation between slaves is high. Worst case would be no correlation between master and slave and full correlation between slaves. Figure 13 shows worst cases for Master-Slave with three, four, and five satellites (please note that between the satellites is not realistic, since different satellites must be exactly colocated therefore) Based on these worst-case constellations, an examination on the probability error (we define as the worst reception case of the multisatellite system. It includes , , and for a system with two, three, and four satellites, respectively) for the Master-Slave approach with more than three satellites is possible.

The highest satellite diversity gain () with a multisatellite system is provided if all satellites are uncorrelated (). In this case, the probability of the combined state is the product of the individual state probabilities: with being the number of satellites. In contrast, a multisatellite system has no gain if all satellites are fully correlated (). Then it holds .

Assuming the constellations in Figure 13 the Master-Slave model would provide uncorrelated slave satellites with . As a consequence, the Master-Slave model has uncorrelated channels and for the combined state results (assuming same state probabilities of the individual satellites). In fact, the multisatellite systems in Figure 13 have only two uncorrelated channels. The probability error of Master-Slave is therefore .

Figure 14 shows the worst-case probability error with dependency on the “bad” state probability of the individual satellites. It represents that the reception for all satellites have the same elevation. As a general tendency, a reduction of the probability (e.g., due to an increased elevation) of the single satellites provides a reduction of . This situation was obtained in Figures 11 and 12 between constellations and .

**(a)**

**(b)**

Figure 15 evaluates the worst-case probability error for different elevation angles of the master and the slaves. It is seen that a high elevation angle of the master provides a lower error .

**(a)**

**(b)**

##### 4.3. Appropriate Constellations for Master-Slave Approach with Three, Four, and Five Satellites

This section concludes appropriate and inappropriate satellite positions for Master-Slave with three satellites and extends it for the multisatellite case.

Figure 16 (left) shows combinations of azimuth separations between master and slave and master and slave for which the modelling error exceeds a certain threshold. Assuming that is an acceptable modelling result, appropriate and inappropriate regions (in terms of ) are found therefore for the Master-Slave approach with three satellites.

Concrete examples are given in Figure 16: by defining a certain position of satellite 2 (it has initially no constraints in terms of appropriate positions) the inappropriate azimuth separations for satellite 3 are detected. Figure 17 presents appropriate zones for satellite 3 for predefined positions of satellite 2 in a polar diagram. Given the example in Figure 16, the inappropriate positions of satellite 3 are —except when satellite 2 is close to the master.

Outgoing from the three-satellite analysis, a quantitative estimation for inappropriate zones for systems with more than three satellites is possible. It is done by decomposing the multisatellite constellation into combinations of each one master and two slaves. For example, from a four-satellite system the combinations master-slave2-slave3, master-slave2-slave4, and master-slave3-slave4 are separately evaluated. The result is that inappropriate positions of each slave depend on the position of another slave. Based on this method, Figure 18 shows exemplary constellations of a four-satellite system with predefined positions of satellites 2 and 3. The result is possible positions of satellite 4, for which the Master-Slave approach provides reliable results. In one of the constellations (third from left) there is no option for satellite 4, since the constellation between S2 and S3 is already inappropriate.

Figure 19 shows constellations of a system with five satellites. Appropriate positions of satellite 5 are resulting from defined positions of satellites 2, 3, and 4.

It can be concluded that with a higher number of satellites the range of appropriate positions decreases.

Please remember that the “inappropriate zones” define only constellations where is below a certain threshold (0.05 was selected exemplarily). For a system with more than three satellites it means that may not be simulated correctly with Master-Slave. It makes no statement about the absolute error . Constellations are even possible where can be neglected, although the Master-Slave constellation was found as inappropriate.

##### 4.4. State Duration Modelling with Master-Slave

In Section 3.1 we presented the Conditional Assembling Method as a first realisation of the Master-Slave approach. It provides an accurate modelling of the state probabilities of the single satellites as well as of the combined satellites for *dual-satellite* reception. However, for the configuration of satellite broadcasting systems with long time interleaving, an accurate modelling of the fading signal over time is crucial.

Figure 20 shows therefore the state duration statistics of the “bad” state from satellite 1 and satellite 2 and of the combined “bad bad” state. The complementary state duration CDFs from measurements are compared with resimulation results of a Master-Slave model and with a dual-satellite semi-Markov model as presented in [7].

**(a)**

**(b)**

**(c)**

It is seen that the Conditional Assembling Method has some limitations for remodelling the duration statistic for the combined “bad bad” state and of the slave satellites’ “bad” state. In contrast, the duration statistic of the master satellite can be exactly remodelled by using complex approaches like dynamic Markov chains.

For comparison, the dual-satellite semi-Markov model with lognormal approximation of the state duration distribution resimulates more accurately the statistics of the slave satellite and of the combined system.

It can be concluded that for state modelling of a dual-satellite system the Master-Slave approach with “Conditional Assembling Method” would not be the first choice.

An extensive validation of different Master-Slave realisations with focus on accurate state duration modelling is therefore the topic of ongoing work.

#### 5. Conclusions

In this paper we presented a new approach for multisatellite state modelling: the Master-Slave approach. The primary concept is that slave satellites are modelled according to an existing master state sequence, whereas the correlation between multiple slaves is omitted. Master-Slave is therefore a generic name for a modelling concept, for which different state models can be applied (such as first-order Markov, semi-Markov, and dynamic Markov models).

As a possible realisation of Master-Slave we described the “Conditional Assembling Method”. It enables an *exact* resimulation of the correlation coefficient, the single satellite state probabilities, and the combined state probabilities between one master and one slave.

Based on this result, a detailed performance evaluation of Master-Slave for a three-satellite system is carried out in terms of state probability modelling. Therefore, measured correlation coefficients between the satellites are compared with recalculated correlation coefficients from Master-Slave with dependency on the elevation angles and the azimuth angle separations of the three-satellite system. Afterwards, the probability of the state “bad bad bad” () resulting from Master-Slave is compared with analytically estimated values from measurements. The difference for from Master-Slave and from measurements is defined as probability error. It was obtained that Master-Slave has a high probability error in case of a low azimuth separation and therefore a high correlation between the slave satellites. An exception is when one of the slaves has a high correlation to the master. Furthermore, a master satellite with a high elevation provides a lower probability error compared to a master with low elevation. It can be concluded that a probability error with Master-Slave can be mitigated by appropriate definition of master and slave satellites.

Based on the analysis with three satellites, a performance analysis of Master-Slave is extended for systems with more than three satellites in terms of state probability modelling. Therefore, appropriate and inappropriate constellations of satellites were estimated, and the probability errors for worst-case constellations with Master-Slave were predicted.

In the last part of this paper we indicated that the Master-Slave implementation of this paper (Conditional Assembling Method) has some limitations according to state duration modelling. An evaluation of state duration modelling with Master-Slave is the topic of ongoing research activities. It depends amongst others on the applied state model, such as semi-Markov and dynamic Markov and requires new concepts for the realisation of Master-Slave. However, the elaboration on state probability modelling in this paper is universally valid for Master-Slave.

#### Appendix

#### State Probabilities and Correlation Coefficients for the Conditional Assembling Method

The equilibrium state probabilities of the master states and of the conditional slave states are given by defines the probability of the slave state in case of the master state is .

*Calculation of Probabilities for Joint States and Slave States (One Master, One Slave)*

From the master probabilities and conditional slave probabilities the combined probability vector is calculated:
Based on the state probabilities of the slave can be calculated finally:

*Calculation of State Probabilities with One Master and Multiple Slaves*

Analogously to the two-satellite case, the combined probability vector for multiple satellites (with satellites) is calculated from the master probabilities and multiple conditional slave probabilities:

*Correlation Coefficient between Master and Slave*

Since only two states per satellite are assumed, the phi coefficient [11] is used to describe the correlation coefficient between two satellites:
is the number of samples where sequence 1 has state and sequence 2 has state . For state modelling this formula can be adapted by using the joint state probabilities of two satellites:

*Correlation Coefficient between Two Slaves*

The correlation coefficient between the slaves depends on the combined state probabilities of the slaves:
The combined state probabilities of the slaves can be calculated from the joint state probabilities of the whole system Master -Slave -Slave . According to (A.4) it holds:
By inserting (A.8) in (A.7) the correlation coefficient between the slaves can be described as a function of the correlation coefficients between master and slave and master and slave :

#### Conflict of Interests

The authors declare that they have no conflict of interests.

#### Acknowledgments

The state analysis in this paper is based on measurement data from the project MiLADY [9]. This project was funded by the ARTES 5.1 Programme of the Telecommunications and Integrated Applications Directorate of the European Space Agency.