State modelling of the land mobile propagation channel with multiple satellites

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 di ﬀ erent 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 coe ﬃ cient between the satellites and the probabilities of single and combined states. Based on this condition, the performance of Master-Slave for three, four, and ﬁve satellites is evaluated in terms of state probability modelling. Therefore, the correlation coe ﬃ cients and the all bad-state probabilities with Master-Slave are compared with the measurements for di ﬀ erent 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.


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][2][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-theart 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 International Journal of Antennas and Propagation 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 (M A , Σ A , and MP) 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.

State Modelling for Single-Satellite and Dual-Satellite Reception
In this section available state modelling approaches for single-and dual-satellite reception are summarised.with N 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].

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 i depends directly on its SDPDF.As with the first-order Markov model, the state transitions are described with the state transition probability p i j , but with i / = j.Assuming a single-satellite model of only two states, the state transition probability is p i j = 1.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].

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 n with p i j = f (n).Therefore, the two-dimensional STPM is extended to a three-dimensional state transition probability tensor (STPT) P trans ∈ R N×N×nmax

0+
, where n max 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].

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 2 × 2 STPM becomes a 4 × 4 STPM.For the dynamic Markov model a 4 × 4 × n max STPT is required for the state series simulation.In case of the semi-Markov approach a 4 × 4 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 N k combined states with N being the number of states per satellite and k being the number of satellites.

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 (4 × 4 elements) from two independent single-satellite STPMs (each with 2 × 2 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.

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 k satellites the Master-Slave approach describes k − 1 satellite links, whereas Master (1) Slave ( 2) Slave (4) Figure 2: Master-Slave approach for multisatellite modelling.Several slave satellites are modelled according to the correlation to one master satellite, while neglecting the correlation between the slave satellites.The Master-Slave method has a reduced complexity compared to the conventional approach, where each individual correlation is described.
the full description requires k • (k − 1)/2 links.As a result, Master-Slave reduces the complexity of a k-satellite model by 2/k.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.

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 Given a master sequence, two independent conditional slave sequences are generated for the cases "master is good" and "master is bad."Different state models (1st-order Markov, semi-Markov, etc.) can be chosen individually therefore.Finally, the slave sequence is composed by piecewise assembling the two conditional sequences according to the master sequence.In this way the correlation coefficient between master and slave, the individual state probabilities, and the combined state probabilities can be accurately described.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 sampleby-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).

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.
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 ρ slaves depends on the individual correlation coefficients For comparison, the results of a first-order Markov model (using a 8 × 8 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 ρ slaves = 0.06).The higher correlation coefficients between master and slaves (ρ 12 = 0.15 and ρ 13 = 0.40) provide some correlation of the slaves as well.In Mapping 1 ρ slaves 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 ρ slaves = 0.40).Due to low correlation coefficients between the master and the slaves (ρ 12 = 0.15 and ρ 23 = 0.06), ρ slaves ≈ 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.

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 (ΔP bbb ), 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.

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 φ 1 , φ 2 , and φ 3 and the azimuth angle separations Δθ 12 , Δθ 13 , and Δθ 23 .Following the Master-Slave approach, only the correlations between master and slave are described (ρ 12 and ρ 13 ), while the correlation between the slaves (ρ 23 ) is neglected.In the following, the consequences of the neglected correlation ρ 23 are investigated.

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 φ 1 = φ 2 = 15 • ; the driving direction can be assumed as uniformly distributed between 0 and 360 • .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 Δθ Based on these results for the dual-satellite case we now analyse the consequences of neglecting the slave correlation ρ 23 in the Master-Slave approach.Following the configuration according to Figure 6, we compare ρ 23 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 Δθ 12 , Δθ 13 , and Δθ 23 .The error is given by Δρ 23 = |ρ 23, measured − ρ 23, Master-Slave |.Clearly, this error is large when the two slaves come close to each other, that is, Δθ 12 ≈ Δθ 13 -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, P bbb .It characterises the probability of service unavailability for the entire system and requires special attention, therefore.According to Figure 9, in the measured scenario P bbb ≈ 0.8 if all three satellites are colocated, and P bbb ≈ 0.6 if two of the three satellites are colocated.Considering the probability error ΔP bbb = |P bbb, measured − P bbb, Master-Slave | (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   the correlation cannot be neglected.Furthermore, the correlation should be considered for antipodal constellations, as can be seen for ΔP bbb, dual-sat + 1 independent .mappings) for the three-satellite Master-Slave approach can be determined.The following is obtained.

Influence of the Elevation
(i) Assuming constant elevation angles (represented by one subfigure), the error ΔP bbb 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., φ 1 > φ 2 and φ 3 ), then the error ΔP bbb for low slave separations is reduced with respect to φ 1 = φ 2 = φ 3 .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 ΔP bbb 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 ΔP bbb for elevation Figure 13: Worst cases for Master-Slave modelling: there is no correlation between master and slave satellites (ρ = 0), but full correlation between the slaves (ρ slaves = 1).
the higher elevation angles provide a lower probability error ΔP bbb .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 φ 1 = 15 • , φ 2 = 45 • , and φ 3 = 45 • and φ 1 = 45 • , φ 2 = 45 • , and φ 3 = 15 • .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.

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 ρ = 1 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 ΔP all bad (we define P all bad as the worst reception case of the multisatellite system.It includes P bb , P bbb , and P bbbb 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 (= P b − P all bad ) with a multisatellite system is provided if all satellites are uncorrelated (ρ = 0).In this case, the probability of the combined state is the product of the individual state probabilities: P all bad = P k b with k being the number of satellites.In contrast, a multisatellite system has no gain if all satellites are fully correlated (ρ = 1).Then it holds P all bad = P b .
Assuming the constellations in Figure 13 the Master-Slave model would provide uncorrelated slave satellites with ρ = 0.As a consequence, the Master-Slave model has k uncorrelated channels and for the combined state results P all bad = P k b (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 ΔP all bad = P 2 b − P k b .Figure 14 shows the worst-case probability error ΔP all bad with dependency on the "bad" state probability P b 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 P b (e.g., due to an increased elevation) of the single satellites provides a reduction of ΔP all bad .This situation was obtained in Figures 11 and 12 between constellations Figure 15 evaluates the worst-case probability error ΔP all bad 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 ΔP all bad .

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 (1) and slave (2) Δθ 12 and master (1) and slave (3) Δθ 13 for which the modelling error ΔP bbb exceeds a certain threshold.Assuming that ΔP bbb ≤ 0.05 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 Δθ 13 = Δθ 12 ± 20 • -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 15: Worst-case probability error of the Master-Slave approach (constellations from Figure 13) with dependency on the bad state probability of the slave satellites, the number of satellites, and the master probability P b, Master = 0.2 (a) and P b, Master = 0.8 (b).A low bad state probability of the master ( = high elevation) is beneficial.
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 ΔP bbb is below a certain threshold (0.05 was selected exemplarily).For a system with more than three satellites it means that ΔP all bad may not be simulated correctly with Master-Slave.It makes no statement about the absolute error ΔP all bad .Constellations are even possible where ΔP all bad can be neglected, although the Master-Slave constellation was found as inappropriate.

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 dualsatellite semi-Markov model as presented in [7].
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.

International Journal of Antennas and Propagation
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 dualsatellite 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.

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" (P bbb ) resulting from Master-Slave is compared with analytically estimated values from measurements.The difference for P bbb 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 worstcase 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.

Calculation of State Probabilities with One Master and
Multiple Slaves.Analogously to the two-satellite case, the combined probability vector p joint for multiple satellites (with k satellites) is calculated from the master probabilities and multiple conditional slave probabilities:

Figure 1 :
Figure 1: Signal components of a two-state LMS channel model.

Figure 4 :
Figure 4: (a) Appropriate mapping of a three-satellite constellation for the Master-Slave approach (Mapping 1)-the correlation between the slaves is low.(b) Inappropriate mapping of the three-satellite constellation (Mapping 2)-the correlation between the slaves is high.

Figure 5 :
Figure 5: Correlation coefficient (a), state probability (b), and mean state duration (c) for the three-satellite constellation mappings (Mapping 1 and Mapping 2) from Figure 4. Data was taken from measurements of 5.4 km length in a suburban environment.Resimulationgives the following results: clearly, the resulting correlation coefficient between the slave satellites deviates for Mapping 2. The joint state probability is modelled accurately in case of Mapping 1, while Mapping 2 shows deviations.The mean state duration shows deviations in all cases (notation of states: g-"good," b-"bad," ggg-"good good good," etc.).

Figure 6 :
Figure 6: Three-satellite constellation for the Master-Slave analysis with elevation angles φ, azimuth angle separations Δθ, and correlation coefficients ρ.

Figure 7 :
Figure 7: State correlation coefficient ρ between two satellites with dependency on the azimuth angle separation Δθ derived from measurements and recalculated with Master-Slave.The Master-Slave approach with Conditional Assembling Method presented in this paper resimulates accurately the correlation coefficients for the two-satellite constellation (master and slave).The elevation angles are φ 1 = φ 2 = 15 • .

Figure 9 :
Figure 9: "bad bad bad" state probabilities (P bbb ) of three satellites for different azimuth angle separations Δθ from measurements and recalculated with Master-Slave.The elevation angles are φ 1 = φ 2 = φ 3 = 15 • , where the "bad" state probability of each individual satellite is P b = 0.77.

Figure 14 :
Figure14: (a) Satellite diversity gain of uncorrelated satellites with dependency on the bad state probability of the individual satellites P b , and the number of satellites k.It represents the improvement of a system with k uncorrelated channels to a system with one channel.(b) Worst-case probability error of the Master-Slave approach (constellations from Figure13) with dependency on the bad state probability of the individual satellites, and the number of satellites (it is k uncorrelated channels minus 2 uncorrelated channels).

8 )
• P S(2)  g|g • P S(3)  g|g + P M b • P S(2)  g|b • P S(3)   g|bP M g • P S(2) g|g • P S(3) b|g + P M b • P S(2) g|b • P S(3) b|b P M g • P S(2) b|g • P S(3) g|g + P M b • P S(2) b|b • P S(3) g|b P M g • P S(2) b|g • P S(3) b|g + P M b • P S(2) b|b • P S(3By inserting (A.8) in (A.7) the correlation coefficient between the slaves ρ states, Slaves can be described as a function of the correlation coefficients between master and slave (2) and master and slave (3):ρ states, Slaves = ρ states, Master, Slave2 • ρ states, Master, Slave3 .(A.9) 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:φ = n 11 • n 22 − n 12 • n 21 (n 1X )(n 2X )(n X1 )(n X2 ) .(A.5)n i j isthe number of samples where sequence 1 has state i and sequence 2 has state j.For state modelling this formula can be adapted by using the joint state probabilities of two satellites: ρ states = P gg • P bb − P gb • P bg P gg + P gb P bb + P bg P gg + P bg P bb + P gb Correlation Coefficient between Two Slaves.The correlation coefficient between the slaves depends on the combined state probabilities P i j of the slaves: ρ states, Slaves = P gg • P bb − P gb • P bg P gg + P gb P bb + P bg P gg + P bg P bb + P gbThe combined state probabilities P i j of the slaves can be calculated from the joint state probabilities of the whole system Master (1)-Slave (2)-Slave (3).According to (A.4) it holds: 3)b|b •, . .., •PS(k)