Major challenge of 4G/5G mobile cellular network is to enhance the carried traffic capacity. One possible solution is offloading of mobile data traffic from the primary network access to the WiFi infrastructure (for example femtocell of LTE). In  and its extended version , the authors show the performance of WiFi Direct Device to Device (D2D) technology and demonstrate how cellular traffic can be effectively offloaded through a WiFi and provides the estimated gains in energy efficiency and capacity from such offloading. Another big problem of outdoor LTE network is its latency, which directly affects its throughput. In LTE eNB maintains a large queue to reduce overflow of burst traffic at the expense of packet delay. In  “Adaptive Receiver-Window Adjustment Algorithm” is used to reduce the latency. Variation of throughput and latency against “modulation and coding scheme” (MCS) and number of users are shown explicitly. Still there is a scope to work on the paper, including other traffic parameters like blocking probability and carried traffic. In , instead of conventional traffic of user data, the paper proposed a M2M traffic model to reduce the signaling traffic and bandwidth utilization. In  the traffic of LTE-WiMAX integrated network is analyzed including the algorithm of traffic simulation. The variation of FTP, HTTP and voice traffic with time is compared explicitly, where the interference is found as the prime factor. Although femtocell offloads traffic of macro cell but the interference of femto and macro cell can deteriorate throughput of the network. The solution, to limit the interference of femtocells due to their uncontrolled placement is found in  .
The redundant bits of channel coding degrade the throughput of a network. The compromise between source code rate and channel code rate to minimize the number of redundant bits generated by channel coding with an acceptable Mean Opinion Score (MOS) is analyzed in . Here the authors proposed a dynamic adaptation algorithm of joint source-channel code rate for enhancing voice transmission over LTE network. Paper  surveys state-of-art mobile data offloading and proposed to use integrated Femto and Wi-Fi to have improved “service performance” and “interference management” with the constraint of non-overlapping femtocells. The traffic model of low dense femtocell network under fading condition is found in  considering the concept of  with variable BW.
The rest of the paper is organized as follows: Section 2 provides previous analysis of traffic of LTE network; Section 3 gives the analytical model of Macro-Femto Cellular Network based on Markov chain and probability tree, MIMO channel under fading condition is combined with traffic parameter and traffic model of variable bandwidth (BW). Section 4 provides results based on the analysis of Section 3 and verifies it by simulation and finally, section 5 concludes the entire analysis.
2. Previous Work
In  the traffic of macro-femto cellular communication system of LTE is modeled by Markov chain. The authors first determined the expression of total call arrival rate of macro and femtocell then separate Markov chains were drawn for macro and femtocell traffic.
The total arrival rate of a macrocell is: ; where λom is the newly originating call at a macrocell, λh,mm→macro to macro handover rate, PBfλof→a call originated in a femtocell and found that all the channels are occupied then the call will use the channel of a macrocell, αPDfλh,ff→after a femto to femto handover a UE (user equipment) found that SINR > τ (the probability of getting SINR > τ is α) but all the channels of the new femtocell are occupied then the channel of a macrocell will be used, (1 − α)λh,ff→the case like before but SINR < τ hence the channel of a macrocell will be used.
The total arrival rate of a femtocell, , where λof→the arrival rate of newly originated call in femtocells, λh,mf→the rate of handover arrival from macro to femtocells, αλh,ff→the femto to femto handover provided that SINR > τ, where τ is the threshold SINR to maintain communication and PD,mβλh,ff→the femto to femto handover provided that received SINR lies between two thresholds τ1 and τ2 and .
Traffic offloading strategy of LTE network is analyzed in  using two-dimensional Markov chain, where the product of session arrival and session duration is considered offered traffic. Here one dimension is Licensed Assisted LTE Access (LLA) traffic and the other is WiFi AP traffic. The authors consider CTMC with probability state: (NW(t), NL(t)), where NW(t) is the number of WiFi sessions and NL(t) is the number of offloaded LTE sessions at time t > 0. The similar traffic model is found in , applied for cognitive radio network (CRN), where arrival rate of PU and SU with their termination rate is considered as the offered traffic. Both models were unable to include the variable bandwidth in the state transition chain.
3. System Model
3.1. Traffic Model of Macro-Femto Cellular Network
In micro-macro cellular system the microcells under a macrocell form the coverage in continuous fashion but the situation is different in LTE with macro-femto combination since femtocells do not form a continuous coverage under a macrocell. Hence the concept of take-back of   is not applicable here, visualized from Figure 1. Let the area of the test macrocell is Am and the
Figure 1. Macro-femto integrated cellular model.
sum of the area of femtocells within the test macro cell is Af. Let the average call arrival rate inside the test cell is λT. The arrival rate of macro and femto cells will be distributed as:
Arrival rate of macrocell is and that of femtocells, In micro-macro cellular system the overflow traffic from microcell is related to macrocell but such concept is not applicable in femto-macro cell hence combined state transition chain of the papers   cannot be applied here.
We consider the separates Markov chain of femto and macrocell of  ; where n is the number of channels (combined data rate/individual user’s rate), N is the total number of femtocells within the test macrocell, λf is the call arrival rate within the femtocells and is the handover arrival rate from surrounding cells as shown in Figure 2.
The blocking probability under a femtocell is derived as  :
, where (1)
Similarly, the call blocking probability under a macrocell is,
, where (2)
and l is the number of channels of the test macrocell, λm is the call arrival rate within the test macrocell and λh is the handover arrival rate from the surrounding cells. The above chain can be improved considering call admission control scheme, where m channels are reserved for handoff calls from the surrounding femtocells. The state transition chain of a macro cell under call admission control is shown in Figure 3.
Figure 2. Markov chain of a test femtocell.
Figure 3. Markov chain of a macrocell under reserved channels for handover calls.
Applying cut equation like  we get in the generalized form as:
The probability of FT is,
The entire sample space satisfies,
Applying normalization condition,
Probability of dropping of newly originating call is Pk and probability of forced termination (FT) is Pk+m. Because of use of low dense femtocells, we assume that two micro and two femtocells are enough to complete a session. In this context we have to determine the probability that an MS crosses the cell border of a femtocell before completing its session as  ,
where , v is the velocity of MS, Rf is the radius of femtocell and is the call terminator rate at femtocell. Similarly, for macrocell, the corresponding probability is expressed as,
We design a probability tree of FT for a MS, i starts its session in a femtocell i as shown in Figure 4 to relate traffic parameters of (1) to (8). The probability of FT, PFTf is derived as,
Figure 4. Probability tree of PFTf.
Another probability tree is designed like Figure 5, when a MS, j starts its session in macro cell I and the corresponding probability of FT is found as:
Solving the linear Equation (9) and Equation (10) we can easily obtain PFTm and PFTf. Now we have to relate the traffic parameters of this section with small scale fading to get realistic scenario of traffic of wireless network. To combat the fading affect we consider MIMO link discussed in next subsection.
3.2. The MIMO Channel Model of LTE
For successful communication both availability of channel of a cell and instantaneous SNR ≥ τ should be ensured; where τ is the threshold value of instantaneous SNR.
The minimum SNR to maintain QoS of the MIMO link is expressed as   :
where is the Nt × Nr MIMO channel matrix, Nt is the number of antenna elements of transmit antenna array, Nr is that of the receiver antenna array and is the sum of diagonal components of the matrix .
The trace of can be expressed as,
where r is the random variable of channel gain between the antenna elements. The detail derivation of Equation (11) and Equation (12) are given in .
From Equation (11) and Equation (12) we have, for STBC.
Let us now concentrate on instantaneous and average SNR γ and γav of the MIMO link. First of all, we analyze the MIMO link under TAS (Transmit Antenna
Figure 5. Probability tree of PFTm.
Selection) technique then for the complete MIMO. The channel gain vector, selecting the kth antenna element at the transmitting end is,
Frobenious norm of is,
The transmit antenna element is selected for maximum value of ;
The received signal vector is,
where is the noise vector at the receiving end and P is the transmitted signal power.
The selection combiner at the receiving we choose,
when MRC is used then the weighting vector ,
The output of MRC at the receiver,
The instantaneous SNR of receiver is,
If nb is the variance of the noise vector then,
The weighting vector under TAS of Equation (17) is the contribution of the kth transmit antenna element, hence under the complete MIMO the Equation (17) will be modified as,
The pdf of the random variable, γ under MRC of the MIMO for Rayleigh fading is like ,
that of the Nakagami-m case of  is,
and the outage probability,
where τ is the threshold SNR and it depends on the number of received antennas and the variance of the channel gain of the antenna elements.
Now the probability of FT, considering the outage probability will be,
The Equation (25) and Equation (26) provides the probability of FT, which depends on the instantaneous SNR, the number of received antenna elements and all traffic parameters will be shown in the result section. In LTE, users deal with traffic of variable BW, hence we have to modify the Markov chain of subsection 3.1. In next subsection we proposed a new state transition chain of variable BW.
3.3. Traffic of Variable Bandwidth
In LTE users can take bandwidth on demand. Here we assume the minimum unit of BW of a user is a pair of RB (resource block). In LTE each frame of 10 ms is divided into 10 subframe of 1 ms, also each sub-frame is divided into two time slots of 0.5 ms each. Each slot contains either six or seven OFDM symbols, depending on the type of channel coding schemes. When normal Cyclic Prefix (CP) is used then a RB contains 7 symbols under each subcarrier of OFDM, under extended CP, 6 symbols under each sub-carrier. Each RB contains 12 subcarriers, hence 12 × 6 = 72 symbols in a RB of the extended CP and 12 × 7 = 84 symbols in a RB of the normal CP.
If a user can take the highest BW of 3 RB then we can represent the traffic by the state transition chain of Figure 6, where RB arrival rate is λ calls/unit time and its termination rate is µ. If k is the highest allowed BW or RB then each probability state will emit k arrivals and maximum possible terminate rate will
Figure 6. Markov chain of LTE user with variable BW.
be kμ. Applying node equation on each node we get 9 linear equations. Using normalization of entire sample space, we get the following relations.
Applying node equation on nodes of Figure 6 we get,
for , m is the last state of the chain.
For entire sampling space,
If normalized probability state of the chain is qx then, the blocking probability of a user uses k RB will be,
Combination of analysis of three subsections: 3.1, 3.2 and 3.3 provide the exact traffic model of a LTE network.
This section provides the analytical results based on the statistical model of Section 3 and verified by simulation. Taking the following traffic parameters: µm = µf = 1.2 calls/min, vf = 10 m/s, Rf = 20 m, vm = 100 m/min, Rm = 1 Km and Bf = 0.02, we obtain the probability of handover to femto and macrocells as: Phf = 0.21 and Phm = 0.05 respectively. Figure 7 shows the profile of probability of FT (both the cases i.e. when a UE starts its session at macro and femto cell) against the blocking probability of macrocell taking the velocity of users at macro and femtocell as the parameters. The probability of FT increases with increase in both Bm and velocity of the users. Again taking, µf = 1.8 calls/min, v = 10 m/min, Rf = 20 m, we get Phf = 0.15. For another set of traffic parameters: µf = 1.2 calls/min, v = 12.5 m/min, Rf = 20 m, we get Phf = 0.25. Considering above two sets of parameters, the variation of FT against arrival rate of femtocell is shown in Figure 8.
Figure 7. Variation of PFT against Bm.
Figure 8. Variation of PFT against λf.
Almost similar profile is found as Figure 7 i.e. the PFT increases with increase in both the λf and Phf.
Next, the probability of FT is varied against the number of channels of femtocell taking the same traffic parameters shown in Figure 9. In this case, the probability of FT decreases with increase in the number of channels of femtocell. All the analytical results are verified by simulation which provides more than 95% confidence level.
The next part of the result section is related to subsection 3.2; where the traffic parameters are combined with fading parameters of MIMO channel to measure the performance of the network. The variation of PFT against the number of channels of femtocell is shown in Figure 10(a) for average SNR of −2 dB and (b)
Figure 9. Variation of probability of FT against the number of channels of femtocell.
Figure 10. Variation of FT against the number of channels of femtocell.
for average SNR of 2 dB under STBC taking the probability of handover of femtocell as a parameter. The PFT decreases with increase in the number of channels of femtocell but increases with increase in the probability of handover of femtocells. The channel parameters are taken as: average SNR = −2 dB and 2 dB, the standard deviation of channel gain σ = 0.83, the parameters of Nakagami-m fading m = 2, Nr = 2. The traffic parameters are taken as: number of femto cells Nf = 12, λf = 9 calls/min, λhf = 4 calls/min, λm = 20 calls/min, λhm = 12 calls/min, µm = 4.2 calls/min and µf = 2 calls/min. The performance is found better under Nakagami-m fading case compared to the Rayleigh fading environment because of m direct link (between transmitter and receiver) of Nakagami-m fading. Under Rayleigh fading environment three is no direct link between transmitter and receiver hence performance is found very poor at low SNR but when the received
Figure 11. Comparison of proposed LTE traffic.
SNR attains above a threshold the performance of Rayleigh and Nakagami-m environment are merged. The difference between performance at low (−2 dB) and high (2 dB) SNR is found distinct in Figure 10(a) and Figure 10(b).
Finally, the call blocking probability of variable BW traffic of LTE and Erlang’s model are compared in Figure 11. Here Blocking probability is varied against call arrival rate taking termination rate as a parameter provided the number of channels n = 8. The blocking probability of variable BW traffic of LTE is found smaller than the case of Erlang’s traffic. The arrival rate of the state transition chain of Figure 6 is λ for any transition, hence each state experiences the total arrival of 3λ, on the other hand the termination rate from probability state Px is 3 xµ. The termination rate of the proposed model is more prominent than the Erlang’s model of M/M/n/K; hence blocking probability of the proposed model is found smaller than the conventional Erlag’s model.
The paper combines the traffic model with the outage probability of wireless link under two different small scale fading cases. A new traffic model of variable BW (variable RB for up or down link) is proposed based on Markov chain and comparison is made with M/M/n/K traffic. In future we will combine the wireless link parameters with our proposed traffic model of variable BW. In this paper we consider the pdf of SNR in determining the outage probability but we have the scope to use the pdf of SIR (signal to interference ratio) or SINR to evaluate the same parameter. In this case SINR will be a random variable as the ratio of another two random variables; hence the concept of mixed random variable can be applied for it. We consider M/M/n/K traffic model throughout the paper but in voice data integrated traffic we can consider M/G/1/K traffic model as well.
 Pyattaev, A., Johnsson, K., Andreev, S. and Koucheryavy, Y. (2013) 3GPP LTE Traffic Offloading onto WiFi Direct. IEEE WCNC Workshop on Mobile Internet: Traffic Modeling, Subscriber Perception Analysis and Traffic Aware Network Design, Shanghai, 7-10 April 2013, 135-140.
 Hagos, D.H. (2015) The Performance of Network-Controlled Mobile Data Offloading from LTE to WiFi Networks. Springer Science + Business Media, New York, 675-694.
 Adesh, N.D. and Renuka, A. (2019) Adaptive Receiver-Window Adjustment for Delay Reduction in LTE Networks. Journal of Computer Networks and Communications, 2019, Article ID:
 Joseph, J., Ruhl, R. and Samuel, H. (2018) Reducing Traffic Congestion for Machine to Machine Type Communication over 4G-LTE Network by Decreasing Total Bytes Transmitted. 17th IEEE International Conference on Trust, Security and Privacy in Computing and Communications/12th IEEE International Conference on Big Data Science and Engineering, New York, 1-3 August 2018, 1720-1721.
 Kankare, M., Asp, A., Sydorov, Y., Niemel?, J. and Valkama, M. (2013) Large-Scale Femtocell Network Deployment and Measurements. 6th International Workshop on Multiple Access Communications, Lithuania, December 2013, 100-112.
 Nguyen, D.-H. and Nguyen, H. (2016) A Dynamic Rate Adaptation Algorithm Using WB E-Model for Voice Traffic over LTE Network. IEEE Conference on Wireless Days, Toulouse, 1-6 March 2016, 1-6.
 Pawar, A.R., Bhardwaj, S.S. and Wandre, S.N. (2015) Mobile Data Offloading Techniques and Related Issues. International Journal of Advanced Research in Computer Engineering & Technology, 4, 1367-1371.
 Ruiz, J.M., Matías, T. and Luna, S. (2015) A Femtocell Location Strategy for Improving Adaptive Traffic Sharing in Heterogeneous LTE Networks. EURASIP Journal on Wireless Communications and Networking, 2015, Article No. 38.
 Chowdhury, M.Z. and Jang, Y.M. (2013) Handover Management in High-Dense Femtocellular Networks. EURASIP Journal on Wireless Communications and Networking, 2013, Article No. 6.
 Markova, E., Moltchanov, D., Gudkova, I., Samouylov, K. and Koucharyavy, Y. (2019) Performance Assessment of QoS-Aware LTE Sessions Offloading onto LAA/WiFi Systems. IEEE Access, 7, 36300-36311.
 Islam, I., Rabbi, F., Khan, R.T. and Akhter, J. (2015) Performance Evaluation of Cognitive Radio Network Based on 2-D Markov Chain. Journal of Telecommunications and Information Technology, 3, 39-44.
 Chaturvedi, A. and Solanki, A. (2012) Analyzing Area Spectral Efficiency and Its Impacts on the Total Spectrum Requirements for IMT Advanced. International Journal of Distributed and Parallel Systems, 3, 141-152.
 Oh, E. (2015) Distributed Cell Configuration Strategy for Macro/Micro Overlaid Hierarchical Cellular Systems. International Journal of Software Engineering and Its Applications, 9, 177-188.
 Jun, K. and Kim, S. (2014) Modeling and Per-formance Analysis for Cell Access and Handoff Schemes in Two-Tier Cellular Networks. Mathematical Problems in Engineering, 2014, Article ID: 860729.
 Hossain, M.A., Islam, M.I. and Amin, M.R. (2016) Load Sensitive Power Saving Technique for 4G Mobile Network under Limited User Traffic. Communications and Network, 8, 79-87.
 Ortigoza-Guerrero, L. and Hamid Aghvami, A. (1999) A Prioritized Handoff Dynamic Channel Allocation Strategy for PCS. IEEE Transaction on Vehicular Technology, 48, 1203-1215.
 Hashem, T. and Islam, M.I. (2014) Performance Analysis of MIMO Link under Fading Channels. 17th International Conference on Computer and Information Technology, Dhaka, December 2014, 498-503.
 Da Costa, D.B., Ferdinand, N.S., Dias, U.S., de Sousa Jr., R.T. and Latva-aho, M. (2016) Secrecy Outage Performance of MIMO Wiretap Channels with Multiple Jamming Signals. Journal of Communications and Information Systems, 31, 30-40.
 Li, Y.Z., Zhang, L., Cimini Jr., L.J. and Zhang, H.L. (2011) Statistical Analysis of MIMO Beamforming with Co-Channel Unequal-Power MIMO Interferers under Path-Loss and Rayleigh Fading. IEEE Transactions on Signal Processing, 59, 3738-3748.
 Hung, C.-C., Chiang, C.-T., Wu, R.-C. and Lee, C. (2012) Positive Effect of Severe Nakagami-m Fading on the Performance of Multiuser TAS/MRC Systems with High Selection Gain. International Journal of Antennas and Propagation, 2012, Article ID: 987631.