04_mathima

10
1 Δευτέρα, 12 Οκτωβρίου 2009 Ροπογεννήτριες: Ροπογεννήτρια συνάρτηση μιας τ.μ. Χ, ορίζεται ως η αναναμενόμενη τιμή: () ( ) ( ) () ( ) , , , , , xu uX X X xu e PX x ή ί M u Ee x u e f x dx X x ή ί διακριτ περ πτωση συνεχ ς περ πτωση = = = = Διδιάστατη: ( ) ( ) 11 2 2 1 2 , 1 2 , xu xu X X M uu Ee + = Πολυδιάστατη: () 1 v i i i xu X M u Ee = = Θεώρημα: Ν.δ.ό. ( ) () 0 X u dM u Ex du = = Απόδειξη: α) Διακριτή Περίπτωση: () ( ) 1 n xu X i M u e PX x = = = ( ) ( ) 1 n xu X i dM u d e PX x du du = = = = ( ) 1 n xu i d e PX x du = = = ( ) 1 n xu i xe PX x = = ( ) ( ) 1 0 0 n xu X i u u dM u xe PX x du = = = = = = ( ) 0 1 0 n x i u xe PX x = = = = ( ) () 1 n i xPX x Ex = = = Γενικά: ( ) ( ) 0 k k X k u dM u Ex du = = β) Συνεχής Περίπτωση: ( ) X M u = () ux e f x dx −∞ = () 2 2 3 3 1 ... 2 3! ux ux ux f x dx −∞ + + + + = () () () 2 2 ... 2 ux f x dx uxf x dx f x dx −∞ −∞ −∞ + + + = () () () 2 2 ... 2 u f x dx u xf x dx xf x dx −∞ −∞ −∞ + + + () () ( ) ( ) 2 3 2 3 1 ... 2 3! ux X u u M u Ee uE x Ex Ex = =+ + + + ροπογεννήτρια.

description

xu xu X X de P X x du = = 2 3 2 3 u f x dx u xf x dx xf xdx ux f x dx uxf x dx f x dx M u ePX x 1 1 2 2 1 2 ( ) xPX x Ex ⋅ = = u u Mu Ee uE x E x E x  = = + + + +   ροπογεννήτρια . M u = = = ∞ ∞ ∞ n xu () 2 2 ∆ευτέρα, 12 Οκτωβρίου 2009 Ροπογεννήτριες: M u u E e + = Πολυδιάστατη : efxdx ∞ ∞ ∞ ∑ ∫ M u E e = n x = ⇒ xe P X x n xu xe P X x ( ) ( ) ( ) ( )   ∑  =     ∑ ∑ ∑ ∑ ∑

Transcript of 04_mathima

Page 1: 04_mathima

1

∆ευτέρα, 12 Οκτωβρίου 2009

Ροπογεννήτριες: Ροπογεννήτρια συνάρτηση µιας τ.µ. Χ, ορίζεται ως η αναναµενόµενη τιµή:

( ) ( )( )

( ) ( )

,

, ,

, ,

xu

uX XX

xu

e P X x ή ί

M u E e x u

e f x dx X x ή ί

διακριτ περ πτωση

συνεχ ς περ πτωση

=

= ∈ = ∈ =

∫ℝ ℝ

∆ιδιάστατη:

( ) ( )1 1 2 2

1 2, 1 2,x u x u

X XM u u E e+=

Πολυδιάστατη:

( ) 1

v

i i

i

x u

XM u E e =

∑ =

Θεώρηµα:

Ν.δ.ό. ( ) ( )

0

X

u

dM uE x

du=

=

Απόδειξη:

α) ∆ιακριτή Περίπτωση:

( ) ( )1

nxu

X

i

M u e P X x=

= = ⇒∑

( ) ( )1

nxuX

i

dM u de P X x

du du =

= = =∑ ( )1

nxu

i

de P X x

du=

= =∑ ( )1

nxu

i

xe P X x=

= ⇒∑

( ) ( )1 00

nxuX

i uu

dM uxe P X x

du = ==

= = =∑ ( )0

1 0

nx

i u

xe P X x⋅

= =

= =∑

( ) ( )1

n

i

x P X x E x=

⋅ = =∑

Γενικά:

( ) ( )0

k

k X

k

u

d M uE x

du=

=

β) Συνεχής Περίπτωση:

( )XM u = ( )uxe f x dx

−∞

=∫ ( )2 2 3 3

1 ...2 3!

u x u xux f x dx

−∞

+ + + + =

( ) ( ) ( )2 2

...2

u xf x dx uxf x dx f x dx

∞ ∞ ∞

−∞ −∞ −∞

+ + + =∫ ∫ ∫

( ) ( ) ( )2

2 ...2

uf x dx u xf x dx x f x dx

∞ ∞ ∞

−∞ −∞ −∞

+ + + ⇒∫ ∫ ∫

( ) ( ) ( ) ( )2 3

2 31 ...2 3!

ux

X

u uM u E e uE x E x E x = = + + + + ροπογεννήτρια.

Page 2: 04_mathima

2

Οπότε ( ) ( ) ( ) ( ) ( )

22 3

00

...2!

x

uu

dM u uE x uE x E x E x

du==

= + + + =

και

( ) ( ) ( ) ( )2

2 3 2

20

0

...x

uu

d M uE x uE x E x

du ==

= + + =

………………………..

( ) ( ) ( ) ( )1

00

...

k

k k kx

ku

u

d M uE x uE x E x

du

+

==

= + + =

Παραδείγµατα

Χ~Poisson(λ)

( ) ux

XM u E e = = ( )0

ux

x

e P X x∞

=

= =∑ 0 !

xux

x

e ex

λ λ∞−

=

⋅ =∑ ( )

0 !

xu

x

ee

x

λλ∞

=

=∑

( )1u

u u ee ee e e eλλ λ λ λ −− −⋅ = = ⇒ ( ) ( )1ue

XM u eλ −

=

( )0

0

ue uX

uu

dM ue e

du

λ λ λ−

==

= ⋅ = 0e eλ λ λ λ− ⋅ =

Χ~∆ιωνυµική(n,p)

( ) ( )0

nxu

X

i

M u e P X x=

= = =∑ ( )0

1n

n xxu x

i

ne p p

x

=

− ⇒

∑ ( ) ( )

0

1n

x n xu

i

npe p

x

=

− =

( )( )1n

upe p+ − = ( )n

upe q+ ⇒

( ) ( )( )( )1n

uXdM u dpe p

du du= + − =

( )( )( ) ( )( ) 1

1 1n n

u u udpe p n pe p pe

du

−+ − = + −

( ) ( )( ) 1

00

1n

u uX

uu

dM un pe p pe

du

==

= + − = ( )( ) 1

0 01n

n pe p pe−

+ − =

( ) 11

nn p p p

−+ − = [ ]np E x=

Page 3: 04_mathima

3

Θεώρηµα:

Για την κανονική µορφή της Ε.Ο.Κ. ισχύουν:

α) Αναµενόµενη τιµή: ( ) ( )i

i

A nE T X

n

∂= ∂

.

β) Συνδιακύµανση: ( ) ( ) ( )2

,i j

i j

A nCov T X T X

n n

∂ = ∂ ∂

.

γ) Ροπογεννήτρια: ( ) ( ) ( ) ( ){ }1

exp exps

T i i

i

M u E u T X A u n A n=

= = + −

∑ � � �

Απόδειξη:

α) ( ), 1f x n dx

−∞

= ⇒∫ ( ) ( ) ( )1

exps

i i

i

nT x A n h x dx

=−∞

− =

∑∫

( ) ( ) ( )1

exp exp 1s

i i

i

A n nT x h x dx

=−∞

− = ⇒

∑∫

( ) ( ) ( )1

exp exps

i i

i

nT x h x dx A n

=−∞

=

∑∫ .

Παραγωγίζοντας ως προς in , παίρνουµε:

( ) ( ) ( ){ }1

exp exps

i i

ii i

nT x h x dx A nn n

=−∞

∂ ∂ = ⇒ ∂ ∂

∑∫

και, αλλάζοντας τη σειρά «παραγώγιση - ολοκλήρωση» παίρνουµε:

( ) ( ) ( ){ } ( )1

exp exps

i i

ii i

A nnT x h x dx A n

n n

=−∞

∂ ∂ = ⇒ ∂ ∂

∑∫

( ) ( ) ( ) ( ) ( )1 1

exp exps s

i i i i

i ii i

A nnT x nT x h x dx A n

n n

= =−∞

∂∂ = ⇒ ∂ ∂

∑ ∑∫

( ) ( ) ( ) ( ) ( ) ( )1

exp exp *s

i i i

i i

A nT x nT x h x dx A n

n

=−∞

∂ = ⇒ ∂

∑∫

( ) ( ) ( ) ( ) ( )1

exp exps

i i i

i i

A nA n T x nT x h x dx

n

=−∞

∂ − = ⇒ ∂

∑∫

( ) ( ) ( ) ( ) ( )1

exp exps

i i i

i i

A nT x nT x A n h x dx

n

=−∞

∂ − = ⇒ ∂

∑∫

( ) ( ) ( ) ( ) ( )1

exps

i i i

i i

A nT x nT x A n h x dx

n

=−∞

∂ − = ⇒ ∂

∑∫ ( )( ) ( )i

i

A nE T x

n

∂=

β) Παραγωγίζουµε την (*) ως προς j

n .

( ) ( ) ( ) ( ) ( )1

exp exps

i i i

ij j i

A nT x nT x h x dx A n

nη η

=−∞

∂∂ ∂ = ⇒ ∂ ∂ ∂

∑∫

Page 4: 04_mathima

4

( ) ( ) ( ) ( ) ( ) ( ) ( ){ }1

expexp exp

s

i i i

ij j i i j

A nA n A nT x nT x h x dx A n

n n n n n

=−∞

∂ ∂ ∂ ∂ ∂ = + ⇒ ∂ ∂ ∂ ∂ ∂ ∑∫

( ) ( ) ( ) ( ) ( ) ( ) ( )( )

2

1

exp exp exps

i i i

ij i j i j

A nA n A nT x nT x h x dx A n A n

n n n n n

=−∞

∂ ∂ ∂ ∂ = + ⇒ ∂ ∂ ∂ ∂ ∂ ∑∫

( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )2

1

exp exps

i j i i

i j i j i

A n A nA nT x T x nT x h x dx A n

n n n n

=−∞

∂ ∂ ∂ = + ⇒ ∂ ∂ ∂ ∂ ∑∫

( ) ( ) ( ) ( ) ( )( ) ( ) ( )2

1

exps

i j i i

i j i j i

A n A n A nT x T x nT x A n h x dx

n n n n

=−∞

∂ ∂ ∂ − = + ⇒ ∂ ∂ ∂ ∂ ∑∫

( ) ( ) ( ) ( ) ( )2

i j i j

j i

A nE T x T x E T x E T x

n n

∂ − = ⇒ ∂ ∂

( ) ( ) ( )2

,i j

j i

A nCov T x T x

n n

∂ = ∂ ∂

διότι, ως γνωστόν: [ ] [ ] [ ] [ ],Cov X Y E X Y E X E Y= ⋅ −

γ) ( )( )

( ) ( )1

1

exp ;

s

i i

i

su T x

T i i

i

M u E e u T x f x n dx=

=−∞

∑ = = =

∑∫

( ) ( ) ( ) ( )1 1

exp exps s

i i i i

i i

u T x nT x A n h x dx

= =−∞

− =

∑ ∑∫

(και, προσθαφαιρώντας το ( )A u n+ )

( ){ } ( ){ } ( ) ( ) ( ) ( )1

exp exp exps

i i i

i

A n A u n u n T x A u n h x dx

=−∞

− + + − + =

∑∫

( ) ( ){ } ( ) ( ) ( ) ( )1

exp exps

i i i

i

A u n A n u n T x A u n h x dx

=−∞

+ − + − + ⇒

∑∫

( ) ( ) ( ){ }expT

M u A u n A n= + −

διότι: ( ) ( ) ( ) ( )1

exp 1s

i i i

i

u n T x A u n h x dx

=−∞

+ − + =

∑∫

Συνοπτικά:

( )( ) ( )i

i

A nE T x

n

∂=

∂ ( ) ( ) ( )2

,i j

j i

A nCov T x T x

n n

∂ = ∂ ∂

( ) ( ) ( ){ }expT

M u A u n A n= + −

και για 1s = (ειδική περίπτωση)

( )( ) ( )AE T x

ηη

∂=

∂ ( ) ( ) ( ) ( )2

2,

ACov T x T x V T x

ηη

∂= = ∂

( ) ( ) ( ){ }expT

M u A u Aη η= + −

Page 5: 04_mathima

5

Ασκήσεις

Έστω ~X Poisson

( );!

x

f x ex

λ λλ −= = ( ) ( ){ } 1

exp exp log!

x

xλ λ− = ( ){ } 1

exp log!

x

xλ λ− =

{ } 1exp log

!x

xλ λ− . Άρα ανήκει ΕΟΚ µε:

( ) logn λ λ= , ( )T x x= , ( )B λ λ= − , ( ) 1

!h x

x=

Κανονική µορφή:

( ) logn nλ λ= = ⇒ neλ =

( ) { } 1; exp log

!f x x

xλ λ λ= − ⇒ ( ) { } 1

; exp!

f x x ex

ηη η= − , δηλαδή ( )A eηη .

( ) ( ) ( )A nE T x E x

n

∂= = = ∂

( )n

ne

en

λ∂

= =∂

( ) ( ) ( ) ( )2

2,

A nCov T x T x V T x

n

∂= = = ∂

( )2

2

n

ne

en

λ∂

= =∂

( )TM u = ( ) ( ){ }exp A u n A n+ − = { }exp u n ne e

+ − = ( ){ }exp 1n u

e e − =

( ){ }exp 1u

eλ −

Άσκηση: ( )~ ,X Binomial v p• ν γνωστό, p άγνωστο.

( ) ( ), 1v xx

vf x p p p

p

− = − =

( )11

x

vv pp

x p

− = −

( )exp ln ln 11

v px v p

x p

+ − −

άρα ανήκει στην ΕΟΚ, µε ( )v

h xx

=

, ( )T X x= , ( ) ln1

pn p

p=

− και

( ) ( )ln 1B p v p= − −

Μετατρέπουµε σε κανονική µορφή: ( ) ( ) ( ){ } ( )exp n p T x A n h x− − , θέτοντας:

( ) ln1

pn p n

p= = ⇒

− 1

npe

p= ⇒

− n np e pe= − ⇒ n n

p pe e+ = ⇒

( )1 n np e e+ = ⇒

1

n

n

ep

e=

+

Page 6: 04_mathima

6

( ), exp ln 11

n

n

v ef x n xn v

x e

= + − = +

1

exp ln1

n n

n

v e exn v

x e

+ −+ = +

1exp ln

1n

vxn v

x e

+ = + ( ){ }exp ln 1 n

vxn v e

x

− +

, δηλαδή

( ) ( )ln 1 nA n v e= +

( ) ( ) ( )E T X E X A nn

∂= = = ∂

( )( )ln 1 nv e

n

∂+ =

1

1

n

nv e

e=

+ vp

( ) ( ) ( ),Cov T X T X V X= = ( )2

2A n

n

∂=

1

n

n

ev

n e

∂= ∂ +

( )( )2

1

1

n n n n

n

e e e ev

e

+ −=

+

( )( ) ( )2 2

1

1 1

n n n n

n n

e e e ev

e e

+ − = + +

2

1 1

n n

n n

e ev

e e

− = + +

( )2v p p− = ( )1vp p−

( ) ( ) ( ){ }expT

M u A u n A n= + − = ( ) ( ){ }exp ln 1 ln 1u n n

v e v e++ − + =

( ) ( ){ }exp ln 1 ln 1v v

u n ne e

++ − + = 1

exp ln1

vu n

n

e

e

+ + = +

1

1

vu n

n

e

e

+ += +

11

11

v

u pe

p

p

p

+ −=

+ −

1

1

vup pe

p p

− += − +

( )vuq pe+ .

Page 7: 04_mathima

7

Πολυδιάστατη Κατανοµή Θεώρηµα

Αν 1 2, ,..., vX X X τ.δ., όπου η κατανοµή κάθε iX ανήκει στη µονοδιάστατη s-

παραµετρική Ε.Ο.Κ. µε ( )1 2, ,..., sT T T T= , τότε το ( )1 2, ,..., vX X X X= ανήκει σε

πολυδιάστατη s-παραµετρική Ε.Ο.Κ. µε διάνυσµα συναρτήσεως

( )1 2* *, *,..., *sT T T T= , όπου ( )1

* , 1,2,..., , 1,2,...,v

i i j

j

T T X i s j v=

= = =∑ (ν-διάστατη

κατανοµή µε s παραµέτρους.

*τ.δ. σηµαίνει ανεξάρτητες και ισόνοµες.

Απόδειξη:

Έστω ( ),jXf x θ η σ.π.π. της τ.µ. , 1,2,...,

jX i v= και επίσης ότι f EOK∈ , θα

έχουµε ότι ( ) ( ) ( ) ( ) ( ) ( )1

, exp , 1j X j

s

X i i j j j f

i

f x n T x h x x Sθ β θ θ=

= ∀ ∈

∑ .

Έστω ( ),Xf x θ η από κοινού σ.π.π. του τυχαίου διανύσµατος

( )1 2, ,..., vX X X X= .

1 2

...X X X Xv

f f f fS S S S= × × ×

( ) ( )1

, ,v

X X j

j

f x f xθ θ=

=∏ , λόγω της ανεξαρτησίας των 1 2, ,..., vX X X .

Και από την (1) ( ) ( ) ( ) ( ) ( )11

, expv s

X i i j j

ij

f x n T x h xθ β θ θ==

= =

∑∏

( ) ( ) ( ) ( )1 1 1

expvv s

v

i i j j

j i j

n T x h xβ θ θ= = =

=

∑ ∑ ∏

( ) ( ) ( ) ( )1 1 1

expvs v

v

i i j j

i j j

n T x h xβ θ θ= = =

=

∑ ∑ ∏

( ) ( ) ( ) ( )* * *

1

exps

i i

i

n T x h xβ θ θ=

∑ ,

όπου: ( ) ( )*

1

v

i i j

j

T x T x=

=∑ , ( ) ( )* v

β θ β θ= και ( ) ( )*

1

v

j

j

h x h x=

=∏ ,

άρα η κατανοµή του X είναι ν- διάστατη, s- παραµετρική και ανήκει στην

ΕΟΚ.

Page 8: 04_mathima

8

Παράδειγµα ( ) ( ) ( )2

1 2, , 0,θ θ θ µ σ= = ∈Θ = × ∞ℝ .

Έστω 1 2, ,..., vX X X τ.δ. από κανονική κατανοµή ( )2,N µ σ , µε 2,µ σ άγνωστα,

Ν.δ.ό. η από κοινού σ.π.π. του τ.δ. ( )1 2, ,..., vX X X X= ανήκει στη ν-διάστατη

και διπαραµετρική ΕΟΚ και να υπολογιστούν τα εξής:

2

1

v

j

j

E x=

∑ ,

1

v

j

j

E x=

∑ και 2

1 1

,v v

j j

j j

Cov x x= =

∑ ∑ .

Απόδειξη:

κατά τα γνωστά θα έχουµε: ( ) ( )2

2

22

1; , exp ,

22

j

j f

xf x x S

µµ σ

σπσ

− = − ∈ =

2 22

2 2 2

1 2 1exp ln

2 2 2 22

x xµ µσ

σ σ σπ

= − − + − =

2 22

2 2 2

2 1 1exp ln

2 2 2 2 2

x xµ µσ

σ σ σ π

= − − − =

22 2

2 2 2

1 1 1exp ln

2 2 2 2j jx x

µ µσ

σ σ σ π

= − + − −

,

22 2

2 2 2

1 1 1exp ln

2 2 2 2j jx x

µ µσ

σ σ σ π

= − + − −

( )2

1 2

1,

2n µ σ

σ= − , ( ) 2

1 j jT x x= , ( )2

2 2,n

µµ σ

σ= , ( )2 j jT x x= ,

( )2

2 2

2

1, exp ln

2 2

µβ µ σ σ

σ

= − −

, ( ) 1

2jh x

π=

Από το θεώρηµα έχουµε ότι:

( ) ( ) ( ) ( ) ( )2 * 2 2 * *

1

; , , exp ,s

X i i

i

f x n T x h xµ σ β µ σ µ σ=

=

∑ ,

όπου:

( ) ( )* 2 2, ,

v

β µ σ β µ σ = = 2

2

2

1exp ln

2 2

v

µσ

σ

− − =

22

2

1exp ln

2 2v v

µσ

σ − −

( ) ( )* 2

1 1

1 1

v v

j j

j j

T x T x x= =

= =∑ ∑ ,

( ) ( )*

2 2

1 1

v v

j j

j j

T x T x x= =

= =∑ ∑

και ( ) ( )*

1

v

j

j

h x h x=

= =∏ 1

2

v

π

=

( ) 22v

π−

Page 9: 04_mathima

9

Μετατροπή σε κανονική µορφή:

Κατά τα γνωστά:

( )2

1 12

1,

2n nµ σ

σ= − = ⇒

2

1

1

2nσ = −

( )2

2 22,n n

µµ σ

σ= = ⇒ 2

2nµ σ= ⇒ 2

1

1

2n

= − ⇒

2

12

n

nµ = −

( ) ( ) ( )( )2

1 2 1 2 1 2, , , ,A n n B n n n nµ σ= = 2

2

2

1ln

2 2v v

µσ

σ+ =

2

2

1

1

1

2 1 1ln

2 212

2

n

nv v

n

n

− + − = −

2

2

2

1

1

1

1 14ln

1 2 2

n

nv v

n

n

+ − =

− ( )

2

1 212

1

ln 24 2

vn n vn

n− − − =

( )2

21

1

ln 24 2

vn vn

n− − − =

2

2

2

2

1ln 2

1 2 24

2

vv

µσ

σσ

− − − − = −

2

4

2

2

1ln

4 2

vv

µσ

σσ

− =

( )2

2

2ln

4 2

v vµσ

σ+

2

1

v

j

j

E x=

=

∑ ( )*

1E T x = ( )1 2

1

,A n n

n

∂=

∂ ( )

2

21

1 1

ln 24 2

vn vn

n n

∂− − − = ∂

( )2

2

2

1 1

12

4 2 2

vn v

n n− − =

2

2

2

1 14 2

vn v

n n− =

2

2

2

22

11 2422

vv

µσ

σσ

− =

−−

2

42

4

1

v

v

µσ σ

σ

+ =

2 2v vµ σ+

1

v

j

j

E x=

=

∑ ( )*

2E T x = ( )1 2

2

,A n n

n

∂=

∂ ( )

2

21

2 1

ln 24 2

vn vn

n n

∂− − − = ∂

2

1

24

vn

n− =

2

12

vn

n− =

2

2

12

2

vµσ

σ

− =

Page 10: 04_mathima

10

2

1 1

,v v

j j

j j

Cov x x= =

=

∑ ∑ ( ) ( )* *

1 2,Cov T x T x = ( )2

1 2

1 2

,A n n

n n

∂=

∂ ∂

( )1 2

1 2

,A n n

n n

∂ ∂= ∂ ∂

2

1 12

vn

n n

∂− = ∂

2

2

12

vn

n=

2

2

2

12

2

vµσ

σ

=

2

4

12

4

vµσ

σ

= 2

4

1

2

vµσ

σ

= 4

2

2vµ σσ

= 22vµσ