Service Outage Based Power and Rate Allocation

Size: px
Start display at page:

Download "Service Outage Based Power and Rate Allocation"

Transcription

1 1 Service Outage Based Power and Rate Allocation Jiangong Luo, Lang Lin, Roy Yates, and Predrag Spasojević Abstract Tis paper combines te concepts of ergodic capacity and capacity versus outage for fading cannels, and explores variable rate transmissions under a service outage constraint in a block flat fading cannel model. A service outage occurs wen te transmission rate is below a given basic rate. We solve te problem of maximizing te expected rate subject to te average power constraint and te service outage probability constraint. Wen te problem is feasible, te optimum power policy is sown to be a combination of water filling and cannel inversion allocation, were te outage occurs at a set of cannel states below a certain tresold. Te service outage approac resolves te conflicting objectives of ig average rate and low outage probability. Index Terms Adaptive transmission, block-fading cannel, ergodic capacity, outage capacity, service outage I. INTRODUCTION Wireless communication cannels vary wit time due to mobility of users and canges in te environment. For a time varying cannel, dynamic allocation of resources suc as power, rate, and bandwidt can yield improved performance over fixed allocation strategies. Indeed, adaptive tecniques are employed in EDGE [1], GPRS [2], and HDR [3], and are proposed as standards for next generation cellular systems. Te system performance criterion is usually application specific; terefore, different classes of applications will benefit from specific adaptive transmission scemes. In order to differentiate real-time service from non real-time service, tree capacity measures ave been defined in te literature: ergodic capacity [4], delay limited capacity [5], and capacity versus outage probability [9 11]. A compreensive survey of tese concepts can be found in [6]. Te ergodic capacity [4] was developed for non real-time data services. It determines te maximum acievable rate averaged over all fading states. Te corresponding optimum power allocation is te well known water filling allocation [7, 8]. Te ergodic capacity may not be relevant for real-time applications in a slow fading environment, were substantial delay can occur wen averaging over all fading states. Delay limited capacity [5] and te capacity versus outage probability [9 11] were developed for constant-rate real-time applications. Te delay limited capacity specifies te igest acievable rate wit a decoding delay independent of fading correlation structures [5]. Te outage capacity in te capacity versus outage probability problem determines te ɛ-acievable rate [12] of te M-block fading cannel. Te corresponding optimum power allocation was derived in [1] for M parallel flat fading blocks (frequency diversity or space diversity), and in [11] for M consecutive flat fading blocks (time diversity). Te zero-outage capacity in [1, 11] is also referred to as te delay-limited capacity. Toug te outage capacity studies te capacity for constantrate real-time applications, te constant-rate assumption may not be essential for many real-time applications. For applications wit simultaneous voice and data transmissions, for example, as soon as a basic rate r o for te voice service as been guaranteed, any excess rate can be used to transmit data in a best effort fasion. For some video or audio applications, te source rate can be adapted according to te fading cannel conditions to provide multiple quality of service levels. Typically tese applications require a nonzero basic service rate r o to acieve a minimum acceptable service quality. Motivated by tese variable-rate real-time applications, we study variable rate transmission scemes subject to a basic rate requirement in a slow fading environment. By allowing variable rate transmissions, te variation of te fading cannel can be exploited to acieve an average rate iger tan te outage capacity. By imposing a basic rate requirement, te system can be guaranteed to operate properly. Since infinite average power is needed to acieve any nonzero rate at all times in a Rayleig fading cannel, we impose a probabilistic basic service rate requirement, tat we call a service outage constraint. Te service is said to be in an outage wen te information rate is smaller tan te basic service rate r o. Service quality is acceptable as long as te probability of te service outage is less tan ɛ, a parameter indicating te outage tolerance of te application. Unlike te information outage in te capacity versus outage probability problem [1, 11], te bits transmitted during te service outage may still be valuable in tat tey will be transmitted reliably and will contribute to te average rate. For variable-rate systems, te expected rate determines ow muc rate can be transmitted on te average and is a meaningful measure of system performance. Terefore, in tis paper, te allocation problem is to find te power and rate allocation tat maximize te expected rate subject to te service outage constraint and te average power constraint. Under te assumptions of a block flat fading AWGN cannel model and perfect cannel state information at te transmitter, we verify tat te outage sould occur at bad cannel states below a certain tresold. Te resulting optimum power allocation is sown to be a combination of cannel inversion and water filling wen te allocation problem is feasible. Te service outage approac strikes good balance between te average rate and te outage probability. Tis approac as been generalized to te case of code words spanning multiple blocks in [13]. Altoug a continuous fading distribution is assumed in tis paper, te results can be extended into te case of discrete fading distributions by employing a probabilistic power allocation, as examplified by te policies in [1]. Altoug our problem as been motivated by real time applications, it also caracterizes coverage versus capacity tradeoffs. In particular, mobility in cellular systems results in can-

2 2 nel variations due to canges in distance attenuation. An important objective of a cellular system is to provide a basic service rate over as muc of te service area as possible. In tis case, te service outage constraint caracterizes te spatial coverage requirement of te system. Te objective is ten to maximize te average rate over all geograpic locations subject to meeting te service outage constraint. Te remainder of tis paper is organized as follows. In Section II, te system model and te optimization problem are presented. In Section III, te optimum allocation policy is derived. In Section IV, a supporting teorem for te optimum allocation policy is proved. Furter discussion of te optimum solution is presented in Section V. II. THE ALLOCATION PROBLEM In tis work, we employ te block flat fading Gaussian cannel (BF-AWGN) model [9]. In te BF-AWGN cannel, a block of N symbols experiences te same cannel state, wic is constant over te wole block, but may vary from block to block. Witin eac block we ave te time-invariant Gaussian cannel y = x + n. Here x is te cannel input, y is te cannel output, n is wite Gaussian noise wit variance σ 2, and is te cannel state. Let f() denote te probability density function of te cannel state and F () denote te corresponding cumulative distribution function (CDF). Here, we consider te case were F () is a continuous function and te power allocation is a deterministic function of cannel state. We make te following assumptions: Te cannel state information is known perfectly at bot transmitter and receiver. One codeword spans one fading block and te block size N is large. Te fading process is ergodic over te time scale of te application. As pointed out in [1], it makes sense to study te BF-AWGN cannel as N, since for typical practical systems N is fairly large and outage is te dominant error event wen using an actual code. Let p() denote te transmission power at cannel state. Ten te maximum acievable rate at eac block is te capacity of Gaussian cannel wit received power p(), and is denoted as R[p()], were R[P ] = 1 2 log 2 (1 + Pσ ) 2. (1) Terefore, te resource allocation problem requires finding only te optimum power allocation p (). Under te assumption of te ergodicity of fading process, te time average rate of te system can be caracterized by te expected rate. Terefore, given te average power P av, te basic service rate r o, and te allowable service outage probability ɛ, we wis to maximize te expected code rate, as follows: R = max p() E {R[p()]} (2) subject to: E {p()} P av (2a) p() Pr{R[p()] < r o } ɛ. (2b) (2c) In te absence of te service outage constraint (2c), R would be te ergodic capacity for te fading cannel, and te well known water filling allocation [7, 8] would be te corresponding optimum power assignment. III. OPTIMUM POWER AND RATE ALLOCATION In tis section, we derive an optimum power allocation p () for problem (2). Te difficulty in deriving p () is primarily due to te probabilistic nature of te constraint (2c). Here, we sow ow an optimum power allocation can be derived based on a problem analogous to (2) wit a deterministic constraint on te assigned rate. Given a basic service rate r o and a power policy p(), te service set is defined as H s (p()) = { R[p()] r o }, and te outage set is H o (p()) = { R[p()] < r o }. Our approac will be to sow tat tere is an optimal solution to problem (2) wit a particular form of a service set. Prior to sowing tis, we solve te following subproblem in wic it is required tat te service set contains an arbitrary set H a. R (H a ) = max p() E {R[p()]} (3) subject to: E {p()} P av (3a) p() (3b) R[p()] r o H a. (3c) Let p (, H a ) denote an optimum solution to problem (3). Terefore, p (, H a ) acieves te igest average rate among all te scemes wose service set contains H a. Problem (3) does not necessarily ave a solution for a given (P av, r o, H a ). Constraint (3c) implies tat a feasible allocation p() must satisfy p() σ2 (2 2ro 1) H a. (4) Tis implies tat te minimum average power needed to meet te constraint (3c) for a given (r o, H a ) is σ 2 (2 2ro 1) P min (r o, H a ) = f() d. (5) H a Consequently, problem (3) as a solution only if P av P min (r o, H a ). Wen P av = P min (r o, H a ) te corresponding power allocation is te on-off cannel inversion policy σ 2 (2 2ro 1) p (, H a ) = H a, (6) oterwise. Wen P av > P min (r o, H a ) te corresponding power allocation is given by te following teorem. We use te notation (x) + = max(x, ). Teorem 1 Wen P av > P min (r o, H a ) te optimum solution for problem (3) is: σ 2 (2 2ro 1) H a { } 2 2ro, p (, H a ) = ( 1 σ oterwise. ) (7)

3 3 were is te solution of E {p (, H a )} = P av. Teorem 1 follows from standard variational arguments; te proof appears in te Appendix. Note tat wen P av = P min (r o, H a ) te resulting power allocation (6) can be viewed as a limiting case of expression (7) as. Te power allocation p (, H a ) is a combination of cannel inversion and water filling allocations. To obtain a ig average rate, we would like to allocate power in te form of te water filling allocation, wile to meet te service constraint (3c), we must allocate power no less tan te cannel inversion allocation witin te set H a. Te solution p (, H a ) balances tese two factors. To caracterize te solution to te optimization problem (2), we define ɛ as te solution to F ( ɛ ) = ɛ. Te tresold ɛ partitions te cannels into a set H ɛ = { ɛ } of good cannels and te complementary set H ɛ = { < ɛ } of bad cannels. In te following, we sow tat te solution of (3) wit H a = H ɛ, specifically p () = p (, H ɛ ), is an optimum solution of problem (2). In order to prove tis, we define te partial ordering and sow a number of preliminary results. Definition 1 H 1 H 2 if 1 < 2 for all 1 H 1 and 2 H 2. Teorem 2 Problem (2) as an optimal solution p () wit te outage set H o (p ()) and te service set H s (p ()) satisfying H o (p ()) H s (p ()). Teorem 2 sows tat tere exists an optimum power allocation suc tat te outage occurs wen te cannel state is worse tan a particular tresold. Proof of Teorem 2 involves a somewat complicated two-step construction and is deferred to Section IV. Using Teorem 2 and te fact tat Pr{H s (p ())} 1 ɛ by constraint (2c), it is easy to sow te following corollary. Corollary 1 Problem (2) as an optimum solution p () suc tat H ɛ H s (p ()). Now we can prove p () = p (, H ɛ ) by sowing tat R = R (H ɛ ). Wit H a = H ɛ in te outage constraint (3c), te service set of p (, H ɛ ) must contain H ɛ. Tus p (, H ɛ ) satisfies te outage constraint (2c) and is a feasible power allocation sceme for problem (2), implying R (H ɛ ) R. On te oter and, Corollary 1 implies tat problem (2) as an optimal solution p () acieving an average rate of R tat satisfies constraint (3c) wit H a = H ɛ. Tus, p () is a feasible power allocation sceme for problem (3) and R R (H ɛ ). Consequently, R = R (H ɛ ). In conclusion, an optimum solution is p () = p (, H ɛ ) and te following conclusions apply to problem (2). Problem (2) is feasible if only if (P av, r o, ɛ) satisfies P av P min (r o, H ɛ ) = σ 2 (2 2ro 1) ɛ Wen P av = P min (r o, H ɛ ) we ave σ 2 (2 2ro 1) p () = ɛ, < ɛ. f() d, (8) (9) Wen P av > P min (r o, H ɛ ), we can apply Teorem 1 wit H a = H ɛ yielding an optimum solution to problem (2) of te form p () = σ 2 (2 2ro 1) ɛ min{ ɛ, 22ro }, ( 1 σ oterwise, ) (1) were is te solution of E {p ()} = P av. As P av approaces P min (r o, H ɛ ), and te power allocation (1) will reduce to te on-off cannel inversion allocation (9). IV. OPTIMUM SERVICE SETS In tis section, we will prove Teorem 2, wic implies tat we can find an optimal solution wose service set H s (p ()) includes te good cannel states H ɛ. Our approac will be to sow tat given an arbitrary feasible power allocation sceme ˆp(), we can always construct a better sceme p () wic satisfies H o (p ()) H s (p ()). Tis implies tat tere is an optimum power allocation sceme p () wit H o (p ()) H s (p ()). Let Ĥ s denote te service set and ˆR te average rate for te policy ˆp(). Feasibility of ˆp() implies E {ˆp()} P av and Pr{Ĥs} 1 ɛ. We use a two-step construction. First, we construct p () from ˆp() by setting H a = Ĥs in problem (3), yielding te solution p () = p (, Ĥs) σ 2 (2 2ro 1) = Ĥs { } <, ( 22ro 1 σ oterwise. ) (11) were is te solution of E {p (, Ĥs)} = P av. Here in te case of P av = P min (r o, Ĥs), p () can be expressed by (11) as. Clearly, p () is feasible and acieves a iger average rate tan ˆp(). Second, we construct p () by decomposing p () into a water filling component and a residual power component. Given, we define te following functions over te wole cannel state space: ( 1 p wf () = σ2 1 +, (12) ) ( σ p 2 + (2 2ro 1) res() = p wf()). (13) Te function p wf () is a water filling allocation over te wole cannel space. Te function p res() is te nonnegative difference of cannel inversion and water filling allocations. From (11), we observe tat water filling alone meets te service condition R[p ()] r o over te set of cannel states H wf = { 2 2ro}. (14)

4 4 ^ R(p()) r p () p () service set (a) residual power residual power (b) (c) 2r 2 2r 2 water filling power water filling power Fig. 1. (a) Rate allocation R[ˆp()] for policy ˆp(), (b) te improved policy p () given by (16) wit water filling p wf () and residual power 1 ( H inv )p res (), (c) te new policy p () given by (17) wit water filling p wf () and residual power 1 ( H inv )p res (). In particular, p res() = for H wf wile residual power p res () > is needed to meet te service condition over te cannel inversion set H inv = Ĥs\H wf. (15) Tus, in terms of te indicator function 1 (x) suc tat 1 (x) = 1 wen x is true, and oterwise, p () can be rewritten in te form p () = p wf() + 1 ( H inv)p res(). (16) Here, we call 1 ( H inv )p res() te residual power allocation for p (). As sown in Figure 1, p () can be viewed as a twolayer allocation: te first layer is te water filling allocation over te wole cannel space and te second layer is te residual power allocation over H inv. Based on p (), we construct p () by preserving te first layer water filling allocation and redistributing te residual power. Intuitively, te best allocation sceme for te residual power is to allocate it to te good cannel states. Since p res() is strictly positive witin < 22ro, we will allocate te residual power over te interval [ b, 2 2ro ] were b is cosen to consume te total residual power. As sown in Figure 1, we ave were p () = p wf () + 1 ( H inv )p res (). (17) H inv = { b < 2 2ro}, (18) and b is te solution to p res ()f() d = H inv H inv p res ()f() d. (19) Note tat (16), (17), and (19) imply tat p () as te same total power as p (). Let R and R denote te average rates for p () and p (), respectively. Te following lemma gives us te properties of p (). Lemma 1 Te power sceme p () as te following properties: (a) E {p ()} = E {p ()} = P av (b) H o (p ()) H s (p ()) (c) R R (d) Pr{H s (p ())} Pr{H s (p ())}. Proof of Lemma 1 is given in Appendix. Hence, we summarize te proof: 1) Start wit arbitrary ˆp() wit average rate ˆR and service set Ĥs. 2) Set H a = Ĥs and solve (3) yielding p () wit rate R ˆR and service set H s (p ()) containing Ĥs. 3) Decompose p () into water filling p wf () and residual power components 1 ( H inv )p res (). 4) Fix te water filling component p wf () and reallocate te residual power to generate p (). Te power policy p () satisfies Pr{H s (p ())} Pr{H s (p ())} and R R. Hence, p () is a better power allocation sceme tan p () for problem (2). We can conclude tat from any feasible ˆp() we can obtain a better power allocation p () in wic H o (p ()) H s (p ()) olds. Tis implies tat problem (2) as an optimum solution p () satisfying H o (p ()) H s (p ()). V. PROPERTIES OF THE OPTIMUM POLICY In Section III, we derived te optimum allocation sceme for problem (2). In tis section, we will discuss tis optimum solution, and sow ow problem (2) in tis paper is related to te capacity versus outage probability problem. Te optimum power allocation sceme (1) includes a combination of cannel inversion and water filling. For a given probability distribution f(), te optimum solution belongs to one of te following possible types depending on te value of (P av, r o, ɛ): 1 I Wen P av = P min (r o, H ɛ ), p () includes no transmission for < ɛ and cannel inversion for ɛ. II Wen P av > P min (r o, H ɛ ) but ɛ, p () includes no transmission for < ɛ, cannel inversion for ɛ <, and water filling for 22ro 22ro. III Wen P av is sufficiently ig suc tat ɛ 2 2ro < < ɛ, p () includes no transmission for <, water filling for < ɛ, cannel inversion for ɛ < 2 2ro, and water filling for 2 2ro. 1 In te case of r o = or ɛ = 1, te solution types II and III will degenerate into solution type IV, wic is te pure water filling allocation.

5 5 p () I R[p ()] I ergodic capacity service outage approac wit r =.5 bits/symbol service outage approac wit r =.7 bits/symbol outage capacity r o 1.6 ε ε 1.4 II r o II R av (bits/symbol) ε 2 2r o ε 2 2r o.8 III III.6 r o P av (db) IV ε 2 2r o IV ε 2 2r o Fig. 3. Comparison of service outage approac wit oter capacity notions in te Rayleig fading cannel, for a fixed ɛ =.1.. ε r o ε.35.3 r =.5 bits/symbol r =.7 bits/symbol ε=.1 Fig. 2. For optimum solution types I-IV, power policies are given on te left and corresponding rate allocation are on te rigt. IV Wen P av is ig enoug for ɛ2 2ro, p () is just te water filling allocation. Tese four types of power allocation scemes as well as te corresponding rate allocations are depicted in Figure 2. For solution types I, II, and III, te optimum service set is H s (p ()) = H ɛ and te resulting outage probability is ɛ, wile for type IV solution H ɛ H s (p ()) and te resulting outage probability is less tan ɛ. Type I solution is te on-off cannel inversion allocation. In tis case, we ave just enoug average power to satisfy te service outage constraint. Wen we ave extra power beyond P min (r o, H ɛ ), we can allocate te power in a more efficient way to obtain a iger average rate and, at te same time, to meet te service outage constraint. Wen P av is sufficiently ig for te water filling allocation to satisfy te service outage constraint, ten it must also be te optimum solution for problem (2). Tus, for a given pair (r o, ɛ), te optimum power allocation sceme gradually canges from te onoff cannel inversion allocation to te water filling allocation as P av increases. Now we examine te connection of te service outage problem wit te outage capacity in [1] and te ergodic capacity in [4]. Te outage capacity C ɛ (P av ) in [1] specifies te maximum supportable rate for a given average power P av wit outage probability ɛ, wic implies tat te basic service rate in tis work must satisfy r o C ɛ (P av ). It is easy to see tat te above condition is equivalent to te feasibility condition (8), tat is, P av P min (r o, H ɛ ). Furtermore, we can see tat te resulting average rate R canges from r o (1 ɛ) to te ergodic capacity wit increasing P av for a given (r o, ɛ). In Fig 3, te expected rate acieved by te service outage approac is plotted against te ergodic capacity and te outage capacity in Rayleig fading cannel wit normalized mean for cannel gain and normalized noise variance. As we can see, for a given outage prob- Outage probability ε P av (db) Fig. 4. Outage probability acieved by te water filling allocation for different basic rate r. ability ɛ =.1, te outage capacity as nearly a 5 db loss in average power compared to te ergodic capacity for a given rate. Between te outage capacity and te ergodic capacity, a number of service outage approaces wit different basic rates exist. Te outage probability for different r acieved by te water filling allocation is also plotted against te service outage approac wit a given ɛ =.1 in Fig 4. It can be observed tat, for a range of P av, te service outage approac acieves a rate very close to te ergodic capacity, and at te same time significantly reduces te outage probability. Hence, te service outage approac strikes good balance between average rate and outage probability. APPENDIX Proof: Teorem 1 Wen P av P min (r o, H a ), problem (3) is feasible and can be, equivalently, translated into te following

6 6 problem: max E {R[p()]} (2) p() subject to E {p()} = P av (2a) p() p() σ2 (2 2ro 1) H a. (2b) (2b) Tis is a standard variational optimization problem [14]. Te objective function is concave on p() and te constraints are linear functions of p(). Ten p (, H a ) is te optimum solution iff it satisfies te te Kun-Tucker conditions [15]. Using a Lagrange multiplier g(p(),, ) = 2log(2)σ 2 >, we define [ ] R[p()] 2 log(2)σ 2 p() f(). (21) Let H denote te set wit non-boundary points as { } H = H a p (, H a ) > σ2 (2 2ro 1) { / H a p (, H a ) > }. (22) It is easy to verify tat wen H te p (, H a ) satisfies dg(p (, H a ),, ) dp (, H a ) oterwise te p (, H a ) satisfies dg(p (, H a ),, ) dp (, H a ) =, (23). (24) Terefore, te p (, H a ) is te optimum solution. Proof: Lemma 1 Power scemes p () and p () differ in te allocation of te residual power. In order to sow tat p () allocates te residual power in a better way tan p (), we define te following power efficiency function for p res() over its strictly positive space. Definition 2 Te power efficiency function η() for p res() is η() = r o R[p wf ()] p res () < 22ro. (25) Te power efficiency function indicates te rate increment corresponding to a unit power assigned from p res (). We ave te following property for η(). Proposition 1 Te power efficiency function η() is a strictly increasing function of over te interval < 2 2ro. Proof: Proposition 1 We consider te cases and separately. For, we ave p wf () = and η() = wic is an increasing function of. r o σ 2 (2 2ro 1), (26) For, (12), (13), and (25) imply η() = r o (1/2) log 2 (/ ( ) ). (27) σ 2 2 2ro 1 We define u() = r o 1/2 log 2 (/ ), so tat η() = ˆη(u()) were ˆη(u) = u σ 2 2 2u (28) 1 It is straigtforward to verify tat ˆη(u) is a strictly decreasing function of u for u. Since u() is a strictly decreasing function of and u() wen 22ro, it follows tat η() is an increasing function of for 2 2ro. We also employ te following proposition for te proof of Lemma 1. Proposition 2 For disjoint sets Ψ and Ψ, let f(x) be an arbitrary function suc tat f(x ) f(x ) for all x Ψ and x Ψ. For any nonnegative function g(x) satisfying Ψ g(x)dx = Ψ g(x)dx, we ave Ψ f(x)g(x)dx Ψ f(x)g(x)dx. Wit tese preliminaries, we now verify te claims of Lemma 1. (a) Equations (16), (17), and (19) imply E {p ()} = E {p ()} = P av. (b) From equations (17) and (18), te service and outage sets of p () are H s (p ()) = { b } and H o (p ()) = { < b } respectively. Terefore, H o (p ()) H s (p ()). (c) Let Ψ = H inv H inv so tat Ψ = H inv \Ψ and Ψ = H inv \Ψ are two disjoint sets and nonempty wen p () p (). Te average rate of p () can be expressed as R = R[p wf ()]f() d + (r o R[p wf()])f() d.(29) H inv Te rate contribution of te water filling component is R wf = R[p wf()]f() d. (3) Since H inv = Ψ Ψ, Definition 2 for te efficiency function η() allows us to write R = R wf + η()p res ()f() d (31) H inv = R wf + η()p res ()f() d Ψ + η()p res ()f() d. (32) Ψ Similarly, H inv = Ψ Ψ, so te average rate for p () can be expressed as R = R wf + η()p res()f() d Ψ + η()p res()f() d. (33) Ψ

7 7 Tus, R R = η()p res()f() d Ψ η()p res()f() d. (34) Ψ Note tat te construction of H inv implies Ψ Ψ. Tat is, for any Ψ and Ψ. By Proposition 1, η() is a strictly increasing function of for < 22ro. Tus, η( ) η( ). Furtermore, equation (19) implies ()f() d = ()f() d. (35) Ψ p res Ψ p res Terefore, te conditions of Proposition 2 are satisfied and we ave R R. (d) From equations (11), (15), (17) and (18), te service sets H s (p ()) and H s (p ()) are disjoint unions given by H s (p ()) = H wf H inv = H wf Ψ Ψ, (36) H s (p ()) = H wf H inv = H wf Ψ Ψ. (37) Tis implies Pr{H s (p ())} Pr{H s (p ())} (38) = Pr{Ψ } Pr{Ψ } (39) 1 = Ψ p res () p res ()f() d 1 Ψ p res ()p res ()f() d. (4) From equations (12) and (13), we observe tat 1/p res() is a increasing function of. Since Ψ Ψ, Proposition 2 implies Pr{H s (p ())} Pr{H s (p ())}. REFERENCES [1] S. Nanda, K. Balacandran, and S. Kumar, Adaptation tecniques in wireless packet data services, IEEE Communications Magazine, vol. 38, no. 1, pp , Jan. 2. [2] ETSI, General Packet Radio Service(GPRS); Mobile Station (MS)- Base Station System(BSS) Interface; Radio Link Control/Medium Access Control (RLC/MAC) Protocol (GSM 4.6), Tec. Spec., vol. 8, no. 1, Nov [3] P. Bender, P. Black, M. Grob, R. Padovani, N. Sindusayana, and A. Viterbi, CDMA/HDR: A Bandwidt-Efficient Hig-Speed Wireless Data Service for Nomadic Users, IEEE Communications Magazine, vol. 38, no. 7, pp. 7 77, July 2. [4] A. J. Goldsmit and P. Varaiya, Capacity of fading cannels wit cannel side information, IEEE Transactions on Information Teory, vol. 43, no. 6, pp , Nov [5] S. V. Hanly and D. N. C. Tse, Multiaccess fading cannels: part II: Delay-limited capacities, IEEE Transactions on Information Teory, vol. 44, no. 7, pp , Nov [6] E. Biglieri, J. Proakis, and S. Samai, Fading cannels: Informationteoretic and communications aspects, IEEE Transactions on Information Teory, vol. 44, no. 6, pp , Oct [7] R. Gallager, Information Teory and Reliable Communication, Jon Wiley and Sons, [8] T. Cover and J. Tomas, Elements of Information Teory, Jon Wiley Sons, Inc., [9] L. H. Ozarow, S. Samai, and A. D. Wyner, Information teoretic considerations for cellular mobile radio, IEEE Transactions on Veicular Tecnology, vol. 43, no. 2, pp , May [1] G. Caire, G. Taricco, and E. Biglieri, Optimum power control over fading cannels, IEEE Transactions on Information Teory, vol. 45, no. 5, pp , July [11] R. Negi, M. Carikar, and J. Cioffi, Minimum outage transmission over fading cannel wit delay constraint, in Proc. of te IEEE International Communication Conference(ICC 99), pp , June 2. [12] S. Verdú, A general formula for cannel capacity, IEEE Transactions on Information Teory, vol. 4, no. 4, pp , July [13] J. Luo, R. Yates, and P. Spasojević, Service outage based capacity and power allocation in parallel fading cannels, in IEEE International Symposium on Information Teory, pp. 18, June 22. [14] J. C. Clegg, Calculus of Variations, Jon Wiley Sons, Inc., [15] M. S. Bazaraa and C. M. Setty, Nonlinear programming teory and algoritms, Jon Wiley and Sons, LIST OF FIGURES 1 (a) Rate allocation R[ˆp()] for policy ˆp(), (b) te improved policy p () given by (16) wit water filling p wf () and residual power 1 ( H inv )p res (), (c) te new policy p () given by (17) wit water filling p wf () and residual power 1 ( H inv )p res () For optimum solution types I-IV, power policies are given on te left and corresponding rate allocation are on te rigt Comparison of service outage approac wit oter capacity notions in te Rayleig fading cannel, for a fixed ɛ = Outage probability acieved by te water filling allocation for different basic rate r Jiangong Luo received er B.S. and M.S. degrees in Electrical Engineering from University of Science and Tecnology of Cina, Hefei, Cina. Since 1997 se as been a P.D. candidate at WINLAB, Rutgers - te State University of New Jersey. Lang Lin (S 97) received is B.S. and M.Pil. degrees in electrical engineering from Beijing University, Cina, and te Hong Kong University of Science and Tecnology, respectively. Since 1997 e as been a P.D. candidate at WINLAB, Rutgers - te State University of New Jersey, New Jersey. Roy Yates received te B.S.E. degree in 1983 from Princeton University, and te S.M. and P.D. degrees in 1986 and 199 from M.I.T., all in Electrical Engineering. Since 199, e as been wit te Wireless Information Networks Laboratory (WINLAB) and te ECE department at Rutgers University. Presently, e serves as an Associate Director of WINLAB and Professor in te ECE Dept. He is a co-autor (wit David Goodman) of te text Probability and Stocastic Processes: A Friendly Introduction for Electrical and Computer Engineers, Jon Wiley and Sons, His researc interests include power control, interference suppression, and media access protocols for wireless communications systems. Predrag Spasojević (M ) received te Diploma of Engineering degree from te Scool of Electrical Engineering, University of Sarajevo, in 199; and Master of Science and Doctor of Pilosopy degrees in electrical engineering from Texas

8 A&M University, College Station, Texas, in 1992 and 1999, respectively. From 2 to 21, e as been wit WINLAB, Rutgers University, as Lucent Post-Doctoral Fellow, were e is currently an Assistant Professor in te Department of Electrical & Computer Engineering. His researc interests are in te general areas of communication teory and signal processing. 8

Multi-User Communication: Capacity, Duality, and Cooperation. Nihar Jindal Stanford University Dept. of Electrical Engineering

Multi-User Communication: Capacity, Duality, and Cooperation. Nihar Jindal Stanford University Dept. of Electrical Engineering Multi-User Communication: Capacity, Duality, and Cooperation Niar Jindal Stanford University Dept. of Electrical Engineering February 3, 004 Wireless Communication Vision Cellular Networks Wireless LAN

More information

DETERMINING the information theoretic capacity of

DETERMINING the information theoretic capacity of IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 6, NO. 7, JULY 007 369 Transactions Letters Outage Capacity and Optimal Power Allocation for Multiple Time-Scale Parallel Fading Channels Subhrakanti

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of

More information

Characterization of Relay Channels Using the Bhattacharyya Parameter

Characterization of Relay Channels Using the Bhattacharyya Parameter Caracterization of Relay Cannels Using te Battacaryya Parameter Josepine P. K. Cu, Andrew W. Ecford, and Raviraj S. Adve Dept. of Electrical and Computer Engineering, University of Toronto, Toronto, Ontario,

More information

Service Outage Based Power and Rate Allocation for Parallel Fading Channels

Service Outage Based Power and Rate Allocation for Parallel Fading Channels 1 Service Outage Based Power and Rate Allocation for Parallel Fading Channels Jianghong Luo, Member, IEEE, Roy Yates, Member, IEEE, and Predrag Spasojević, Member, IEEE Abstract The service outage based

More information

MIMO decorrelation for visible light communication based on angle optimization

MIMO decorrelation for visible light communication based on angle optimization MIMO decorrelation for visible ligt communication based on angle optimization Haiyong Zang, and Yijun Zu Citation: AIP Conference Proceedings 80, 09005 (07); View online: ttps://doi.org/0.03/.4977399 View

More information

lecture 26: Richardson extrapolation

lecture 26: Richardson extrapolation 43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)

More information

Bayesian ML Sequence Detection for ISI Channels

Bayesian ML Sequence Detection for ISI Channels Bayesian ML Sequence Detection for ISI Cannels Jill K. Nelson Department of Electrical and Computer Engineering George Mason University Fairfax, VA 030 Email: jnelson@gmu.edu Andrew C. Singer Department

More information

Dedicated to the 70th birthday of Professor Lin Qun

Dedicated to the 70th birthday of Professor Lin Qun Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,

More information

Copyright c 2008 Kevin Long

Copyright c 2008 Kevin Long Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula

More information

NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL. Georgeta Budura

NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL. Georgeta Budura NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL Georgeta Budura Politenica University of Timisoara, Faculty of Electronics and Telecommunications, Comm. Dep., georgeta.budura@etc.utt.ro Abstract:

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

Combining functions: algebraic methods

Combining functions: algebraic methods Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)

More information

Differentiation in higher dimensions

Differentiation in higher dimensions Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends

More information

Energy efficiency in cognitive radio network: study of cooperative sensing using different channel sensing methods

Energy efficiency in cognitive radio network: study of cooperative sensing using different channel sensing methods Graduate Teses and Dissertations Graduate College 2015 Energy efficiency in cognitive radio network: study of cooperative sensing using different cannel sensing metods Cenxuan Cui Iowa State University

More information

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department

More information

. If lim. x 2 x 1. f(x+h) f(x)

. If lim. x 2 x 1. f(x+h) f(x) Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc

More information

IEOR 165 Lecture 10 Distribution Estimation

IEOR 165 Lecture 10 Distribution Estimation IEOR 165 Lecture 10 Distribution Estimation 1 Motivating Problem Consider a situation were we ave iid data x i from some unknown distribution. One problem of interest is estimating te distribution tat

More information

3. THE EXCHANGE ECONOMY

3. THE EXCHANGE ECONOMY Essential Microeconomics -1-3. THE EXCHNGE ECONOMY Pareto efficient allocations 2 Edgewort box analysis 5 Market clearing prices 13 Walrasian Equilibrium 16 Equilibrium and Efficiency 22 First welfare

More information

Math 161 (33) - Final exam

Math 161 (33) - Final exam Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.

More information

EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS

EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Statistica Sinica 24 2014, 395-414 doi:ttp://dx.doi.org/10.5705/ss.2012.064 EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Jun Sao 1,2 and Seng Wang 3 1 East Cina Normal University,

More information

Order of Accuracy. ũ h u Ch p, (1)

Order of Accuracy. ũ h u Ch p, (1) Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical

More information

Trust Degree Based Beamforming for Multi-Antenna Cooperative Communication Systems

Trust Degree Based Beamforming for Multi-Antenna Cooperative Communication Systems Trust Degree Based Beamforming for Multi-Antenna Cooperative Communication Systems Mojtaba Vaezi, Hazer Inaltekin, Wonjae Sin, H. Vincent Poor, and Junsan Zang* Department of Electrical Engineering, Princeton

More information

Complexity of Decoding Positive-Rate Reed-Solomon Codes

Complexity of Decoding Positive-Rate Reed-Solomon Codes Complexity of Decoding Positive-Rate Reed-Solomon Codes Qi Ceng 1 and Daqing Wan 1 Scool of Computer Science Te University of Oklaoma Norman, OK73019 Email: qceng@cs.ou.edu Department of Matematics University

More information

Chapter 5 FINITE DIFFERENCE METHOD (FDM)

Chapter 5 FINITE DIFFERENCE METHOD (FDM) MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential

More information

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example, NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing

More information

A SHORT INTRODUCTION TO BANACH LATTICES AND

A SHORT INTRODUCTION TO BANACH LATTICES AND CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,

More information

On convexity of polynomial paths and generalized majorizations

On convexity of polynomial paths and generalized majorizations On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal

More information

MANY scientific and engineering problems can be

MANY scientific and engineering problems can be A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial

More information

HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS

HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3

More information

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems 5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

Analysis of an SEIR Epidemic Model with a General Feedback Vaccination Law

Analysis of an SEIR Epidemic Model with a General Feedback Vaccination Law Proceedings of te World Congress on Engineering 5 Vol WCE 5 July - 3 5 London U.K. Analysis of an ER Epidemic Model wit a General Feedback Vaccination Law M. De la en. Alonso-Quesada A.beas and R. Nistal

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x

More information

Chapters 19 & 20 Heat and the First Law of Thermodynamics

Chapters 19 & 20 Heat and the First Law of Thermodynamics Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,

More information

How to Find the Derivative of a Function: Calculus 1

How to Find the Derivative of a Function: Calculus 1 Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te

More information

MVT and Rolle s Theorem

MVT and Rolle s Theorem AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state

More information

Optimal Power Control in Decentralized Gaussian Multiple Access Channels

Optimal Power Control in Decentralized Gaussian Multiple Access Channels 1 Optimal Power Control in Decentralized Gaussian Multiple Access Channels Kamal Singh Department of Electrical Engineering Indian Institute of Technology Bombay. arxiv:1711.08272v1 [eess.sp] 21 Nov 2017

More information

The derivative function

The derivative function Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative

More information

Exercises for numerical differentiation. Øyvind Ryan

Exercises for numerical differentiation. Øyvind Ryan Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can

More information

Optimal Cross Layer Scheduling of Transmissions over a Fading Multiaccess Channel

Optimal Cross Layer Scheduling of Transmissions over a Fading Multiaccess Channel Optimal Cross Layer Sceduling of Transmissions over a Fading Multiaccess Cannel Munis Goyal, Anurag Kumar, Vinod Sarma Department of Electrical Communication Engg Indian Institute of Science, Bangalore,

More information

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as

More information

3.2 THE FUNDAMENTAL WELFARE THEOREMS

3.2 THE FUNDAMENTAL WELFARE THEOREMS Essential Microeconomics -1-3.2 THE FUNDMENTL WELFRE THEOREMS Walrasian Equilibrium 2 First welfare teorem 3 Second welfare teorem (conve, differentiable economy) 12 Te omotetic preference 2 2 economy

More information

1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).

1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t). . Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use

More information

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator. Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative

More information

Click here to see an animation of the derivative

Click here to see an animation of the derivative Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,

More information

Adaptive Neural Filters with Fixed Weights

Adaptive Neural Filters with Fixed Weights Adaptive Neural Filters wit Fixed Weigts James T. Lo and Justin Nave Department of Matematics and Statistics University of Maryland Baltimore County Baltimore, MD 150, U.S.A. e-mail: jameslo@umbc.edu Abstract

More information

ch (for some fixed positive number c) reaching c

ch (for some fixed positive number c) reaching c GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan

More information

Taylor Series and the Mean Value Theorem of Derivatives

Taylor Series and the Mean Value Theorem of Derivatives 1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential

More information

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225 THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:

More information

Poisson Equation in Sobolev Spaces

Poisson Equation in Sobolev Spaces Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on

More information

On the Identifiability of the Post-Nonlinear Causal Model

On the Identifiability of the Post-Nonlinear Causal Model UAI 9 ZHANG & HYVARINEN 647 On te Identifiability of te Post-Nonlinear Causal Model Kun Zang Dept. of Computer Science and HIIT University of Helsinki Finland Aapo Hyvärinen Dept. of Computer Science,

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat

More information

Learning based super-resolution land cover mapping

Learning based super-resolution land cover mapping earning based super-resolution land cover mapping Feng ing, Yiang Zang, Giles M. Foody IEEE Fellow, Xiaodong Xiuua Zang, Siming Fang, Wenbo Yun Du is work was supported in part by te National Basic Researc

More information

Effect of the Dependent Paths in Linear Hull

Effect of the Dependent Paths in Linear Hull 1 Effect of te Dependent Pats in Linear Hull Zenli Dai, Meiqin Wang, Yue Sun Scool of Matematics, Sandong University, Jinan, 250100, Cina Key Laboratory of Cryptologic Tecnology and Information Security,

More information

CORRELATION TEST OF RESIDUAL ERRORS IN FREQUENCY DOMAIN SYSTEM IDENTIFICATION

CORRELATION TEST OF RESIDUAL ERRORS IN FREQUENCY DOMAIN SYSTEM IDENTIFICATION IAC Symposium on System Identification, SYSID 006 Marc 9-3 006, Newcastle, Australia CORRELATION TEST O RESIDUAL ERRORS IN REQUENCY DOMAIN SYSTEM IDENTIICATION István Kollár *, Ri Pintelon **, Joan Scouens

More information

Differential Calculus (The basics) Prepared by Mr. C. Hull

Differential Calculus (The basics) Prepared by Mr. C. Hull Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit

More information

Characterization of reducible hexagons and fast decomposition of elementary benzenoid graphs

Characterization of reducible hexagons and fast decomposition of elementary benzenoid graphs Discrete Applied Matematics 156 (2008) 1711 1724 www.elsevier.com/locate/dam Caracterization of reducible exagons and fast decomposition of elementary benzenoid graps Andrej Taranenko, Aleksander Vesel

More information

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation

More information

Some Review Problems for First Midterm Mathematics 1300, Calculus 1

Some Review Problems for First Midterm Mathematics 1300, Calculus 1 Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,

More information

Linearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization

Linearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization Linearized Primal-Dual Metods for Linear Inverse Problems wit Total Variation Regularization and Finite Element Discretization WENYI TIAN XIAOMING YUAN September 2, 26 Abstract. Linear inverse problems

More information

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0 3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,

More information

3.4 Worksheet: Proof of the Chain Rule NAME

3.4 Worksheet: Proof of the Chain Rule NAME Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are

More information

Higher Derivatives. Differentiable Functions

Higher Derivatives. Differentiable Functions Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.

More information

The Priestley-Chao Estimator

The Priestley-Chao Estimator Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are

More information

3.1 Extreme Values of a Function

3.1 Extreme Values of a Function .1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find

More information

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a? Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is

More information

Preconditioning in H(div) and Applications

Preconditioning in H(div) and Applications 1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition

More information

HT TURBULENT NATURAL CONVECTION IN A DIFFERENTIALLY HEATED VERTICAL CHANNEL. Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008

HT TURBULENT NATURAL CONVECTION IN A DIFFERENTIALLY HEATED VERTICAL CHANNEL. Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008 Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008 August 10-14, 2008, Jacksonville, Florida USA Proceedings of HT2008 2008 ASME Summer Heat Transfer Conference August 10-14, 2008, Jacksonville,

More information

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to

More information

HOMEWORK HELP 2 FOR MATH 151

HOMEWORK HELP 2 FOR MATH 151 HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,

More information

Precalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!

Precalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here! Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.

More information

Continuous Stochastic Processes

Continuous Stochastic Processes Continuous Stocastic Processes Te term stocastic is often applied to penomena tat vary in time, wile te word random is reserved for penomena tat vary in space. Apart from tis distinction, te modelling

More information

Equilibrium and Pareto Efficiency in an exchange economy

Equilibrium and Pareto Efficiency in an exchange economy Microeconomic Teory -1- Equilibrium and efficiency Equilibrium and Pareto Efficiency in an excange economy 1. Efficient economies 2 2. Gains from excange 6 3. Edgewort-ox analysis 15 4. Properties of a

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

Online Learning: Bandit Setting

Online Learning: Bandit Setting Online Learning: Bandit Setting Daniel asabi Summer 04 Last Update: October 0, 06 Introduction [TODO Bandits. Stocastic setting Suppose tere exists unknown distributions ν,..., ν, suc tat te loss at eac

More information

Continuity and Differentiability of the Trigonometric Functions

Continuity and Differentiability of the Trigonometric Functions [Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te

More information

Continuity and Differentiability Worksheet

Continuity and Differentiability Worksheet Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;

More information

5.1 We will begin this section with the definition of a rational expression. We

5.1 We will begin this section with the definition of a rational expression. We Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go

More information

LECTURE 14 NUMERICAL INTEGRATION. Find

LECTURE 14 NUMERICAL INTEGRATION. Find LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use

More information

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra

More information

2.3 Product and Quotient Rules

2.3 Product and Quotient Rules .3. PRODUCT AND QUOTIENT RULES 75.3 Product and Quotient Rules.3.1 Product rule Suppose tat f and g are two di erentiable functions. Ten ( g (x)) 0 = f 0 (x) g (x) + g 0 (x) See.3.5 on page 77 for a proof.

More information

Tail Conditional Expectations for Extended Exponential Dispersion Models

Tail Conditional Expectations for Extended Exponential Dispersion Models American Researc Journal of Matematics Original Article ISSN 378-704 Volume 1 Issue 4 015 Tail Conditional Expectations for Extended Exponential Dispersion Models Ye (Zoe) Ye Qiang Wu and Don Hong 1 Program

More information

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY (Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative

More information

Functions of the Complex Variable z

Functions of the Complex Variable z Capter 2 Functions of te Complex Variable z Introduction We wis to examine te notion of a function of z were z is a complex variable. To be sure, a complex variable can be viewed as noting but a pair of

More information

Bob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk

Bob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of

More information

New Streamfunction Approach for Magnetohydrodynamics

New Streamfunction Approach for Magnetohydrodynamics New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite

More information

7.1 Using Antiderivatives to find Area

7.1 Using Antiderivatives to find Area 7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between

More information

Complexity of Decoding Positive-Rate Primitive Reed-Solomon Codes

Complexity of Decoding Positive-Rate Primitive Reed-Solomon Codes 1 Complexity of Decoding Positive-Rate Primitive Reed-Solomon Codes Qi Ceng and Daqing Wan Abstract It as been proved tat te maximum likeliood decoding problem of Reed-Solomon codes is NP-ard. However,

More information

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,

More information

Subdifferentials of convex functions

Subdifferentials of convex functions Subdifferentials of convex functions Jordan Bell jordan.bell@gmail.com Department of Matematics, University of Toronto April 21, 2014 Wenever we speak about a vector space in tis note we mean a vector

More information

Section 3.1: Derivatives of Polynomials and Exponential Functions

Section 3.1: Derivatives of Polynomials and Exponential Functions Section 3.1: Derivatives of Polynomials and Exponential Functions In previous sections we developed te concept of te derivative and derivative function. Te only issue wit our definition owever is tat it

More information

Flavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.

Flavius Guiaş. X(t + h) = X(t) + F (X(s)) ds. Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical

More information

Minimum degree thresholds for bipartite graph tiling

Minimum degree thresholds for bipartite graph tiling Minimum degree tresolds for bipartite grap tiling Albert Bus Yi Zao Department of Matematics and Statistics Georgia State University January 17, 011 Abstract Given a bipartite grap H and a positive integer

More information

Chapter 2 Performance Analysis of Call-Handling Processes in Buffered Cellular Wireless Networks

Chapter 2 Performance Analysis of Call-Handling Processes in Buffered Cellular Wireless Networks Capter 2 Performance Analysis of Call-Handling Processes in Buffered Cellular Wireless Networks In tis capter effective numerical computational procedures to calculate QoS (Quality of Service) metrics

More information

The Laplace equation, cylindrically or spherically symmetric case

The Laplace equation, cylindrically or spherically symmetric case Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,

More information

2.11 That s So Derivative

2.11 That s So Derivative 2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point

More information

Efficient algorithms for for clone items detection

Efficient algorithms for for clone items detection Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire

More information