Packet Video 99. A Corruption Model for Motion Compensated Video Subject to Bit Errors. Gustavo de los Reyes Amy R. Reibman Shih-Fu Chang

Size: px
Start display at page:

Download "Packet Video 99. A Corruption Model for Motion Compensated Video Subject to Bit Errors. Gustavo de los Reyes Amy R. Reibman Shih-Fu Chang"

Transcription

1 A Corruptio Model for Motio Compesated Video Subject to Bit Errors Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag AT&T Labs AT&T Labs-Research Columbia Uiversity Abstract: We preset a model of how corruptio propagates i a block-based-ecoded video decoder whe the bitstream is subjected to bit errors. Our Markov model takes ito accout both corruptio itroduced via bit errors ad corruptio propagated via motio compesatio. The model icorporates the impact of spatial ad temporal source resiliece o limitig the propagatio of corruptio withi a frame ad across time. For a give video sequece, the model requires source-specific iputs characterizig motio ad spatial eergy. It also requires ecoder-specific iputs characterizig resiliece slice legth ad the proportio of I-blocks per frame. Oe fial parameter eeded by the model is the average umber of bits i a compressed block. We apply the model to the specific example of determiig the optimal allocatio of bit rate amog spatial resiliece, temporal resiliece, ad source rate i a trascoder that ijects resiliece ito the compressed bitstream.. Itroductio Desigig video codig systems for oisy chaels is a difficult problem that has bee receivig much attetio due to the icrease i the use of wireless devices. These efforts have resulted i recet video codig stadards for low-bitrate commuicatios that iclude methods for improvig spatial ad temporal source resiliece i additio to the usual chael codig methods [2]. Spatial ad temporal localizatio techiques limit the propagatio of errors withi a frame or to subsequet frames. Despite the umber of algorithms that have bee proposed for improvig the resiliece of video over oisy chael [4, 6], few aalytical models have bee proposed to support these ad other algorithms. I this paper, we derive a Markov model that describes how corruptio, iitially caused by bit errors durig trasport, propagates i videos that are coded usig motio compesatio. Our model takes ito accout the effects of spatial ad temporal resiliece as a meas for limitig the propagatio of the corruptio. Budagavi [] preseted a model for error propagatio usig a restrictive motio model where errors affect oly future blocks i the same spatial locatio as the origial error. I this paper, we allow urestricted block-based motio compesatio as defied i H.263 [2] ad similar block-based ecodig stadards. Hece the model is applicable to a wide rage of compressed video sequeces. I additio, little work has bee doe to combie the various methods of resiliece ad provide guidace as to how best to combie the methods give a fixed bitrate budget. I [3], we described a trascoder that ijects combied spatial ad temporal source resiliece ito the compressed bitstream at the wireless base statio or mobile switch. I this paper, we use our Markov model to derive operatioal rate distortio fuctios for spatial ad temporal source resiliece, ad use these to determie the optimal bitrate allocatio amog source rate, spatial resiliece, ad temporal resiliece. We use software simulatio to verify the results obtaied from our model. 2. Resiliece I our model, we focus o techiques for spatial ad temporal localizatio. These methods prevet sigal errors caused by a bit error or packet loss from propagatig withi a video frame or to subsequet frames. Spatial error propagatio ca occur whe the decoder loses sychroizatio while decodig the variable legth codes. I this paper, we limit spatial error propagatio by isertig additioal sychroizatio headers to reduce the umber of blocks i a slice. With a shorter slice, the decoder re-sychroizes more quickly, resultig i less lost data. For temporal localizatio, we trasmit more itra-frames (I-frames) or I-blocks. These frames or blocks are coded to be temporally idepedet Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. of

2 of previous frames or blocks ad are used as refereces for subsequet temporal predictio. More frequet I-blocks reduce the duratio of error propagatio. The temporal ad spatial localizatio techiques improve resiliece by icreasig the overall bitrate. Therefore, we also use rate-reductio techiques to recover some bit-rate to allocate for improvig resiliece. Rate reductio techiques fall ito two geeral classes: requatizig coefficiets ad discardig coefficiets. I this paper, we take the simple approach of discardig coefficiets to reduce the rate i a give frame by usig a zoal mask, where the size of the mask is selected to provide the required rate reductio. It is the cotetio for bit rate betwee resiliece ad source rate that leads us to use rate distortio theory as a guide for selectig optimal bit rates as discussed below. 3. Overview of Corruptio Propagatio Model I this sectio we describe the motivatio for usig a Markov model, ad we describe the basic compoets of the model. Blocks withi a frame ca be corrupted for two mai reasos: bit errors occur durig trasport, or corrupted iformatio from the previous frame is used durig motio compesatio. Therefore, the umber of corrupted blocks i the curret frame oly depeds o the iput errors for the curret frame, the motio vectors i the curret frame, ad the umber of corrupted blocks i the previous frame. This fits aturally ito the structure of a Markov model. The overall corruptio propagatio model calculates the umber of corrupted blocks i each video frame which is the sum of the umber of blocks corrupted due to bit errors ad the umber of blocks corrupted due to motio compesatio mius the itersectio betwee the two types of corrupted blocks. The umber of corrupted blocks, X, withi frame,, is the: ( X X ) X X, mc + X, loss, mc, loss = (Eq. ) where mc idicates corruptio due to motio compesatio, ad loss idicates corruptio due to bit errors. For blocks that are lost durig trasport due to bit errors, X,loss, we derive a statistical model for the expected umber of blocks lost as a fuctio of the spatial resiliece (slice legth) ad the bit error rate (BER) i the chael. For blocks that are corrupted durig motio compesatio, X,mc, we derive a Markov model that icludes the probability that a block is corrupted by referecig a block that has previously bee corrupted. Withi each frame, the probability that a block is corrupted by motio compesatio is treated as a Beroulli process with probability of success, p. This probability depeds heavily o the motio i the video ad the corruptio i the previous frame. Sice the corruptio i the previous frame icludes blocks that were lost to bit errors, the Markov model icludes the effects of the spatial resiliece model. We will show that this Markov model also icludes the effects of temporal resiliece. To fold the spatial model ito the Markov model, we simplify Equatio ad use the expected value of X,loss,. Equatio the becomes. X = X, mc X where r = F + E [ X, loss ] ( X, mc E[ X, loss ]) X E[ X ] = re[ X ] We estimate the itersectio as,,, ad F mc is the total umber of loss loss blocksi a frame. (Eq. 2) I the ext sectios we describes the spatial model ad the Markov model i more detail. 4. Spatial Resiliece Model I this sectio we derive a statistical model of the video quality as a fuctio of the spatial resiliece (as measured by slice legth) ad the chael coditios (as measured by mea bit error rate Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. 2 of

3 (BER)). This model takes ito accout the corruptio propagatio withi a frame caused by bit errors durig trasport. This spatial resiliece model is oe part of the overall error propagatio model. 4. Derivatio of Spatial Resiliece Model The importat parameters i this model are slice legth, BER, ad the umber of bits per block. The slice legth is measured i macroblocks as described i [2]. I this paper, we use the more geeric term blocks rather tha the specific term macroblocks. However, durig implemetatio, we use the specific defiitio of macroblocks i [2]. For this aalysis, we oly use mea BER as the idicator of wireless chael coditios ad assume bit errors are idepedet i time. We assume that oce the decoder detects a error, it discards the block cotaiig the error ad all subsequet blocks util the ext slice header. We further assume that every bit error results i a error that is detected by the decoder. This overestimates the umber of lost blocks but does ot characterize the distortio caused by errors that caot be detected from video sytax. For the spatial resiliece model, we use a probabilistic model for the expected block losses. The probability that block m ad all subsequet blocks i the slice are lost is equal to the probability that o bit errors occur i blocks to (m-) ad at least oe bit error occurs i block m. Therefore, if we assume a costat block legth i bits, the expected umber of blocks lost per slice is: [ ] ib+ c ( i+ ) b+ c {( legth i)( P ) ( P ) } X legth, loss E = e e slice i= where X, loss are the blocks lost, i is the block positio, legth is the legth of the slice i blocks, Pe is the probability of error (or BER), c is the header legth, ad b is the costat legth of the block size i bits. (Eq. 3) 4.2 Effect of Spatial Resiliece The result of the spatial resiliece model is show i Figure, which shows how addig more slice headers (reducig the slice legth) leads to a lower Mea Square Error (MSE) per video frame, especially at higher BER. This figure also shows strog o-liearities across both slice legth ad BER. Oe of our goals is to take advatage of these o-liearities to select a operatig poit for spatial resiliece that produces the best video quality i the presece of bit errors. 3 MSE (per video frame) Slice legth (blocks) Log(BER) -3-2 Figure. Effect of spatial resiliece ad wireless chael coditios o block loss. Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. 3 of

4 5. Overall Corruptio Propagatio Model I this sectio we derive the Markov model that describes the propagatio of corruptio due to motio compesatio. Also, we iclude the results of the spatial model i the motio compesatio model ad derive a overall corruptio-propagatio model. 5. Corruptio due to Motio Compesatio We begi this derivatio by defiig a trasitio matrix that govers the propagatio of corruptio i frame + durig motio compesatio give a umber of corrupted blocks i frame. As discussed previously, our model treats the probability of corruptio due to motio compesatio as a Beroulli process with probability of success, p. The motio compesatio trasitio matrix is the defied by a biomial distributio: k j k j P( i, jmc ) = P{ X +, mc = jmc X = i} = p mcq mc j mc where X is the umber of corrupted blocks i frame, (Eq. 4) p is the probability a block is corrupted via motio compesatio, q = p, k is the umber of blocks i a frame, ad i, jmc =,,, k. Note that this equatio describes how motio compesatio propagates corruptio from the previous (referece) frame. This equatio does ot specify how the corruptio occurred i the previous frame. That will be described i the ext sectios. 5.2 Corruptio due to Motio Compesatio ad Bit Errors Recall that Equatio 2 describes how we model corruptio due joitly to bit errors, X,loss, ad motio compesatio, X,mc. By substitutig Equatio 2 ito Equatio 4 ad solvig for X +, we derive the trasitio matrix that govers the overall propagatio of errors from frame to frame. k j k j P( i, j) = P{ X + = j X = i} = p q j where X is the umber of corrupted blocks i frame, (Eq. 5) p is the probability a block is corrupted via motio compesatio, q = p, k is the umber of blocks i a frame, ad i, j =,,, k. Recall that P(i, j) is ideed a fuctio of i because p is a fuctio of the corruptio i the referece frame,. We will derive a equatio for p shortly. Figure 2 shows schematically how the spatial model ad the motio-compesated Markov model are used together to describe the overall corruptio propagatio. The large squares idicate video frames, ad the smaller squares idicate lost blocks withi the frames. The letters, M, ad L are used to idicate et corrupted blocks, corruptio due to motio compesatio, ad losses due to chael bit errors, respectively. I the iitial video frame, losses are itroduced due to bit errors. These losses are propagated via motio compesatio as described i Equatio 4. Accordig to Equatio 2, ew losses are itroduced i the ext frame due to bit errors, ad itersectig losses are removed. This results i the et corruptio due to both bit-error losses ad motio compesatio. However, by usig the trasitio matrix of Equatio 5, the et result from Frame to Frame ca be calculated directly. The the process repeats for the followig frame. Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. 4 of

5 X X =E[X,loss ] X,mc X,mc + E[X,loss ] = X,mc + E[X,loss ] -X,mc E[X,loss ] M M ML ML M M M M L M L M M M Frame Frame E[X,loss ] Motio Compesatio P(i, j mc ) E[X,loss ] Motio Compesatio P(i, j mc ) P(i, j) Figure 2. Schematic propagatio of corruptio via bit errors ad motio compesatio. Now that we have the trasitio matrix, we ca determie the expected umber of corrupted blocks due to motio compesatio ad bit errors i each video frame,, from the -step trasitio probabilities: P { X = j X = i} ( i, j) = P (Eq. 6) The expected umber of corrupted blocks due to motio compesatio i each frame i the video is: E [ X X = i] = jp{ X = j X = i} = j=,,, k jp j=,,, k where i is ( i, j) the umber of corrupted blocks i the first frame,. (Eq. 7) From Equatio 7, oe ca see how the model describes the error propagatio frame by frame. The corruptio i each frame,, is just a fuctio of the -step trasitio probabilities for that frame. Now we determie the probability, p, that a block is corrupted through motio compesatio so that we ca compute the trasitio matrix, P(i, j) of Equatio 5. A block i the curret frame will be corrupted through motio compesatio if it refereces corrupted iformatio (pixels) i the previous frame. Sice we are workig with blocks as the uit of measure rather tha pixels, this is equivalet to sayig that a block i the curret frame will be corrupted if it refereces a corrupted block, ad it refereces corrupted pixels withi that corrupted block. To compute these probabilities, we group the blocks ito four categories characterized by their motio vectors (MVs). The first category is for blocks that have MVs that describe o motio. The secod category is for blocks that have MVs that describe horizotal or vertical motio. The third category is for blocks that have MVs that describe combied horizotal ad vertical motio. The last category is for I-blocks that have o associated MVs. For each category of block, we compute the probability, p i, that a block of that type is corrupted via motio compesatio. The overall probability of corruptio for a block give a certai amout of corruptio i the previous frame is the: p = p P( block ) (Eq. 8) i i=,2,3,4 i Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. 5 of

6 where i is oe of the four possible categories i = is a block of o motio, i = 2 is a block of horizotal or vertical motio, i = 3 is a block of combied horizotal ad vertical motio, ad i = 4 is a I-block. For I-blocks, the probability of a I-block beig corrupted by motio compesatio is clearly zero. O the other had, for a block with o motio, the probability of corruptio is oe if it refereces a corrupted block. For these blocks, the probability of corruptio is simply the probability that a block i the refereced frame is corrupted, which we estimate by the proportio of blocks i the previous frame that are corrupted, r. Calculatig the probability of corruptio for the other two categories of blocks is somewhat more ivolved. These blocks ca referece either two or four blocks i the previous frame, ad various combiatios of corrupted ad ot corrupted blocks ca exist for these refereced blocks. Furthermore, to be corrupted, a block must still referece the corrupted iformatio i a corrupted block. For these two categories of blocks, the probability that a block refereces oe or more corrupt blocks, p ref, ca be described by a biomial distributio: l l j l j l pref = pcorr ( pcorr ) = ( pcorr ) j= j where pcorr is the probability a block i the referece frame is corrupted, ad l is the umber of blocks that a curret block refereces : l = for o motio, l = 2 for vertical or horizotal motio, ad l = 4 for combied vertical ad horizotal motio.. 9) As before, we estimate the probability that a block i the refereced frame is corrupted by the proportio of blocks i the refereced frame that are corrupted, p corr = r. Thus far, we have estimated the probability that a block i the curret frame refereces a corrupted block. Now we estimate the probability that a block i the curret frame refereces corrupted iformatio give that it refereces a corrupted block. For each block, we assume a uiform distributio of corruptio from miimum to complete corruptio. The, the fial equatio for the probability that a block is corrupted through motio compesatio, p, is ow: (Eq p = a rp( block ad ) + a 2 2 [ ( r) ] P( block ) + a ( r) 2 4 [ ] P( block P( blocki ) = =,2,3,4 (Eq. ) i 3 3 ) where block a ( = ), a 2 i are the categories of blocks defied above, ad, ad a blocks of category block, 3 are the probability of referecig corrupt iformatio i block, ad block 2 3, respectively. 5.3 Effect of Temporal Resiliece I this sectio we show how our model icorporates the effect of usig additioal I-blocks for temporal resiliece. I Equatio, sice the probability of all four categories of blocks must sum to oe, as the proportio of I-blocks, P(block 4 ), icreases, the probability of receivig the other blocks, P(block ) + P(block 2 ) + P(block 3 ), decreases. This decreases p, the probability that a block is corrupted via motio compesatio. This behavior is show i Figure 3, which shows the ivariat (steady state) Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. 6 of

7 distributio of the Markov model with a BER = 5x -4, ad two differet probability of I blocks, P(block 4 ). O average, there are fewer corrupted blocks as P(block 4 ) icreases. P(X = j) [Ivariat (steady state) distributio] P(I-block)=.5, BER = 5x -4 P(I-block)=.25, BER = 5x j (Number of blocks corrupted) Figure 3. Ivariat distributio of Markov model as a fuctio of P(block 4 )=probability of I blocks. 6. Optimal Allocatio of Bit Rate I the rest of the paper we show oe way that we use the corruptio-propagatio models: to optimally allocate resiliece. Give that we ca allocate bit rate over three sources source rate, spatial resiliece, ad temporal resiliece we would like to allocate the bit rate optimally to miimize distortio i the video. Applyig rate distortio theory to our resiliece problem, we wat to fid a vector of rates which miimizes a overall distortio measure subject to a rate costrait, B. Solvig this costraied optimizatio problem is equivalet to solvig the ucostraied optimizatio problem with a objective fuctio, H: H = D( L) + λr( L) where L is the vector of resources, λ is the Lagragia multiplier, D is the distortio that results from applyig the resources, ad R is the rate that results from applyig the resources. (Eq. ) The objective fuctio ca be writte as: [ Di ( Li ) + λri ( Li )] H (Eq. 2) = N i= Solvig this set of simultaeous equatios leads to: D DN λ = = = R RN (Eq. 3) subject to R( L) B That is, give that there is at least oe solutio to the problem, the solutios exist at the poits where the slopes of the distortio-rate fuctios D i (R i ) are equal, ad the sum of the rates R i are less tha or equal to the available bit rate. I our case, the vector of resources icludes the spatial resiliece, the temporal Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. 7 of

8 resiliece, ad the source rate. We have a additioal costrait i that the resources vary i discrete amouts. For example, slice legths vary i iteger values from oe to 22. Therefore, the bit rates used by the resources also vary i discrete amouts. This discreteess affects the optimal solutio. We cosider the optimal solutio to be the oe where the total bit rate comes closest to the available bit rate of the chael. Also, we defie the rate distortio fuctios as operatioal fuctios that are geerated from the Markov models or from simulatios. 7. Trascoder I [3] we exteded the use of trascodig beyod rate reductio [5], to iclude the additio of resiliece iformatio i the bitstream. The trascoder operatio is show i Figure 4. The trascoder decodes the icomig bitstream to the degree required to add resiliece. Resiliece is added by the trascoder, ad, if ecessary, the bit-rate is reduced. The the resiliet bitstream is re-quatized (for temporal resiliece oly) ad variable legth re-ecoded. Icomig ecoded bitstream Decode video, as required Add spatial ad/or temporal resiliece Reduce rate of bitstream Figure 4. High-level trascoder operatio. Resiliet bitstream at ear-origial rate. Re-quatize ad re-ecode, as required The spatial localizatio techique of addig more slice headers requires oly simple bitstream parsig ad arithmetic operatios. A variable-legth decoder parses the bitstream ad iserts additioal slice start codes where ecessary. I additio, motio vectors (MVs) are decoded, a ew differetial MV is calculated from referece MVs oly withi the curret slice, ad the result is variable legth ecoded. The trascoder requires more complicated processig to icrease the temporal resiliece by usig more frequet I-blocks. The icomig bitstream is fully decoded, icludig motio compesatio, to create the ew I-blocks. The coefficiets of the ew I-blocks are the DCT-coded, re-quatized, scaed, ad variable-legth ecoded. Because the trascoder ca modify the resiliece to match the prevailig chael coditios, a method to determie the optimal resiliece is ivaluable. 8. Results The operatioal rate distortio fuctio (ORDF) for source rate was geerated from simulatio, while the ORDF for spatial ad temporal resiliece were geerated from our aalytical models. I additio, the models estimate block losses ad corruptio that must be coverted to distortio to produce the ORDF. This coversio is a fuctio of the video cotet, so we estimated the distortio by simulatig losses i the decoded video. For each of the ORDF i the figures below, we curve fitted the data ad the differetiated the fit curves to get smooth slopes as a fuctios of bit rate. The figures below show the origial ORDF with the superimposed curve fit o oe plot, ad the correspodig slope fuctios o the adjacet plot. Figures 5-7 show the operatioal rate distortio fuctios ad their slopes at a BER of 5x -4 for source rate, spatial resiliece (slice legth), ad temporal resiliece (I-block proportio), respectively. The y-axis for the slope of the spatial resiliece ORDF (Figure 6) has bee scaled to highlight the area of iterest. For source rate, the x-axes are the ratio of the resultig bit rate whe coefficiets are dropped to the origial bit rate with all coefficiets. For spatial resiliece ad temporal resiliece, the x-axes are the ratio of the resultig bit rate whe resiliece is added to the baselie resiliece (mius ). Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. 8 of

9 Distortio (MSE) Bit rate decrease due to Dropped Coefficets Slope d(distortio)/d(rate) Bit rate decrease due to Dropped Coefficiets Figure 5. Operatioal rate distortio fuctio ad slope fuctio for source rate (from simulatio). 4 Distortio (MSE) Slope d(distortio)/d(rate) Bit rate icrease due to Slices (relative to GOB) Bit rate icrease due to Slices Figure 6. Operatioal rate distortio fuctio ad slope fuctio for spatial resiliece (aalytical model). 8-5 Distortio (MSE) Slope d(distortio)/d(rate) Bit rate icrease due to I-blocks Bit rate icrease due to I-blocks Figure 7. Operatioal rate distortio fuctio ad slope fuctio for temporal resiliece (Markov model). The operatig poit is selected by tryig to match slopes such that the total bit rate is less tha or equal to the bit rate of the origial video. However, the possible solutios are further costraied because the resources vary i discrete amouts. I additio, the effect of MSE due to source rate is perceptually differet tha the effect of MSE due to resiliece. Decreasig the source rate teds to blur the video, while decreasig spatial ad temporal resiliece teds to icrease the umber of severe artifacts. For these severe artifacts, several blocks of the video may chage color, or parts of the video become mismatched. Therefore, a certai amout of judgemet must be applied whe selectig the optimal solutio. Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. 9 of

10 I [3] we used a trial-ad-error approach to select the operatig poit, while here we described a optimal method for selectig the operatig poit. The optimal operatig poit for a BER of 5x -4 correspods to a slice legth of 4 blocks, a I-block proportio equal to each block beig refreshed every 5 frames, ad a source rate such that a maximum of coefficiets are retaied i each block. This results i a Peak Sigal to Noise Ratio (PSNR) of 8. db. The results preseted i [3] for the same BER resulted i a PSNR of 7.4 db for a slice legth of 5 blocks, a I-block proportio equal to each block beig refreshed every frames, ad a source rate such that a maximum 5 coefficiets are retaied i each block. The result for the baselie video with omial resiliece is a PSNR of 3.5 db. For omial resiliece we use a slice legth of 22 blocks ad a sigle I-frame followed by all P-frames. We also computed the optimal resiliece for a BER of -4. The operatig poit with the optimal method is a slice legth of blocks, a I-block proportio equal to each block beig refreshed every 7 frames, ad a source rate such that a maximum 6 coefficiets are retaied i each block. This results i a Peak Sigal to Noise Ratio (PSNR) of 2.4 db. This is slightly less tha the PSNR of 2.8 db preseted i [3] for the same BER. However, both methods lead to much better results tha the baselie video with omial resiliece ad a PSNR of 8.5 db. Therefore, we ca achieve early the same results with the optimal approach as compared to our previous trial-ad-error approach. This is a great beefit sice the optimal approach depeds more o aalytical models ad less o the specific video sequece. This makes the optimal approach better suited to ear-real-time calculatios that ca be used i a resiliece trascoder. 9. Summary ad Coclusios We developed a aalytical model for the propagatio of corruptio as a fuctio of spatial ad temporal resiliece ad BER. For a give video sequece, the model requires source-specific iputs characterizig motio ad spatial eergy. The spatial eergy of the video is used to estimate the distortio (MSE) that results from block losses ad corruptio. For a system implemetatio, a likely place to estimate this distortio is at a scee chage. The eergy cotet per frame betwee scees chages will be reasoably steady, ad methods exist for determiig scee chages i real time. The model also takes as iput the ecoder-specific parameters of slice legth ad I-block proportio per frame. We apply the model to the specific example of determiig the optimal allocatio of bit rate amog spatial resiliece, temporal resiliece, ad source rate i a trascoder that ijects resiliece ito the compressed bitstream. The optimal approach leads to early the same results tha previously reported whe usig a trial-ad-error approach.. Refereces [] Budagavi, M. ad J. D. Gibso. Error Propagatio i Motio Compesated Video over Wireless Chaels, IEEE Iteratioal Coferece o Image Processig, October 997. [2] Draft ITU-T Recommedatio H.263, Video Codig for Low Bitrate Commuicatio, Sept [3] De los Reyes, G., A. R. Reibma, J. C. Chuag, ad S.-F. Chag, Video Trascodig for Resiliece i Wireless Chaels, IEEE Iteratioal Coferece o Image Processig, October 998. [4] Farber. N., E. Steiback, ad B. Girod, Robust H.263 Compatible Video Trasmissio over Wireless Chaels, Iteratioal Picture Codig Symposium, Melboure, Australia, March 996. [5] Keesma, G., R. Hellighuize, F. Hoeksema, G. Heidma, Trascodig of MPEG Bitstreams, Sigal Processig: Image Commuicatios, 8 (996), pp [6] Liao, J. Y. ad J. D. Villaseor, Adaptive Itra Update for Video Codig over Noisy Chaels. IEEE Iteratioal Coferece o Image Processig, September 996. Gustavo de los Reyes Amy R. Reibma Shih-Fu Chag 8/26/99 p. of

Time-Domain Representations of LTI Systems

Time-Domain Representations of LTI Systems 2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable

More information

Run-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE

Run-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE Geeral e Image Coder Structure Motio Video (s 1,s 2,t) or (s 1,s 2 ) Natural Image Samplig A form of data compressio; usually lossless, but ca be lossy Redudacy Removal Lossless compressio: predictive

More information

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution EEL5: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we begi our mathematical treatmet of discrete-time s. As show i Figure, a discrete-time operates or trasforms some iput sequece x [

More information

Reliability and Queueing

Reliability and Queueing Copyright 999 Uiversity of Califoria Reliability ad Queueig by David G. Messerschmitt Supplemetary sectio for Uderstadig Networked Applicatios: A First Course, Morga Kaufma, 999. Copyright otice: Permissio

More information

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows

More information

CHAPTER 10 INFINITE SEQUENCES AND SERIES

CHAPTER 10 INFINITE SEQUENCES AND SERIES CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece

More information

NUMERICAL METHODS FOR SOLVING EQUATIONS

NUMERICAL METHODS FOR SOLVING EQUATIONS Mathematics Revisio Guides Numerical Methods for Solvig Equatios Page 1 of 11 M.K. HOME TUITION Mathematics Revisio Guides Level: GCSE Higher Tier NUMERICAL METHODS FOR SOLVING EQUATIONS Versio:. Date:

More information

Cooperative Communication Fundamentals & Coding Techniques

Cooperative Communication Fundamentals & Coding Techniques 3 th ICACT Tutorial Cooperative commuicatio fudametals & codig techiques Cooperative Commuicatio Fudametals & Codig Techiques 0..4 Electroics ad Telecommuicatio Research Istitute Kiug Jug 3 th ICACT Tutorial

More information

Error Resilience Analysis of Multi-Hypothesis Motion Compensated Prediction for Video Coding

Error Resilience Analysis of Multi-Hypothesis Motion Compensated Prediction for Video Coding 1. Research Team Error Resiliece Aalysis of Multi-Hypothesis Motio Compesated Predictio for Video Codig Project Leader: Other Faculty: Post Doc(s): Graduate Studets: Idustrial Parter(s): C.-C. Jay Kuo,

More information

ASYMPTOTIC CLOSED-LOOP DESIGN OF ERROR RESILIENT PREDICTIVE COMPRESSION SYSTEMS. Sina Zamani, Tejaswi Nanjundaswamy, Kenneth Rose

ASYMPTOTIC CLOSED-LOOP DESIGN OF ERROR RESILIENT PREDICTIVE COMPRESSION SYSTEMS. Sina Zamani, Tejaswi Nanjundaswamy, Kenneth Rose ASYMPTOTIC OSED-LOOP DESIGN OF ERROR RESILIENT PREDICTIVE COMPRESSION SYSTEMS Sia Zamai, Tejaswi Najudaswamy, Keeth Rose Departmet of Electrical ad Computer Egieerig, Uiversity of Califoria Sata Barbara,

More information

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) = AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,

More information

Encoding-Assisted Temporal Direct Mode Decision for B Pictures in H.264/AVC

Encoding-Assisted Temporal Direct Mode Decision for B Pictures in H.264/AVC Ecodig-Assisted Temporal Direct Mode Decisio for B Pictures i H.64/AVC Ya-Neg Fag ad Yiyi Li Hui-Jae Hsieh Departmet of Commuicatio Egieerig Departmet of Iformatio Maagemet Natioal Cetral Uiversity, Taiwa

More information

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

GUIDELINES ON REPRESENTATIVE SAMPLING

GUIDELINES ON REPRESENTATIVE SAMPLING DRUGS WORKING GROUP VALIDATION OF THE GUIDELINES ON REPRESENTATIVE SAMPLING DOCUMENT TYPE : REF. CODE: ISSUE NO: ISSUE DATE: VALIDATION REPORT DWG-SGL-001 002 08 DECEMBER 2012 Ref code: DWG-SGL-001 Issue

More information

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at

More information

Vector Quantization: a Limiting Case of EM

Vector Quantization: a Limiting Case of EM . Itroductio & defiitios Assume that you are give a data set X = { x j }, j { 2,,, }, of d -dimesioal vectors. The vector quatizatio (VQ) problem requires that we fid a set of prototype vectors Z = { z

More information

Sensitivity Analysis of Daubechies 4 Wavelet Coefficients for Reduction of Reconstructed Image Error

Sensitivity Analysis of Daubechies 4 Wavelet Coefficients for Reduction of Reconstructed Image Error Proceedigs of the 6th WSEAS Iteratioal Coferece o SIGNAL PROCESSING, Dallas, Texas, USA, March -4, 7 67 Sesitivity Aalysis of Daubechies 4 Wavelet Coefficiets for Reductio of Recostructed Image Error DEVINDER

More information

ANALYSIS OF EXPERIMENTAL ERRORS

ANALYSIS OF EXPERIMENTAL ERRORS ANALYSIS OF EXPERIMENTAL ERRORS All physical measuremets ecoutered i the verificatio of physics theories ad cocepts are subject to ucertaities that deped o the measurig istrumets used ad the coditios uder

More information

Definitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients.

Definitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients. Defiitios ad Theorems Remember the scalar form of the liear programmig problem, Miimize, Subject to, f(x) = c i x i a 1i x i = b 1 a mi x i = b m x i 0 i = 1,2,, where x are the decisio variables. c, b,

More information

FIR Filter Design: Part II

FIR Filter Design: Part II EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple sigle-peak

More information

SNAP Centre Workshop. Basic Algebraic Manipulation

SNAP Centre Workshop. Basic Algebraic Manipulation SNAP Cetre Workshop Basic Algebraic Maipulatio 8 Simplifyig Algebraic Expressios Whe a expressio is writte i the most compact maer possible, it is cosidered to be simplified. Not Simplified: x(x + 4x)

More information

ECE 564/645 - Digital Communication Systems (Spring 2014) Final Exam Friday, May 2nd, 8:00-10:00am, Marston 220

ECE 564/645 - Digital Communication Systems (Spring 2014) Final Exam Friday, May 2nd, 8:00-10:00am, Marston 220 ECE 564/645 - Digital Commuicatio Systems (Sprig 014) Fial Exam Friday, May d, 8:00-10:00am, Marsto 0 Overview The exam cosists of four (or five) problems for 100 (or 10) poits. The poits for each part

More information

UC Berkeley CS 170: Efficient Algorithms and Intractable Problems Handout 17 Lecturer: David Wagner April 3, Notes 17 for CS 170

UC Berkeley CS 170: Efficient Algorithms and Intractable Problems Handout 17 Lecturer: David Wagner April 3, Notes 17 for CS 170 UC Berkeley CS 170: Efficiet Algorithms ad Itractable Problems Hadout 17 Lecturer: David Wager April 3, 2003 Notes 17 for CS 170 1 The Lempel-Ziv algorithm There is a sese i which the Huffma codig was

More information

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig

More information

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015 ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],

More information

Linear Regression Demystified

Linear Regression Demystified Liear Regressio Demystified Liear regressio is a importat subject i statistics. I elemetary statistics courses, formulae related to liear regressio are ofte stated without derivatio. This ote iteds to

More information

Filter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and

Filter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and Filter bas Separately, the lowpass ad highpass filters are ot ivertible T removes the highest frequecy / ad removes the lowest frequecy Together these filters separate the sigal ito low-frequecy ad high-frequecy

More information

OPTIMAL ALGORITHMS -- SUPPLEMENTAL NOTES

OPTIMAL ALGORITHMS -- SUPPLEMENTAL NOTES OPTIMAL ALGORITHMS -- SUPPLEMENTAL NOTES Peter M. Maurer Why Hashig is θ(). As i biary search, hashig assumes that keys are stored i a array which is idexed by a iteger. However, hashig attempts to bypass

More information

The Random Walk For Dummies

The Random Walk For Dummies The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli

More information

Optimally Sparse SVMs

Optimally Sparse SVMs A. Proof of Lemma 3. We here prove a lower boud o the umber of support vectors to achieve geeralizatio bouds of the form which we cosider. Importatly, this result holds ot oly for liear classifiers, but

More information

Module 5 EMBEDDED WAVELET CODING. Version 2 ECE IIT, Kharagpur

Module 5 EMBEDDED WAVELET CODING. Version 2 ECE IIT, Kharagpur Module 5 EMBEDDED WAVELET CODING Versio ECE IIT, Kharagpur Lesso 4 SPIHT algorithm Versio ECE IIT, Kharagpur Istructioal Objectives At the ed of this lesso, the studets should be able to:. State the limitatios

More information

Chapter 7: The z-transform. Chih-Wei Liu

Chapter 7: The z-transform. Chih-Wei Liu Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability

More information

Summary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram.

Summary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram. Key Cocepts: 1) Sketchig of scatter diagram The scatter diagram of bivariate (i.e. cotaiig two variables) data ca be easily obtaied usig GC. Studets are advised to refer to lecture otes for the GC operatios

More information

6.867 Machine learning, lecture 7 (Jaakkola) 1

6.867 Machine learning, lecture 7 (Jaakkola) 1 6.867 Machie learig, lecture 7 (Jaakkola) 1 Lecture topics: Kerel form of liear regressio Kerels, examples, costructio, properties Liear regressio ad kerels Cosider a slightly simpler model where we omit

More information

U8L1: Sec Equations of Lines in R 2

U8L1: Sec Equations of Lines in R 2 MCVU U8L: Sec. 8.9. Equatios of Lies i R Review of Equatios of a Straight Lie (-D) Cosider the lie passig through A (-,) with slope, as show i the diagram below. I poit slope form, the equatio of the lie

More information

OPTIMAL PIECEWISE UNIFORM VECTOR QUANTIZATION OF THE MEMORYLESS LAPLACIAN SOURCE

OPTIMAL PIECEWISE UNIFORM VECTOR QUANTIZATION OF THE MEMORYLESS LAPLACIAN SOURCE Joural of ELECTRICAL EGIEERIG, VOL. 56, O. 7-8, 2005, 200 204 OPTIMAL PIECEWISE UIFORM VECTOR QUATIZATIO OF THE MEMORYLESS LAPLACIA SOURCE Zora H. Perić Veljo Lj. Staović Alesadra Z. Jovaović Srdja M.

More information

Castiel, Supernatural, Season 6, Episode 18

Castiel, Supernatural, Season 6, Episode 18 13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio

More information

Machine Learning for Data Science (CS 4786)

Machine Learning for Data Science (CS 4786) Machie Learig for Data Sciece CS 4786) Lecture & 3: Pricipal Compoet Aalysis The text i black outlies high level ideas. The text i blue provides simple mathematical details to derive or get to the algorithm

More information

subject to A 1 x + A 2 y b x j 0, j = 1,,n 1 y j = 0 or 1, j = 1,,n 2

subject to A 1 x + A 2 y b x j 0, j = 1,,n 1 y j = 0 or 1, j = 1,,n 2 Additioal Brach ad Boud Algorithms 0-1 Mixed-Iteger Liear Programmig The brach ad boud algorithm described i the previous sectios ca be used to solve virtually all optimizatio problems cotaiig iteger variables,

More information

Session 5. (1) Principal component analysis and Karhunen-Loève transformation

Session 5. (1) Principal component analysis and Karhunen-Loève transformation 200 Autum semester Patter Iformatio Processig Topic 2 Image compressio by orthogoal trasformatio Sessio 5 () Pricipal compoet aalysis ad Karhue-Loève trasformatio Topic 2 of this course explais the image

More information

Lecture 10: Universal coding and prediction

Lecture 10: Universal coding and prediction 0-704: Iformatio Processig ad Learig Sprig 0 Lecture 0: Uiversal codig ad predictio Lecturer: Aarti Sigh Scribes: Georg M. Goerg Disclaimer: These otes have ot bee subjected to the usual scrutiy reserved

More information

Optimization Methods MIT 2.098/6.255/ Final exam

Optimization Methods MIT 2.098/6.255/ Final exam Optimizatio Methods MIT 2.098/6.255/15.093 Fial exam Date Give: December 19th, 2006 P1. [30 pts] Classify the followig statemets as true or false. All aswers must be well-justified, either through a short

More information

Principle Of Superposition

Principle Of Superposition ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give

More information

Information Theory Model for Radiation

Information Theory Model for Radiation Joural of Applied Mathematics ad Physics, 26, 4, 6-66 Published Olie August 26 i SciRes. http://www.scirp.org/joural/jamp http://dx.doi.org/.426/jamp.26.487 Iformatio Theory Model for Radiatio Philipp

More information

Scheduling under Uncertainty using MILP Sensitivity Analysis

Scheduling under Uncertainty using MILP Sensitivity Analysis Schedulig uder Ucertaity usig MILP Sesitivity Aalysis M. Ierapetritou ad Zheya Jia Departmet of Chemical & Biochemical Egieerig Rutgers, the State Uiversity of New Jersey Piscataway, NJ Abstract The aim

More information

( ) = p and P( i = b) = q.

( ) = p and P( i = b) = q. MATH 540 Radom Walks Part 1 A radom walk X is special stochastic process that measures the height (or value) of a particle that radomly moves upward or dowward certai fixed amouts o each uit icremet of

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

Infinite Sequences and Series

Infinite Sequences and Series Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet

More information

Kinetics of Complex Reactions

Kinetics of Complex Reactions Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet

More information

Chapter 7 z-transform

Chapter 7 z-transform Chapter 7 -Trasform Itroductio Trasform Uilateral Trasform Properties Uilateral Trasform Iversio of Uilateral Trasform Determiig the Frequecy Respose from Poles ad Zeros Itroductio Role i Discrete-Time

More information

Double Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution

Double Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution Iteratioal Mathematical Forum, Vol., 3, o. 3, 3-53 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.9/imf.3.335 Double Stage Shrikage Estimator of Two Parameters Geeralized Expoetial Distributio Alaa M.

More information

The Riemann Zeta Function

The Riemann Zeta Function Physics 6A Witer 6 The Riema Zeta Fuctio I this ote, I will sketch some of the mai properties of the Riema zeta fuctio, ζ(x). For x >, we defie ζ(x) =, x >. () x = For x, this sum diverges. However, we

More information

Section 5.1 The Basics of Counting

Section 5.1 The Basics of Counting 1 Sectio 5.1 The Basics of Coutig Combiatorics, the study of arragemets of objects, is a importat part of discrete mathematics. I this chapter, we will lear basic techiques of coutig which has a lot of

More information

Random Models. Tusheng Zhang. February 14, 2013

Random Models. Tusheng Zhang. February 14, 2013 Radom Models Tusheg Zhag February 14, 013 1 Radom Walks Let me describe the model. Radom walks are used to describe the motio of a movig particle (object). Suppose that a particle (object) moves alog the

More information

w (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ.

w (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ. 2 5. Weighted umber of late jobs 5.1. Release dates ad due dates: maximimizig the weight of o-time jobs Oce we add release dates, miimizig the umber of late jobs becomes a sigificatly harder problem. For

More information

CSE 202 Homework 1 Matthias Springer, A Yes, there does always exist a perfect matching without a strong instability.

CSE 202 Homework 1 Matthias Springer, A Yes, there does always exist a perfect matching without a strong instability. CSE 0 Homework 1 Matthias Spriger, A9950078 1 Problem 1 Notatio a b meas that a is matched to b. a < b c meas that b likes c more tha a. Equality idicates a tie. Strog istability Yes, there does always

More information

1 Introduction to reducing variance in Monte Carlo simulations

1 Introduction to reducing variance in Monte Carlo simulations Copyright c 010 by Karl Sigma 1 Itroductio to reducig variace i Mote Carlo simulatios 11 Review of cofidece itervals for estimatig a mea I statistics, we estimate a ukow mea µ = E(X) of a distributio by

More information

The Maximum-Likelihood Decoding Performance of Error-Correcting Codes

The Maximum-Likelihood Decoding Performance of Error-Correcting Codes The Maximum-Lielihood Decodig Performace of Error-Correctig Codes Hery D. Pfister ECE Departmet Texas A&M Uiversity August 27th, 2007 (rev. 0) November 2st, 203 (rev. ) Performace of Codes. Notatio X,

More information

Increasing timing capacity using packet coloring

Increasing timing capacity using packet coloring 003 Coferece o Iformatio Scieces ad Systems, The Johs Hopkis Uiversity, March 4, 003 Icreasig timig capacity usig packet colorig Xi Liu ad R Srikat[] Coordiated Sciece Laboratory Uiversity of Illiois e-mail:

More information

10-701/ Machine Learning Mid-term Exam Solution

10-701/ Machine Learning Mid-term Exam Solution 0-70/5-78 Machie Learig Mid-term Exam Solutio Your Name: Your Adrew ID: True or False (Give oe setece explaatio) (20%). (F) For a cotiuous radom variable x ad its probability distributio fuctio p(x), it

More information

7. Modern Techniques. Data Encryption Standard (DES)

7. Modern Techniques. Data Encryption Standard (DES) 7. Moder Techiques. Data Ecryptio Stadard (DES) The objective of this chapter is to illustrate the priciples of moder covetioal ecryptio. For this purpose, we focus o the most widely used covetioal ecryptio

More information

1 of 7 7/16/2009 6:06 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 6. Order Statistics Defiitios Suppose agai that we have a basic radom experimet, ad that X is a real-valued radom variable

More information

ADVANCED SOFTWARE ENGINEERING

ADVANCED SOFTWARE ENGINEERING ADVANCED SOFTWARE ENGINEERING COMP 3705 Exercise Usage-based Testig ad Reliability Versio 1.0-040406 Departmet of Computer Ssciece Sada Narayaappa, Aeliese Adrews Versio 1.1-050405 Departmet of Commuicatio

More information

6.003 Homework #3 Solutions

6.003 Homework #3 Solutions 6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the

More information

Orthogonal Gaussian Filters for Signal Processing

Orthogonal Gaussian Filters for Signal Processing Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios

More information

Math 2784 (or 2794W) University of Connecticut

Math 2784 (or 2794W) University of Connecticut ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really

More information

The target reliability and design working life

The target reliability and design working life Safety ad Security Egieerig IV 161 The target reliability ad desig workig life M. Holický Kloker Istitute, CTU i Prague, Czech Republic Abstract Desig workig life ad target reliability levels recommeded

More information

ELEC1200: A System View of Communications: from Signals to Packets Lecture 3

ELEC1200: A System View of Communications: from Signals to Packets Lecture 3 ELEC2: A System View of Commuicatios: from Sigals to Packets Lecture 3 Commuicatio chaels Discrete time Chael Modelig the chael Liear Time Ivariat Systems Step Respose Respose to sigle bit Respose to geeral

More information

September 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1

September 2012 C1 Note. C1 Notes (Edexcel) Copyright   - For AS, A2 notes and IGCSE / GCSE worksheets 1 September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright

More information

Chapter 9: Numerical Differentiation

Chapter 9: Numerical Differentiation 178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical

More information

Invariability of Remainder Based Reversible Watermarking

Invariability of Remainder Based Reversible Watermarking Joural of Network Itelligece c 16 ISSN 21-8105 (Olie) Taiwa Ubiquitous Iformatio Volume 1, Number 1, February 16 Ivariability of Remaider Based Reversible Watermarkig Shao-Wei Weg School of Iformatio Egieerig

More information

True Nature of Potential Energy of a Hydrogen Atom

True Nature of Potential Energy of a Hydrogen Atom True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial

More information

Introduction to Signals and Systems, Part V: Lecture Summary

Introduction to Signals and Systems, Part V: Lecture Summary EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive

More information

Number Representation

Number Representation Number Represetatio 1 Number System :: The Basics We are accustomed to usig the so-called decimal umber system Te digits :: 0,1,2,3,4,5,6,7,8,9 Every digit positio has a weight which is a power of 10 Base

More information

Disjoint set (Union-Find)

Disjoint set (Union-Find) CS124 Lecture 7 Fall 2018 Disjoit set (Uio-Fid) For Kruskal s algorithm for the miimum spaig tree problem, we foud that we eeded a data structure for maitaiig a collectio of disjoit sets. That is, we eed

More information

CHAPTER I: Vector Spaces

CHAPTER I: Vector Spaces CHAPTER I: Vector Spaces Sectio 1: Itroductio ad Examples This first chapter is largely a review of topics you probably saw i your liear algebra course. So why cover it? (1) Not everyoe remembers everythig

More information

DESCRIPTION OF THE SYSTEM

DESCRIPTION OF THE SYSTEM Sychroous-Serial Iterface for absolute Ecoders SSI 1060 BE 10 / 01 DESCRIPTION OF THE SYSTEM TWK-ELEKTRONIK GmbH D-001 Düsseldorf PB 1006 Heirichstr. Tel +9/11/6067 Fax +9/11/6770 e-mail: ifo@twk.de Page

More information

A statistical method to determine sample size to estimate characteristic value of soil parameters

A statistical method to determine sample size to estimate characteristic value of soil parameters A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig

More information

Polynomial Functions and Their Graphs

Polynomial Functions and Their Graphs Polyomial Fuctios ad Their Graphs I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively

More information

Chapter 6. Advanced Counting Techniques

Chapter 6. Advanced Counting Techniques Chapter 6 Advaced Coutig Techiques 6.: Recurrece Relatios Defiitio: A recurrece relatio for the sequece {a } is a equatio expressig a i terms of oe or more of the previous terms of the sequece: a,a2,a3,,a

More information

Analysis of Algorithms. Introduction. Contents

Analysis of Algorithms. Introduction. Contents Itroductio The focus of this module is mathematical aspects of algorithms. Our mai focus is aalysis of algorithms, which meas evaluatig efficiecy of algorithms by aalytical ad mathematical methods. We

More information

An Introduction to Randomized Algorithms

An Introduction to Randomized Algorithms A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis

More information

A New Solution Method for the Finite-Horizon Discrete-Time EOQ Problem

A New Solution Method for the Finite-Horizon Discrete-Time EOQ Problem This is the Pre-Published Versio. A New Solutio Method for the Fiite-Horizo Discrete-Time EOQ Problem Chug-Lu Li Departmet of Logistics The Hog Kog Polytechic Uiversity Hug Hom, Kowloo, Hog Kog Phoe: +852-2766-7410

More information

( ) ( ), (S3) ( ). (S4)

( ) ( ), (S3) ( ). (S4) Ultrasesitivity i phosphorylatio-dephosphorylatio cycles with little substrate: Supportig Iformatio Bruo M.C. Martis, eter S. Swai 1. Derivatio of the equatios associated with the mai model From the differetial

More information

Chapter 4. Fourier Series

Chapter 4. Fourier Series Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,

More information

CSI 2101 Discrete Structures Winter Homework Assignment #4 (100 points, weight 5%) Due: Thursday, April 5, at 1:00pm (in lecture)

CSI 2101 Discrete Structures Winter Homework Assignment #4 (100 points, weight 5%) Due: Thursday, April 5, at 1:00pm (in lecture) CSI 101 Discrete Structures Witer 01 Prof. Lucia Moura Uiversity of Ottawa Homework Assigmet #4 (100 poits, weight %) Due: Thursday, April, at 1:00pm (i lecture) Program verificatio, Recurrece Relatios

More information

A. Basics of Discrete Fourier Transform

A. Basics of Discrete Fourier Transform A. Basics of Discrete Fourier Trasform A.1. Defiitio of Discrete Fourier Trasform (8.5) A.2. Properties of Discrete Fourier Trasform (8.6) A.3. Spectral Aalysis of Cotiuous-Time Sigals Usig Discrete Fourier

More information

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + 62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of

More information

Complex Number Theory without Imaginary Number (i)

Complex Number Theory without Imaginary Number (i) Ope Access Library Joural Complex Number Theory without Imagiary Number (i Deepak Bhalchadra Gode Directorate of Cesus Operatios, Mumbai, Idia Email: deepakm_4@rediffmail.com Received 6 July 04; revised

More information

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018 CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical

More information

ECONOMIC OPERATION OF POWER SYSTEMS

ECONOMIC OPERATION OF POWER SYSTEMS ECOOMC OEATO OF OWE SYSTEMS TOUCTO Oe of the earliest applicatios of o-lie cetralized cotrol was to provide a cetral facility, to operate ecoomically, several geeratig plats supplyig the loads of the system.

More information

Feedback in Iterative Algorithms

Feedback in Iterative Algorithms Feedback i Iterative Algorithms Charles Byre (Charles Byre@uml.edu), Departmet of Mathematical Scieces, Uiversity of Massachusetts Lowell, Lowell, MA 01854 October 17, 2005 Abstract Whe the oegative system

More information

ENGI 9420 Engineering Analysis Assignment 3 Solutions

ENGI 9420 Engineering Analysis Assignment 3 Solutions ENGI 9 Egieerig Aalysis Assigmet Solutios Fall [Series solutio of ODEs, matri algebra; umerical methods; Chapters, ad ]. Fid a power series solutio about =, as far as the term i 7, to the ordiary differetial

More information

Quadratic Functions. Before we start looking at polynomials, we should know some common terminology.

Quadratic Functions. Before we start looking at polynomials, we should know some common terminology. Quadratic Fuctios I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively i mathematical

More information

Ray-triangle intersection

Ray-triangle intersection Ray-triagle itersectio ria urless October 2006 I this hadout, we explore the steps eeded to compute the itersectio of a ray with a triagle, ad the to compute the barycetric coordiates of that itersectio.

More information

Section 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations

Section 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations Differece Equatios to Differetial Equatios Sectio. Calculus: Areas Ad Tagets The study of calculus begis with questios about chage. What happes to the velocity of a swigig pedulum as its positio chages?

More information

Numerical Methods in Fourier Series Applications

Numerical Methods in Fourier Series Applications Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets

More information

Chapter 6 Principles of Data Reduction

Chapter 6 Principles of Data Reduction Chapter 6 for BST 695: Special Topics i Statistical Theory. Kui Zhag, 0 Chapter 6 Priciples of Data Reductio Sectio 6. Itroductio Goal: To summarize or reduce the data X, X,, X to get iformatio about a

More information

Chapter 2: Numerical Methods

Chapter 2: Numerical Methods Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,

More information

Chapter 6 Infinite Series

Chapter 6 Infinite Series Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat

More information