Trưng ði hc Nông nghip Hà Ni – Giáo trình Vn trù hc
………………………………..........150
q
00
π
0
+ q
10
π
1
+ q
20
π
2
+... = 0,
q
01
π
0
+ q
11
π
1
+ q
21
π
2
+... = 0,
q
02
π
0
+ q
12
π
1
+ q
22
π
2
+... = 0,
....
Do tính ch
t
ñ$
c bi
t, nh
ư
ñ
ã phân tích
*
trên, c
a ma tr
n c
ư
ng
ñ
Q c
a quá trình
sinh
t
#
, h
trên
ñư
c vi
t m
t cách t
ư
ng minh h
ơ
n nh
ư
sau:
−λ
0
π
0
+
µ
1
π
1
+... = 0,
λ
0
π
0
(
λ
1
+
µ
1
)
π
1
+
µ
2
π
2
+... = 0,
λ
1
π
1
(
λ
2
+
µ
2
)
π
2
+
µ
3
π
3
+... = 0,
...
T
"
ñ
ây
ñ.
ng tìm
ñư
c
π
n+1
= (
λ
n
/
µ
n+1
)
π
n
,
n = 1, 2, 3,...
ñ
ñ
i t
i công th
c tính
π
i
,
i.
π
1
= (
λ
0
/
µ
1
)
π
0
,
π
2
= (
λ
1
/
µ
2
)
π
1
= (
λ
1
λ
0
/
µ
2
µ
1
)
π
0
,
π
3
= (
λ
2
/
µ
3
)
π
2
= (
λ
2
λ
1
/
µ
3
µ
2
)
π
1
= (
λ
2
λ
1
λ
0
/
µ
3
µ
2
µ
1
)
π
0
,
...
π
n+1
= (
λ
n
/
µ
n+1
)
π
n
=... = (
λ
n
λ
n1...
λ
0
/
µ
n+1
µ
n...
µ
1
)
π
0
,
...
V
i
ñ
i
u ki
n
i
i 0
1,
=
π =
cu
i cùng ta có:
π
0
= 1/(1 +
=0k
(λ
k
λ
k1...
λ
0
/µ
k+1
µ
k...
µ
1
)).
ð$
c bi
t khi µ
n
= 0, n thì quá trình sinht
#
tr
*
thành quá trình sinh thu
n khi
t
(
pure birth process
). Quá trình sinh thu
n khi
t v
i λ
n
= λ
quá trình Poátxông v
i
tham s
λ
.
d 2: Gi s# dòng khách hàng ñn mua * mt văn phòng n vi M qu y
phc v dòng Poátxông vi tham s λ = 6 khách hàng/1 phút (ñiu này cũng
nghĩa khách hàng ñn phòng bán vi c thi ñim ñn tuân theo lut phân phi
mũ vi tham s λ = 6).
Ngoài ra, còn bit nguyên tc phc v là FCFS (First come first served) và thi gian
phc v ti m&i qu y có lut phân phi mũ vi kì vng 1/3 (phút).
C n tr li hai câu hi sau ñây:
S qu y hàng ti thiu là bao nhiêu ñ hàng ch không tr* nên dài vô hn?
Gi s# N
t
s khách hàng ñang ch hay ñang ñưc phc v ti thi ñim t. Chn
M = 4 mt khách hàng s, ch ñ ñưc phc v nu N
t
4, ch vi xác sut 0,5 nu
N
t
= 5 và s, b ñi nu N
t
= 6. Hãy xác ñnh phân phi d"ng ca quá trình này?
Trưc ht, trong ví d này chúng ta mt quá trình sinht# vi không gian trng
thái S = {S
0
, S
1
, S
2
,..., S
n
,...}, trong ñó S
n
trng thái trong văn phòng có n khách hàng.
Các cưng ñ chuyn là λ
k
= 6 vi k = 0, 1, 2,... còn µ
k
= 3k vi k M và µ
k
= 3M vi
k > M. ðiu này do bin cc tiu ca các bin ngu nhiên vi pn phi mũ ñc
lp cũng có phân phi mũ vi tham s bng tng các tham s ca các pn phi
mũ tương ng. Tht vy, gi s# X = Min {X
1
, X
2
} vi X
1
và X
2
tuân theo phân
phi mũ ñc lp vi c tham s µ
1
và µ
2
, th thì P(X t) = P(X
1
t) P(X
2
t) =
1 2 1 2
t t ( )t
e e e
µ −µ µ +µ
=
. Do ñó X phân phi mũ vi tham s µ = µ
1
+ µ
2
.
Do λ
k
/µ
k+1
= 6/3M < 1 (khi k M) nên vi M 3 thì:
=0k
k k 1 0 k 1 k 1
( ... / ... )
+
λ λ λ µ µ µ
< .
B*i vy hàng ñi s, không dài vô hn (nu trái li, khi chu&i phân thì
0 1 2
... 0,
nên s khách trong hàng ñi s, d n ti mt s h(u hn khi t
vi xác sut bng 0).
Trong câu hi th hai, ta λ
0
= λ
1
= λ
2
= λ
3
= λ
4
= 6, λ
5
= 3. Theo công thc tính
5
0
k 0
1/(1
=
π = +
k k 1 0 k 1 k 1
( ... / ... ))
+
λ λ λ µ µ µ
ta có ngay π
0
= 12/89. T
"
ñ
ó tính ra
π
1
= 24/89,
π
2
= 24/89, π
3
= 16/89, π
4
= 8/89, π
5
= 4/89 và π
6
= 1/89.
3. MÔ PH$NG XÍCH MARKOV
3.1. Mô ph%ng xích Markov th(i gian r(i rc
Phương pháp 1
Xích Markov ri rc thu n nht còn th ñưc hiu X
0
, X
1
, X
2
,... Gi s#
không gian trng thái S gm h(u hn trng thái: S = {0, 1, 2,..., N} ma trn xác
sut chuyn trng thái ñã ñưc bit là P = [p
ij
]
N
×
N
. Chúng ta s, phng xích Markov
ri rc và thu n nht thông qua ví d ñã trình bày * các mc 1.2 và 2.1 ca chương này.
Ta có phân phi ban ñ u là:
X
0
1 2 3
Π
(0)
π
1(0)
= 0,2 π
2(0)
= 0,5 π
3(0)
= 0,3
Trưng ði hc Nông nghip Hà Ni – Giáo trình Vn trù hc
………………………………..........152
ð phng X
0
ta áp dng phương pháp phng phân phi ri rc ñã hc *
chương III. Trên y tính, ta phát sinh ra mt s ngu nhiên r = RANDOM[0,1) theo
lut phân phi ñu U[0,1) trong [0,1). Nu r 0,2 ta ly X
0
= 1; nu 0,2 < r 0,7 thì ta
ly X
0
= 2 ; còn nu r > 0,7 thì ñ$t X
0
= 3. Căn c kt qu phng X
0
, ta phng
X
1
da trên ma trn xác sut chuyn trng thái:
P =
083,0
07,0
8,0
067,0
9,0
1,0
85,0
03,0
1,0
.
Gi
s
#
ñ
ã bi
t X
0
= 2, lúc
ñ
ó ta c
n ph
ng bi
n ng
u nhiên X
1
c
ă
n c
phân ph
i
sau:
X
1
1 2 3
Xác sut tương ng
p
21
= 0,07 p
22
= 0,9 p
23
= 0,03
ð
i
u này có th
ñư
c th
c hi
n t
ươ
ng t
nh
ư
khi mô ph
ng X
0
. C
n chú ý r
ng,
trong hàng th
hai c
a b
ng trên ta phân ph
i xác su
t có
ñ
i
u ki
n c
a X
1
v
i
ñ
i
u
ki
n X
0
= 2. Các b
ư
c ti
p theo ph
ng X
2
, X
3
,...
ñư
c ti
n hành t
ươ
ng t
(cho t
i
X
500
ch
-
ng h
n).
L
$
p l
i quy trình y b
t
ñ
u t
"
X
0
cho m
t s
b
ư
c l
$
p L
ñ
l
n (ch
-
ng h
n 1000
l
n), ta s
,
có m
t b
1000 s
li
u cho X
500
. T
"
ñ
ó, có th
tìm
ñư
c b
ng phân ph
i t
n
su
t (còn g
i xác su
t th
c nghi
m) c
a X
500
qua thí nghi
m ph
ng trên
ñ
ây
ñ
i
v
i X
500
. Nh
ư
v
y, ta tìm
ñư
c véc t
ơ
phân ph
i (xác su
t th
c nghi
m)
Π
(500)
. Cu
i
cùng, chúng ta có k
t qu
tìm g
n
ñ
úng phân ph
i d
"
ng là:
Π
Π
(500)
.
Chú ý:
Trong d
trên
ñ
ây, ta th
y th
dùng ph
ng
ñ
tìm phân ph
i d
"
ng. Tuy
nhiên, m
c
ñ
ích ch
y
u c
a ph
ươ
ng pháp 1 là nh
m mô ph
ng các xích Markov r
i r
c
thu
n nh
t, là các quá trình có th
x
y ra trong các h
th
ng ph
c t
p.
Khi không gian tr
ng thái S g
m m
t s
l
n các tr
ng thái thì ph
ươ
ng pháp
ph
ng trên yêu c
u th
i gian ch
y máy tính khá l
n.
ð
kh
c ph
c
ñ
i
u này, chúng ta
xem xét ph
ươ
ng pháp 2 sau
ñ
ây.
Phương pháp 2
Xét m
t h
th
ng k
ĩ
thu
t
ñư
c bi
u di
.
n b
*
i xích Markov r
i r
c thu
n nh
t {X
t
}, t
= 0, 1, 2,... v
i không gian tr
ng thái S N tr
ng thái (N khá l
n) ma tr
n chuy
n
tr
ng thái P = [p
ij
]
N
×
N
. Xét th
i
ñ
i
m n, t
i th
i
ñ
i
m này gi
s
#
ñ
ã mô ph
ng
ñư
c X
n
=
s. Ta s
,
mô ph
ng th
i gian T
n
là th
i gian t
i l
n nh
y ti
p theo s
m nh
t mà X
t+Tn
s.
Do xích Markov r
i r
c n T
n
ch
th
nh
n các giá tr
1, 2,...
ð$
t p = p
ss
, d
.
th
y
T
n
phân ph
i hình h
c nh
ư
sau:
T
n
1 2 ... k ...
c sut
tương ng 1p (1p)p ... (1p)p
k1
...
ph
ng phân ph
i này ta tìm
ñư
c giá tr
T
n
.
Còn X
n+Tn
phân ph
i xác su
t
nh
ư
sau:
X
n+Tn
1 2 ... s ... N
Xác sut
tương ng p
s1
/(1p
ss
) p
s2
/(1p
ss
) ... 0 ... p
sN
/(1p
ss
)
Cách ph
ng này s
,
ti
t ki
m h
ơ
n th
i gian ch
y máy tính (khi N khá l
n),
nh
ư
ng vi
c l
p trình s
,
ph
c t
p h
ơ
n ít nhi
u.
Xét ví d
nh
ư
ñ
ã trình bày trên, n
u dùng ph
ươ
ng pháp 2, m
t cách hoàn toàn t
ươ
ng
t
, chúng ta c
ũ
ng tìm
ñư
c phân ph
i d
"
ng
Π
(*)
Π
(500)
.
3.2. Mô ph%ng xích Markov th(i gian liên tc
Xét xích Markov thi gian liên tc {X(t)}
t
[0,
)
. Gi s# rng xích ñi vào trng thái i
ti thi ñim o ñó, ch-ng hn thi ñim 0 không ri khi trng thái y cho ñn
thi ñim s. Lúc ñó, do tính không nhca quá trình Markov, xác sut ñ xích vn
tip tc * nguyên trng thái ñó cho ti thi ñim (t + s) s, là:
P{(T
i
> s + t )/(Ti > s)} = P{T
i
> t}
trong ñó T
i
thi gian quá trình d"ng li * trng thái i. D. thy, nu T
i
phân
phi mũ vi hàm phân phi F(T
i
< τ) = 1 - e
−λτ
thì ñ-ng thc trên ñưc tho mãn. ðiu
ngưc li cũng có th chng minh ñưc. Vy T
i
có phân phi mũ.
T" nhn xét trên, ta th ñưa ra mt ñnh nghĩa khác cho xích Markov thi gian
liên tc. Xích Markov thi gian liên tc là mt quá trình ngu nhiên có các tính cht sau
m&i khi nó ñi vào trng thái i:
Lưng thi gian T
i
xích d"ng li ti trng thái i trưc khi chuyn sang trng
thái khác là mt bin ngu nhiên vi phân phi mũ có tham s v
i
(hay có kì vng 1/v
i
).
Mt khi quá trình ri khi trng thái i, nó s, ñi vào trng thái j nào ñó (ñc lp vi
T
i
) vi các xác sut p
ij
tho mãn
 

= =
.
Vy ñphng xích Markov thi gian liên tc, chúng ta c n mô phng dãy τ
0
, τ
1
,
τ
2
,... (các lưng thi gian τ
r
xích d"ng li ti trng thái J
r
trưc khi nó chuyn sang trng
thái khác) dãy J
0
, J
1
, J
2
,... (các trng thái xích chuyn ñn). ð phát sinh τ
r,
như
trên ñã nói, ta c n bit tham s v
Jr
ca phân phi mũ tương ng. Còn ñ phát sinh trng
thái xích Markov chuyn ñn J
r
r, chúng ta có bng phân phi xác sut sau:
Trng thái ñn 1 2 ... i ... N
Trưng ði hc Nông nghip Hà Ni – Giáo trình Vn trù hc
………………………………..........154
c sut tương ng p
i1
p
i2
... 0 ... p
iN
Trong bng trên, i =J
r
-
1
là trng thái ca xích ti bưc r 1 (vi c xác sut p
ij
tho
mãn
 

= =
).
ð thc hin mô phng xích Markov thi gian liên tc, có th s# dng s liu ca
d ñã xét trong mc 2.4 hay 2.5.
BÀI TP CHƯƠNG V
1. Ch s tiêu th ñin là mt lưng ngu nhiên phân phi ti thi ñim ban ñ u
như sau:
X
0
Dưi 50 s 50 ti 100 100 ti 150 Trên 150
T l % 5% 40% 40% 15%
Bit ma trn xác sut chuyn trng thái là:
P =
05,0
02,0
05,0
85,0
05,0
03,0
85,0
10,0
10,0
90,0
08,0
05,0
80,0
05,0
02,0
0
Hãy gi
i thích ý ngh
ĩ
a c
a ma tr
n P.
Tìm phân ph
i d
"
ng c
a xích Markov th
i gian r
i r
c trên
ñ
ây cho bi
t ý
ngh
ĩ
a c
a k
t qu
thu
ñư
c.
2.
M
t ch
trang tr
i tr
ng hoa hàng n
ă
m th
c hi
n phân tích thành ph
n
ñ
t c
a
trang tr
i. K
t qu
phân tích
ñư
a ra
ñ
t thu
c vào m
t trong ba tr
ng thái: t
t, bình
th
ư
ng x
u. Các kh
o t th
ng cho bi
t các ma tr
n xác su
t chuy
n tr
ng thái
(sau m
t n
ă
m) trong các tr
ư
ng h
p không bón phân và có bón phân (h
(
u c
ơ
) nh
ư
sau:
1
0,2 0,5 0,3
P 0 0,5 0,5
0 0 1
=
2
0,3 0,6 0,1
P 0,1 0,6 0,3
0,05 0,4 0,55
=
.
Các ma tr
n l
i nhu
n/n
ă
m (
ñơ
n v
tính là 10 ngàn USD) t
ươ
ng
ng v
i các ma tr
n
xác su
t chuy
n tr
ng thái trên là:
1
7 6 3
R 0 5 1
0 0 1
=
2
6 5 1
R 7 4 0
6 3 2
=
.