14.02.2013 Views

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Appendix I / Utilities / General<br />

_ _<br />

C<br />

:<br />

65539 ⋅ S i −1<br />

≥ 0 ⇒ Si<br />

: = 65539 ⋅ S i −1<br />

⎛<br />

A ⋅ ⎜<br />

⎜<br />

⎝<br />

⎡<br />

12<br />

= ⎢ ∑<br />

i:<br />

= 1<br />

⎢⎣<br />

⎤<br />

−10<br />

( S ⋅ 4.<br />

656613⋅10<br />

) ⎥ − 6 ⎟ + B<br />

i<br />

22.2-5<br />

⎥⎦<br />

Equation 22.2-24<br />

⎞<br />

⎟<br />

⎠<br />

Equation 22.2-25<br />

GAUSS returns a Gaussian distributed random variable (C) from the <strong>pdf</strong>,<br />

p<br />

( C )<br />

:<br />

A ⋅<br />

1 ⎛<br />

⋅ exp ⎜ −<br />

2 ⋅ π ⎜<br />

⎝<br />

( C − B )<br />

= 2<br />

2 ⋅ A<br />

2<br />

⎞<br />

⎟<br />

⎠<br />

Equation 22.2-26<br />

The exponential represents a scaling so that the cumulative probability, the<br />

integral of p(C) over C ∈ [-∞,∞], is 1; the peak of p(C) being 1/√2π at (B).<br />

The following special values are useful when analysing stochastic filters,<br />

( C ∈ [ −1<br />

, 1 ] ) : 0.<br />

6837<br />

p =<br />

( C ∈ [ − 2 , 2 ] ) : 0.<br />

9545<br />

p =<br />

( C ∈ [ − 3 , 3 ] ) : 0.<br />

9973<br />

p =<br />

When initialising SIGNV the seed from GAUSS is used,<br />

0<br />

( 1 ) : S0<br />

S =<br />

[ 1(<br />

1 ) 999 ] ⇒ S ( k + 1 ) : = S ( k )<br />

k ∈<br />

0<br />

12<br />

Equation 22.2-27<br />

Equation 22.2-28<br />

Equation 22.2-29<br />

Equation 22.2-30<br />

Equation 22.2-31<br />

GAUSS uses integer overflow to generate random numbers and must be<br />

compiled with overflow checking disabled. Initialisation of the random<br />

number generator is completed by taking S12(1000) and using it to seed<br />

1000 sequences stored in vector IRAND. During this process the seeds<br />

returned by GAUSS are manipulated as follows,<br />

( 1001 ) : S ( 1000 )<br />

S0 = 12<br />

Equation 22.2-32

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!