Đề 24:
Đánh giá hệ số sẵn sàng của hệ thống có phục hồi
1. Tính hệ số sẵn sàng :
Tlv
Kss =
(Tlv +Tph)
2. Vẽ đường cong làm việc của hệ thống
Cho biết :
λ
, 1/giây 0.0001 0.0002 0.0003
µ
,1/ giây 0.001 0.002 0.003
Phần 1 . Giới thiệu phương pháp mô hình hoá hệ ngẫu nhiên

Ở giai đoạn nghiên cứu và thiết kế hệ thống bằng máy tính ,
người ta thường ứng dụng phương pháp thử thống kê - còn có
tên là phương pháp Monte-Carlo .Phương pháp này sử dụng các
số ngẫu nhiên có phân bố xác suất nhất định.Những mô hình dựa
trên phương pháp này được gọi là mô hình xác suất.
Thực chất của phương pháp mô hình xác suất là xây dựng
trên máy tính hệ thống S với các quan hệ nội tại của nó. Đầu vào
của hệ có tác động mang tính ngẫu nhiên như số lượng các sự
kiện xảy, thời gian giữa các sự kiện hoặc tác động của môi trường
xung quanh E. Trên cơ sở đó phân tích các tín hiệu đầu ra người
ta nhận được dáng điệu phản ứng của hệ thống. Đó là phương
pháp mô phỏng. Mỗi một lần thực hiện phép thử người ta thu
được một lời giải chứa đựng thông tin về dáng điệu của hệ thống
S. Nếu số phép thử N đủ lớn thì kết quả thu được bằng cách lấy
trung bình theo xác xuất sẽ ổn định và đạt độ chính xác cần thiết.
Trong hệ ngẫu nhiên có nhiều biến ngẫu nhiên. Các đặc trưng
quan trọng của biến ngẫu nhiên là hàm mật độ xác xuất, hàm

với a≤x≤ b
0 với giá trị còn lại của x
1
b
−
a
f(x)=
với a≤x≤b
0 với giá trị còn lại của x
trường hợp riêng khi a=0 , b=1 ta có phân bố đều trong (0,1)
*Phân bố mũ expo (β)

Có hàm phân bố và hàm mật độ là:
1
a
F(x)=
e
− x
/ a
với x>=0
0 với x<=0
1-
e
−
f(x)=
0
x / a
với x>=0
với x<=0
2

0
=0.2152
X
0
2
=0.04631104
⇒

X
1
=0.6311
X
1
2
=0.39828721
⇒
X
2
=0.8287
Thuật toán nhân

Xn+1= phần lẻ của λ*xn trong đó λ là số có dạng ( 8t-3) hoặc
(8t+3) trong đó t là số nguyên dương
X
0
=0.37843
Lấy λ=37 thì
X
1
=0.00191

≈
(0,1)
- Lấy X=a+(b-a)U
Phần2. Phân tích thuật giải

1.Thuật giải

1.Gieo số ngẫu nhiên Ui ~ U(0,1) để tính thời gian làm việc
2.Thời gian làm việc :
Tlv = -
λ

ln(Ui), Ui ~ U(0,1)
3.Gieo số ngẫu nhiên Ui ~ U(0,1) để tính thời gian phục hồi
4.Thời gian phục hồi :
Tph=-
µ

ln(Ui) , Ui ~ U(0,1)
5.Đếm thời gian làm việc : Tlàm việc =
∑
T
lvi
6.Đếm thời gian phục hồi: T phục hồi =
∑
T
phi
7.Tính hệ số sẵn sàng :
Tlv
Kss =

+


In kÕt qu¶
KÕt thóc
Phần 3. Lập trình Visual Basic
Option Explicit
Const e = 2.71828
Dim maxx, maxy, d, dt, t(10), w(10), n, f, g, Tlv, Tph, c, b As Double
Dim kx, Kss, ex, ky, Plv(10), Clv(10000) As Double
Dim muy, lamda, Tphuchoi, Tlamviec, Tongthoigiansua,
Tongthoigiankiemtra, u As Double
Dim Thoigianlamviec(1 To 10) As Double
Dim Thoigianphuchoi(1 To 10) As Double
Dim i As Long 'su dung trong cac vong lap
Dim i1 As Long 'nt
Dim i2 As Long 'nt
Dim biendem As Long ' biendem co tac dung lien ket mang cac
'thong so da tinh toan de truyen so lieu vao chuong trinh
'duoc thuc hien trong cac timer
Dim a As Double
'bien darun giup kiem tra chuong trinh da chay lan nao chua
Dim Dacapnhatthongso As Boolean
Dim darun As Boolean
Private Sub Cmdexit_Click()
Dim a As Double
a = MsgBox("Do you want to exit?", vbQuestion + vbYesNo, " Exit?")
If a = 6 Then
End
End If

Label16.Caption = ""
'***********************************************
darun = True
If Not Dacapnhatthongso Then
capnhatthongso
End If
Khoitaotudau
Tinhtoan
Timer1.Enabled = True
End Sub
Private Sub Command1_Click()
End Sub
Private Sub Form_Load()
Khoitaotudau
End Sub
Private Sub khoitaothongso()
Text1.Text = "0.0001"
Text2.Text = "0.001"
lamda = 0.0001
muy = 0.001
Dacapnhatthongso = False
darun = False
End Sub
Private Sub Khoitaotudau()
Picture3.Cls
Timer1.Enabled = False
Timer2.Enabled = False
Timer3.Enabled = False
'Khoi dong lai vi tri thanh progress
Picphantram.Move -Picprogress.Width, 0, Picprogress.Width,

Tongthoigiansua = 0
Tongthoigiankiemtra = 0
For i = 1 To 10
Tongthoigiansua = Tongthoigiansua + Thoigianphuchoi(i)
Tongthoigiankiemtra = Tongthoigiankiemtra + Thoigianlamviec(i)
Next i
c = (Tongthoigiansua + Tongthoigiankiemtra) / 1000
'Tinh he so san sang
Kss = (Tongthoigiankiemtra) / (Tongthoigiansua +
Tongthoigiankiemtra)
End Sub
Private Sub Image2_Click()
End Sub
Private Sub mnuchay_Click()
'chay moi thi cac thong so tinh tong thoi gian da thay doi, vi vay ma
'cac nhan thong so cu phai duoc xoa di
Label15.Caption = ""
Label16.Caption = ""
'***********************************************
darun = True
If Not Dacapnhatthongso Then
capnhatthongso
End If
Khoitaotudau
Tinhtoan
Timer1.Enabled = True
End Sub
Private Sub mnudanhsachnho_Click()
Dim an As Integer
an = MsgBox("Nguyen Quoc Thang" & Chr(13) & "Nguyen Van

Khoitaotudau
End Sub
Private Sub Text1_Change()
Dacapnhatthongso = False
darun = False
End Sub
Private Sub Text2_Change()
Dacapnhatthongso = False
darun = False
End Sub
Private Sub Timer1_Timer()
'tang bien dem len 1 de goi so lieu moi trong mang so lieu da tinh
biendem = biendem + 1
'hien thi cac thong so tren cac nhan tu mang da tinh toan
Label8.Caption = ""
Label9.Caption = ""
Timer1.Enabled = False
Label17.Caption = " WORKING"
'Hien thi thoi gian lam viec len man hinh
Label8.Caption = " " & CStr(Round(Thoigianlamviec(biendem), 2)) & "
h"
Tlv = Thoigianlamviec(biendem)
Tph = Thoigianphuchoi(biendem)
'Tinhh gia tri exp(-lamda*t)
Plv(0) = 0
Plv(biendem) = Exp(-lamda * Tlv)
'Khoi dong timer2 de bat dau qua trinh kiem tra oto
Timer2.Enabled = True
Image1.Visible = True
Image2.Visible = False

Picture3.Print Format(Round((i + 1) / 10, 2))
Next
'Danh so diem'0'
Picture3.CurrentX = 400
Picture3.CurrentY = maxy - 200
Picture3.Print Format(0)
'Ve cac vach doc tren truc nam ngang va danh so
dt = c * 100
For i = 1 To 10
Picture3.DrawStyle = 0
Picture3.Line (500 + i * (kx - 500) / 10, maxy - 240)-(500 + i * (kx -
500) / 10, maxy - 160), QBColor(5)
Picture3.CurrentX = 250 + i * (kx - 500) / 10
Picture3.CurrentY = maxy - 200
Picture3.Print Format(i * dt, "0")
Next
've ghi chu
Picture3.Line (ex - 1500, 100)-(ex - 1000, 100), QBColor(2)
Picture3.Print " WORK "
Picture3.Line (ex - 1500, 300)-(ex - 1000, 300), QBColor(12)
Picture3.Print " RESTORE "
Picture3.CurrentX = 2000
Picture3.CurrentY = 100
Picture3.Print "STATUS OF OPERATION CHART"
'Ve thoi gian lam viec theo ham cong e mu
d = (kx - 500) / 1000
Tlv = Thoigianlamviec(biendem)
Tph = Thoigianphuchoi(biendem)
t(0) = 0
t(0) = 0

Picture3.Line (500 + t(biendem) + w(biendem) - (Tph / c) * d, ky - 10 *
ky / 12 + (1 - Plv(biendem)) * 10 * ky / 12)-(500 + t(biendem) +
w(biendem), ky - 10 * ky / 12), QBColor(12)
'Ve duong dong
Picture3.DrawStyle = 2
Picture3.Line -(500 + t(biendem) + w(biendem), ky + 40)
End Sub
Private Sub Timer3_Timer()
'Bien i1, i2 co tac dung la cac bien dem thoi gian tao ra
'cac ham tre thoi gian theo y muon
i1 = i1 + 1
i2 = i2 + 1
If i2 >= 10 Then
If Picphantram.Left < 0 Then
Label17.Caption = " RESTORING "
Picphantram.Move Picphantram.Left + 10
i1 = 0
End If
'sau khi phuc hoi xong dua ra thong bao DA PHUC HOI
If Picphantram.Left >= 0 And i1 = 0 Then
Label17.Caption = " RESTORED ALREADY "
End If
'***************************************
End If
If i1 >= 100 Then
i2 = 0
If biendem < 10 Then
Label17.Caption = " ARE YOU READY?"
Timer3.Enabled = False
Timer1.Enabled = True

Tổng thời gian phục hồi : 7784.83h
Hệ số sẵn sàng : 0.93
3.Nhận xét bình luận kết quả mô phỏng :
Mô hình mô phỏng hoàn toàn phù hợp với yêu cầu của bài toán.
Em đã hoàn thành bài tập môn mô hình hóa với đề tài : “Đánh
giá hệ số sẵn sàng của hệ thống có phục hồi”. Một lần nữa em xin
chân thành cảm ơn cô giáo Phạm Thị Hồng Anh đã hướng dẫn giúp
em thực hiện tốt bài tập này.
TÀI LIỆU THAM KHẢO
1.Mô hình hóa hệ thống và mô phỏng. Gs.Ts.Nguyễn công
Hiền
Ts.Nguyễn Phạm Thục Anh
NXB Khoa học kỹ thuật HN
2.Giáo trình xác suất thống kê. Tống Đình Quỳ
NXB Giáo dục HN
3.Lập trình cơ sở dữ liệu 6.0. Nguyễn thị ngọc Mai
NXB Thống Kê
4.Tự học lập trình Visual Basic trong 21 Ngày
NXB Giáo Dục