The Basics interpolation Model fitting

Simulating the Ionosphere Anna Scaife Cavendish Astrophysics University of Cambridge.

MCCT SKADS Summer School on Third Generation Calibration

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Outline

1

The Basics

2

interpolation

3

Model fitting

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Z–Jones Z =

 iφ e iono 0

0



eiφiono

φiono = −25λTEC TEC = ....

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Z–Jones Z =

 iφ e iono 0

0



eiφiono

φiono = −25λTEC TEC = ....

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Z–Jones Z =

 iφ e iono 0

0



eiφiono

φiono = −25λTEC TEC = ....

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

TEC 1 TECU = 1016 m−2 Z

ne d` + const

TEC = 0

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

It seems that we should be able to predict the behaviour of the ionosphere... ...but, in fact we are limited to interpolation of scattered measurements. Counter-intuitively, simulating the ionosphere is in fact an inverse problem. Fortunately there are lots of ways of approaching inverse problems.

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Interpolation

Let’s look at interpolating between measurements... DEFINITION: Approximating measurements at intermediate scales/positions from scattered measurements. We have sparse measurements i.e. We have under-sampled data

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Interpolation

Let’s look at interpolating between measurements... DEFINITION: Approximating measurements at intermediate scales/positions from scattered measurements. We have sparse measurements i.e. We have under-sampled data

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Interpolation schemes

1

Kriging

2

Triangulation based

3

Natural neighbour

4

Splining

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Kriging

Semivariance γ(h) =

1 (f (x) − f (x + h))2 2

A typical model:  3 h s( 2 ( a ) − 12 ( ha )3 ), 0 ≤ h ≤ a γ= s, h>a

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Triangulation Given three measured points we can interpolate to any point within the triangle using: Delaunay triangulation f (x, y ) =

3 X

φi (x, y )fi

i=1

φi (x) is our basis function In a simple case we can use linear equations: f (x, y ) = c1 x + c2 y + c3 We can just solve Ac = f, where f = (f1 , f2 , f3 )T and A = {(xi , yi , 1)}1≤i≤3 . Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Natural neighbour Voronoi tesselation Voronoi tesselation divides the data into cells defined by the positions of the measurements. We use the interpolation point to define a new Voronoi cell. The value of this cell can be evaluated as the weighted sum of the contributions from its overlapping cells.

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Splining

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Bayes Theorem p(D|M, Θ) = posterior =

L(M, Θ|D)π(Θ) E

likelihood × prior Evidence

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Bayes Theorem p(D|M, Θ) = posterior =

L(M, Θ|D)π(Θ) E

likelihood × prior Evidence

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Maximum likelihood

We maximise the likelihood of the DATA w.r.t the MODEL. When we perform a χ2 test we are in fact calculating a Gaussian ML. N(Di , σ) ∝ exp

Anna Scaife

−(Di − Mi )2 2σ 2

Z–Jones

The Basics interpolation Model fitting

Maximum A Posteriori If we know something about our parameters then we can utilise that prior information. Example Say we are fitting a spectral index... S = Aν −α

ln L = −0.5

(di − mi )2 σi2

+ ln π

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Maximum A Posteriori If we know something about our parameters then we can utilise that prior information. Example Say we are fitting a spectral index... S = Aν −α

ln L = −0.5

(di − mi )2 σi2

+ ln π

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Model comparison

What if we have more than one model...? The Evidence: Z E = L(Θ)π(Θ)d D Θ The model selection ratio: R=

Pr(D|H1 )Pr(H1 ) Z1 Pr(H1 ) Pr(H1 |D) = = Pr(H2 |D) Pr(D|H1 )Pr(H1 ) Z2 Pr(H2 )

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Worked Example Let’s look at Kriging again... Kriging variance γ(h) =

1 (f (x) − f (x + h))2 2

TEC(x0 ) =

n X

λi TEC(xi )

i=1

∆=2

n X

λi γ(xi , x0 ) −

i=1

n n X X i=1 j=1

Γij λ = Γi0 Anna Scaife

Z–Jones

λi λj γ(xi , x0 )

The Basics interpolation Model fitting

Worked Example Let’s look at Kriging again... Kriging variance γ(h) =

1 (f (x) − f (x + h))2 2

TEC(x0 ) =

n X

λi TEC(xi )

i=1

∆=2

n X

λi γ(xi , x0 ) −

i=1

n n X X i=1 j=1

Γij λ = Γi0 Anna Scaife

Z–Jones

λi λj γ(xi , x0 )

The Basics interpolation Model fitting

Worked Example Let’s look at Kriging again... Kriging variance γ(h) =

1 (f (x) − f (x + h))2 2

TEC(x0 ) =

n X

λi TEC(xi )

i=1

∆=2

n X

λi γ(xi , x0 ) −

i=1

n n X X i=1 j=1

Γij λ = Γi0 Anna Scaife

Z–Jones

λi λj γ(xi , x0 )

The Basics interpolation Model fitting

We need a model for γ We can calculate γ directly from the data: γ(xi − xj ) =

1 (TEC(xi ) − TEC(xj ))2 2

But we need an analytic form for γ So we have to pick a model... Typical models would be: γ1 = a + b ∗ h γ2 = a + b ∗ hα

Anna Scaife

Z–Jones

The Basics interpolation Model fitting

Anna Scaife

Z–Jones

S. Dodelson. Modern Cosmology. Elsevier, 2003. W. Hu & S. Dodelson Cosmic Microwave Background Anisotropies Ann. Rev. Astron. and Astrophys., 2002 R. Subrahmanyan & R. Ekers CMB observations using the SKA SKA Memo Series, 26

Anna Scaife

Z–Jones

## Simulating the Ionosphere - GitHub

Sep 30, 2009 - DEFINITION: Approximating measurements at intermediate scales/positions from scattered measurements. We have sparse measurements.

#### Recommend Documents

Simulating a two dimensional particle in a square quantum ... - GitHub
5.3.12 void runCuda(cudaGraphicsResource **resource) . . . . . 17 ... the probabilities of the position and the energy of the particle at each state. ..... 2PDCurses is an alternative suggested by many http://pdcurses.sourceforge.net/. The.

Simulating History: The Problem of Contingency - CiteSeerX
outcome measures that affect the security and prosperity of the world community. The problem of contingency for testing hypotheses about the causes of certain.

Simulating Reach Motions
technique to fit polynomial equations to the angular displacements of each ... A description of current research on human motion simulations at the University of.

N00014-95-10223 to John Anderson at Carnegie Mellon University. ..... In fact, even when given explicit instructions that the computer could not ...... Computer Science and Engineering, Oregon Graduate Institute of Science & Technology.

Simulating Computational Societies
and implementors, this can be used to model, analyse and explain the .... model and the software developed to facilitate experimentation with agent so-.