Experimental identification of low order model
|
|
- Laura Conley
- 6 years ago
- Views:
Transcription
1 METTI METTI5 Tutorial T o Experimetal idetificatio of low order model Jea-Luc Battaglia, Laboratory IM, Departmet TREFLE. Abstract The goal of this tutorial is to apply the system idetificatio techique i order to obtai a accurate direct model devoted to measuremets iversio. This tutorial is closely related to Lesso L8. A simple experimet will be used i order to give the basic ideas of the optimal experimet i system idetificatio approach. It will be particularly emphasized o the choice of the excitatio sequece. Two methods will be used: the correlatio method ad the parametric method. A direct model will thus be obtaied ad it will be aalyzed i terms of reliability ad accuracy. As a coclusio, it will be poited out the advatages of this approach with respect to the classical oe based o the resolutio of the heat diffusio equatio. System descriptio. Real life cotext I high ethalpy plasma flows for aerospace applicatios but also i high power pulsed laser physics as well as i laser surface hardeig moitorig, very high heat fluxes o the order of several MW/m have to be measured (see Figure ). Oly trasiet measuremet techiques have bee developed so far. Recetly, fast trasiet heat flux measuremets have bee coducted usig a ovel calibratio approach. I priciple, these sesors ca stay a very short time, of the order of millisecods, i the harsh eviromet i order ot to reach meltig temperature. The trasiet respose of the sesor recorded durig this short time is used to estimate the heat flux.
2 Plasma torche MW/m² Pe fluxmeter Fluxmètre coique Figure : heat flux measuremet i a plasma flow usig a fluxmeter. Oe basic priciple of such a heat flux sesor is to measure the temperature, usually with a thermocouple embedded iside a appropriate material which is part of the sesor, ad to estimate the heat flux from iverse heat coductio calculatio. Cosequetly, highest reliability i terms of measuremet accuracy ad local resolutio withi the sesor is reached whe the distace betwee the temperature measuremet ad the heated surface is small (see Figure ). Solvig the iverse heat coductio problem requires a model, so-called direct model (DM) of the heat trasfer from the heated surface of the sesor to the locatio of the thermocouple iside the sesor. Model preseted i the literature geerally assume that the sesor behaves like a semi-ifiite wall. Figure : fluxmeter descriptio Cosiderig a particular sesor ad liear heat trasfer, i.e. costat material properties for the measuremet time, the Duhamel s theorem makes sure that the DM ca be viewed as the impulse respose. This is the trasiet temperature of the thermocouple due to a heat flux o the form of a Dirac fuctio. Assumig the sesor behaves as a semi ifiite wall, the impulse respose is aalytically expressed accordig to the thermal properties of the medium as well as the locatio of the thermocouple iside the sesor. However, i real cofiguratio, the heat flux sesors ivolve several types of materials arraged i a more or less complicated way. The kowledge of their thermophysical properties as well as the thermal cotact resistaces at the iterfaces have to be kow if oe wats to use the fiite elemet method to solve a detailed DM. Furthermore, it appears that the thermocouple rise time ca be greater tha the samplig time iterval durig the acquisitio process. This meas that the thermocouple juctio is ot at a uiform temperature for each acquisitio time. I other words,
3 there exists a temperature gradiet i the juctio ad oe must also cosider heat diffusio iside the juctio. Obviously, it is ot allowed applyig the electrothermal coversio betwee the voltage drop ad the temperature of the thermocouple sice it rests o the assumptio of a uiform temperature of the juctio. I order to overcome these problems, the basic idea rests o the system idetificatio of the heat flux sesor. It cosists i calibratig the sesor by applyig a measurable trasiet heat flux i the time domai of iterest usig a modulated laser source. Give that the idetified system has to be accurate for all the swept frequecies, it is searched a heat flux waveform whose Power Spectral Desity is comparable to that of a white oise. This calibratio thus ivolves the realizatio of a specific experimetal setup i laboratory which allows applyig the heat flux o the form of a radom or Pseudo Radom Biary Sigal (PRBS) ad measurig it precisely as well as the voltage drop at the thermocouple. This approach does ot require kowig the thermophysical properties of materials as well as the exact locatio of the thermocouple. Also the kowledge of the thermocouple iertia is circumveted sice it is take ito accout withi the calibratio. From a theoretical poit of view, the idetified system is the direct model whe solvig the iverse heat coductio problem. It meas that the idetified system ad the voltage drop measuremet at the sesor durig the use of the sesor for a give applicatio is sufficiet to estimate the heat flux. This approach does ot deped o a particular sesor geometry which the allows maufacturig also particular sesor geometries for particular applicatios. Major drawback of this approach is that the calibratio must be realized i the same coditios tha those ecoutered o the process durig the use of the heat flux sesor otherwise a liear behavior has to be assumed. I other words, oe must reproduces i the laboratory the same boudary coditios i terms of trasiet ad magitude of heat flux.. System idetificatio hardware The sesor cosists of a cylidrical copper tube, where a thermocouple is itegrated. The tip is of spherical form. The thermocouple is of type K with a juctio diameter of.8 mm. As specified by the maufacturer, the rise time for the thermocouple is about msec. However, owadays the heat flux sesor has to be used o comparable time duratio, but with. msec samplig iterval. The Figure 3 shows the schematic view of the scale experimetal setup for these calibratio measuremets. Sice the ull-poit calorimeter is used to measure very high heat flux desities up to MW/m², the laser pulse eergy has to be high i order to achieve a sufficietly resolvable sigal of the thermocouple. The laser i use is a laser diode that ca provide W at 98 m wavelegth. It is kow i system idetificatio, that the best results are achieved whe laser pulses of variable legth are geerated i order to better distiguish each of the characteristic times of the system: the diffusio i the tip, the respose of the thermocouple ad the ifluece of the iterfaces. The laser is drive usig a fuctio geerator that ca geerate a Pseudo Radom Biary Sigal (PRBS) for example. O the Figure 4, the Power Spectral Desity (PSD) of such a sigal is compared to that of a
4 Dirac fuctio ad that of the Heaviside fuctio. It is obviously foud that the best represetatio of the frequecies is give with a perfect Dirac fuctio which is very difficult to implemet experimetally, particularly i terms of reproducibility. Cocerig the two others, the best oe is the PRBS because its PSD is superior to the Heaviside fuctio PSD i the studied frequecy rage. Fast Photodiode for laser trasiet behaviour Optical device for laser beam desig IR Pyrometerfor temperature measuremeto the heated area laser diode W Camera for measurigthe heated area Figure 3: system idetificatio hardware at scale Power spectral desity 5 5 PSD of PRBS PSD of a Heaviside fuctio PSD of a Dirac fuctio 4 6 Frequecy (Hz) Figure 4: Power Spectral Desity of a Dirac fuctio, a Heaviside fuctio ad a PRBS. A very small part (less tha %) of the heat flux from the laser is recorded usig a fast photodiode,. All data are recorded usig a fast oscilloscope. Moreover, the thermocouple data is amplified with a costat gai of 5. The icidet heat flux desity for the laser is deduced from the calibratio curve of the maufacturer based o blackbody absorptio. The sesor is made off copper which has bee oxidized i a furace at 4 C durig four hours i order to achieve high surface catalycity ad to reach high absorptivity. Its emissivity has the bee measured to.7 at the laser wavelegth.
5 The oise measuremet at the thermocouple is recorded without heatig the sesor by the laser ad results (samplig time is µsec) ad the oise histogram taht shows that the oise has a Gaussia distributio are preseted i Figure 5. O the other had, the computatio of the oise auto-correlatio fuctio is represeted i Figure 6. This fuctio is close to from the secod poit ad is thus equivalet to a Dirac fuctio. I coclusio, all the assumptios cocerig oise measuremet have bee checked ad the applicatio of the least square algorithm i order to estimate the parameters of the o iteger system is fully cosistet. Measuremet oise (V) Time (sec) Number of poits Sigal (V) Figure 5: Noise measuremet ad oise distributio fuctio Number of shifts Figure 6: Autocorrelatio fuctio applied to the oise measuremet. The ucertaity o the emissivity measuremet is about 6%. Cocerig the laser heat flux measured with a fa cooled broad bad sesor, the costructor gives a calibratio certificate ad a relative precisio of about %. Fially, the radius of the beam laser is measured maually ad the ucertaity is estimated to be 5%. I coclusio, the ucertaities of the measured heat flux desity absorbed by
6 the sesor is approximately 7% which is high compared to the ucertaities ivolved by the parameter estimatio method. For this tutorial it has bee realized a experimetal setup at scale /4. Ideed, the laser diode is 5W maximum power. It is drive through Natioal Istrumet card uder Labview software. The other parts of the experimet remai idetical to those of the experimet scale..3 Heat trasfer model i the fluxmeter Thermal properties of the fi are λ for the thermal coductivity, α for the thermal diffusivity, ρ for the desity ad C p for the specific heat. Heat losses betwee the fi ad the surroudig fluid are characterized by the heat exchage coefficiet h. The sectio of the thi is deoted S ad the circumferece is deoted s. If the temperature depeds oly o the logitudial directio, (Bi = h d/λ <<.) carryig out a heat flux balace o a slice of the fi of width x, we get: dt φe φs φ p = ρ Cp S x () dt φe φ p φs with: x x+ x φ s φ e d T ( x, t) = λ S () dx d T ( x + x, t) = λ S (3) dx ( (, ) ) φ p = h s x T x t T (4) Substitutig these relatios ito () we obtai: ( + ) d T x x, t d T x, t λ S h s x( T ( x, t) T ) = dx dx d T x, t ρ Cp S x dt (5) that is, Puttig ' d T x, t d, T x t λ S x h s x ( T ( x, t) T ) = ρ C p S x (6) dx dt T x = T x T, the above expressio becomes:
7 Applyig the Laplace trasform to with '(, ) = ρ (7) d T ' x, t d T ' x, t λ S h st x t C p S dx dt T ', we obtai: k ( x p) k θ ( x, p) = (8) d θ, dx ρ C S p + h s p h s = = + (9) λ S α λ S p The solutio is: Whe x we fid A=. The: θ ( x, p) = Aexp( k x) + B exp( k x) () ( ) θ = L T, t T = B () ad thus: The heat dissipated by the semi-ifiite fi is: that is, ψ θ ( x, p) = θ exp( k x) () d θ ( x, p) ( p) = λ S (3) θ dx x= ψ ( p) ( x, p) = (4) λ S k Relatio (4) shows that thermal impedace of the fi is: Z a exp( k x) =, k λ S k p h s α λ S = + (5) ψ θ Z a Figure 7: thermal impedace of the fi.
8 Correlatio ad spectral techiques. Theoretical backgroud Based o the liearity assumptio, the measured temperature y ( t ) is related to the heat flux ( t ) from the covolutio product: = ( τ ) ( τ ) dτ + y t h t e t (6) h ( t ) is the impulse respose ad e( t ) deotes the measuremet error. Now, let us multiply the two members of this equality by ( τ ) It is the obtaied: t ad itegrate from t= to ifiity. y ( t) ( t τ ) dt = h ( t τ ) ( τ ) ( t τ ) dt dτ + ( t τ ) e( t) dt (7) Oe recogizes the cross ad auto correlatio fuctios, C y ( τ ), C e ( τ ) ad is thus writte o the followig form: By choosig the heat flux t y e C τ. Relatio (7) C ( τ ) = h ( t τ ) C ( τ ) dτ + C ( τ ) (8) as a white oise lead to: C ( τ ) = δ ( τ ) (9) Fially, if it is admitted that the oise measuremet is ot correlated to the heat flux ( Ce = ), relatio (8) is summarized to: C h τ = τ () y Oe sees that the impulse respose is equal to the cross correlatio fuctio betwee the temperature of the sesor ad the heat flux. This approach is very sesitive to oise measuremet magitude, so oe would rather use the power spectral desity (PSD) istead of the correlatio fuctio. I practice, it cosists i applyig the Fourier trasform o the cross correlatio ad auto correlatio fuctios, i. e.: ad FFT Cy FFT τ = h t τ C τ dτ = Y f Φ f = S y f ()
9 Y ( f ) ad ( f ) FFT C = FFT t d = Φ f = S f () ( τ ) ( τ ) ( τ ) τ Φ are the Fourier trasforms of the temperature ad the heat flux respectively as well as S ( f ) ad S ( f ) y are the auto ad cross PSD. The, by applyig the Fourier trasform o relatio (8) it is immediately obtaied: S f = H f S f + S f (3) y e Fially, assumig that the oise measuremet is ot correlated with the heat flux ( S ( f ) expressio of the trasfer fuctio is: H ( f ) S ( f ) ( f ) e = ), the y = (4) Sice the legth of the experimet is set to a fixed value τ, the real iput sigal is: I this relatio, the heat flux leads to: S ( t) ( t) ( t) Π = Π τ (5) Π t = whe t τ ad elsewhere. The applyig the Fourier trasform o τ ( f ) ( f ) ( π τ f ) si τ π τ f Φ Π = Φ (6) It appears that the Fourier trasform of the heat flux is covoluted by the sius cardial fuctio. Usually, the heat flux is pre widowed by a specific fuctio g ( t) the fuctio ( t) Π as: τ ( t) ( t) g ( t) For example, it is ofte used of the Haig widow defied by: τ which decreases the ifluece of Π = τ (7) g τ ( t) π t =.5 cos τ (8) It is used a improved estimatio of S ( f ) ad S ( f ) y proposed by Welch. The method cosists i dividig the time series data ito possible overlappig segmets, computig the auto ad cross power spectral desities ad averagig the estimates.. Applicatio Durig the tutorial we will geerate kids of heat flux sequece:
10 .. liear swept-frequecy cosie sigal We will cosider first a heat flux o the form of a liear swept-frequecy cosie sigal: The frequecy varies liearly with time as: ( t) ( π f ( t) t) = si (9) f f f t f t = +, t t t (3) f is the istataeous frequecy at time, ad f is the istataeous frequecy at time t (see Figure 8) Figure 8: example of liear swept-frequecy cosie sigal with f =.Hz ad f = Hz. It will be used the Welch techique. The swept-frequecy cosie heat flux waveform has two major iterestig features. The first is that the offset must be easily removed from the experimetal heat flux i order to fully satisfy the relatio (9). The secod feature is that the explored frequetial domai, defied from the sesitivity aalysis, is perfectly swept. Furthermore, the use of auto ad cross power spectral desity fuctios allows defiig the so called coherece fuctio as: C y ( f ) y y y f S = (3) S f S f This fuctio ca be viewed as the correlatio coefficiet betwee the temperature ad the heat flux ad lies betwee ad. If it is at a certai frequecy, the there is perfect correlatio betwee the two sigals at this frequecy. I other words, there is cosequetly o oise iterferig at this frequecy, what lead to: ( ) S f = S f C f (3) e y y y
11 .. PRBS sigal I a secod stage we will cosider the heat flux sequece as a Pseudo Radom Biary Sigal (PRBS). White oise is the term give to completely radom upredictable oise, such as the hiss you hear o a utued radio. It has the property of havig compoets at every frequecy. A pseudo-radom biary sequece (PRBS) ca also have this property, but is etirely predictable. A PRBS is rather like a log recurrig decimal umber- it looks radom if you examie a short piece of the sequece, but it actually repeats itself every m bits. Of course, the larger m is, the more radom it looks. You ca geerate a PRBS with a shift register ad a XOR gate. Coectig the outputs of two stages of the shift register to the XOR gate, ad the feedig the result back ito the iput of the shift register will geerate a PRBS of some sort. Some combiatios of outputs produce loger PRBSs tha others- the logest oes are called m-sequeces (where m meas maximum legth ). A biary sequece (BS) is a sequece of N bits, a j for j =,,...,N, i.e. m oes ad N m zeros. A BS is pseudo-radom (PRBS) if its autocorrelatio fuctio: has oly two values: C v N = a j a j+ v (33) j= C v m, if v mod N = m c,otherwise (34) where c = m N (35) is called the duty cycle of the PRBS Figure 9: example of a PRBS sequece (X axis is the umber of samples).
12 A PRBS is radom i a sese that the value of a a j elemet is idepedet of the values of ay of the other elemets, similar to real radom sequeces. It is 'pseudo' because it is determiistic ad after N elemets it starts to repeat itself, ulike real radom sequeces, such as sequeces geerated by radioactive decay or by white oise. The PRBS is more geeral tha the -sequece, which is a special pseudo-radom biary sequece of bits geerated as the output of a liear shift register. A -sequece always has a / duty cycle ad its umber of elemets N = k. 3 Parametric idetificatio 3. AR models, theoretical backgroud As expressed by relatio (5) heat trasfer i a fi is modelled as: θ exp( k x) ( x, p) = ψ ( p), with λ S k k p h s α λ S = + (36) We ca write that: p h s h s λ S k = + = p + (37) α λ S λ S α h s Usig the series it is foud: O the other had oe has: h s λ S h s S k = p + = p λ S α h s λ S α h s e z = ( )! λ = ( )! (38) z =, z (39)! Replacig k i the expoetial fuctio ad this former with its series, it is foud that o obtai a equivalet expressio of relatio (36) o the followig form: Where α ad + α s θm ( s) = β s Φ ( s) (4) = = β have complex but aalytical expressios. Give to the iitial coditio (temperature is zero at each poit of the domai) ad usig the property: it thus appear that relatio (4) is equivalet to: d f t k d f L = s F ( s) s (4) dt k= dt
13 ( t) + d T M, t d α = β (4) = dt = dt Usig the discrete form of the derivatives a equivalet form of relatio (4) that lead to express the temperature at time k t from the heat flux ad the temperature at previous times as: T = b + b + + b a T a T L L, k (,, N ) k k k k k b k k b k a k a I this equatio the cotiuous time t have bee discretized i N samplig itervals of legth as T N t =, where T is the fial time. Parameters ( a, a,, a, b, b, b ) a b = L (43) t, such L L are real umbers. The iteger k correspods to the lag time that depeds o the distace betwee the output ad the iput. This laggig effect oly depeds o the distace betwee the tip of the isert ad the poit M. Sice the white oise term e( k ) defied by ek = Tk Yk (44) eters as a direct error i the differece equatio, the model (43) is ofte called a equatio error model. By substitutig equatio (44) i equatio (43) oe obtais where is the regressio matrix ad is the vector of ukows. Y = H Θ + e (45) k k k [ Y Y Y ] H = L L (46) k k k k a k k k k k k b [ a a a b b ] = L L (47) Θ a b By performig N successive measuremets, relatio (45) ca be writte o a matrix form as where Y = ΗΘ + E (48) [ Y Y Y ] T Y = d d + L N (49) [ e e e ] T E = d d + L N (5)
14 Η H H M H d d + = N (5) ad Resolutio of equatio (48) i the least square sese leads to The choice of = [ a, b, k] d = max a +, b + k + (5) ˆ T T Θ = Η Η Η Y (53) Λ is made by collectig i a matrix all the values of Λ to be ivestigated ad lookig o the value of the Aikake criterio defied by + N Ψ = V, = a + b + (54) N where is the total umber of estimated parameters ad V is the loss fuctio defied by V = (55) N ek k= Stadard errors of the estimates are calculated from the covariace matrix of ˆΘ. If the assumptios of additive, zeros mea, costat variace is expressed as A estimate of the variace σ, deoted σ ad ucorrelated errors are verified, the covariace matrix cov ( ˆ T ) = σ s, is Θ H H (56) s = N E E (57) T From equatio (53), by substitutig the estimated vector ˆΘ i equatio (48), oe fids ˆ T T Thereby, asymptotic estimatio of the mea of ˆΘ is Θ = Θ + H H H E (58) ˆ T T = + k k k k E E E e Θ Θ H H H (59) From this last relatio it is clear that the estimatio is ubiased if e k is ucorrelated with E[ e k ] =. Thereby, the typical whiteess test is to determie the covariace estimate k or if
15 N ek r k r N r= Idepedece betwee e ad q is tested usig the sample covariace R = e e (6) R N e k er k r N r= = (6) It may be oted that whe usig a model as equatio (43), the least square procedure automatically makes the correlatio betwee e k ad k δ zero for δ = k, k +, L, k + b. I practice, relatios (6) ad (6) are computed usig the Fourier Trasform of the covolutio product. 3. Applicatio The method will be applied durig the tutorial startig from the respose to a PRBS sequece for the heat flux. 3.3 NI models, theoretical backgroud At very short times, heat losses are egligible face to the heat diffusio i the fi. Therefore: Sice: p h s p ad k α λ S α (6) θ ( x, p) = ( k x) λ S k =! ψ ( p), with k p h s α λ S = + (63) Replacig k i this relatio leads to: θ x p / α! / = ( x, p) = ψ ( p) λ ρ C S p p / (64) That ca be also writte as: Sice we have show i L8 course that relatio: x λ ρ C S p θ x p p ψ p p (, ) = (65) / / / = α! ν ν d f t ν k d f L s F = ( s) s (66) ν dt k = dt remais true eve if ν is a real umber or more geerally complex, we ca express relatio (65) i the time domai as:
16 α D = β (67) { T } D { } m t t = With: x α = λ ρ Cp S ad β = (68) / α! It appears thus tha model (4) will ot be accurate eough to describe the trasiet behaviour at the short times. As said i lesso 8, a optimal structure of a low order model for heat trasfer problem by diffusio must be of the followig form: D m = = { T ( t) } = D { ( t) } α β (69) 3.4 Applicatio Durig the tutorial, we will apply the NISI method to estimate parameters i relatio (69) i order to fit temperature measuremets at the sesor startig from a PRBS of the heat flux.
3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationThe 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 informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationA. 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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departmet of Electrical Egieerig ad Computer Sciece 6.34 Discrete Time Sigal Processig Fall 24 BACKGROUND EXAM September 3, 24. Full Name: Note: This exam is closed
More informationChapter 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 informationThermal characterization of an insulating material through a tri-layer transient method
Advaced Sprig School «Thermal Measuremets & Iverse techiques», Domaie de Fraço, Biarritz, March -6 05 T6 Thermal characterizatio of a isulatig material through a tri-layer trasiet method V. Félix, Y. Jaot,
More informationFrequency Response of FIR Filters
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we itroduce the idea of the frequecy respose of LTI systems, ad focus specifically o the frequecy respose of FIR filters.. Steady-state
More informationStatistical Analysis on Uncertainty for Autocorrelated Measurements and its Applications to Key Comparisons
Statistical Aalysis o Ucertaity for Autocorrelated Measuremets ad its Applicatios to Key Comparisos Nie Fa Zhag Natioal Istitute of Stadards ad Techology Gaithersburg, MD 0899, USA Outlies. Itroductio.
More information6.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 informationFIR 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 informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationOrthogonal 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 informationA 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 informationChapter 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 informationTime-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 informationThe Discrete Fourier Transform
The iscrete Fourier Trasform The discrete-time Fourier trasform (TFT) of a sequece is a cotiuous fuctio of!, ad repeats with period. I practice we usually wat to obtai the Fourier compoets usig digital
More informationLinear time invariant systems
Liear time ivariat systems Alejadro Ribeiro Dept. of Electrical ad Systems Egieerig Uiversity of Pesylvaia aribeiro@seas.upe.edu http://www.seas.upe.edu/users/~aribeiro/ February 25, 2016 Sigal ad Iformatio
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationTHE KALMAN FILTER RAUL ROJAS
THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider
More informationDiscrete Mathematics for CS Spring 2008 David Wagner Note 22
CS 70 Discrete Mathematics for CS Sprig 2008 David Wager Note 22 I.I.D. Radom Variables Estimatig the bias of a coi Questio: We wat to estimate the proportio p of Democrats i the US populatio, by takig
More informationChapter 10 Advanced Topics in Random Processes
ery Stark ad Joh W. Woods, Probability, Statistics, ad Radom Variables for Egieers, 4th ed., Pearso Educatio Ic.,. ISBN 978--3-33-6 Chapter Advaced opics i Radom Processes Sectios. Mea-Square (m.s.) Calculus
More informationDiscrete-Time Signals and Systems. Signals and Systems. Digital Signals. Discrete-Time Signals. Operations on Sequences: Basic Operations
-6.3 Digital Sigal Processig ad Filterig..8 Discrete-ime Sigals ad Systems ime-domai Represetatios of Discrete-ime Sigals ad Systems ime-domai represetatio of a discrete-time sigal as a sequece of umbers
More informationChapter 6 Sampling Distributions
Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationLecture 8 & Tutorial 1: Experimental idention of low order model
Eurother Advaced School Metti 5 Roscoff Jue 13-18, 211 Lecture 8 & utorial 1: Experietal idetio of low order odel Jea-Luc Battaglia, Laboratory I2M, Departet REFLE, I2M Bordeaux Measureet iversio 35 3
More informationTHE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS
R775 Philips Res. Repts 26,414-423, 1971' THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS by H. W. HANNEMAN Abstract Usig the law of propagatio of errors, approximated
More informationModule 18 Discrete Time Signals and Z-Transforms Objective: Introduction : Description: Discrete Time Signal representation
Module 8 Discrete Time Sigals ad Z-Trasforms Objective:To uderstad represetig discrete time sigals, apply z trasform for aalyzigdiscrete time sigals ad to uderstad the relatio to Fourier trasform Itroductio
More informationFilter 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 informationChapter 11 Output Analysis for a Single Model. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation
Chapter Output Aalysis for a Sigle Model Baks, Carso, Nelso & Nicol Discrete-Evet System Simulatio Error Estimatio If {,, } are ot statistically idepedet, the S / is a biased estimator of the true variace.
More informationStatistical Noise Models and Diagnostics
L. Yaroslavsky: Advaced Image Processig Lab: A Tutorial, EUSIPCO2 LECTURE 2 Statistical oise Models ad Diagostics 2. Statistical models of radom iterfereces: (i) Additive sigal idepedet oise model: r =
More informationMCT242: Electronic Instrumentation Lecture 2: Instrumentation Definitions
Faculty of Egieerig MCT242: Electroic Istrumetatio Lecture 2: Istrumetatio Defiitios Overview Measuremet Error Accuracy Precisio ad Mea Resolutio Mea Variace ad Stadard deviatio Fiesse Sesitivity Rage
More informationEECE 301 Signals & Systems
EECE 301 Sigals & Systems Prof. Mark Fowler Note Set #8 D-T Covolutio: The Tool for Fidig the Zero-State Respose Readig Assigmet: Sectio 2.1-2.2 of Kame ad Heck 1/14 Course Flow Diagram The arrows here
More informationProbability, Expectation Value and Uncertainty
Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5
CS434a/54a: Patter Recogitio Prof. Olga Veksler Lecture 5 Today Itroductio to parameter estimatio Two methods for parameter estimatio Maimum Likelihood Estimatio Bayesia Estimatio Itroducto Bayesia Decisio
More informationCEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering
CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio
More informationProblem Set 4 Due Oct, 12
EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios
More informationSignal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform
Sigal Processig i Mechatroics Summer semester, 1 Lecture 3, Covolutio, Fourier Series ad Fourier rasform Dr. Zhu K.P. AIS, UM 1 1. Covolutio Covolutio Descriptio of LI Systems he mai premise is that the
More informationMeasurement uncertainty of the sound absorption
Measuremet ucertaity of the soud absorptio coefficiet Aa Izewska Buildig Research Istitute, Filtrowa Str., 00-6 Warsaw, Polad a.izewska@itb.pl 6887 The stadard ISO/IEC 705:005 o the competece of testig
More informationFrequency Domain Filtering
Frequecy Domai Filterig Raga Rodrigo October 19, 2010 Outlie Cotets 1 Itroductio 1 2 Fourier Represetatio of Fiite-Duratio Sequeces: The Discrete Fourier Trasform 1 3 The 2-D Discrete Fourier Trasform
More informationSignals & Systems Chapter3
Sigals & Systems Chapter3 1.2 Discrete-Time (D-T) Sigals Electroic systems do most of the processig of a sigal usig a computer. A computer ca t directly process a C-T sigal but istead eeds a stream of
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationIntroduction 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 informationExpectation and Variance of a random variable
Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio
More informationEE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course
Sigal-EE Postal Correspodece Course 1 SAMPLE STUDY MATERIAL Electrical Egieerig EE / EEE Postal Correspodece Course GATE, IES & PSUs Sigal System Sigal-EE Postal Correspodece Course CONTENTS 1. SIGNAL
More informationADVANCED DIGITAL SIGNAL PROCESSING
ADVANCED DIGITAL SIGNAL PROCESSING PROF. S. C. CHAN (email : sccha@eee.hku.hk, Rm. CYC-702) DISCRETE-TIME SIGNALS AND SYSTEMS MULTI-DIMENSIONAL SIGNALS AND SYSTEMS RANDOM PROCESSES AND APPLICATIONS ADAPTIVE
More informationBayesian Methods: Introduction to Multi-parameter Models
Bayesia Methods: Itroductio to Multi-parameter Models Parameter: θ = ( θ, θ) Give Likelihood p(y θ) ad prior p(θ ), the posterior p proportioal to p(y θ) x p(θ ) Margial posterior ( θ, θ y) is Iterested
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationMechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter
Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,
More informationMonte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem
Australia Joural of Basic Applied Scieces, 5(): 097-05, 0 ISSN 99-878 Mote Carlo Optimizatio to Solve a Two-Dimesioal Iverse Heat Coductio Problem M Ebrahimi Departmet of Mathematics, Karaj Brach, Islamic
More informationTIME-CORRELATION FUNCTIONS
p. 8 TIME-CORRELATION FUNCTIONS Time-correlatio fuctios are a effective way of represetig the dyamics of a system. They provide a statistical descriptio of the time-evolutio of a variable for a esemble
More informationChapter 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 informationa b c d e f g h Supplementary Information
Supplemetary Iformatio a b c d e f g h Supplemetary Figure S STM images show that Dark patters are frequetly preset ad ted to accumulate. (a) mv, pa, m ; (b) mv, pa, m ; (c) mv, pa, m ; (d) mv, pa, m ;
More informationFIR Filter Design: Part I
EEL3: Discrete-Time Sigals ad Systems FIR Filter Desig: Part I. Itroductio FIR Filter Desig: Part I I this set o otes, we cotiue our exploratio o the requecy respose o FIR ilters. First, we cosider some
More information577. Estimation of surface roughness using high frequency vibrations
577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia
More informationChapter 2 The Monte Carlo Method
Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful
More informationResampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.
Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator
More informationOlli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5
Sigals ad Systems Sigals ad Systems Sigals are variables that carry iformatio Systemstake sigals as iputs ad produce sigals as outputs The course deals with the passage of sigals through systems T-6.4
More informationEECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1
EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Trasform The -trasform is a very importat tool i describig ad aalyig digital systems. It offers the techiques for digital filter desig ad frequecy aalysis of digital sigals. Defiitio of -trasform:
More informationECE 901 Lecture 12: Complexity Regularization and the Squared Loss
ECE 90 Lecture : Complexity Regularizatio ad the Squared Loss R. Nowak 5/7/009 I the previous lectures we made use of the Cheroff/Hoeffdig bouds for our aalysis of classifier errors. Hoeffdig s iequality
More informationTEACHER CERTIFICATION STUDY GUIDE
COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationWe 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 informationProperties and Hypothesis Testing
Chapter 3 Properties ad Hypothesis Testig 3.1 Types of data The regressio techiques developed i previous chapters ca be applied to three differet kids of data. 1. Cross-sectioal data. 2. Time series data.
More informationSolution of EECS 315 Final Examination F09
Solutio of EECS 315 Fial Examiatio F9 1. Fid the umerical value of δ ( t + 4ramp( tdt. δ ( t + 4ramp( tdt. Fid the umerical sigal eergy of x E x = x[ ] = δ 3 = 11 = ( = ramp( ( 4 = ramp( 8 = 8 [ ] = (
More informationChapter 2 Systems and Signals
Chapter 2 Systems ad Sigals 1 Itroductio Discrete-Time Sigals: Sequeces Discrete-Time Systems Properties of Liear Time-Ivariat Systems Liear Costat-Coefficiet Differece Equatios Frequecy-Domai Represetatio
More informationUNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014
UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 6C Problem Set 4 Bejami Stahl November 6, 4 BOAS, P. 63, PROBLEM.-5 The Laguerre differetial equatio, x y + ( xy + py =, will be solved
More informationSince X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain
Assigmet 9 Exercise 5.5 Let X biomial, p, where p 0, 1 is ukow. Obtai cofidece itervals for p i two differet ways: a Sice X / p d N0, p1 p], the variace of the limitig distributio depeds oly o p. Use the
More informationInfinite 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 informationSignal Processing in Mechatronics
Sigal Processig i Mechatroics Zhu K.P. AIS, UM. Lecture, Brief itroductio to Sigals ad Systems, Review of Liear Algebra ad Sigal Processig Related Mathematics . Brief Itroductio to Sigals What is sigal
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationBasics of Probability Theory (for Theory of Computation courses)
Basics of Probability Theory (for Theory of Computatio courses) Oded Goldreich Departmet of Computer Sciece Weizma Istitute of Sciece Rehovot, Israel. oded.goldreich@weizma.ac.il November 24, 2008 Preface.
More information1 Inferential Methods for Correlation and Regression Analysis
1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet
More information5. Fractional Hot deck Imputation
5. Fractioal Hot deck Imputatio Itroductio Suppose that we are iterested i estimatig θ EY or eve θ 2 P ry < c where y fy x where x is always observed ad y is subject to missigess. Assume MAR i the sese
More informationAnalysis of Experimental Data
Aalysis of Experimetal Data 6544597.0479 ± 0.000005 g Quatitative Ucertaity Accuracy vs. Precisio Whe we make a measuremet i the laboratory, we eed to kow how good it is. We wat our measuremets to be both
More informationDiscrete-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 informationLinear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy
Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y
More informationChapter 8. DFT : The Discrete Fourier Transform
Chapter 8 DFT : The Discrete Fourier Trasform Roots of Uity Defiitio: A th root of uity is a complex umber x such that x The th roots of uity are: ω, ω,, ω - where ω e π /. Proof: (ω ) (e π / ) (e π )
More informationQuantile regression with multilayer perceptrons.
Quatile regressio with multilayer perceptros. S.-F. Dimby ad J. Rykiewicz Uiversite Paris 1 - SAMM 90 Rue de Tolbiac, 75013 Paris - Frace Abstract. We cosider oliear quatile regressio ivolvig multilayer
More informationKinetics 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 informationGeneralizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations
Geeraliig the DTFT The Trasform M. J. Roberts - All Rights Reserved. Edited by Dr. Robert Akl 1 The forward DTFT is defied by X e jω = x e jω i which = Ω is discrete-time radia frequecy, a real variable.
More informationThe standard deviation of the mean
Physics 6C Fall 20 The stadard deviatio of the mea These otes provide some clarificatio o the distictio betwee the stadard deviatio ad the stadard deviatio of the mea.. The sample mea ad variace Cosider
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationStatisticians use the word population to refer the total number of (potential) observations under consideration
6 Samplig Distributios Statisticias use the word populatio to refer the total umber of (potetial) observatios uder cosideratio The populatio is just the set of all possible outcomes i our sample space
More informationProblem Set 2 Solutions
CS271 Radomess & Computatio, Sprig 2018 Problem Set 2 Solutios Poit totals are i the margi; the maximum total umber of poits was 52. 1. Probabilistic method for domiatig sets 6pts Pick a radom subset S
More informationStatistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.
Statistical Iferece (Chapter 10) Statistical iferece = lear about a populatio based o the iformatio provided by a sample. Populatio: The set of all values of a radom variable X of iterest. Characterized
More informationPractical Spectral Anaysis (continue) (from Boaz Porat s book) Frequency Measurement
Practical Spectral Aaysis (cotiue) (from Boaz Porat s book) Frequecy Measuremet Oe of the most importat applicatios of the DFT is the measuremet of frequecies of periodic sigals (eg., siusoidal sigals),
More informationModule 11: Applications : Linear prediction, Speech Analysis and Speech Enhancement Prof. Eliathamby Ambikairajah Dr. Tharmarajah Thiruvaran School
Module : Applicatios : Liear predictio, Speech Aalysis ad Speech Ehacemet Prof. Eliathamby Ambiairajah Dr. Tharmarajah Thiruvara School of Electrical Egieerig & Telecommuicatios The Uiversity of New South
More informationThe Discrete Fourier Transform
The Discrete Fourier Trasform Complex Fourier Series Represetatio Recall that a Fourier series has the form a 0 + a k cos(kt) + k=1 b k si(kt) This represetatio seems a bit awkward, sice it ivolves two
More informationCS284A: Representations and Algorithms in Molecular Biology
CS284A: Represetatios ad Algorithms i Molecular Biology Scribe Notes o Lectures 3 & 4: Motif Discovery via Eumeratio & Motif Represetatio Usig Positio Weight Matrix Joshua Gervi Based o presetatios by
More informationSECTION 2 Electrostatics
SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial
More informationRandom 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 information11 Correlation and Regression
11 Correlatio Regressio 11.1 Multivariate Data Ofte we look at data where several variables are recorded for the same idividuals or samplig uits. For example, at a coastal weather statio, we might record
More informationECE-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 informationQuantum Computing Lecture 7. Quantum Factoring
Quatum Computig Lecture 7 Quatum Factorig Maris Ozols Quatum factorig A polyomial time quatum algorithm for factorig umbers was published by Peter Shor i 1994. Polyomial time meas that the umber of gates
More informationRandom Walks on Discrete and Continuous Circles. by Jeffrey S. Rosenthal School of Mathematics, University of Minnesota, Minneapolis, MN, U.S.A.
Radom Walks o Discrete ad Cotiuous Circles by Jeffrey S. Rosethal School of Mathematics, Uiversity of Miesota, Mieapolis, MN, U.S.A. 55455 (Appeared i Joural of Applied Probability 30 (1993), 780 789.)
More informationCh3 Discrete Time Fourier Transform
Ch3 Discrete Time Fourier Trasform 3. Show that the DTFT of [] is give by ( k). e k 3. Determie the DTFT of the two sided sigal y [ ],. 3.3 Determie the DTFT of the causal sequece x[ ] A cos( 0 ) [ ],
More informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 016 MODULE : Statistical Iferece Time allowed: Three hours Cadidates should aswer FIVE questios. All questios carry equal marks. The umber
More informationSample Size Estimation in the Proportional Hazards Model for K-sample or Regression Settings Scott S. Emerson, M.D., Ph.D.
ample ie Estimatio i the Proportioal Haards Model for K-sample or Regressio ettigs cott. Emerso, M.D., Ph.D. ample ie Formula for a Normally Distributed tatistic uppose a statistic is kow to be ormally
More information