From deterministic regular waves to a random field. From a determinstic regular wave to a deterministic irregular solution

Size: px
Start display at page:

Download "From deterministic regular waves to a random field. From a determinstic regular wave to a deterministic irregular solution"

Transcription

1 Classiicatio: Iteral Status: Drat z ξ(x,y, w& x w u u& h Particle ositio From determiistic regular waves to a radom ield Sverre Haver, StatoilHydro, Jauary 8 From a determistic regular wave to a determiistic irregular solutio Solutio o liearized goverig equatios (i.e. boudary equatios are alied at the mea ree surace ad terms o order higher tha are eglected: ω ξ ( x, φ si( ωt kx ξ si( ωt kx g Real ocea waves excet swell do ot look like this. They look less regular! Let us assume that the sea surace reeat ater T ad use Fourier. The we ca write the observed surace as sum o siusoidals. ξ ( π acos t + b T π si t T Each comoet is a solutio to the liearized equatios, i.e. the sum is also a solutio to the liearized equatios. we have a solutio that looks more reasoable.

2 A radom solutio By itroducig: b ξ a + b ad θ arcta a The sum ca be writte: ξ ( ξ cos( ωt θ, requecy resolutio (ω i+ ω i is: This solutio is a determistic solutio. How could we establish a more geeral solutio? Let us itroduce the hase as a radom variable beig uiormly distributed betwee ad π, while keeig our origial values or amlitude ξ, i.e.: ( ξ cos( ω t I order to obtai oe realizatio, o (, we geerate dieret hases, θ,,. By geeratig a ew set o hases we will obtai a ew realizatio o the radom rocess, ( a.s.o.. Θ ote: Frequecy resolutio is deied by legth o time widow, T. 3 Δω π T We have a solutio rom which we ca geerate a iiite amout o ossible realizatios. The observed time history we started with is oe o the may ossible realizatios. 4 Probabilistic structure o our solutio I Solutio: ( ξ cos( ωt Θ What is the distributio o ( t? Observatio: t is a sum o may ideedet radom comoets. I o o the comoets domiates over the other, we will accordig to the cetral limit theorem have: ( ξ; ξ ex π ( ( t is Gaussia with zero mea ad variace equal to the sum o the comoet variaces. The variace o a harmoic wave is give by: The variace o t thus becomes: ( ξ Observatio: As ar as we do ot chage the comoet amlitudes durig the time widow, the rocess variace will be costat.

3 5 Descritio o a Gaussia surace rocess The oit distributio o (t i ad (t is a two-dimesioal Gaussia: i ( ξi, ξ ex π ρi ξ i ρi ( ξ ξ ξ i ρ + i Where: (, Cov i ρ i E[ ( t ( t ] r ( τ i, τ t i -t r (τ ideedet o absolute time, i rocess is statioary. A Gaussia stochastic rocess is i ricile described by a oit distributio o its values at a large umber o times, t, t,., t. The arameters o this oit model ca be determied rom r (τ. A statioary Gaussia rocess is i a statistical sese comletely characterized by the autocorrelatio uctio r (τ. Questio: How ca the variace be determied rom the auto-correlatio uctio? 6 Wave sectrum I ractise it is more commo to work with the Fourier trasorm o auto-correlatio uctio which is the wave sectrum. The oe-sided wave sectrum reads: iωτ s ( ω e r ( τ dτ ; ω π r i ( τ r ( τ e ωτ s ( ω dω; τ We kow rom questio o revious slide that: r ( From exressio or the auto-correlatio above, we see that: s ( ω dω Imortat result!! I we have exaded the rocess o legth T i comoet, the variace o the th comoet reads: s ( ω Δω where Δω π T 3

4 7 Direct estimatio o sectrum rom a wave history We kow rom a revious slide that the variace o the th comoet also is give by: ξ Settig the two exressios or the variace equal, ad solve with resect to the wave sectrum, gives: ξ s ( ω Δω The Fourier amlitides, ξ, are tyically calculated rom a time history o ξ(, usig a Fast Fourier Trasorm techique. Simulatig a realizatio o X( corresodig to a kow wave sectrum 8 The radom rocess is give by: ( ξ cos( ωt Θ From revious slide we have: ξ s ( ω Δω ( s ( ω Δω cos( ω t Θ I we ow geerate m sets o radom hases, we get m realizatios o the rocess. They will all look dieret, but wave wave heights ad wave eriods will o course show some similarity. I we calculate the wave sectrum rom these simulated time histories, we will id that they all roduce (as exected the iut sectrum, i.e. a smooth sectrum i iut sectrum is smooth. I we estimate the wave sectrum rom a observed time history o legth T, we will id that the estimated sectrum is o very irregular shae. I our simulatios shall relect the same behaviour, we will also have to itroduce the Fourier amlitude,, as a radom variable. I that case it should be itroduced as a Rayleigh variable. 4

5 9 Stadard models or wave sectrum Fully develoed wid sea: Pierso-Moskowitz (PM s ( { 4. t 4 } 4.35 h t 5 ex 5 s J Growig wid sea: JOSWAP sectrum (J PM s (.35h s t 5 ex.5 4 (.87 lγ γ ex.5 T Pure swell sea: JOSWAP sectrum with γ e.g. -5 Combied sea (Wid sea lus swell: Torsethauge sectrum (T Torsethauge sectrum suerosistio o two Joswa tye sectra Wid Sea I sectrum give i ad oe would reer to have i w, oe ca id the latter by: ω s ( ω s ( π π Swell Sea Torsethauge versus JOSWAP Hs.m, T.s Hs.m, T s 3 5 Sectral Desity (m***s 5 Torsethauge (4 JOSWAP Frequecy (Hz Sectral Desity (m***s Torsethauge(4 JOSWAP Frequecy (Hz Hs.m, T 4.s Sectral Desity (m***s Torsethauge (4 JOSWAP Frequecy (Hz 5

6 Rage o validity or sectral models (Figure is meat to give a qualitative idicatio. Oe should kow which sea states that are o maor imortace whe selectig sectral model. -year cotour or h s ad t Secod order surace rocess The Gaussia surace model is a aroximatio. It is i agreemet with liear wave theory. A secod order model is give by: ξ ( t ξ ( t ( t G + ξ Comariso Gaussia rocess ad secod order sum rocess. I ricile there will also be a dierece requecy comoet. Gaussia rocess Secod order correctio ξ ( Re c m m c + [ h i( ωm + ω t e + h i( ω t m ω e ] m m I roblem uder cosideratio is sesitive to the detail surace elevatio, a o-gaussia model should be adoted. Surace elevatio is straight orward, the corresodig kiematics are more challegig. 6

7 3 Crest heigts o Gaussia ad Scod order rocesses Cumulative robability Rayleigh.3 Jahs ad Wheeler. Forristall. Witerstei Crest height (m Distributio uctio or the 3-hour maximum crest height or a sea state with h s 8.m ad t 7s (**(-4 aual robability sea states. (Forristall Secod order 7

Confidence intervals for proportions

Confidence intervals for proportions Cofidece itervals for roortios Studet Activity 7 8 9 0 2 TI-Nsire Ivestigatio Studet 60 mi Itroductio From revious activity This activity assumes kowledge of the material covered i the activity Distributio

More information

Distribution of Sample Proportions

Distribution of Sample Proportions Distributio of Samle Proortios Probability ad statistics Aswers & Teacher Notes TI-Nsire Ivestigatio Studet 90 mi 7 8 9 10 11 12 Itroductio From revious activity: This activity assumes kowledge of the

More information

Lecture 2: Monte Carlo Simulation

Lecture 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 information

The Discrete-Time Fourier Transform (DTFT)

The Discrete-Time Fourier Transform (DTFT) EEL: Discrete-Time Sigals ad Systems The Discrete-Time Fourier Trasorm (DTFT) The Discrete-Time Fourier Trasorm (DTFT). Itroductio I these otes, we itroduce the discrete-time Fourier trasorm (DTFT) ad

More information

ECE534, Spring 2018: Final Exam

ECE534, Spring 2018: Final Exam ECE534, Srig 2018: Fial Exam Problem 1 Let X N (0, 1) ad Y N (0, 1) be ideedet radom variables. variables V = X + Y ad W = X 2Y. Defie the radom (a) Are V, W joitly Gaussia? Justify your aswer. (b) Comute

More information

Chapter 12 EM algorithms The Expectation-Maximization (EM) algorithm is a maximum likelihood method for models that have hidden variables eg. Gaussian

Chapter 12 EM algorithms The Expectation-Maximization (EM) algorithm is a maximum likelihood method for models that have hidden variables eg. Gaussian Chapter 2 EM algorithms The Expectatio-Maximizatio (EM) algorithm is a maximum likelihood method for models that have hidde variables eg. Gaussia Mixture Models (GMMs), Liear Dyamic Systems (LDSs) ad Hidde

More information

Confidence Intervals

Confidence Intervals Cofidece Itervals Berli Che Deartmet of Comuter Sciece & Iformatio Egieerig Natioal Taiwa Normal Uiversity Referece: 1. W. Navidi. Statistics for Egieerig ad Scietists. Chater 5 & Teachig Material Itroductio

More information

Simulation. Two Rule For Inverting A Distribution Function

Simulation. Two Rule For Inverting A Distribution Function Simulatio Two Rule For Ivertig A Distributio Fuctio Rule 1. If F(x) = u is costat o a iterval [x 1, x 2 ), the the uiform value u is mapped oto x 2 through the iversio process. Rule 2. If there is a jump

More information

NUMERICAL AND THEORETICAL STUDIES

NUMERICAL AND THEORETICAL STUDIES SELF-SIMILARITY OF WIND-WAVE SPECTRA. NUMERICAL AND THEORETICAL STUDIES Sergei I. Baduli (1), Adrei N. Pushkarev (,4), Doald Resio (3), Vladimir E. Zakharov (4,5) (1) P.P.Shirshov Istitute of Oceaology

More information

Nonequilibrium Excess Carriers in Semiconductors

Nonequilibrium Excess Carriers in Semiconductors Lecture 8 Semicoductor Physics VI Noequilibrium Excess Carriers i Semicoductors Noequilibrium coditios. Excess electros i the coductio bad ad excess holes i the valece bad Ambiolar trasort : Excess electros

More information

Chapter 6: BINOMIAL PROBABILITIES

Chapter 6: BINOMIAL PROBABILITIES Charles Bocelet, Probability, Statistics, ad Radom Sigals," Oxford Uiversity Press, 016. ISBN: 978-0-19-00051-0 Chater 6: BINOMIAL PROBABILITIES Sectios 6.1 Basics of the Biomial Distributio 6. Comutig

More information

FIR Filter Design: Part I

FIR 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 information

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated

More information

To make comparisons for two populations, consider whether the samples are independent or dependent.

To make comparisons for two populations, consider whether the samples are independent or dependent. Sociology 54 Testig for differeces betwee two samle meas Cocetually, comarig meas from two differet samles is the same as what we ve doe i oe-samle tests, ecet that ow the hyotheses focus o the arameters

More information

Jacob Hays Amit Pillay James DeFelice 4.1, 4.2, 4.3

Jacob Hays Amit Pillay James DeFelice 4.1, 4.2, 4.3 No-Parametric Techiques Jacob Hays Amit Pillay James DeFelice 4.1, 4.2, 4.3 Parametric vs. No-Parametric Parametric Based o Fuctios (e.g Normal Distributio) Uimodal Oly oe peak Ulikely real data cofies

More information

Estimation of the Mean and the ACVF

Estimation of the Mean and the ACVF Chapter 5 Estimatio of the Mea ad the ACVF A statioary process {X t } is characterized by its mea ad its autocovariace fuctio γ ), ad so by the autocorrelatio fuctio ρ ) I this chapter we preset the estimators

More information

Resampling 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.

Resampling 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 information

Random Variables, Sampling and Estimation

Random 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 information

ENGI 4421 Discrete Probability Distributions Page Discrete Probability Distributions [Navidi sections ; Devore sections

ENGI 4421 Discrete Probability Distributions Page Discrete Probability Distributions [Navidi sections ; Devore sections ENGI 441 Discrete Probability Distributios Page 9-01 Discrete Probability Distributios [Navidi sectios 4.1-4.4; Devore sectios 3.4-3.6] Chater 5 itroduced the cocet of robability mass fuctios for discrete

More information

Songklanakarin Journal of Science and Technology SJST R1 Teerapabolarn

Songklanakarin Journal of Science and Technology SJST R1 Teerapabolarn Soglaaari Joural of Sciece ad Techology SJST--.R Teeraabolar A No-uiform Boud o Biomial Aroimatio to the Beta Biomial Cumulative Distributio Fuctio Joural: Soglaaari Joural of Sciece ad Techology For Review

More information

Circle the single best answer for each multiple choice question. Your choice should be made clearly.

Circle the single best answer for each multiple choice question. Your choice should be made clearly. TEST #1 STA 4853 March 6, 2017 Name: Please read the followig directios. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directios This exam is closed book ad closed otes. There are 32 multiple choice questios.

More information

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

DETERMINATION 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 information

Coping with Insufficient Data: The Case of Household Automobile Holding Modeling by Ryuichi Kitamura and Toshiyuki Yamamoto

Coping with Insufficient Data: The Case of Household Automobile Holding Modeling by Ryuichi Kitamura and Toshiyuki Yamamoto Coig with Isufficiet Data: he Case of ousehold utomobile oldig odelig by Ryuichi Kitamura ad oshiyuki Yamamoto It is ofte the case that tyically available data do ot cotai all the variables that are desired

More information

Hypothesis Testing. Evaluation of Performance of Learned h. Issues. Trade-off Between Bias and Variance

Hypothesis Testing. Evaluation of Performance of Learned h. Issues. Trade-off Between Bias and Variance Hypothesis Testig Empirically evaluatig accuracy of hypotheses: importat activity i ML. Three questios: Give observed accuracy over a sample set, how well does this estimate apply over additioal samples?

More information

Estimation for Complete Data

Estimation for Complete Data Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of

More information

BIOSTATISTICAL METHODS FOR TRANSLATIONAL & CLINICAL RESEARCH

BIOSTATISTICAL METHODS FOR TRANSLATIONAL & CLINICAL RESEARCH BIOSAISICAL MEHODS FOR RANSLAIONAL & CLINICAL RESEARCH Direct Bioassays: REGRESSION APPLICAIONS COMPONENS OF A BIOASSAY he subject is usually a aimal, a huma tissue, or a bacteria culture, he aget is usually

More information

Parameter, Statistic and Random Samples

Parameter, Statistic and Random Samples Parameter, Statistic ad Radom Samples A parameter is a umber that describes the populatio. It is a fixed umber, but i practice we do ot kow its value. A statistic is a fuctio of the sample data, i.e.,

More information

Lecture 7: Properties of Random Samples

Lecture 7: Properties of Random Samples Lecture 7: Properties of Radom Samples 1 Cotiued From Last Class Theorem 1.1. Let X 1, X,...X be a radom sample from a populatio with mea µ ad variace σ

More information

THE KALMAN FILTER RAUL ROJAS

THE 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 information

Homework 3 Solutions

Homework 3 Solutions Math 4506 Sprig 04 Homework 3 Solutios. a The ACF of a MA process has a o-zero value oly at lags, 0, ad. Problem 4.3 from the textbook which you did t do, so I did t expect you to metio this shows that

More information

Nonlinear regression

Nonlinear regression oliear regressio How to aalyse data? How to aalyse data? Plot! How to aalyse data? Plot! Huma brai is oe the most powerfull computatioall tools Works differetly tha a computer What if data have o liear

More information

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated

More information

( ) = is larger than. the variance of X V

( ) = is larger than. the variance of X V Stat 400, sectio 6. Methods of Poit Estimatio otes by Tim Pilachoski A oit estimate of a arameter is a sigle umber that ca be regarded as a sesible value for The selected statistic is called the oit estimator

More information

Topic 18: Composite Hypotheses

Topic 18: Composite Hypotheses Toc 18: November, 211 Simple hypotheses limit us to a decisio betwee oe of two possible states of ature. This limitatio does ot allow us, uder the procedures of hypothesis testig to address the basic questio:

More information

tests 17.1 Simple versus compound

tests 17.1 Simple versus compound PAS204: Lecture 17. tests UMP ad asymtotic I this lecture, we will idetify UMP tests, wherever they exist, for comarig a simle ull hyothesis with a comoud alterative. We also look at costructig tests based

More information

Provläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE

Provläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE TECHNICAL REPORT CISPR 16-4-3 2004 AMENDMENT 1 2006-10 INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE Amedmet 1 Specificatio for radio disturbace ad immuity measurig apparatus ad methods Part 4-3:

More information

Fundamental Concepts: Surfaces and Curves

Fundamental Concepts: Surfaces and Curves UNDAMENTAL CONCEPTS: SURACES AND CURVES CHAPTER udametal Cocepts: Surfaces ad Curves. INTRODUCTION This chapter describes two geometrical objects, vi., surfaces ad curves because the pla a ver importat

More information

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

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

More information

ECE 442. Spring, Lecture - 4

ECE 442. Spring, Lecture - 4 ECE 44 Power Semicoductor Devices ad Itegrated circuits Srig, 6 Uiversity of Illiois at Chicago Lecture - 4 ecombiatio, geeratio, ad cotiuity equatio 1. Geeratio thermal, electrical, otical. ecombiatio

More information

Stat 421-SP2012 Interval Estimation Section

Stat 421-SP2012 Interval Estimation Section Stat 41-SP01 Iterval Estimatio Sectio 11.1-11. We ow uderstad (Chapter 10) how to fid poit estimators of a ukow parameter. o However, a poit estimate does ot provide ay iformatio about the ucertaity (possible

More information

Estimation of Parameters of Johnson s System of Distributions

Estimation of Parameters of Johnson s System of Distributions Joural of Moder Alied tatistical Methods Volume 0 Issue Article 9 --0 Estimatio of Parameters of Johso s ystem of Distributios Florece George Florida Iteratioal iversity, fgeorge@fiu.edu K. M. Ramachadra

More information

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 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 information

1 Inferential Methods for Correlation and Regression Analysis

1 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 information

The Poisson Distribution

The Poisson Distribution MATH 382 The Poisso Distributio Dr. Neal, WKU Oe of the importat distributios i probabilistic modelig is the Poisso Process X t that couts the umber of occurreces over a period of t uits of time. This

More information

BHW #13 1/ Cooper. ENGR 323 Probabilistic Analysis Beautiful Homework # 13

BHW #13 1/ Cooper. ENGR 323 Probabilistic Analysis Beautiful Homework # 13 BHW # /5 ENGR Probabilistic Aalysis Beautiful Homework # Three differet roads feed ito a particular freeway etrace. Suppose that durig a fixed time period, the umber of cars comig from each road oto the

More information

First Year Quantitative Comp Exam Spring, Part I - 203A. f X (x) = 0 otherwise

First Year Quantitative Comp Exam Spring, Part I - 203A. f X (x) = 0 otherwise First Year Quatitative Comp Exam Sprig, 2012 Istructio: There are three parts. Aswer every questio i every part. Questio I-1 Part I - 203A A radom variable X is distributed with the margial desity: >

More information

Last time: Moments of the Poisson distribution from its generating function. Example: Using telescope to measure intensity of an object

Last time: Moments of the Poisson distribution from its generating function. Example: Using telescope to measure intensity of an object 6.3 Stochastic Estimatio ad Cotrol, Fall 004 Lecture 7 Last time: Momets of the Poisso distributio from its geeratig fuctio. Gs () e dg µ e ds dg µ ( s) µ ( s) µ ( s) µ e ds dg X µ ds X s dg dg + ds ds

More information

Diagnosis of Kinematic Vertical Velocity in HYCOM. By George Halliwell, 28 November ( ) = z. v (1)

Diagnosis of Kinematic Vertical Velocity in HYCOM. By George Halliwell, 28 November ( ) = z. v (1) Diagosis of Kiematic Vertical Velocity i HYCOM By George Halliwell 28 ovember 2004 Overview The vertical velocity w i Cartesia coordiates is determied by vertically itegratig the cotiuity equatio dw (

More information

An Improved Algorithm and It s Application to Sinusoid Wave Frequency Estimation

An Improved Algorithm and It s Application to Sinusoid Wave Frequency Estimation Iteratioal Coerece o Iormatio ad Itelliget Computig IPCSI vol.8 () () IACSI Press, Sigapore A Improved Algorithm ad It s Applicatio to Siusoid Wave Frequecy Estimatio iaohog Huag ad Li Zhag College o Iormatio

More information

Estimation for a Class of Generalized State-Space Space Models:

Estimation for a Class of Generalized State-Space Space Models: Estimatio for a Class of Geeralized State-Sace Sace Models: Richard A. Davis ad Gabriel Rodriuez-Yam Colorado State Uiversit htt://www.stat.colostate.edu/~rdavis/lectures Joit work with: William Dusmuir

More information

Carriers in a semiconductor diffuse in a carrier gradient by random thermal motion and scattering from the lattice and impurities.

Carriers in a semiconductor diffuse in a carrier gradient by random thermal motion and scattering from the lattice and impurities. Diffusio of Carriers Wheever there is a cocetratio gradiet of mobile articles, they will diffuse from the regios of high cocetratio to the regios of low cocetratio, due to the radom motio. The diffusio

More information

Exercises and Problems

Exercises and Problems HW Chapter 4: Oe-Dimesioal Quatum Mechaics Coceptual Questios 4.. Five. 4.4.. is idepedet of. a b c mu ( E). a b m( ev 5 ev) c m(6 ev ev) Exercises ad Problems 4.. Model: Model the electro as a particle

More information

Elementary Statistics

Elementary Statistics Elemetary Statistics M. Ghamsary, Ph.D. Sprig 004 Chap 0 Descriptive Statistics Raw Data: Whe data are collected i origial form, they are called raw data. The followig are the scores o the first test of

More information

Nanomaterials for Photovoltaics (v11) 6. Homojunctions

Nanomaterials for Photovoltaics (v11) 6. Homojunctions Naomaterials for Photovoltaics (v11) 1 6. Homojuctios / juctio diode The most imortat device cocet for the coversio of light ito electrical curret is the / juctio diode. We first cosider isolated ad regios

More information

MA Advanced Econometrics: Properties of Least Squares Estimators

MA Advanced Econometrics: Properties of Least Squares Estimators MA Advaced Ecoometrics: Properties of Least Squares Estimators Karl Whela School of Ecoomics, UCD February 5, 20 Karl Whela UCD Least Squares Estimators February 5, 20 / 5 Part I Least Squares: Some Fiite-Sample

More information

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 4 Solutions [Numerical Methods]

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 4 Solutions [Numerical Methods] ENGI 3 Advaced Calculus or Egieerig Facult o Egieerig ad Applied Sciece Problem Set Solutios [Numerical Methods]. Use Simpso s rule with our itervals to estimate I si d a, b, h a si si.889 si 3 si.889

More information

Output Analysis (2, Chapters 10 &11 Law)

Output Analysis (2, Chapters 10 &11 Law) B. Maddah ENMG 6 Simulatio Output Aalysis (, Chapters 10 &11 Law) Comparig alterative system cofiguratio Sice the output of a simulatio is radom, the comparig differet systems via simulatio should be doe

More information

x = Pr ( X (n) βx ) =

x = Pr ( X (n) βx ) = Exercise 93 / page 45 The desity of a variable X i i 1 is fx α α a For α kow let say equal to α α > fx α α x α Pr X i x < x < Usig a Pivotal Quatity: x α 1 < x < α > x α 1 ad We solve i a similar way as

More information

1 Introduction to reducing variance in Monte Carlo simulations

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

More information

ECE 145A / 218 C, notes set 13: Very Short Summary of Noise

ECE 145A / 218 C, notes set 13: Very Short Summary of Noise class otes, M. odwell, copyrighted 009 C 45A / 8 C, otes set 3: Very hort ummary o Noise Mark odwell Uiversity o Calioria, ata Barbara rodwell@ece.ucsb.edu 805-893-344, 805-893-36 ax Backgroud / tet class

More information

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 5 STATISTICS II. Mea ad stadard error of sample data. Biomial distributio. Normal distributio 4. Samplig 5. Cofidece itervals

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.

MASSACHUSETTS 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 information

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

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

More information

Time series models 2007

Time series models 2007 Norwegia Uiversity of Sciece ad Techology Departmet of Mathematical Scieces Solutios to problem sheet 1, 2007 Exercise 1.1 a Let Sc = E[Y c 2 ]. The This gives Sc = EY 2 2cEY + c 2 ds dc = 2EY + 2c = 0

More information

Dimension of a Maximum Volume

Dimension of a Maximum Volume Dimesio of a Maximum Volume Robert Kreczer Deartmet of Mathematics ad Comutig Uiversity of Wiscosi-Steves Poit Steves Poit, WI 54481 Phoe: (715) 346-3754 Email: rkrecze@uwsmail.uws.edu 1. INTRODUCTION.

More information

IX Stochastic Chemical Kinetics

IX Stochastic Chemical Kinetics IX Stochastic Chemical Kietics Surrisi this hae whe we take the discreteess o molecule umber seriously, abadoi the otio that chemical cocetratios may be treated as cotiuous variables. Here we will show

More information

Phys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12

Phys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12 Physics Departmet, Yarmouk Uiversity, Irbid Jorda Phys. Mathematical Physics Dr. Nidal M. Ershaidat Doc. Fourier Series Deiitio A Fourier series is a expasio o a periodic uctio (x) i terms o a iiite sum

More information

Hybridized Heredity In Support Vector Machine

Hybridized Heredity In Support Vector Machine Hybridized Heredity I Suort Vector Machie May 2015 Hybridized Heredity I Suort Vector Machie Timothy Idowu Yougmi Park Uiversity of Wiscosi-Madiso idowu@stat.wisc.edu yougmi@stat.wisc.edu May 2015 Abstract

More information

Entropy Rates and Asymptotic Equipartition

Entropy Rates and Asymptotic Equipartition Chapter 29 Etropy Rates ad Asymptotic Equipartitio Sectio 29. itroduces the etropy rate the asymptotic etropy per time-step of a stochastic process ad shows that it is well-defied; ad similarly for iformatio,

More information

Machine Learning Brett Bernstein

Machine Learning Brett Bernstein Machie Learig Brett Berstei Week 2 Lecture: Cocept Check Exercises Starred problems are optioal. Excess Risk Decompositio 1. Let X = Y = {1, 2,..., 10}, A = {1,..., 10, 11} ad suppose the data distributio

More information

Computing Confidence Intervals for Sample Data

Computing Confidence Intervals for Sample Data Computig Cofidece Itervals for Sample Data Topics Use of Statistics Sources of errors Accuracy, precisio, resolutio A mathematical model of errors Cofidece itervals For meas For variaces For proportios

More information

Joint Probability Distributions and Random Samples. Jointly Distributed Random Variables. Chapter { }

Joint Probability Distributions and Random Samples. Jointly Distributed Random Variables. Chapter { } UCLA STAT A Applied Probability & Statistics for Egieers Istructor: Ivo Diov, Asst. Prof. I Statistics ad Neurology Teachig Assistat: Neda Farziia, UCLA Statistics Uiversity of Califoria, Los Ageles, Sprig

More information

The Blood Testing Problem

The Blood Testing Problem Teachig Cotemorary Mathematics Coferece Feb 6-7, 00 The Blood Testig Problem Suose that you have a large oulatio that you wish to test for a certai characteristic i their blood or urie (for examle, testig

More information

Classification of DT signals

Classification of DT signals Comlex exoetial A discrete time sigal may be comlex valued I digital commuicatios comlex sigals arise aturally A comlex sigal may be rereseted i two forms: jarg { z( ) } { } z ( ) = Re { z ( )} + jim {

More information

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

Signal 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 information

WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT

WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? Harold G. Loomis Hoolulu, HI ABSTRACT Most coastal locatios have few if ay records of tsuami wave heights obtaied over various time periods. Still

More information

Package pgnorm. November 24, 2015

Package pgnorm. November 24, 2015 Tye Package Title The -Geeralized Normal Distributio Versio 2.0 Date 2015-11-23 Author Package gorm November 24, 2015 Maitaier Evaluatio of the df ad the cdf of the uivariate,

More information

5. Fast NLMS-OCF Algorithm

5. Fast NLMS-OCF Algorithm 5. Fast LMS-OCF Algorithm The LMS-OCF algorithm preseted i Chapter, which relies o Gram-Schmidt orthogoalizatio, has a compleity O ( M ). The square-law depedece o computatioal requiremets o the umber

More information

Stochastic Structural Dynamics. Lecture-28. Monte Carlo simulation approach-4

Stochastic Structural Dynamics. Lecture-28. Monte Carlo simulation approach-4 Stochastic Structural Dyamics Lecture-8 Mote Carlo simulatio approach-4 Dr C S Maohar Departmet of Civil Egieerig Professor of Structural Egieerig Idia Istitute of Sciece Bagalore 56 Idia maohar@civil.iisc.eret.i

More information

2 One-Dimensional Elasticity

2 One-Dimensional Elasticity Oe-Dimesioal Elasticity There are two tyes of oe-dimesioal roblems, the elastostatic roblem ad the elastodyamic roblem. The elastostatic roblem gives rise to a secod order differetial equatio i dislacemet

More information

COMPUTING FOURIER SERIES

COMPUTING FOURIER SERIES COMPUTING FOURIER SERIES Overview We have see i revious otes how we ca use the fact that si ad cos rereset comlete orthogoal fuctios over the iterval [-,] to allow us to determie the coefficiets of a Fourier

More information

Lecture 9: September 19

Lecture 9: September 19 36-700: Probability ad Mathematical Statistics I Fall 206 Lecturer: Siva Balakrisha Lecture 9: September 9 9. Review ad Outlie Last class we discussed: Statistical estimatio broadly Pot estimatio Bias-Variace

More information

The Poisson Process *

The Poisson Process * OpeStax-CNX module: m11255 1 The Poisso Process * Do Johso This work is produced by OpeStax-CNX ad licesed uder the Creative Commos Attributio Licese 1.0 Some sigals have o waveform. Cosider the measuremet

More information

This exam contains 19 pages (including this cover page) and 10 questions. A Formulae sheet is provided with the exam.

This exam contains 19 pages (including this cover page) and 10 questions. A Formulae sheet is provided with the exam. Probability ad Statistics FS 07 Secod Sessio Exam 09.0.08 Time Limit: 80 Miutes Name: Studet ID: This exam cotais 9 pages (icludig this cover page) ad 0 questios. A Formulae sheet is provided with the

More information

January 25, 2017 INTRODUCTION TO MATHEMATICAL STATISTICS

January 25, 2017 INTRODUCTION TO MATHEMATICAL STATISTICS Jauary 25, 207 INTRODUCTION TO MATHEMATICAL STATISTICS Abstract. A basic itroductio to statistics assumig kowledge of probability theory.. Probability I a typical udergraduate problem i probability, we

More information

Expectation and Variance of a random variable

Expectation 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 information

Measurement uncertainty of the sound absorption

Measurement 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 information

Chapter 12 Sound Waves

Chapter 12 Sound Waves Chapter 2 Soud Waves We study the properties ad detectio o a particular type o wave soud waves. A speaker geerates soud. The desity o the air chages as the wave propagates. The rage o requecies that ca

More information

2.2 One-dimensional Elastodynamics

2.2 One-dimensional Elastodynamics Sectio.. Oe-dimesioal Elastodyamics I rigid body dyamics, it is assumed that whe a force is alied to oe oit of a object, every other oit i the object is set i motio simultaeously. O the other had, i static

More information

Machine Learning. Logistic Regression -- generative verses discriminative classifier. Le Song /15-781, Spring 2008

Machine Learning. Logistic Regression -- generative verses discriminative classifier. Le Song /15-781, Spring 2008 Machie Learig 070/578 Srig 008 Logistic Regressio geerative verses discrimiative classifier Le Sog Lecture 5 Setember 4 0 Based o slides from Eric Xig CMU Readig: Cha. 3..34 CB Geerative vs. Discrimiative

More information

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare Spectral Aalysis This week i lab Your ext experimet Homework is to prepare Next classes: 3/26 ad 3/28 Aero Testig, Fracture Toughess Testig Read the Experimets 5 ad 7 sectios of the course maual Spectral

More information

Problem Set 4 Due Oct, 12

Problem 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 information

ANALYSIS OF EXPERIMENTAL ERRORS

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

More information

The Hong Kong University of Science & Technology ISOM551 Introductory Statistics for Business Assignment 3 Suggested Solution

The Hong Kong University of Science & Technology ISOM551 Introductory Statistics for Business Assignment 3 Suggested Solution The Hog Kog Uiversity of ciece & Techology IOM55 Itroductory tatistics for Busiess Assigmet 3 uggested olutio Note All values of statistics i Q ad Q4 are obtaied by Excel. Qa. Let be the robability that

More information

( ) = ( ) + ( ), One Degree of Freedom, Harmonically Excited Vibrations. 1 Forced Harmonic Vibration. t dies out with time under each of.

( ) = ( ) + ( ), One Degree of Freedom, Harmonically Excited Vibrations. 1 Forced Harmonic Vibration. t dies out with time under each of. Oe Degree of Freedom, Harmoically Excited Vibratios Forced Harmoic Vibratio A mechaical syem is said to udergo forced vibratio wheever exteral eergy is sulied to the syem durig vibratio Exteral eergy ca

More information

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

CEE 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 information

DS 100: Principles and Techniques of Data Science Date: April 13, Discussion #10

DS 100: Principles and Techniques of Data Science Date: April 13, Discussion #10 DS 00: Priciples ad Techiques of Data Sciece Date: April 3, 208 Name: Hypothesis Testig Discussio #0. Defie these terms below as they relate to hypothesis testig. a) Data Geeratio Model: Solutio: A set

More information

2(25) Mean / average / expected value of a stochastic variable X: Variance of a stochastic variable X: 1(25)

2(25) Mean / average / expected value of a stochastic variable X: Variance of a stochastic variable X: 1(25) Lecture 5: Codig of Aalog Sources Samplig ad Quatizatio Images ad souds are ot origially digital! The are cotiuous sigals i space/time as well as amplitude Typical model of a aalog source: A statioary

More information

John H. J. Einmahl Tilburg University, NL. Juan Juan Cai Tilburg University, NL

John H. J. Einmahl Tilburg University, NL. Juan Juan Cai Tilburg University, NL Estimatio of the margial exected shortfall Laures de Haa, Poitiers, 202 Estimatio of the margial exected shortfall Jua Jua Cai Tilburg iversity, NL Laures de Haa Erasmus iversity Rotterdam, NL iversity

More information

Statistical Noise Models and Diagnostics

Statistical 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 information