Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the network - layer performance analysis seeks to provide estimates on the delays experienced by traffic traversing the elements of a network , as well as the corresponding buffer requirements . for wireless networks , a question of interest is how the stochastic properties of wireless channels impact delay and backlog performance . wireless channels are characterized by rapid variation of channel quality caused by the mobility and location of communicating devices . this is due to _ fading _ , which is the deviation in the attenuation experienced by the transmitted signal when traversing a wireless channel . the term _ fading channel _ is used to refer to a channel that experiences such effects . in this paper we explore the network - layer performance of a multi - hop network where each link is represented by a fading channel . we model the wireless network by tandem queues with variable capacity servers , where each server expresses the random capacity of a fading channel . we ignore the impact of coding by assuming that transmission rates over the fading channels are equal to their information - theoretic capacity limit , @xmath0 , which is expressed as a function of the instantaneous signal - to - noise ratio ( snr ) at the receiver , @xmath1 , by @xmath2 , where @xmath3 is the channel bandwidth ( in _ hz _ ) . numerous models are available to describe the gain of fading channels depending on the type of fading ( slow or fast ) , and the environment ( e.g. , urban or rural ) . the instantaneous , information - theoretic channel capacity of a fading channel can be represented as the logarithm of @xmath1 by ( see chp .  14.2 in @xcite ) @xmath4 where @xmath5 is a constant and the function @xmath6 is used to characterize the fading channel . we are interested in finding bounds on the end - to - end delay and on buffer requirements for a cascade of fading channels , with store - and - forward processing at each channel . the analysis in this paper takes a system - theoretic stochastic network calculus approach @xcite , which describes the network properties using a @xmath7 dioid algebra . arrival and departure processes at a network element are described by bivariate stochastic processes @xmath8 and @xmath9 , respectively , denoting the cumulative arrivals ( departures ) in the time interval @xmath10 . a network element is characterized by the service process @xmath11 , denoting the available service in @xmath10 . the input - output relationship at the network element is governed by @xmath12 where the convolution operator ` @xmath13 ' is defined as @xmath14 . if network traffic passes through a tandem of @xmath15 network elements with service processes @xmath16 , the service of the network as a whole can be expressed by the convolution @xmath17 . the stochastic properties of fading channels present a formidable challenge for a network - layer analysis since the service processes corresponding to the channel capacity of common fading channel models such as rician , rayleigh , or nakagami-@xmath18 , require to take a logarithm of their distributions . as discussed in the next section , researchers frequently turn to higher - layer abstractions to model fading channels . widely used abstraction are the two - state gilbert - elliott model and its extensions to finite - state markov channels ( fsmc ) @xcite . fsmc models simplify the analysis to a degree that the network model becomes tractable , at least at a single node . extensions to multi - hop settings encounter a rapidly growing state space . as of today , a general multihop analysis that is applicable to models of fading channels , such as rician , rayleigh , or nakagami-@xmath18 , remains open .    in this paper , we pursue a novel approach to the analysis of multi - hop wireless networks . we develop a calculus for wireless networks that can be applied to fading channel models from the wireless communication literature to provide network - layer performance bounds . we view the network - layer model with arrival , departure and service processes as residing in a _ bit domain _ , where traffic and service is measured in bits . we view the fading channel models used in wireless communications as residing in an alternate domain , which we call the _ snr domain _ , where channel properties are expressed in terms of the distribution of the signal - to - noise ratio at the receiver . we then derive a method to compute performance bounds from these traffic and service characterizations . a key observation in our work is that service elements in the snr domain obey the laws of a dioid algebra . we devise a suitable dioid , referred to as _  algebra _ , where the minimum takes the role of the standard addition , and the second operation is the usual multiplication , and use it for analysis in the snr domain . in this domain multi - hop descriptions of fading channels become tractable . in particular , we find that a cascade of fading channels can be expressed in terms of a convolution in the new algebra of the constituting channels . the key to our analysis is that we derive performance bounds entirely in the snr domain . observing that the bit and snr domains are linked by the exponential function , we transfer arrival and departure processes from the bit to the snr domain . then , we derive backlog and delay bounds in the transfer domain using the  algebra . the results are mapped back to the original bit domain to finally give us the desired performance bounds . our derivations in the snr domain require the computation of products and quotients of random variables . here , we take advantage of the _ mellin transform _ to facilitate otherwise cumbersome calculations . then , the computational problem is reduced to finding the mellin transform for service and traffic processes . the main contribution of this paper is the development of a framework for studying the impact of channel gain models on the network - layer performance of wireless networks . for the purposes of this paper , the snr domain is used solely as a transfer domain that enables us to solve an otherwise intractable mathematical problem . on the other hand , the ability to map quantities that appear in network - layer models and concepts found in a physical - layer analysis may prove useful in a broader context , e.g. , for studying cross - layer performance issues in wireless communications . moreover , the  algebra and the mellin transform form a tool set that can be applied more generally in wireless communications for studying the channel gain of cascades of fading channels . as the first paper on the  network calculus algebra , our paper only considers simple network scenarios and makes numerous convenient assumptions ( which are made explicit in sec . [ sec : model ] ) . there is room for significant future work on extensions of the model and a relaxation of the presented assumptions . the remainder of the paper is organized as follows . in sec . [ sec : related ] we discuss related work . we describe the system model in sec .  [ sec : model ] , where we also motivate the use of the snr domain . in sec . [ sec : calculus ] we present the  algebra and derive performance bounds . in sec . [ sec : rayleigh ] we apply the analysis to a cascade of rayleigh channels , and present numerical examples . in sec . [ sec : rayleigh ] , we investigate a network with cross traffic at each node . we discuss brief conclusions in sec . [ sec : conclusions ] . analytical approaches for network - layer performance analysis of wireless networks include queueing theory , effective bandwidth and , more recently , network calculus . since the service processes corresponding to the channel capacity of common fading channel models such as rician , rayleigh , or nakagami-@xmath18 , require to take a logarithm of their distributions , researchers often turn to higher - layer abstractions to model fading channels , which lend themselves more easily to an analysis . a widely used abstraction is the two - state channel model developed by gilbert @xcite and elliott @xcite , and subsequent extensions to a finite - state markov channel ( fsmc ) @xcite . markov channel models are well suited to express the time correlation of fading channel samples . we refer to @xcite for a survey of the development and applications of fsmc models . zorzi et al . @xcite evaluated the accuracy of first - order markov channel models of fading channels , where the next channel sample depends only on the current state of the markov process , and higher order processes that can capture memory extending further back in the process history . the authors found that a first - order markov model is a good approximation of the fading channel , and that using higher order markov processes does not significantly improve the accuracy of the model . queue - based channel ( qbc ) @xcite is an alternative model for fading channels , which models a binary additive noise channel with memory based on a finite queue . here , a queue with size @xmath19 contains the last @xmath19 noise symbols , and the noise process is an @xmath20-order markov chain . the model was found to provide a better approximation to the rayleigh and rician slow fading channels compared to the gilbert - elliot model @xcite . an extension of the qbc model , called weighted qbc @xcite permit queue cells ( i.e. , channel samples ) to contribute with different weights to the noise process . queueing theoretic studies of fading channels generally apply approximations to reduce the complexity of multi - hop models . le , nguyen , and hossain @xcite apply a decomposition approximation to analyze the loss probability and average delay of a multi - hop wireless network with slotted transmissions for a batch bernoulli arrival process , and with independent cross traffic at each node . the wireless link is assumed to employ adaptive modulation and coding with multiple modes , where each mode corresponds to a given link rate , which is chosen based on the snr of the channel . the channel state of a link is assumed to be stationary , and channel states in consecutive time slots are independent . another decomposition approximation is presented by le and hossain @xcite , who consider a multi - hop tandem network with a batch arrival process and multi - rate transmissions , to develop a routing scheme that can meet given delay and loss requirements . the analysis obtains end - to - end loss rates and delays with a decomposition analysis , and feeds the results as metrics to the routing algorithm . bisnik and abouzeid @xcite model a multi - hop wireless network as an open network of g / g/1 queueing systems . using diffusion approximation , they obtain closed - form expressions for average end - to - end delays . ishizaki and hwang @xcite studied the impact of multiuser diversity assisted packet scheduling on the packet delay performance in a wireless network with nakagami - m fading channels . the network is modeled by a multi - queue system , where each channel is described by an fsmc . under assumptions of stationarity , homogeneity , and independence of the channel processes , they approximate the tail distribution of packet delays , and compare it to that of a round - robin scheduler . the results indicate that the delay performance of multiuser diversity assisted scheduling algorithms is not necessarily superior to that of round - robin scheduling . since the application of queuing theoretical methods to study cascade of fading channels requires many assumptions on arrival and service distributions , and simplifications of the model , the use of classical queueing theoretic methods for the performance analysis of multi - hop wireless networks has been put into question @xcite . an effective bandwidth @xcite analysis seeks to develop ( asymptotic ) bounds on performance metrics , e.g. , an exponential decay of the backlog . wu and negi @xcite have adapted an effective bandwidth analysis to the analysis of fading channels . they introduce the concept of effective capacity , which characterizes a wireless channel by a log - moment generating function ( log - mgf ) of the channel capacity . they obtain an asymptotic approximation of the delay bound violation probability of a rayleigh fading channel . due to the difficulty of computing the moment generating function ( mgf ) of the rayleigh distribution , they assume non - correlated distributions with low snr and estimate channel parameters from measurements . the work has been extended to correlated rayleigh and correlated nakagami-@xmath18 channels , and to cascades of fading channels @xcite . a closely related concept is the effective channel capacity presented by li et al . @xcite , which describes the available channel capacity by a first order markov chain and computes the log - mgf of the underlying markov process . taking advantage of methods developed in @xcite , they compute statistical delay bounds for nakagami - m fading channel . hassan , krunz , and matta @xcite use an effective bandwidth analysis to study delay and loss performance at a single wireless link , which is modeled by an fsmc . for fluid on - off traffic and fifo buffering , they obtain a closed form expression for the effective bandwidth required to guarantee bounds on delay and packet loss . there is a collection of recent works that apply stochastic network calculus methods @xcite to wireless networks with fading channels . the stochastic network calculus is closely related to the effective bandwidth theory , in that it seeks to develop bounds on performance metrics under assumptions also found in the effective bandwidth literature . different from effective bandwidth literature , stochastic network calculus methods seek to develop non - asymptotic bounds . an attractive element of a network calculus analysis is that sometimes it is possible to extend a single node analysis to a tandem of nodes , using the convolution operation seen in the introduction . fidler @xcite presents a network calculus methodology for a two - state fsmc model of a single - hop fading channel . he applies the mgf network calculus , which was suggested in the problem sets of chapter  7 in @xcite , and which has been developed in @xcite . the mgf network calculus takes its name from the extensive use of moment generating functions in the derivation of performance bounds . mahmood , rizk , and jiang @xcite apply the mgf network calculus to mimo channels and derive delay bounds for periodic traffic sources . @xcite also use an mgf network calculus to study the performance of two - hop relay networks . a similar methodology is applied by mahmod , vehkaper , and jiang @xcite to compute the throughput of a multi - user ds - cdma system with delay constraints . in the works above that apply the mgf network calculus , models for a cascade of multiple fading channels become complex , so that multi - node results for networks with more than two nodes have not been obtained . the  network calculus developed in this paper uses similar descriptions and assumptions for traffic and service as the mgf network calculus . by performing computations in a transfer domain , where fading channel models take a simpler form , we are able to compute multi - node service descriptions for an arbitrarily large number of nodes . the mgf network calculus assumes that arrivals and service at each node are independent . these assumptions can be relaxed using statistical envelope descriptions for traffic ( effective envelopes ) and service ( statistical service curve ) @xcite . jiang and emstad @xcite have applied an approach with envelopes to a fading channel where the wireless channel is characterized by two stochastic processes : an ideal service process and an impairment process , where the impairment process captures effects due to fading , noise , and cross traffic . verticale and giacomazzi @xcite obtain a closed form expression for the variance of a service curve , which describes the available service by a markov chain . this is used for the analysis of an fsmc model of a rayleigh fading channel . for computing the bounds for markovian arrivals , they apply the bounded - variance network calculus introduced in @xcite , which is an extension of the central limit theorem methods by choe and shroff @xcite to multi - hop paths . verticale @xcite has applied the same methodology to constant bit rate traffic . ciucu , pan , and hohlfeld @xcite and ciucu @xcite present non - asymptotic ( i.e. , finite number of hops ) , closed - form expressions for the delay and throughput distributions for multi - hop wireless networks . as many of the works discussed above , the fading channel is modeled by an abstraction that uses a link layer model of the transmission channel . here , the channel is assumed to be governed by a slotted - aloha system in half - duplex mode . the model of this channel is a two - state on - off server , where a node can transmit ( i.e. , is in the on state ) only when other all nodes in the interference range are not transmitting . there is also a literature on physical - layer performance metrics of fading channels in multi - hop wireless networks . hasna and alouini @xcite have presented a framework for evaluating the end - to - end outage probability of a multi - hop wireless relay network with independent , non - regenerative relays , i.e. , amplify - and - forward ( af ) , over nakagami fading channels . they provide a closed - form expression for the mgf of the reciprocal of the equivalent end - to - end snr for independent nakagami fading channels . for the same af relay network , tsiftsis @xcite obtains a closed - form bound for the average error probability . this bound is reportedly tight at low snr , but may become loose for higher snr values and for more severe fading environment , e.g. , rayleigh fading . amarasuriya , tellambura , and ardakani @xcite present an alternative bound to @xcite on the end - to - end snr in multi - hop af relay networks . they derive the distribution function and the mgf for i.i.d . nakagami - m fading and for independent , but non - identically distributed rayleigh fading . the works above study physical layer performance bounds of channel - assisted , amplify - and - forward relaying over a multi - hop fading channels . they do not consider buffering or traffic burstiness , and are not concerned with network performance metrics addressed in this paper . delay and backlog analysis and optimization of multihop wireless networks remains an open research problem @xcite . we consider a wireless @xmath15-node tandem network as shown in fig . [ fig : tandem ] , where each node is modeled by a server with an infinite buffer . we are interested in the performance experienced by a ( through ) flow that traverses the entire network and may encounter cross traffic at each node . one can think of the cross traffic at a node as the aggregate of all traffic traversing the node that does not belong to the through flow . the service given to the through flow at a node is a random process , which is governed by the instantaneous channel capacity as well as the cross traffic at the node . we consider a fluid - flow traffic model where the flow is infinitely divisible . we will work in a discrete - time domain @xmath21 , where @xmath22 is the set of integers and @xmath23 is length of the time unit . setting @xmath24 allows us to replace @xmath25 by @xmath26 , which we interpret as the index of a time slot . we assume that the system is started with empty queues at time @xmath27 . different nodes and different traffic flows will be distinguished by subscripts . the cumulative arrivals to , the service offered by , and the departures from the node are represented by random processes @xmath28 , @xmath29 , and @xmath30 that will be described more precisely below , with @xmath31 for @xmath32 . we denote by @xmath33 and @xmath34 the arrivals to and the departures from the tandem network . throughout , we assume that arrival and service processes satisfy stationary bounds . consider for the moment a single node . dropping subscripts , we write @xmath35 for the cumulative arrivals and departures , respectively , at the node in the time interval @xmath36 , where @xmath37 denotes the arrivals and @xmath38 the departures in the @xmath26-th time slot . due to causality , we have @xmath39 . the processes lie in the set @xmath40 of non - negative bivariate functions @xmath41 that are increasing in the second argument and vanish unless @xmath42 . the backlog at time @xmath43 is given by @xmath44 and the delay at the node is given by @xmath45    the service of the node in the time interval @xmath36 is given by a random process @xmath46 , such that eq . holds for every arrival process @xmath47 and the corresponding departure process @xmath48 . this service description with bivariate functions is referred to as _ dynamic server_. initially defined for non - random service  @xcite , dynamic servers have been extended to random processes in @xcite . the above model is a typical network - layer model . where traffic is measured in bits , and service is measured in bits per second . we thus refer to this model of arrivals , departures , and service as residing in a _ bit domain_.    the network calculus exploits that networks which satisfy the input - output relation of eq .   with equality can be viewed as linear systems in a dioid algebra @xcite . in the ( @xmath49 ) dioid , the minimum and addition take the place of the standard addition and multiplication operations . the network calculus is based on the fact that @xmath50 is again a dioid  @xcite . note that the min - plus convolution , which provides the second operation in the dioid , is not commutative in @xmath40 .      to compute a service model for a wireless channel , we assume that the channel state information is sampled at equal time intervals @xmath23 . with @xmath51 , let @xmath52 denote the instantaneous signal - to - noise ratio observed at the receiver in the @xmath26-th sampling epoch . then , @xmath52 is a nonnegative random variable that has the probability distribution of the underlying fading model . we assume that the random variables @xmath52 are independent and identically distributed . this assumption is justified when @xmath23 is longer than the channel coherence time . otherwise , the assumption will give optimistic bounds . we emphasize that the network calculus in this paper applies to settings without independence , however , the derivation of performance bounds will proceed differently . using eq . , the instantaneous service offered by the channel in the @xmath26-th slot is given by @xmath53 . and the corresponding service process is given by @xmath54 where we haven chosen units such that the constant in eq . takes the value @xmath55 .    the service description in eq . requires us to work with the logarithm of fading distributions , which presents a non - trivial technical difficulty via the usual network calculus or queueing theory . on the other hand , observe that the exponential @xmath56 is described more simply by @xmath57 this motivates the development of a system model that allows us to exploit the more tractable service representation in eq .  . in this alternative model , arrivals , departures , and service reside in a different domain , where we can work directly with the distribution functions of the fading channel gain and the corresponding snr at the receiver . we now proceed by mapping the network model from fig . [ fig : tandem ] into a transfer domain , which we refer to as _ snr domain_. we will seek to derive performance bounds in the transfer domain , and then map the results to the bit domain to obtain network - layer bounds for backlog and delays . the relationship of the network models in bit domain and snr domain is illustrated in fig . [ fig : domains ] .    in the previous subsection , we constructed the service process for a wireless link in the snr domain in eq . as @xmath58 by analogy , we describe the arrivals and departures in the snr domain by @xmath59 throughout this paper , we use calligraphic upper - case letters to represent processes that characterize traffic or service as a function of the instantaneous snr in the sense of eq .  . due to the monotonicity of the exponential function , @xmath60 and @xmath61 are increasing in @xmath62 , and satisfy the causality property @xmath63 . the backlog process is accordingly described by @xmath64 since time is not affected by this transformation , the delay is given by @xmath65 to interpret these processes in the transfer domain , let @xmath66 be the instantaneous channel snr required to transmit @xmath37 in a single time slot , assuming transmission at the rate of the capacity limit . the arrival process in the snr domain can then be expressed in terms of these variables as @xmath67 here , we are treating channel quality expressed in terms of the instantaneous snr as a commodity . an arrival in a time unit represents a workload , where @xmath68 expresses the amount of resources that will be consumed by the workload . the backlog can similarly be expressed in terms of the instantaneous snr as @xmath69 with the interpretation that a node with backlog @xmath70 at time @xmath62 requires full use of the channel capacity for @xmath71 time units to clear the backlog . most importantly , the concept of the dynamic server translates to the snr domain . in a network system , the service process in the bit domain satisfies eq . if and only if the process @xmath72 in the snr domain satisfies @xmath73 we refer to a network element that satisfies eq .   for any sample path as _ dynamic snr server_. in this general setting , we not require that @xmath74 takes the form in eq . , in particular , @xmath75 does not be equal to @xmath76 . traffic aggregation in the snr domain is expressed in terms of a product . when @xmath19 flows have arrivals at a node with arrival processes denoted by @xmath77 , then the total arrival , @xmath78 , are given for any @xmath79 by @xmath80 if we let @xmath81 and @xmath82 denote the corresponding processes in the snr domain , we see that @xmath83    with the above definitions , the usual network description by a dioid algebra in the bit domain can be expressed in the snr domain by a dioid algebra on @xmath40 where the second operator is a multiplication . this enables the development of the  network calculus in sec . [ sec : calculus ] . we observe that the exponential function defines a one - to - one correspondence between arrival and departure processes in the bit and snr domains . the physical arrival , departure , service , and backlog processes can be recovered from their counterparts in the snr domain by taking a logarithm ( see fig .  [ fig : domains ] ) . this section contains our main contribution : an analytical framework for statistical end - to - end performance bounds for a network , where service is expressed in terms of fading distributions residing in the snr domain .    by an _ snr process _ we mean a bivariate process @xmath84 taking values in @xmath85 that is increasing in the second argument , with @xmath86 for all @xmath62 . the space of snr processes will be denoted by @xmath87 . for any pair of snr processes @xmath84 and @xmath88 , set @xmath89 and @xmath90 we refer to ` @xmath91 ' and ` @xmath92 ' as the the  convolution and  deconvolution operators , respectively .    the arrival , departure , and service processes constructed in the previous section are snr processes . with the  convolution , we can express the defining property of a dynamic snr server from eq . as @xmath93 for every pair of snr arrival and departure processes @xmath94 and @xmath95 . we note that , in fact , for any system description in the bit domain by the ( @xmath49 ) and the ( @xmath96 ) dioid algebras there exists a corresponding characterization in the snr domain using ( @xmath97 ) and ( @xmath98 )  dioids . the following argument confirm that the properties of a dioid are satisfied . [ lm : mx - dioid ] @xmath99 is a dioid . we show that ( @xmath97 ) satisfies the dioid axioms . for @xmath100 :    1 .   _ commutativity of @xmath101 : _ @xmath102 . associativity of @xmath101 : _ @xmath103 . idempotency of @xmath101 : _ @xmath104 . associativity of @xmath105 : _ @xmath106 . _ distributivity of @xmath105 : _ @xmath107 . existence of a null element : _ the null element is @xmath108 . since @xmath109 . _ absorption of the null element : _ _ existence of a unity element : _ the unit of multiplication is 1 , since @xmath111 . @xmath112 is a dioid . given bivariate functions @xmath113 . since the @xmath101 operation is a pointwise minimum , properties of the @xmath101 operation , that is , properties ( 1)(3 ) from the proof of lemma  [ lm : mx - dioid ] , follow from the @xmath99 dioid . for the remaining properties we have    * _ associativity of @xmath91 : _ @xmath114 * _ distributivity of @xmath91 over finite sums : _ @xmath115 * _ existence of a null element : _ the null element is @xmath116 for all values of @xmath117 and @xmath62 . * _ absorption of the null element : _ @xmath118 note that functions in @xmath87 are strictly positive by definition . * _ existence of a unity element : _ the unity element is @xmath119 , where @xmath120 this gives @xmath121      the existing network calculus in the bit domain allows for the concatenation of tandem service elements using the  convolution ( see page  1 ) . as an immediate consequence , single node performance bounds are extended to a multi - hop setting . we now show establish the corresponding result in the  network calculus . specifically , the concatenation of dynamic snr servers is again a dynamic snr server . we will prove the result for a tandem network of two nodes , as shown in fig . [ fig : concat ] . [ thm : conc ] let @xmath122 and @xmath123 be two dynamic snr servers in tandem as shown in fig . [ fig : concat ] . then , the service offered by the tandem of nodes is given by the dynamic snr server @xmath124 with @xmath125    using eq . , the departure process @xmath126 can be written as @xmath127   \end{aligned}\ ] ]    the extension to networks with more than two nodes follows by iteratively applying lemma  [ thm : conc ] . hence , the dynamic network snr server with @xmath15  nodes in tandem is given by @xmath128        performance bounds in the  network calculus are computed with the  deconvolution operator . this is analogous to role of the  deconvolution in the existing  network calculus . the bounds are expressed in the following lemma .    [ lem : obbound ] given a system with snr arrival process @xmath94 and dynamic snr server @xmath75 . * output burstiness . the snr departure process is bounded by @xmath129 . * backlog bound . the snr backlog process is bounded by @xmath130 . * delay bound . the delay process is bounded by @xmath131 . for the output bound , we fix @xmath117 and @xmath62 with @xmath132 and derive @xmath133 where we used the inequality @xmath134 in the second line .    for any fixed sample path , fix an arbitrary @xmath135 . the bound on the backlog is derived by @xmath136 where we used @xmath137 in the second step . recall that the delay is invariant under the transform of domains , that is , @xmath138 . by definition of the delay in eq . , a delay bound @xmath139 satisfies @xmath140 where we used the inequality @xmath141 in the second line .    with an algebraic description for network performance bounds in the snr domain in hand , we now turn to the problem of computing the bounds . the concise ( and familiar ) expressions from the previous section for the network service and performance bounds in the snr domain hide the difficulty of computing the expressions . in fact , all expressions of the  network calculus contain products or quotients of random variables . the _ mellin transform _ @xcite facilitates such computations , particularly when the arrival and service processes are independent . the mellin transform of a nonnegative random variable @xmath142 is defined by @xmath143\ ] ] for any complex number @xmath144 such that this expected value exists .    among its many properties , we will exploit that the mellin transform of a product of two independent random variables @xmath142 and @xmath145 equals the product of their mellin transforms , @xmath146 similarly , the mellin transform of the quotient of independent random variables is given by @xmath147    we will evaluate the mellin transform only for @xmath148 , where it is always well - defined ( though it may take the value @xmath108 ) . when @xmath149 , the mellin transform is order - preserving , i.e. , for any pair of random variables @xmath150 with @xmath151 we have @xmath152 for all @xmath144 . when @xmath153 , the order is reversed . hence bounds on the distribution of a random variable @xmath142 generally imply bounds on its mellin transform . a more subtle question is how to obtain bounds on the distribution of a random variable from its mellin transform . here , the complex inversion formula is not helpful . instead , we will use the moment bounds @xmath154 for all @xmath155 and @xmath156 . for bivariate random processes @xmath84 , we will write @xmath157 .    in our calculus , we work with the mellin transform of  convolutions and deconvolutions , which not only involves products and quotients , but also requires to compute infimums and supremums . the computation of the exact mellin transform for these operations is generally not feasible . we therefore resort to bounds , as specified in the next lemma . [ lem : mtconvdecon ] let @xmath84 and @xmath88 be two independent nonnegative bivariate random processes . for @xmath153 , the mellin transform of the  convolution @xmath158 is bounded by @xmath159 for @xmath149 , the mellin transform of the  deconvolution @xmath160 is bounded by @xmath161    note that the function @xmath162 is increasing for @xmath149 and decreasing for @xmath153 . for @xmath153 , the convolution is estimated by @xmath163\\ & \hspace{-2 cm } = e \bigl [ \ \sup_{\tau\le u \le t }   \ {   ( \x(\tau , u))^{s-1 } \cdot   ( \y ( u , t))^{s-1 } \ }   \bigr]\\   & \hspace{-2 cm } \le \sum_{u=\tau}^{t } e \big [ (   \x(\tau , u )   ) ^{s-1 }   \big ] \cdot e \big [   (     \y ( u , t )   ) ^{s-1 } \big]\,.\end{aligned}\ ] ] in the last line , we have used the non - negativity of @xmath164 and @xmath165 to replace the supremum with a sum , and their independence to evaluate the expectation of the products . .   follows by inserting the definition of the mellin transform . the deconvolution is similarly estimated for @xmath149 by @xmath166 \\ & \hspace{-2 cm } =     e \bigl [ \ \sup_{0   \le u \le \tau }   \bigl\ { \bigl ( \x ( u , t)\bigr)^{s-1 } \cdot \bigl(\y ( u,\tau ) \bigr)^{1-s } \bigr\ } \bigr]\\ & \hspace{-2 cm } \le \sum_{u=0}^{\tau }   e \bigl [ \bigl (   \x ( u , t)\bigr)^{s-1}\bigr ] \cdot e\bigl [ \bigl(\y ( u,\tau)\bigr)^{1-s}\bigr]\,,\end{aligned}\ ] ] and eq .   follows from the definition of the mellin transform .    as an application of lemmas  [ thm : conc ] and  [ lem : mtconvdecon ] , we compute a bound on the mellin transform of the service process for a cascade of fading channels . we also make the idealizing assumption that the channels are independent . [ cor : cascade ] consider a cascade of @xmath15 independent , identically distributed fading channels , where the service process for each channel is given by eq . , with i.i.d . random variables @xmath52 . let @xmath124 denote the snr service process for the entire cascade . then , the mellin transform for this process satisfies @xmath167 for all @xmath153 . we use the server concatenation formula in eq .   to represent the service of the cascade as @xmath168 . by lemma  [ lem : mtconvdecon ] , its mellin transform satisfies for @xmath153 @xmath169 where the sum runs over all sequences @xmath170 with @xmath171 and @xmath172 . the assumptions on the service processes of the individual channels imply that each product evaluates to the same function @xmath173 where @xmath1 is a random variable that has the same distribution as the @xmath52 . since the number of summands in eq . is given by @xmath174 , the claim follows . we next obtain network - level performance bounds for the bit domain . this involves a transformation from the snr domain to the bit domain via the relationship in fig . [ fig : domains ] , which provides the translation of the abstract metrics @xmath175 and @xmath176 into processes @xmath48 and @xmath177 , which , along with @xmath3 , are concrete measures for traffic burstiness , buffer requirements , and delay . [ thm : blbound ] given a system where arrivals are described by a bivariate process @xmath8 , and the available service is given by a dynamic server @xmath11 . let @xmath94 and @xmath75 be the corresponding snr processes . fix @xmath178 and define , for @xmath156 , @xmath179 then , we have the following probabilistic performance bounds . @xmath180 output burstiness : @xmath181 , where @xmath182 @xmath180 backlog : @xmath183 , where @xmath184 @xmath180 delay : @xmath185 , where @xmath186 is the smallest number satisfying @xmath187    for the bound on the distribution of the output burstiness , we start from the inequality @xmath188 . it follows from the moment bound and lemma  [ lem : mtconvdecon ] that , for any choice of @xmath189 and all @xmath156 @xmath190 to obtain the claim , we set the right hand side equal to @xmath191 , solve for @xmath192 , and optimize over the value of @xmath156 to obtain @xmath193 . the proof of the backlog bound proceeds in the same way , starting from the inequality @xmath194 . the delay bound is slightly more subtle . fix @xmath195 . using lemma [ lem : obbound ] and the moment bound with @xmath196 , we obtain that @xmath197 for every @xmath156 . by lemma  [ lem : mtconvdecon ] , the mellin transform @xmath198 satisfies a bound that agrees with the function @xmath199 , except that the upper limit in the summation that defines @xmath199 would have to be replaced by @xmath200 . to obtain a sharper estimate , we use instead eq .   from the proof of lemma [ lem : obbound ] . the resulting bound is that @xmath201 satisfies @xmath202 where we have used that the supremum in the definition of @xmath203 extends only up to @xmath204 , and then repeated the proof eq .  . the claim follows by optimizing over @xmath144 . corresponding bounds as in theorem  [ thm : blbound ] can be obtained using the algebra and the network calculus with moment - generating functions @xcite . the significance of theorem  [ thm : blbound ] is that it permits the application of the the network calculus , where traffic is characterized in the bit domain , and service is naturally expressed in the snr domain . this will become evident in the next section , where the theorem  gives us concise bounds for delays and backlog of multi - hop networks with rayleigh fading channels . we now apply the techniques developed in the two previous sections to a network of rayleigh channels . consider the dynamic snr server description for a rayleigh fading channel , as constructed in sec . [ sec : model].b . we use eq .  , with the function @xmath6 given by @xmath205 where @xmath206 is the average snr of the channel and @xmath207 is the fading gain . for rayleigh fading , @xmath207 is a rayleigh random variable with probability density @xmath208 . in a physical system , @xmath209 , where @xmath210 and @xmath211 are the received signal power and the ( additive white gaussian ) noise power at the receiver , respectively . then , @xmath212 is exponentially distributed , and the mellin transform of @xmath6 is given by @xmath213 where @xmath214 is the upper incomplete gamma function . using the assumption that the @xmath52 are independent and identically distributed , we obtain for the mellin transform of the dynamic server @xmath215    for the arrival process , we use a characterization due to chang @xcite , where the moment - generating function of the cumulative arrival process in the bit domain is bounded by @xmath216 \le   \rho(s)\cdot(t-\tau ) + \sigma(s)\end{aligned}\ ] ] for some @xmath156 . in general , @xmath217 and @xmath218 are nonnegative increasing functions of @xmath144 that may become infinite when @xmath144 is large . this traffic class , referred to as  bounded arrivals , is broad enough to include markov - modulated arrival processes . the corresponding class of snr arrival processes is defined by the condition that @xmath219 for some @xmath149 . we consider the transmission of  bounded arrivals on a rayleigh fading channel . to obtain single - hop performance bounds , we apply theorem  [ thm : blbound ] with the expressions for the mellin transforms for the snr service and arrival processes from eqs .   and . for the function @xmath220 from the statement of the theorem , we compute @xmath221_+}^ { \infty }    \big (   \underbrace{e^{s   \rho(s ) }    e^{1/\bar\gamma }   \bar{\gamma}^{-s } \gamma(1-s,\bar\gamma^{-1})}_{\deq v(s ) } \big ) ^ { u } \,,\end{aligned}\ ] ] where @xmath222_+$ ] is the maximum of @xmath223 and @xmath224 . the sum converges when @xmath225 , which can be interpreted as a stability condition . inserting the result into theorem  [ thm : blbound ] , we obtain for the output burstiness the probabilistic bound @xmath226 the backlog bound is obtained by setting @xmath227 , @xmath228 the delay bound is the smallest number @xmath186 such that @xmath229    we also derive end - to - end bounds for a cascade of @xmath15 rayleigh channels with  bounded arrivals , using the same parameters as before . let @xmath124 be the service process for the entire cascade . by corollary  [ cor : cascade ] , its mellin transform satisfies for @xmath230 @xmath231 we will use again theorem  [ thm : blbound ] . for @xmath230 , we compute @xmath232 where @xmath233 is as defined above , and where we applied the combinatorial identity @xmath234 for any @xmath235 and for @xmath236 . inserting @xmath237 into theorem  [ thm : blbound ] gives for the end - to - end output bound , denoted by @xmath238 , the value @xmath239 note that for @xmath240 , this bound agrees with the previous bound for a single node . in the same way , we derive the probabilistic backlog bound @xmath241 for the delay bound , we estimate for @xmath242 @xmath243 here , the first term in the minimum is obtained by extending the summation down to @xmath244 , and the second term results from the inequality @xmath245 for @xmath246 . in both cases , the resulting sum can be evaluated with eq .  . the delay bound @xmath186 is determined according to theorem  [ thm : blbound ] by setting the right hand side of eq . equal to @xmath191 , solving for @xmath139 , and minimizing over @xmath144 . because of the complexity of the bound in eq . , the last two steps can only be performed numerically . it is apparent that the complexity of computing end - to - end bounds is no different than bounds for a single channel . more importantly , observe that the end - to - end bounds scale linearly in the number of nodes @xmath15 . relaxing the strong independence assumptions on the channel properties would give different scaling properties .    ) vs. average channel snr ( @xmath247 ) for multihop rayleigh fading channels with @xmath248 ,  bounded traffic with @xmath249  kb and @xmath250  kbps , and @xmath251  khz.,width=480 ]    ) vs. arrival rate ( @xmath217 ) for multi - hop rayleigh fading channels with @xmath248 ,  bounded traffic with @xmath249  kb , @xmath252  db , and @xmath251  khz.,width=480 ]    ) vs. arrival rate ( @xmath217 ) for multi - hop rayleigh fading channels with @xmath248 ,  bounded traffic with @xmath249  kb , @xmath252  db , and @xmath251  khz.,width=480 ]      we next present numerical results , where we assume rayleigh channels with a transmission bandwidth of @xmath253  khz . for traffic , we use  bounded arrivals with default values @xmath254  kb and @xmath255  kbps , i.e. , the bounds on rate and bursts are deterministic . hence , the source of randomness in the examples results only from randomness of the channels . we use a violation probability of @xmath256 .    in fig . [ fig : backlog - vs - snr - multihop ] we show the end - to - end backlog for cascades of @xmath15  rayleigh channels , as a function of the average snr of each channel . even though the backlog bounds increase only linearly in the number of nodes , the per - node requirements  at least for the last node of the cascade  must satisfy the end - to - end bounds , since it can not be assumed that backlog is equally distributed across the nodes . when the snr of the nodes is sufficiently high , the backlog remains low even for a large number of hops . we observe that the channel becomes saturated for @xmath257  db . when the number of channels is small , the backlog increases sharply in the vicinity of @xmath257 , but remains low everywhere else .    in figs . [ fig : backlog - vs - rho - multihop ] and  [ fig : backlog - vs - rho2-multihop ] we present , for values of the average snr of @xmath258  db ( fig . [ fig : backlog - vs - rho - multihop ] ) and @xmath259  db ( fig . [ fig : backlog - vs - rho2-multihop ] ) , how the end - to - end backlog increases as a function of the transmission rate , for different network sizes . we observe that the maximum achievable rate that does not result in a ` blow - up ' of the backlog decreases as the number of nodes is increased . we observe that the blow - up occurs earlier when the channel has a lower snr .    ) vs. average channel snr ( @xmath247 ) for multi - hop rayleigh fading channels for @xmath260 and @xmath261 with buffer size @xmath262  kb ,  bounded traffic with @xmath249  kb and @xmath263 or @xmath264 kbps , and @xmath251  khz.,width=480 ]    ) vs. average channel snr ( @xmath247 ) for multi - hop rayleigh fading channels for @xmath260 and @xmath261 with buffer size @xmath262  kb ,  bounded traffic with @xmath249  kb and @xmath263 or @xmath264 kbps , and @xmath251  khz.,width=480 ]    ) vs. end - to - end delay bound for multi - hop rayleigh fading channels for @xmath260 , @xmath265  db ,  bounded traffic with @xmath249  kb and @xmath250 kbps and @xmath251  khz.,width=480 ]    ) vs. end - to - end delay bound for multi - hop rayleigh fading channels for @xmath260 , @xmath265  db ,  bounded traffic with @xmath249  kb and @xmath250 kbps and @xmath251  khz.,width=480 ]    ) vs. average channel snr ( @xmath247 ) for multi - hop rayleigh fading channels for @xmath266 , @xmath267  ms ,  bounded traffic with @xmath268  kb and @xmath263  kbps and @xmath251  khz.,width=480 ]    suppose that buffer sizes are set to satisfy the end - to - end backlog . for a fixed buffer size @xmath269 , we can then use the probability @xmath270 as an estimate of the probability of dropped traffic , and refer to it as the _ loss probability_. in figs . [ fig : loss - vs - snr - multihop ] and  [ fig : loss - vs - snr2-multihop ] , we depict the loss probability as a function of the average channel snr , for values of @xmath271  kb ( fig . [ fig : loss - vs - snr - multihop ] ) and @xmath272  kb ( fig . [ fig : loss - vs - snr2-multihop ] ) , for traffic with a rate of @xmath273  and @xmath264  kbps , and for @xmath260 , and  @xmath261 nodes . the figure shows the minimum snr needed to support a given loss probability is very sensitive to the number of network nodes . in figs . [ fig : viol - vs - delay - multihop ] and  [ fig : viol - vs - delay - multihop2 ] we show the violation probability for given end - to - end delay bounds , where we compare the delays at a single node ( @xmath240 ) with a multi - hop network ( @xmath266 ) for different average channel snr values , where we use eqs .   and  . the traffic parameters are @xmath249  kb for the burst , and @xmath263  kbps ( in fig . [ fig : viol - vs - delay - multihop ] ) @xmath250  kbps ( in fig . [ fig : viol - vs - delay - multihop2 ] ) . the figures illustrates that at sufficiently high snr values , low delays are achieved even when traffic traverses 10  links . when the snr is decreased , we can observe how the delay performance deteriorates in the multi - hop scenario . note that the curves for the violation probability are essentially piecewise linear with two segments . this is caused by the different exponential decay rates , which follows from eq .  . depending on where in the equation the minimum occurs , we obtain a faster or slower decay . next we investigate how the delay violation probability for given end - to - end delay bounds is impacted by the snr of the channel . we consider a network with @xmath266  nodes , with the same parameters as before . the traffic parameters are @xmath249  kb and @xmath263  kbps . we consider end - to - end delay bounds of @xmath267  ms . [ fig : viol - vs - snr - multihop ] presents the results . an interesting observation is that the snr required to support a given violation probability for a short delay bound @xmath274  ms , the snr requirement increase fast for low violation probabilities . delay bounds @xmath275  ms and higher can be supported with low violation probabilities even when the average snr is small . the graphs presented here can be used in the planning of a multi - hop wireless network where predefined qos bounds are desired for a given transmission rate . since the average snr depends largely on the path loss , which , in turn , is a function of the transmission radius , the graphs could help with determining the maximum distance between nodes to support a desired transmission rate and qos . consider a scenario in fig .  [ fig : leftover - snr ] where a through flow arriving to a fading channel shares the available bandwidth with other flows . we will refer to the traffic from these other flows as _ cross traffic_. we use @xmath276 and @xmath277 to denote the snr arrival processes of the through flow and the cross traffic , respectively , and let @xmath278 and @xmath279 denote the corresponding departure processes . in the snr domain , cross traffic can be viewed as reducing the channel capacity of the through flow by generating interference .    the following lemma states that , in the snr domain , the service available to a through flow that experiences cross traffic at a channel can be expressed by a dynamic snr server .    [ lem : leftover ] consider a channel with a through flow and cross traffic as shown in fig . [ fig : leftover - snr ] . assume that the channel provides a dynamic snr server to the aggregate of through flow and cross traffic , with service process @xmath75 , i.e , @xmath280 then @xmath281 is a dynamic snr server satisfying for all @xmath135 that @xmath282 if , moreover , the service to the cross flows satisfies the upper bound @xmath283 , then @xmath284    for any a sample path , and any @xmath195 , we have @xmath285 let @xmath286 be the point where the infimum is assumed . dividing by @xmath287 , we obtain @xmath288 where we used that @xmath289 by causality . this the first claim in eq .  . to see the second claim , assume that @xmath283 . then @xmath290 , and therefore @xmath291 combining this with the first part of the proof , we obtain @xmath292 proving the second claim . note that @xmath293 need not be monotone in @xmath62 , i.e. , it may not lie in @xmath87 . we next estimate the service process available to the through flow across a cascade of channels . assuming that @xmath294 and @xmath74 are independent , we obtain for the mellin transform of the service process at a single node @xmath295 these service descriptions may be convolved , using lemmas  [ thm : conc ] and  [ lem : mtconvdecon ] to obtain a bound for the mellin transform of the service provided by a cascade of fading channels to a flow that experiences cross traffic at each node . [ cor : cascade - leftover ] consider a network as in fig . [ fig : tandem ] , where a through flow experiences cross traffic at each fading channel . let the snr service process at each channel be given by eq .  . assume that the arrival process of the cross traffic satisfies eq .   with parameters @xmath296 . assume that arrivals from through flow and cross traffic , as well as the service processes at each channel are independent . then , the service provided by the network to the through flow satisfies @xmath297 , where the mellin transform of @xmath298 satisfies , for @xmath153 , @xmath299    for a single node ( @xmath240 ) , we estimate for @xmath230 and @xmath153 @xmath300 where we have used that @xmath301 for @xmath153 . by lemma  [ thm : conc ] , the service of the network is given by the  convolution , @xmath302 . we use lemma  [ lem : mtconvdecon ] to bound its mellin transform by @xmath303 where the sum runs over all sequences @xmath170 with @xmath171 and @xmath172 . collecting terms , we compute , as in the proof of corollary  [ cor : cascade ] , that each product evaluates to the same term , @xmath304 since the number of summands is @xmath174 , this proves the claim . we now consider rayleigh fading channels , and assume that both through flow and cross traffic are  bounded traffic with parameters @xmath305 and @xmath306 for the through flow , and @xmath307 and @xmath308 for the cross traffic . we next give end - to - end performance bounds , using theorem  [ thm : blbound ] . for the function @xmath309 , we compute as in sec . [ sec : rayleigh].b for @xmath310 , @xmath311 where @xmath312 this computation is valid under the stability condition that @xmath313 . thus , we obtain that the output burstiness at the network egress gives @xmath314    \bigg \ } \ , . \nonumber    \end{aligned}\ ] ] for the end - to - end backlog of the through flow , we obtain @xmath315     \bigg \ } \ , . \nonumber    \end{aligned}\ ] ] similarly , for the delay bound , we estimate for @xmath242 @xmath316    ) vs. arrival rate ( @xmath306 ) for multi - hop rayleigh fading channels with @xmath248 ,  bounded traffic with @xmath317  kb ,  bounded cross traffic with @xmath318  kb and @xmath319 , @xmath252  db , and @xmath251  khz.,width=480 ]    ) vs. arrival rate ( @xmath306 ) for multi - hop rayleigh fading channels with @xmath248 ,  bounded traffic with @xmath317  kb ,  bounded cross traffic with @xmath318  kb and @xmath319 , @xmath252  db , and @xmath251  khz.,width=480 ]    ) vs. arrival rate ( @xmath306 ) for multi - hop rayleigh fading channels with @xmath248 ,  bounded traffic with @xmath317  kb ,  bounded cross traffic with @xmath318  kb and @xmath319 , @xmath252  db , and @xmath251  khz.,width=480 ]      we now present numerical examples for rayleigh fading channels . we use the same traffic and channel model as in the numerical examples in sec . [ sec : numerical ] . the cross traffic is  bounded traffic with parameters @xmath307 and @xmath308 . as before , we assume that cross traffic is deterministically bounded by using fixed values for @xmath307 and @xmath308 . consequently , there is no statistical multiplexing between through and cross traffic . throughout , the violation probability is set to @xmath248 .    in fig . [ fig : backlogvssnr_wct ] we show the end - to - end backlog bound @xmath320 as a function of the average channel snr @xmath206 . the through flow has fixed parameters @xmath317  kb and @xmath321  kbps , and the cross traffic has parameters @xmath318  kb and @xmath319 . we consider networks with @xmath266 and @xmath322 nodes . the graph illustrates how the bursts of the cross traffic contribute to the backlog bound . there is an additive component for each traversed channel , which explains the difference for the backlog with @xmath323 and @xmath324  channels . the minimum required snr needed for stability appears less sensitive to the number of channels .    in figs . [ fig : backlogvsrho_wct1 ] and  [ fig : backlogvsrho_wct2 ] we again evaluate the end - to - end backlog bound @xmath325 . here , we keep the channel snr values constant at @xmath258 ( fig . [ fig : backlogvsrho_wct1 ] ) and @xmath259 ( fig . [ fig : backlogvsrho_wct2 ] ) . we vary the rate of the through flow @xmath306 on the x axis , and plot graphs for different values of @xmath319 . we again consider @xmath266 and @xmath322 . the outcomes are as seen in the previous figure . we observe the effects of the additive component contributed by each node . also , we observe , at least for lower transmission rates , that the stability of end - to - end backlogs is not sensitive to the number of traversed channels . we have developed a novel network calculus that can incorporate fading channel distributions , without the need for secondary models , such as fsmc . we use the calculus to compute statistical bounds on delay and backlog of multi - hop wireless networks with fading channels . we took a fresh point of view , where the descriptions of the arrivals and the fading channels reside in different domains , referred to as bit domain and snr domain . we found that by mapping arrival processes to the snr domain , an end - to - end analysis with fading channels becomes tractable . an important discovery was that arrivals and service in the snr domain obey the laws of a  dioid algebra . the analytical framework developed in this paper appears suitable to study a broad class of fading channels and their impact on the network - layer performance in wireless networks . even though we computed numerical examples for very simplified networks , in particular , we made strong independence assumptions for the fading channels , our  network calculus is applicable to networks where these assumptions are relaxed . generalizing our framework and obtaining a more profound understanding of the dioid algebra and computational methods in the snr domain is the subject of future research . f.  ciucu , p.  hui , and o.  hohlfeld . non - asymptotic throughput and delay distributions in multi - hop wireless networks . in _ proc . 48th annual allerton conf . on communication , control , and computing _ , pages 662669 , sep . 2010 . y.  jiang and p.  j. emstad . analysis of stochastic service guarantees in communication networks : a server model . in _ proc . iwqos , springer lecture notes in computer science 3552 _ , pages 233245 , june 2005 . f.  kelly . notes on effective bandwidths . in _ stochastic networks : theory and applications . ( editors : f.p . kelly , s. zachary and i.b . ziedins ) royal statistical society lecture notes series , 4 _ , pages 141168 . oxford university press , 1996 .
Abstract:

Answer: a fundamental problem for the delay and backlog analysis across multi - hop paths in wireless networks is how to account for the random properties of the wireless channel . since the usual statistical models for radio signals in a propagation environment do not lend themselves easily to a description of the available service rate , the performance analysis of wireless networks has resorted to higher - layer abstractions , e.g. , using markov chain models . in this work , we propose a network calculus that can incorporate common statistical models of fading channels and obtain statistical bounds on delay and backlog across multiple nodes . we conduct the analysis in a transfer domain , which we refer to as the _ snr domain _ , where the service process at a link is characterized by the instantaneous signal - to - noise ratio at the receiver . we discover that , in the transfer domain , the network model is governed by a dioid algebra , which we refer to as  algebra . using this algebra we derive the desired delay and backlog bounds . an application of the analysis is demonstrated for a simple multi - hop network with rayleigh fading channels .