PowerPoint - MONTECARpepca/lectures/pdf/montecarlo.pdf · 9]runrygqmsrgohgÖoh'lwrvwl xqÀnwhuk...

44
0RQWH&DUORPHWRG\ {-RVHI3HOLNgQ .69,0))8.3UDKD HPDLO-RVHI3HOLNDQ#PIIFXQLF] :::KWWSFJJPVPIIFXQLF]aSHSFD

Transcript of PowerPoint - MONTECARpepca/lectures/pdf/montecarlo.pdf · 9]runrygqmsrgohgÖoh'lwrvwl xqÀnwhuk...

0RQWH�&DUOR�PHWRG\

{�����������-RVHI�3HOLNgQ

.69,�0))�8.�3UDKD

H�PDLO��-RVHI�3HOLNDQ#PII�FXQL�F]

:::��KWWS���FJJ�PV�PII�FXQL�F]�aSHSFD�

0RQWH�&DUOR�LQWHJUDFH

( )I f x dx= ∫0

1

( )I fprim = ξ

2GKDGRYDQ¬�LQWHJUgO�

3ÎHGSRNODG� ( ) ( )f x L∈ 2 01,

-H�OL�ξ�QgKRGQk�³mVOR�V�GLVWULEXFm�5�������SDN�I�ξ��MHW]Y��SULPgUQm�RGKDG�LQWHJUgOX�

2GKDG�MH�QHVWUDQQ¬��QHERÔ�

( ) ( )E I f x dx Iprim = =∫0

1

3ULPgUQm�RGKDG�XU³��LQWHJUgOX

ξ

I�[�

I�ξ�

� �

5R]SW\O�RGKDGX

( ) ( ) ( )V I f x I dx f x dx Iprim prim= = − = −∫ ∫σ2 2

0

12

0

12

0ÀÎmWNHP�NYDOLW\�RGKDGX�MH�MHKR�UR]SW\O�QHER�VWDQGDUGQm�RGFK\OND��

3ÎL�Y¬SR³WX�MHGLQkKR�Y]RUNX�MH�UR]SW\O�Y¬VOHGNXSÎmOL§�YHON¬�

�SUR�QHVWUDQQ¬�RGKDG�

6HNXQGgUQm�RGKDG�XU³��LQWHJUgOX

ξ�

I���1

� �

ξ�

I�[��1I���1

� �

ξ�

I���1

� �

6HNXQGgUQm�RGKDG

( ) ( )I f x dx

f xN

dx Ii

N

ii

N

= = =∫ ∫∑ ∑= =0

1

0

1

1 1

5R]OR©HQm�LQWHJUgOX�QD�VRX³HW�1�³OHQÖ�

6HNXQGgUQm�RGKDG�LQWHJUgOX�

6HNXQGgUQm�RGKDG�MH�WDNk�QHVWUDQQ¬�

( )I IN

fi primi

N

ii

N

sec = == =∑ ∑

1 1

1 ξ

5R]SW\O�VHNXQGgUQmKR�RGKDGX

( )

( )

σ

σ

sec ... ...2

1

2

0

1

1

0

12

2

0

12

2

1

1 1

=

− =

= − =

=

=∑∫∫

Nf x dx dx I

Nf x dx

NI

N

ii

N

N

prim ����VWG��FK\ED�MH�√1�NUgW�PHQ§m������NRQYHUJHQFH���√1�

9]RUNRYgQm�SR�³gVWHFK

X �SÎL�Y¬EÀUX�PQR©LQ\�QH]gYLVO¬FK�Y]RUNÖ�VHVWHMQRX�KXVWRWRX�SUDYGÀSRGREQRVWL�GRFKg]m�NHVKOXNRYgQm� ]E\WH³QÀ�YHON¬�UR]SW\O�RGKDGX

·�Y]RUNRYgQm�SR�³gVWHFK���VWUDWLILHG�VDPSOLQJ��� SRWOD³XMH�VKOXNRYgQm� UHGXNXMH�UR]SW\O�RGKDGX

·�LQWHUYDO�VH�UR]GÀOm�QD�³gVWL��NWHUk�VH�RGKDGXMmVDPRVWDWQÀ

9]RUNRYgQm�SR�³gVWHFK

ξ�

I�[�

I�ξL�

� �ξ� ξ� ξ�

9]RUNRYgQm�SR�³gVWHFK

( ) ( )I f x dx f x dx IAi

N

ii

N

i

= = =∫ ∫∑ ∑= =0

1

1 1

5R]GÀOHQm�LQWHUYDOX�������QD�1�³gVWm�$L�

2GKDG�LQWHJUgOX�

( ) ( )I IN

f f Astrat i primi

N

ii

N

i i= = ∈= =∑ ∑

1 1

1 ξ ξ,

5R]SW\O�Y]RUNRYgQm�SR�³gVWHFK

( )

( )

σ

σ

strati

i i

Ai

N

ii

N

f xN

Ndx I

Nf x dx I

i

22

2

1

2

0

12

1

21

=

=

= − ≤

∫∑

∫ ∑

=

=sec

����UR]SW\O�QHPÖ©H�E¬W�YÀW§m�QH©����X�VHNXQGgUQmKR�RGKDGX�

5R]NODG�LQWHUYDOX�QD�³gVWL

X �XQLIRUPQm�UR]NODG�LQWHUYDOX������� SÎLUR]HQg�PHWRGD�SUR�]FHOD�QH]QgPRX�IXQNFL�I

X �]QgPH�OL�DOHVSRÈ�SÎLEOL©QÀ�SUÖEÀK�IXQNFH�I�VQD©mPH�VH�R�WDNRY¬�UR]NODG��DE\�E\O�UR]SW\OIXQNFH�QD�VXELQWHUYDOHFK�FR�QHMPHQ§m

X �UR]NODG�G�UR]PÀUQkKR�LQWHUYDOX�YHGH�QD�1G

Y¬SR³WÖ� rVSRUQÀM§m�PHWRGRX�MH�Y]RUNRYgQm��1�YÀ©m�

9]RUNRYgQm�SRGOH�GÖOH©LWRVWL

X �QÀNWHUk�³gVWL�Y]RUNRYDQkKR�LQWHUYDOX�MVRXGÖOH©LWÀM§m��SURWR©H�]GH�Pg�I�YÀW§m�KRGQRWX� Y]RUN\�]�WÀFKWR�REODVWm�PDMm�YÀW§m�YOLY�QD�Y¬VOHGHN

·�Y]RUNRYgQm�SRGOH�GÖOH©LWRVWL���LPSRUWDQFHVDPSOLQJ���XPLVÔXMH�Y]RUN\�SÎHGQRVWQÀ�GRWDNRY¬FK�REODVWm� Y]RUNRYgQm�MH�IRUPgOQÀ�Îm]HQR�IXQNFm�S�[�����KXVWRWRX�SUDYGÀSRGREQRVWL�QD�GDQkP�LQWHUYDOX

·�PHQ§m�UR]SW\O�SÎL�]DFKRYgQm�QHVWUDQQRVWL

9]RUNRYgQm�SRGOH�GÖOH©LWRVWL

ξ�

I�[�

� �

S�[�

ξ�ξ� ξ�ξ� ξ�

9]RUNRYgQm�SRGOH�GÖOH©LWRVWL

( ) ( )( ) ( )I f x dx

f xp x

p x dx= =∫ ∫0

1

0

1�SUDYD�RGKDGRYDQkKR�LQWHJUgOX�

0g�OL�QgKRGQg�SURPÀQQg�ξ�UR]GÀOHQm�V�KXVWRWRX�S�[��RGKDGXMHPH�LQWHJUgO�,�Y¬UD]HP�

( )( )I

fpimp =

ξξ �WHQWR�RGKDG�MH�QHVWUDQQ¬�

5R]SW\O�Y]RUNRYgQm�SRGOH�GÖOH©�

( )( ) ( )

( )( )

σ impf xp x

p x dx I

f xp x

dx I

22

0

12

2

0

12

=

− =

= −

3RNXG�VH�SUÖEÀK�KXVWRW\�S�[��SRGREg�LQWHJURYDQkIXQNFL�I�[���RGKDGXMHPH�LQWHJUgO�IXQNFH��NWHUgPg�PHQ§m�UR]SW\O�QH©�I�[��

9ODVWQRVWL�KXVWRW\�S�[�

X �S�[��≥����D��S�[��> ���WDP��NGH��I�[��≠��

X �∫�S�[��G[� ��

X �O]H�HIHNWLYQÀ�JHQHURYDW�Y]RUN\�V�GDQRXKXVWRWRX�SUDYGÀSRGREQRVWL� O]H�VSR³mWDW�SÎmVOX§QRX�GLVWULEX³Qm�IXQNFL�3�[��D]LQYHUWRYDW�ML��3���[��

( ) ( )P x p t dtx

= ∫0

3UDNWLFNg�LPSOHPHQWDFH

( ) ( )( ) ( ) ( )( )( )( )I f x dx f P t

dP tdt

dtf P t

p P tdt= = =∫ ∫∫ −

− −

−0

11

1 1

10

1

0

1

0mVWR�SÎmPkKR�Y¬EÀUX�QgKRGQk�SURPÀQQk�V�KXVWRWRXSUDYGÀSRGREQRVWL�S�[��EHUHPH��τ��V�URYQRPÀUQ¬PUR]GÀOHQmP�SUDYGÀSRGREQRVWL�D�WUDQVIRUPXMHPH�ML�

2GKDG�Pg�WHG\�WYDU�

( )ξ τ= −P 1

( )( )( )( )I

f P

p Pimp =−

1

1

τ

τ

.RPELQRYDQk�RGKDG\

I�[�

� �

S��[� S��[�

.RPELQDFH�QÀNROLND�RGKDGÖ

( ) ( )( )I w

fpcomb i i

i

i ii

n

==∑ ξ

ξξ

1

3ÎHGSRNOgGgPH�Q�QgKRGQ¬FK�SURPÀQQ¬FK�ξ������ξQV�KXVWRWDPL�SUDYGÀSRGREQRVWL�S��[������SQ�[��.RPELQRYDQ¬�RGKDG�LQWHJUgOX�EXGH�PmW�WYDU�

NGH�ZL�[��MVRX�QH]gSRUQk�YgKRYk�IXQNFH�

1HVWUDQQRVW�NRPELQ��RGKDGX

( ) ( ) ( )( ) ( )

( ) ( ) ( )

E I w xf x

p xp x dx

w x f x dx f x dx

comb i ii

i ii i i

i

n

ii

n

=

=

=

∫∑

∑∫ ∫

=

=

0

1

1

10

1

0

1

3RGPmQND�SURYgKRYk�IXQNFH�

( )∀ ∑ ==

x w xii

n

: 11

5R]SW\O�NRPELQRYDQkKRRGKDGX

( ) ( )( ) ( )

( ) ( )( ) ( )

( )( ) ( ) ( ) ( )

σcomb

i ii

i ii i i

i ii

i ii i i

i

n

i

ii

n

ii

n

w xf x

p xp x dx

w xf xp x

p x dx

w xp x

f x dx w x f x dx

2

2

0

1

0

1 21

2

10

1

0

1 2

1

=

=

=

∫∑

∑∫ ∫∑

=

= =

$ULWPHWLFN¬�SUÖPÀU��PD[LPXP

( )w xni = 1 ( )

( )In

fpaverage

i

i ii

n

==∑1

1

ξξ

( ) ( ) ( ){ }w x

p x p xi

i j jpro

jinde=

=

1

0

max

( ) ( ){ }( ) ( )( )I p p

fpi i j j i

i

i ii

n

max max ? := ==∑ ξ ξ

ξξ

1

0

9\URYQDQg�KHXULVWLND

( ) ( )( )

w xp x

p xi

i

jj

n=

=∑ 1

( )( )

If

pbal

i

j ij

ni

n

=

== ∑∑ ξ

ξ1

1

( )( )

( )( )

( )σbal

ii

ni

jj

ni

nf x

p xdx

p x

p xf x dx2

2

1 10

12

10

1

= −

= =

=∑ ∑∫∑∫

( )σ σcomb bal n I2 2 1 21≥ − − ⋅

0RFQLQQg�KHXULVWLND

( ) ( )( )

w xp x

p xi

i

jj

n=

=∑β

β1

( )( )

( )Ip

pfpower

i i

j ij

ni

n

i=−

== ∑∑

β

β

ξ

ξξ

1

11

=REHFQÀQm�

β ������Y\URYQDQg���β �∞����PD[LPgOQm�KHXULVWLND�

7UDQVIRUPDFH�LQWHJUDQGX

( ) ( ) ( )I f x dx w x f x dxii

n

= = ⋅∫ ∫∑=0

1

0

1

1

( )( )( )( ) ( )( )I

w P t

p P tf P t dt

i i

i ii

n

i=−

−=

−∫∑1

10

1

1

1

,QWHUSUHWDFH�NRPELQD³QmKR�RGKDGX�MDNRWUDQVIRUPDFH�LQWHJURYDQk�IXQNFH�

.RPELQDFH�RGKDGÖ�SRGOH�GÖOH©LWRVWL�

-HGHQ�³OHQ�NRPELQ��RGKDGX

I�[�

� �

S��[� S��[�

$ULWPHWLFN¬�SUÖPÀU

� �

( )( )( )( )0 5

11

1 11

.f P t

p P t

0D[LPXP

� �

( ) ( )[ ] ( )( )( )( )p t p t

f P t

p P t1 2

11

1 11

1 0>−

−? :

9\URYQDQg�KHXULVWLND

� �

( )( )( )( ) ( )( )f P t

p P t p P t

11

1 11

2 11

− −+

0RFQLQQg�KHXULVWLND�SUR�β �

� �

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

p P t f P t

p P t p P t

1 11

11

12

11

22

11

− −

− −

+

ÍmGmFm�IXQNFH

( ) ( ) ( )[ ] ( )

( ) ( )[ ] ( ) ( )[ ]

I f x dx f x g x dx g x dx

f x g x dx J f x g x J dx

= = − + =

= − + = − +

∫ ∫ ∫

∫ ∫

0

1

0

1

0

1

0

1

0

1

( ) ( )I f g Jcon = − +ξ ξ

)XQNFH�J�[���NWHUg�DSUR[LPXMH�LQWHJUDQG�DGRNg©HPH�ML�DQDO\WLFN\�]LQWHJURYDW�

1HVWUDQQ¬�RGKDG�

7UDQVIRUPDFH�ÎmGmFm�IXQNFm

I�[�

� �

J�[�

I�[��J�[�

ÍH§HQm�LQWHJUgOQmFK�URYQLF

( ) ( ) ( ) ( )f x g x K x y f y dy= + ⋅∫ ,0

1

)UHGKROPRYD�LQWHJUgOQm�URYQLFH�GUXKkKR�W\SX�

QH]QgPg ]QgPk�IXQNFH

·�PHWRG\�NRQH³Q¬FK�SUYNÖ��Y¬SR³HW�FHOk�IXQNFH�

·�PHWRG\�0RQWH�&DUOR��ORNgOQm�Y¬SR³HW�

5HNXU]LYQm�0RQWH�&DUOR�RGKDG

( ) ( ) ( )( ) ( )

( ) ( )( ) ( ) ( )

( ) ( )

( ) ( )( ) ( ) ( )

( )( )

( ) ( )

f x g xK xp

f

g xK xp

gKp

f

g xK xp

gK xp

Kp

g

r r

r

= + ⋅ =

= + ⋅ + ⋅

= + + +

,

, ,

, , ,..

ξξ

ξ

ξξ

ξξ ξ

ξξ

ξξ

ξξ

ξξ ξ

ξξ

1

1 11

1

1 11

1 2

2 22

1

1 11

1

1 1

1 2

2 22

3UDYRX�VWUDQX�URYQLFH�RGKDGXML�VWRFKDVWLFN\V�ÎmGmFmPL�KXVWRWDPL�SUDYGÀSRGREQRVWL�SL�[��

5HNXU]LYQm�0RQWH�&DUOR�RGKDG

( ) ( )( ) ( )f x

K

pg x

rj j

j jj

i

ii=

∞=−

==∏∑ ξ ξ

ξξ ξ1

100

,,

I J I= + T

^ξ���ξ���ξ�������`�VH�QD]¬Yg�0DUNRYRYVN¬�ÎHWÀ]HF�MHVWOL©H�SUDYGÀSRGREQRVW�SL�[��]gYLVm�SRX]H�QD�ξL��

�VSRUQÀM§m�IXQNFLRQgOQm�]gSLV�

ÍH§HQm��1HXPDQQRYD�ÎDGD�� I J J J= + + +T T2 ..

5XVNg�UXOHWD

X �SÎL�RGKDGX�QHNRQH³Qk�1HXPDQQRY\�ÎDG\�VHPÖ©H�VSR³mWDW�MHQ�NRQH³Q¬�³gVWH³Q¬�VRX³HW� SHYQÀ�GDQg�GkOND�SRVORXSQRVWL�]DYgGm�GR�RGKDGXV\VWHPDWLFNRX�FK\EX

·�YKRGQÀM§m�MH�PHWRGD�QgKRGQkKR�XNRQ³HQmY¬SR³WX��W]Y��UXVNg�UXOHWD� RGKDG�]ÖVWgYg�QHVWUDQQ¬

·�WHRUHWLFN\�VH�Gg�WHQWR�SRVWXS�DSOLNRYDW�L�QDY¬SR³HW�MHGQRWOLYkKR�LQWHJUgOX

5XVNg�UXOHWD�SUR�MHGHQ�LQWHJUgO

( ) ( )I f x dx f dt PPtP

P

= = < ≤∫ ∫0

11

0

0 1

( )If P

RussP P pro

jinak=

<

1

0

ξ ξ

1HVWUDQQ¬�RGKDG�V�MHGQmP�QgKRGQ¬P�Y]RUNHP�

7UDQVIRUPDFH�LQWHJUgOX�

^ξ���ξ�������ξN`�MH�NRQH³Qg�QgKRGQg�SURFKg]ND�SURWR©H�RGKDG��⟨I�ξN�⟩ ���.D©G¬�Y]RUHN�ξL�MH�Y\EmUgQ�V�SUDYGÀSRGREQRVWm�3L

D�V�KXVWRWRX�UR]GÀOHQm�SL�[��3RNXG�QgKRGQg�SURPÀQQg��τL���!�3L����FHO¬�SURFHV�VH]DVWDYm��MLQDN�VH�Y\JHQHUXMH�ξL����D�SÎLGg�VH�GDO§m�³OHQ��

5XVNg�UXOHWD�SUR�LQWHJU��URYQLFH

( ) ( )( ) ( )f x

K

P pg x

Russ rj j

j j jj

i

i

k

i,

,,=

=−

==∏∑ ξ ξ

ξξ ξ1

100

( )K x y dy,0

1

1∫ <

9ROED�SUDYGÀSRGREQRVWm

( ) ( ) ( )P K y dy p x

K xPi i ii

i= =−

−∫ ξξ

1

0

11, ,,

7HKG\�O]H�MgGUR�.�SRX©mW�NH�NRVWUXNFL�W]Y�VXENULWLFNkKR�UR]GÀOHQm�SUDYGÀSRGREQRVWL�

9H�I\]LNgOQmFK�DSOLNDFmFK�MH�³DVWR�

2GKDG�SDN�EXGH�PmW�WYDU� ( ) ( )f x gsubcrit i

i

k

==∑ ξ

1

2GKDG�SÎm§Wm�XGgORVWL

( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( ) ( )

f x g x h x

h x K x y f y dy

K x y g y dy K x y h y dy

= +

= ⋅ =

= ⋅ + ⋅

∫ ∫

,

, ,

0

1

0

1

0

1

3ÎHGFKR]m�RGKDG�PmYg�YHON¬�UR]SW\O��PgOR�V³mWDQFÖMH�QHQXORY¬FK���/HS§m�Y¬VOHGN\�GgYg�PHWRGDRGKDGXMmFm�³OHQ�J�[��R�MHGHQ�VWXSHÈ�SÎHVQÀML�

2GKDG�SÎm§Wm�XGgORVWL

X �SUYQm�LQWHJUgO�VH�RGKDGXMH�]D�SRPRFL�SUDYGÀ�SRGREQRVWL�V�KXVWRWRX�SRGREQRX�IXQNFL�J�[�� QgKRGQg�SURPÀQQg��ζL��V�KXVWRWRX��SL�[�

X �GUXK¬�LQWHJUgO�VH�RGKDGXMH�SRPRFmVXENULWLFNk�KXVWRW\�SUDYGÀSRGREQRVWL�MgGUD�.� QgKRGQg�SURPÀQQg��ξL��V�KXVWRWRX��.�ξL���[��3L

( ) ( ) ( )( ) ( )h x

K x gp

hnextev nextev

= +,ζ ζ

ζξ1 1

1 11

2GKDG�SÎm§Wm�XGgORVWL

2GKDG�LQWHJUgOQm�URYQLFH�

2GKDG�IXQNFH�K�

( ) ( ) ( )( )h x

K gpnextev

i i i

i ii

k

= −

=∑ ξ ζ ζ

ζ1

1

,

( ) ( ) ( ) ( )( )f x g x

K gpnextev

i i i

i ii

k

= + −

=∑ ξ ζ ζ

ζ1

1

,

.RQHF

'DO§m�LQIRUPDFH�

Q (��/DIRUWXQH��0DWKHPDWLFDO�0RGHOV�DQG

0RQWH�&DUOR�$OJRULWKPV�IRU�3K\VLFDOO\�%DVHG

5HQGHULQJ��3K'�WKHVLV��.8�/HXYHQ�������

Q 0��.DORV��3��:KLWORFN��0RQWH�&DUOR�0HWKRGV�-RKQ�:LOH\��6RQV��������������

Q $��*ODVVQHU��3ULQFLSOHV�RI�'LJLWDO�,PDJH

6\QWKHVLV��0RUJDQ�.DXIPDQQ���������������