Optimal estimation of clock values and trends from finite data
|
|
- Ralph York
- 5 years ago
- Views:
Transcription
1 Opimal esimaion of clock values and rends from finie daa Charles Greenhall Je Propulsion Laboraory, California Insiue of Technology Pasadena, CA USA Absrac We show how o solve wo problems of opimal linear esimaion from a finie se of phase daa. Clock noise is modeled as a sochasic process wih saionary dh incremens. The covariance properies of such a process are conained in he generalized auocovariance funcion (GACV). We se up wo principles for opimal esimaion; hese principles lead o a se of linear equaions for he regression coefficiens and some auxiliary parameers. The mean square errors of he esimaors are easily calculaed. The mehod can be used o check he resuls of oher mehods and o find good subopimal esimaors based on a small subse of he available daa. I. INTRODUCTION Suppose ha he phase residual x () of a clock (or difference of clocks) is given a a finie se of imes, T = { 1,..., n }, which don have o be equally spaced. We consider wo esimaion arges: 1) he phase value x ( ) a some / T ; 2) he overall rend coefficien of x (). Depending on he model for x () and wha else is known, his can mean he long-erm average value of phase or of is dh ime derivaive: frequency, drif rae, or aging rae (d =1, 2, or P 3). For each arge, we wan o calculae he linear esimaor n i=1 a ix ( i ) ha is opimal in he mean-square sense while saisfying an invariance condiion ha will be explained below. To carry ou his program, we need a sochasic model for x (). The general model used here is he class of sochasic processes wih saionary dh incremens, he subjec of a monograph of Yaglom [1]. These include saionary processes, indefinie inegrals of saionary processes, and all he powerlaw processes familiar o he ime and frequency field: flicker PM, whie FM, flicker FM, and so on. Yaglom showed how o solve problems of opimal predicion and filering for hese processes from heir values on unbounded or bounded ime inervals. For pure power-law processes, oher auhors have derived he mean square error (MSE) of he predicor ha is based on he infinie pas [2] [4]; Boulanger and Douglas [5] calculaed he MSE of wo-poin linear exrapolaion. Vernoe e al. [6] calculaed predicors and heir MSEs based on exrapolaion of leas-squares linear and quadraic fis of equally spaced finie daa ses. For finie-sae clock models, a recursive opimal predicor and is MSE, based on all discreeime pas measuremens, can be calculaed from a Kalman filer [7]. The MSE of various subopimal drif-rae esimaors has been calculaed for power-law noises [8] [11]. Here we show how o calculae he regression coefficiens a i and he MSE of he opimal linear esimaors of boh arges, using sysems of linear equaions ha generalize he equaions of orhogonal projecion for saionary x (). II. CLOCK NOISE MODELS A real-valued, mean-square coninuous 1 sochasic process x () is said o have saionary dh incremens (d 1) iffor each τ he process d µ d d τ x () = ( 1) k x ( kτ) (1) k k=0 is saionary. I is convenien o le SI (d) denoe he class of all such processes, and o le SI (0) denoe he saionary processes. Then SI (d) SI (d +1), and we define he degree of x () as he leas d such ha x SI (d). Everyhing we know abou a process x SI (d) is wrapped up in he dh incremens (1), which do no change if we add a polynomial of degree d 1 o x (). In his sense, a process in SI (d) is ambiguous. Any use of hese processes mus ake accoun of his ambiguiy. For any x SI (d), Yaglom esablished a nonnegaive specral densiy funcion 2 S x (f) ha exends he noion of specral densiy for saionary processes. Here, we are using he wo-sided, even version: S x ( f) =S x (f). Ford 1, S x (f) can diverge as f 0 bu obeys he resricions Z Z f 2d S x (f) df <, S x (f) df <. (2) f <1 f >1 The process also has an average rend coefficien, denoed by c d, ha can be defined as he infinie-ime average of he saionary process d 1x (): c d = lim Z 2 1 d 1x () d, he limi being aken in he mean-square sense. In his reamen, c d can be a random variable. We ofen wan o ge rid of c d, and here are wo ways o do i. Firs, if we know c d,hen we can consider he process x 0 () =x () c d d /d!, whichis also in SI (d) bu has rend coefficien zero. Second, if we don know c d, hen we can rea x () as a member of SI (d +1); as such, is rend coefficien c d+1 is always zero. For example, a saionary process wih an unknown mean can be reaed as 1 This means ha E[x (u) x ()] 2 0 as u. 2 Acually, a measure on he puncured frequency axis f 6= /05/$ IEEE. 377
2 amemberofsi (1). If a random walk of phase (whie FM), which is in SI (1), has an unknown slope (frequency) added o i, we can rea i as a member of SI (2). Inhiswaywe can ge resuls ha are invarian o he unknown rend. Yaglom s heory of SI (d) was based on he specral densiy. For saionary processes, we also have he auocovariance (ACV) funcion which saisfies s x () = Z e i2πf S x (f) df, (3) s x ( u) =Ex () x (u). (4) The esimaion mehods given here for SI (d) are based on a generalized ACV (GACV) funcion [12][13], also wrien as s x (), ha can be obained from S x (f) and c d by a generalized Fourier inegral as follows: Z " # 2d 1 s x () = e i2πf (i2πf) k S x (f) df f 1 k! k=0 Z + e i2πf S x (f) df + E c 2 2 d d. (5) f >1 (2d)! One may add a polynomial of degree 2d 1 o s x () wihou changing he value of any formula in which s x () is properly used. Wih his ambiguiy undersood, Table I gives he GACV for power-law componens of clock noise, specified by he onesided specral densiy of frequency, S y + (f) =2(2πf) 2 S x (f). For any saionary phase componen, he GACV is he ACV. The flicker PM enry is obained by passing pure 1/f noise hrough a moving-average filer of widh τ o saisfy he second condiion in (2). The rend coefficiens are zero. The degree of a sum of noises is he maximum degree of he summands, and he GACV of he sum of orhogonal noises is he sum of he GACVs. By (5) we may and will assume ha he GACV is an even funcion: s x ( ) =s x (). TABLE I GENERALIZED AUTOCOVARIANCE OF POWER-LAW NOISES Name degree S + y (f) s x () Flicker PM 1 h 1 f sin2 (πfτ) (πfτ) 2 h 1 3 4π 2 2 ln τ, =0 (lowpass) h 1 4π 2 ln, À τ Whie FM 1 h 0 h 0 4 Flicker FM 2 h 1 f 1 h 1 2 ln 2 Random walk FM 2 h 2 f 2 h 2 π Flicker walk FM 3 h 3 f 3 h 3π 2 4 ln 6 Random run FM 3 h 4 f 4 h 4π 4 5 Now we have o say how he GACV is used. In he following discussion, a () and b () denoe funcions ha are zero excep on some unspecified finie se of imes; by using his noaion 30 we can freely perform and combine sums over wihou worrying abou he range of he summaions. Such a funcion a () is said o saisfy he momen condiion of order d [6][14] if a () j =0, j =0,...,d 1. (6) In oher words, if x () in P a () x () is replaced by a polynomial p () of degree d 1, he resul is zero. The coefficiens of he dh incremen (1) saisfy his condiion. All he covariances needed for he esimaion problems can be calculaed by he following heorem. Theorem 1: Le s x () be he GACV of x (), a process wih saionary dh incremens, where d 1. Ifa () and b () saisfy he momen condiion of order d, hen Ã!Ã! E a () x () b () x () = a () b (u) s x ( u). (7),u For saionary processes, (7) follows from (4), and no momen condiions are needed. For d 1, we are only allowed o ake he covariance of linear combinaions of x () whose coefficiens saisfy he momen condiion. According o he heorem, we may do so as if (4) were rue. In realiy, s x () is no an ACV, and (4) does no hold. Neverheless, he enire formula (7) is correc, even hough he corresponding erms of he expansions of is lef and righ sides are no equal. P We define M d (x) as he se of random variables a () x () whose coefficiens a () saisfy he momen condiion of order d. ThenM d (x) is a linear subspace, and Theorem 1 ells us how o calculae he covariance of wo members of M d (x). I can be shown ha any member of M d (x) is a mean-square limi of linear combinaions of dh incremens of x (); for his reason, if c d =0hen he members of M d (x) have zero expecaion. We also define M d (x, T ) as hose members of M d (x) whose coefficiens are suppored on he finie se T. III. CLOCK PREDICTION Le x () have saionary dh incremens. The esimaion arge is x ( ), and he esimaors are of form ˆx ( )= a () x (), (8) where T = { 1,..., n } wih n d. (Alhough his problem is called predicion, here is no need o insis ha > max i.). We now have o make he problem invarian o he ambiguiy of x () wih respec o polynomials p () of degree d 1. We do so by insising ha he esimaion error x ( ) ˆx ( ) should no change if x () is replaced by x ()+p (). Thus he expression x ( )+p ( ) a ()[x ()+p ()] 378
3 should no depend on p (). Consequenly, a () p () =p ( ) (9) for any polynomial p () of degree d 1; he esimaor predics p ( ) perfecly from he values p (), T. For d =2, if we add a consan phase and frequency o x (), he esimae ˆx ( ) should auomaically adjus iself so ha he error does no change. An equivalen saemen of his invariance condiion is ha x ( ) ˆx ( ) M d (x). (10) We may also express (9) as a () j = j, j =0,...,d 1. (11) In marix form, Ga = g, where a = a ( 1 ) a ( n ) T, (12) n G =.., (13) d 1 1 n d 1 g = 1 d 1 T. (14) A random variable ˆx ( ) of form (8) ha saisfies he invariance condiion is called a linear invarian esimaor (LIE) of order d of x ( ) based on x (), T. Since n d, alieoforderd exiss, namely, he value a of he inerpolaing polynomial deermined by x ( 1 ),...,x( d ).For d =2, his LIE is he linear exrapolaor of x ( 1 ) and x ( 2 ) o, which was reaed by Boulanger and Douglas [5]. Some of he GACVs of Table I appear in heir formulas. Supposing ha 1 < 2 <, one can ofen opimize 2 1 o give an MSE ha is accepably close o he minimum MSE ha one can ge from LIEs based on he infinie pas 2 [4]. The linear and quadraic exrapolaors discussed by Vernoe e al. [6] are LIEs of order 2 and 3, respecively. The ACVs abulaed in [6] are relaed o he GACVs used here. Of all he order-d LIEs of x ( ), we wan he one wih he smalles MSE, called he bes LIE (BLIE) of order d. These of LIEs, call i {LIE}, is deermined by he inhomogeneous equaions (11) for a (). The difference of any wo LIEs is a random variable P b () x () whose coefficiens saisfy he corresponding homogenous equaions, ha is, i belongs o he linear subspace M d (x, T ). Thus, {LIE} is an affine se, a shifed version of M d (x, T ) ha does no pass hrough he origin. To find he closes poin of {LIE} o x ( ) in he mean-square sense, we drop a perpendicular from x ( ) o {LIE}. Thus,ˆx ( ) has o be a member of {LIE} and also has o saisfy he orhogonaliy condiion which means ha x ( ) ˆx ( ) M d (x, T ), (15) E[x ( ) ˆx ( )] Y =0whenever Y M d (x, T ). (16) Figure 1 should make he geomery clear. I can be proved from basic facs abou orhogonal projecions ha he BLIE ˆx ( ) exiss and is unique. Fig. 1. Md { LIE} ( x, T) x( * ) ˆx( * ) 0 Geomeric inerpreaion of he condiions for a BLIE Le x (T )= x ( 1 ) x ( n ) T T, b = b1 b n. (17) Then (16) can be rewrien as E x ( ) a T x (T ) h i x (T ) T b =0whenever Gb =0. (18) Boh facors on he lef side of (18) belong o M d (x); by Theorem 1, he expecaion can be evaluaed as if (4) were rue. Doing so gives he condiion r T a T R b =0whenever Gb =0, (19) where r is a column vecor and R asymmericn n marix formed from he GACV of x (): r = s x ( 1 ) s x ( n ) T, (20) R =[s x ( i j ):i, j =1,...,n]. (21) In urn, condiion (19) says ha he row vecor r T a T R is orhogonal o all vecors b ha are orhogonal o he rows of G. Therefore, r T a T R belongs o he row space of G, ha is, r Ra = G T θ for some d-vecor θ = T. θ 0 θ d 1 We now have he sysem of equaions Ga = g, Ra + G T θ = r, (22) which consiue n + d equaions in he n + d unknowns a, θ. They generalize he Yule-Walker equaions Ra = r (23) for saionary x (); in ha case, R is a genuine covariance marix, and ˆx ( ) is he orhogonal projecion of x ( ) on he unresriced subspace generaed by x (), T.Ford 1, R is symmeric if s x () is even, bu usually has boh posiive and negaive eigenvalues. Afer solving (22) for a and θ, we can calculae he mean square error of ˆx ( ) as follows: MSE = E [x ( ) ˆx ( )] 2 =E x ( ) a T x (T ) h i x ( ) x (T ) T a = s x (0) r T a a T (r Ra) 379
4 by Theorem 1. Then, by (22), MSE = s x (0) r T a g T θ. (24) We menion hree mehods for solving (22). Firs, if n = d, hen G is nonsingular, and we can solve (22) for a andhenfor θ. Second, if R is nonsingular 3, hen we can wrie a soluion in he following form (R 1 mehod): Λ = GR 1 G T, θ = Λ 1 GR 1 r g, a = R 1 r G T θ. Third, we can se up he sysem (22) as he single marix equaion R G T a r =, (25) G 0 θ g and ell Malab 4 o solve i in one operaion ( brue force mehod). Alhough Malab ofen says ha he big marix is badly scaled or nearly singular, he soluion seems o work anyway in he cases he auhor has ried. Before carrying ou his soluion, we should reduce he rend coefficien o zero as explained earlier, eiher subracing he rend from x () if c d is known, or increasing d by1if c d is unknown. Then ˆx ( ) is unbiased for x ( ) because x ( ) ˆx ( ) M d+1 (x), all of whose members have zero expecaion. The penaly for increasing d is a greaer MSE for he BLIE, because he se of LIEs shrinks as he order increases. A. Examples 1. Model: whie FM, h 0 =1, d =1, average frequency c 1 =0, T = {0, 1,..., 10}, =5.Wegeheexpeced resul: ˆx (5) = x (0), MSE = 1 2 h 0 = The same example, bu wih an unknown frequency (phase slope) added. We calculae he BLIE of order 2 o ge a predicor ha is invarian o he added frequency. The resul is he linear exrapolaor ˆx (5) = 3 2 x (0) 1 2x ( 10), wih MSE = This is a simple demonsraion of he penaly paid for a lack of knowledge of he rend. Bu, as T reaches farher and farher ino he pas, he esimaor recovers he unknown frequency, and he MSE ends o 1 2 h Model: whie FM (h 0 =2) + random walk FM (h 2 = ), d =2. Figure 2 shows normalized rms predicion error, MSE/τ, ofx ( + τ) for hree predicors: a) linear exrapolaion from x (),x( τ), whose rms error is 2τ imes Allan deviaion for τ; b) opimal predicion based on he enire discree pas x ( 10n), n =0, 1,..., calculaed in closed form by solving for he saionary covariance marix of akalmanfiler [7]; c) opimal predicion based on x ( 10n), n =0,...,5, calculaed by he mehod given here. 4. Model: flicker FM, h 1 =1, d =2.For =8,he squares in figure 3 show he BLIE regression coefficiens for he 33-poin predicion se T = { 32, 31,..., 0}. Because 3 The auhor does no know any condiions for R o be nonsingular. 4 Copyrigh by The MahWorks, Inc. rms predicion error / τ σ y (τ) pas = pas = predicion inerval, τ Fig. 2. Normalized rms error of predicors of a sum of whie FM and random walk FM. he coefficiens are small off he se { 32, 31, 1, 0}, we also ry he laer se for predicion (filled circles). Figure 4 shows MSE/τ for boh T ses and for = τ = 1, 2, 4,...,256. The lower horizonal line is for opimal predicion from he coninuous pas, 0. The upper horizonal line is for simple linear exrapolaion from T = { τ,0} as before. We observe only a small error penaly for using he 4- poin se insead of he 33-poin se. As τ becomes large, he linear exrapolaor becomes beer han he oher predicors, which see only 32 unis ino he pas of his long-memory process rms err = rms err = predicion inerval = Fig. 3. Regression coefficiens for predicing flicker FM IV. TREND ESTIMATION Le x () have saionary dh incremens (d 0), wih an unknown rend coefficien c d.then d x () =c d d! + x 0 (), (26) where x 0 () has rend coefficien zero. Again here will be invariance and orhogonaliy condiions o deermine he 380
5 rms predicion error /τ σ y (τ) (,0] { 32,...,0} { 32, 31, 1,0} polynomial deermined by x ( 1 ),...,x( d+1 ).Moreover,he difference of any wo LIEs is in he subspace M d+1 (x, T ), he se of P b () x () such ha b () saisfies he momen condiion of order d +1. As wih he predicion problem, a LIE ĉ d is he BLIE if i saisfies he orhogonaliy condiion By (28), (29), and (30), c d ĉ d M d+1 (x, T ). (33) c d ĉ d = a () x 0 (), (34) Fig. 4. predicion inerval,τ Normalized rms errors of predicors of flicker FM opimal esimaor of c d having he form ĉ d = a () x (). (27) The invariance condiion is deermined by he requiremen ha he esimaion error c d ĉ d be unchanged if a polynomial of degree d is added o x (). Thus he error is invarian o he rend c d d /d! iself as well as o polynomials p () of degree d 1 added o x (). Afer replacing (26) by x () = c d d /d!+x 0 ()+p(), we obain " # c d ĉ d = c d 1 1 a () d d! a ()[x 0 ()+p()]. (28) If c d ĉ d is invarian o c d and p (), hen a () j =0, j =0,...,d 1, (29) a () d = d!, (30) and he esimaor gives he righ answer c d if x () is a polynomial of degree d. A random variable of form (27) is said o be a LIE of c d if (29) and (30) hold. We wrie his in he marix form Ga = g again, bu wih differen definiions of G and g: n G =.., (31) d 1 d n g = 0 0 d! T. (32) Assume ha T has a leas d +1 poins. Then a LIE of c d exiss, namely, he coefficien of d /d! in he inerpolaing which belongs o M d (x 0 ) by (29). Also, M d+1 (x, T ) = M d+1 (x 0,T) M d (x 0 ), because every member of M d+1 (x) has coefficiens ha kill he rend. Therefore, (33) can be wrien as E a T x 0 (T ) h i x 0 (T ) T b =0whenever Gb =0, (35) where x 0 (T ) is a column vecor like x (T ). Because boh facors in he lef side of (35) belong o M d (x 0 ),wemayuse Theorem 1 for x 0 () o evaluae his expression, giving he condiion a T Rb =0whenever Gb =0, (36) where R is defined by (21) wih s x () replaced by s x0 (). By he same argumen as before, Ra = G T θ for some (d +1)-vecor θ = T.Wearriveahesysem θ 0 θ d of equaions Ga = g, Ra + G T θ =0, (37) which is like (22), excep ha now here are n+d+1 equaions and unknowns, he definiions of G and g are differen, and r has become 0. The same soluion mehods are available. By (33), ĉ d is unbiased for c d because all he members of M d+1 (x) have mean zero. By (34), Theorem 1, and (37), we can calculae he MSE of ĉ d by MSE = E (c d ĉ d ) 2 =E a T x 0 (T ) h i x 0 (T ) T a = a T Ra = g T θ = d!θ d. (38) A. Examples 1. Model: x () =whie FM + c 1 (unknown frequency), h 0 = 1, d = 1, T = {0, 1,...,10}. Resul: ĉ 1 = 1 10 [x (10) x (0)], MSE= Model: x () = x 0 () c 2 2 where x 0 () is whie FM, flicker FM, or random walk FM, T = {0, 1,..., 10}. The regression coefficiens for ĉ 2 are shown in Fig. 5. Even hough whie FM is in SI(1), we have o rea i as a member of SI(2) o exrac a quadraic rend, independenly of any linear rend ha may also be presen. For whie FM, ĉ 2 is also he slope of he leas-squares linear fi o he frequency daa. 381
6 Whie FM Flicker FM Rand. walk FM Fig. 5. Regression coefficiens for esimaing quadraic drif rae in he presence of whie FM, flicker FM, or random walk FM. V. ADDITIONAL RESULTS A. Predicion from equally spaced daa For fixed d, iakeso n 3 operaions o solve he Yule- Walker equaions (23) (when d =0) or heir generalizaion (22) by general linear equaion solving mehods, where n is he number of elemens in T.IfT is an equally spaced se of imes, however, hen R is a Toepliz marix. For saionary x (), he Levinson-Durbin algorihm [15], which is a loop on n, calculaes he regression coefficiens and he MSE in O n 2 operaions. The auhor has been able o exend his algorihm o he case d 1 while keeping he O n 2 propery for fixed d. In cases ha have been ried, he resuls of he wo general algorihms (R 1 andbrueforce)andheexendedlevinson algorihm agree wihin roundoff error. [3] P. Tavella and M. Goa, Uncerainy and predicion of clock errors in space and ground applicaions, 14h European Frequency and Time Forum, pp 77 81, Torino, [4] L.-G. Bernier, Use of he Allan deviaion and linear predicion for he deerminaion of he uncerainy on ime calibraions agains prediced imescales, IEEE Trans. Insrum. Meas., vol 52, no 2, pp , April [5] J.-S. Boulanger and R. J. Douglas, Traceabiliy of frequency conrol sysems, 1997 IEEE Inernaional Frequency Conrol Symposium, pp , Orlando. [6] F. Vernoe, J. Delpore, M. Brune, and T. Tournier, Uncerainies of drif coefficiens and exrapolaion errors: applicaion o clock error predicion, Merologia, vol 38, pp , [7] J.A.Davis,C.A.Greenhall,andR.Boudjemaa, Thedevelopmenof a Kalman filer clock predicor, 19h European Frequency and Time Forum, Besançon, [8] M. A. Weiss and C. Hackman, Confidence on he hree-poin esimaor of frequency drif, 24h Annual Precise Time and Time Inerval (PTTI) Applicaions and Planning Meeing, pp , [9] V. A. Logachev and G. P. Pashev, Esimaion of linear frequency drif coefficien of frequency sandards, 1996 IEEE Inernaional Frequency Conrol Symposium, pp [10] G. Wei, Esimaions of frequency and is drif rae, IEEE Trans. Insrum. Meas., vol 46, no 1, pp 79 82, Feb [11] C. A. Greenhall, A frequency-drif esimaor and is removal from modified Allan variance, 1997 IEEE Inernaional Frequency Conrol Symposium, pp [12] C. A. Greenhall, A srucure funcion represenaion heorem wih applicaions o frequency sabiliy esimaion, IEEE Trans. Insrum. Meas., vol IM-32, no 2, pp , June [13] C. A. Greenhall, The generalized auocovariance: a ool for clock noise saisics, The Tracking and Mission Operaions Progress Repor , pp 1 32, Je Propulsion Laboraory, May 1999, hp://mo.jpl.nasa.gov/mo/progress_repor/42-137/137j.pdf [14] F. Vernoe, Applicaion of he momen condiion o noise simulaion and o sabiliy analysis, IEEE Trans. UFFC, vol 49, no 4, pp , April [15] P. J. Brockwell and R. A. Davis, Time Series: Theory and Mehods, Second Ediion, 5.2, Springer, B. Trend from symmeric daa Suppose ha T is symmeric abou some poin, which we assume o be zero; hen T = T. I can be shown ha he opimal rend coefficien esimaor ĉ d has coefficiens ha are even (a ( ) =a ()) ifd is even, and odd (a ( ) = a ()) if d is odd. In eiher case one can se up equaions like (37) for a (), 0, and a smaller auxiliary vecor θ. Thus, he dimension n + d +1 of he sysem (37) can be reduced by approximaely a facor of wo. This work was carried ou a he Je Propulsion Laboraory, California Insiue of Technology, under a conrac wih he Naional Aeronauics and Space Adminisraion. REFERENCES [1] A. M. Yaglom, Correlaion heory of processes wih random saionary nh incremens, American Mahemaical Sociey Translaions, series 2, vol 8, pp , [2] J. A. Barnes and D. W. Allan, An approach o he predicion of coordinaed universal ime, Frequency, pp 1 8, Nov Dec
Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationLecture Notes 2. The Hilbert Space Approach to Time Series
Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationReferences are appeared in the last slide. Last update: (1393/08/19)
SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be
More informationNature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.
Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike
More informationPhysics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution
Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his
More information4.6 One Dimensional Kinematics and Integration
4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of
More informationAn recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes
WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More information12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =
1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of
More informationSOLUTIONS TO ECE 3084
SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no
More informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
More informationUnit Root Time Series. Univariate random walk
Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he
More informationInstitute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler
MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationDYNAMIC ECONOMETRIC MODELS vol NICHOLAS COPERNICUS UNIVERSITY - TORUŃ Józef Stawicki and Joanna Górka Nicholas Copernicus University
DYNAMIC ECONOMETRIC MODELS vol.. - NICHOLAS COPERNICUS UNIVERSITY - TORUŃ 996 Józef Sawicki and Joanna Górka Nicholas Copernicus Universiy ARMA represenaion for a sum of auoregressive processes In he ime
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationMathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation
Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationMATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.
MATH 3B: MIDTERM REVIEW JOE HUGHES. Inegraion by Pars. Evaluae 3 e. Soluion: Firs make he subsiuion u =. Then =, hence 3 e = e = ue u Now inegrae by pars o ge ue u = ue u e u + C and subsiue he definiion
More informationWavelet Methods for Time Series Analysis. What is a Wavelet? Part I: Introduction to Wavelets and Wavelet Transforms. sines & cosines are big waves
Wavele Mehods for Time Series Analysis Par I: Inroducion o Waveles and Wavele Transforms waveles are analysis ools for ime series and images as a subjec, waveles are relaively new (983 o presen) a synhesis
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
More informationACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.
ACE 56 Fall 5 Lecure 8: The Simple Linear Regression Model: R, Reporing he Resuls and Predicion by Professor Sco H. Irwin Required Readings: Griffihs, Hill and Judge. "Explaining Variaion in he Dependen
More informationInventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions
Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.
More informationThe Arcsine Distribution
The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we
More informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
More informationIntroduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.
Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationThe Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite
American Journal of Operaions Research, 08, 8, 8-9 hp://wwwscirporg/journal/ajor ISSN Online: 60-8849 ISSN Prin: 60-8830 The Opimal Sopping Time for Selling an Asse When I Is Uncerain Wheher he Price Process
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN
Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationA Bayesian Approach to Spectral Analysis
Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2
More informationMath 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.
1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be
More informationMath 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:
Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial
More informationCHAPTER 2: Mathematics for Microeconomics
CHAPTER : Mahemaics for Microeconomics The problems in his chaper are primarily mahemaical. They are inended o give sudens some pracice wih he conceps inroduced in Chaper, bu he problems in hemselves offer
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More informationBook Corrections for Optimal Estimation of Dynamic Systems, 2 nd Edition
Boo Correcions for Opimal Esimaion of Dynamic Sysems, nd Ediion John L. Crassidis and John L. Junins November 17, 017 Chaper 1 This documen provides correcions for he boo: Crassidis, J.L., and Junins,
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationEmpirical Process Theory
Empirical Process heory 4.384 ime Series Analysis, Fall 27 Reciaion by Paul Schrimpf Supplemenary o lecures given by Anna Mikusheva Ocober 7, 28 Reciaion 7 Empirical Process heory Le x be a real-valued
More informationTwo Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017
Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
More informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationA generalization of the Burg s algorithm to periodically correlated time series
A generalizaion of he Burg s algorihm o periodically correlaed ime series Georgi N. Boshnakov Insiue of Mahemaics, Bulgarian Academy of Sciences ABSTRACT In his paper periodically correlaed processes are
More informationOnline Appendix to Solution Methods for Models with Rare Disasters
Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,
More informationThis document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC
his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationProperties of Autocorrelated Processes Economics 30331
Properies of Auocorrelaed Processes Economics 3033 Bill Evans Fall 05 Suppose we have ime series daa series labeled as where =,,3, T (he final period) Some examples are he dail closing price of he S&500,
More informationEE363 homework 1 solutions
EE363 Prof. S. Boyd EE363 homework 1 soluions 1. LQR for a riple accumulaor. We consider he sysem x +1 = Ax + Bu, y = Cx, wih 1 1 A = 1 1, B =, C = [ 1 ]. 1 1 This sysem has ransfer funcion H(z) = (z 1)
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationRANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY
ECO 504 Spring 2006 Chris Sims RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY 1. INTRODUCTION Lagrange muliplier mehods are sandard fare in elemenary calculus courses, and hey play a cenral role in economic
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More informationES.1803 Topic 22 Notes Jeremy Orloff
ES.83 Topic Noes Jeremy Orloff Fourier series inroducion: coninued. Goals. Be able o compue he Fourier coefficiens of even or odd periodic funcion using he simplified formulas.. Be able o wrie and graph
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More information= ( ) ) or a system of differential equations with continuous parametrization (T = R
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationOperators related to the Jacobi setting, for all admissible parameter values
Operaors relaed o he Jacobi seing, for all admissible parameer values Peer Sjögren Universiy of Gohenburg Join work wih A. Nowak and T. Szarek Alba, June 2013 () 1 / 18 Le Pn α,β be he classical Jacobi
More informationnon -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.
LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationPade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol
Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne
More information( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationCash Flow Valuation Mode Lin Discrete Time
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics
More informationMATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence
MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More informationVectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1
Vecorauoregressive Model and Coinegraion Analysis Par V Time Series Analysis Dr. Sevap Kesel 1 Vecorauoregression Vecor auoregression (VAR) is an economeric model used o capure he evoluion and he inerdependencies
More informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More informationADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91
ddiional Problems 9 n inverse relaionship exiss beween he ime-domain and freuency-domain descripions of a signal. Whenever an operaion is performed on he waveform of a signal in he ime domain, a corresponding
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More information6.003 Homework #8 Solutions
6.003 Homework #8 Soluions Problems. Fourier Series Deermine he Fourier series coefficiens a k for x () shown below. x ()= x ( + 0) 0 a 0 = 0 a k = e /0 sin(/0) for k 0 a k = π x()e k d = 0 0 π e 0 k d
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationProperties Of Solutions To A Generalized Liénard Equation With Forcing Term
Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick
More informationOscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,
More information