BÀI TẬP LỚN
MÔ HÌNH HÓA HỆ THỐNG NGẪU NHIÊN
1.Tên đề tài.
Mô phỏng trạm xe bus :
- Sinh viên đi từ ký túc xá đến trường bằng xe bus,mỗi xe chứa được 60
sinh viên .Thời gian đi đến trường bắt đầu từ 6h30 đến 7h30,sinh viên đi đến
trạm được mô tả bằng dòng tối giản với cường độ λ=1sv/s.Cứ sau 15’ thì có
một chuyến xe bus đến trường .Nếu số sinh viên chờ xe < 60 sinh viên thì xe
vẫn chạy đúng giờ.Nếu số sinh viên chờ xe > 60 thì số sinh viên thừa sẽ phải
chờ chuyến sau.Sinh viên than phiền rằng không đủ xe để đi đến trường nên
sau 7h30 vẫn còn nhiều sinh viên bị kẹt lại không kịp đến trường đúng giờ.
Để đảm bảo không có sinh viên bị chậm giờ học thì thời gian giữa các
chuyến xe phải là bao nhiêu?
2.Phân tích đề tài và tóm tắt phương pháp mô hình hóa.
- Đề tài đưa ra trên cơ sở phân tích hoạt động của trạm xe bus.
- Xây dựng phương pháp mô hình hóa hệ ngẫu nhiên để tìm ra phương án tối
ưu nhất cho hoạt động của các chuyến xe giúp việc vận chuyển sinh viên đạt
hiệu quả ,không có sinh viên nào bị kẹt tại trạm xe bus.
Ta xây dựng thuật toán mô phỏng qua các bước :
• Xây dựng mô hình sinh viên đi đến trạm xe bus.Khoảng cách giữa
các trạm xe bus là t
i
=-
λ
1
ln(U
i
), U
i

≈

e
x
λ
λ
−
nếu x
∈
(0,1)
0 phần còn lại
+)Hàm phân phối:
0 nếu x< 0
F(x)=
[ ]
∑
=
−
x
i
i
i
e
0
!
λ
λ
nếu x<=0
+)Thông số λ>0
+)kỳ vọng toán M(x)= λ
+)Phương sai:
2

Điều kiện đầu:thời gian t=0,c=0,s=0
Điều kiện mô phỏng : nhập giá trị λ>0 và Txe>0
Điều kiện ngừng mô phỏng :khi t=3600(s)
Ta sẻ sử dụng các biến sau để lập trình mô phỏng cho hệ:
t : tổng thời gian đi đến bến xe bus của sinh viên
Txe : chu kỳ thời gian của xe bus
Tx : tổng thời gian mà xe bus đã đi(có giá tri đầu là: tx=Txe)
a : giá trị cường độ lamda(λ) (sẻ được nhập khi mô phỏng)
dk :để điều khiển dừng khi hết thời gian mô phỏng (giá trị ban đầu dk=1)
c: số sinh viên đang chờ ở bến xe bus
s : số chuyến xe bus đã chạy
7.Chương trình mô phỏng trên VB
Option Explicit
Dim i, j, k, kh, tx, txe, txe1, phut, times, c, s, dk As Integer
Dim t, a, u As Double
Dim ast, bst, at, bt As String
Dim X(160), Y(160) As Integer
'***thu tuc tinh****
'***khai bao bien***
Public Sub tinh()
Const e = 2.7182
a = CDbl(txtCDSV.Text * 100)
u = Rnd
t = t - 100 * Log(u) / (a * Log(e))
If t < 3600 Then
c = c + 1
Else: t = 3600
End If
End Sub
'******thu tuc ve Form 2********

'*******nut lenh hien thi Form2*****************
Private Sub cmdDB_Click()
txt1.Text = 6
txt2.Text = 30
t = 0
s = 0
c = 0
If txttxe.Text = "" Then
MsgBox "Nhap tg giua cac chuyen xe", vbOK
Else
txe = CDbl(txttxe.Text) * 60
End If
tx = txe
dk = 1
pt2.Visible = True
Timer1.Enabled = True
Timer2.Enabled = True
End Sub
Private Sub cmddung_Click()
Timer1.Enabled = False
Timer2.Enabled = False
pt2.Left = -1920
pt2.Visible = False
End Sub
'**********nut lenh ket thuc chuong trinh**********
Private Sub Command1_Click()
Beep
End
End Sub
'*******thu tuc dieu khien chuong trinh bang bo timerl********

Pic4(i).Visible = False
Next
Picture1((kh - 1) / 3).Visible = True
End If
If kh Mod 3 = 0 Then
For i = 0 To (kh / 3) - 1
Picture1(i).Visible = True
Pic2(i).Visible = True
Pic4(i).Visible = True
Next
For i = (kh - 1) / 3 To 44
Picture1(i).Visible = False
Pic2(i).Visible = False
Pic4(i).Visible = False
Next
End If
If kh Mod 3 = 2 Then
For i = 0 To (kh - 2) / 3 - 1
Picture1(i).Visible = True
Pic2(i).Visible = True
Pic4(i).Visible = True
Next
For i = (kh - 1) / 3 To 44
Picture1(i).Visible = False
Pic2(i).Visible = False
Pic4(i).Visible = False
Next
Picture1((kh - 2) / 3).Visible = True
Pic2((kh - 2) / 3).Visible = True
End If

Timer1.Interval = 40
If dk = 1 Then
If t >= (txe - 10) Then
imI.Visible = True
imI.Left = imI.Left + 150
End If
If imI.Left >= 12000 Then
imI.Left = 0
End If
If (t < tx) Then
tinh
veduong
'hien thi thoi gian hien tai
h = 7 + Int((t - 1800) / 3600)
If t < 1800 Then
p = 30 + Int(t / 60)
Else: p = Int((t - 1800) / 60)
End If
List4.Clear
List4.AddItem " " & h
List5.Clear
List5.AddItem " " & p
List3.Clear
List3.AddItem " " & c
If t >= tx Then
Timer1.Interval = 1000
imI.Left = 0
veduong 've SV dang den ben
If c > 60 Then
c = c - 60

10.Nhận xét.
Qua kết quả mô phỏng ta nhận thấy để sinh viên đến trường đúng giờ và ko
có sinh viên phải đợi nên chọn chu kỳ chạy của xe bus là : Txe=7 phút 30 giây

Tài liêu tham khảo:
[1] Bài giảng môn mô hình hóa hệ thống.
[2] Bài giảng Visual Basic
Khoa Công Nghệ Thông Tin-Trường Đại Học Hàng Hải-2008