Mô phỏng cuộc chạy thi giửa 2 tàu chiến gọi là Tàu 1 và Tàu 2.doc.DOC

bài 2.
1.Đề bài:
Mô phỏng cuộc chạy thi giửa 2 tàu chiến gọi là Tàu 1 và Tàu 2.
Nội dung nh sau:
Có 2 tàu chiến gọi là Tàu1 và Tàu 2 xuất phát tại cùng một căn cứ điểm và
đi đến đích trong một lộ trình dài Y (km). Trong quá trình đi thì các tàu phải dừng
lại tiếp nhiên liệu dọc đờng. Do có dựa vào trang thiết bị hiện đại và đợc chế tạo
với công nghệ cao hơn nên Tàu 1 có vận tốc lớn hơn Tàu 2 tuy nhiên cũng vì lý do
đó mà Tàu1 có hệ số sẵn sàng xuất phát bé thua hệ số sẵn sàng xuất phát của Tàu
2. Trên đờng đi cứ mỗi lần dừng lại để tiếp nhiên liệu các tàu lại tiếp tục xuất phát
với hệ số sẵn sàng lần lợc là hssst1 và hssst2. Quãng đờng đi đợc của cả hai tàu đ-
ợc cộng dồn.
Hãy thiết lập một bài toán mô phỏng với các hệ số sẵn sàng là một số ngẫu
nhiên sao cho Tàu 2 có xác suất về đích lớn hơn Tàu 1.
2. Phân tích đề bài.
Bài toán này thuộc loại bài toán mô hình hoá hệ ngẫu nhiên với biến ngẫu
nhiên ở đây là hệ số sẵn sàng xuất phát của các tàu. Trong thuật giải này thì vận
tốc của các tàu có thể đợc thể hiện qua quãng đờng mà nó di chuyển đợc trong
một khoảng thời gian nhất định cho trớc. Nh vậy khi thực hiện gieo một số ngẫu
nhiên thì ta phải so sánh với các hệ số sẵn sàng của cả hai tàu để ra quyết định
tàu nào xuất phát. Thực hiện gieo các số ngẫu nhiên để thử là nội dung của ph-
ơng pháp mô hình hoá hệ ngẫu nhiên.
Nội dung của phơng pháp đợc trình bày ở trong bài 1.
3. Các biến ngẫu nhiên
Nh trên đã trình bày sự xuất phát của các tàu đợc đánh giá qua một số
ngẫu nhiên đó là hệ số xuất phát. Nh vậy biến ngẫu nhiên ở đây chính là hệ số
xuất phát. Nếu nh hệ số sẵn sang của các tàu đợc lấy các giá trị trong khoảng
[0,1] thì số ngẫu nhiên là một số có phân bố U(0,1).
4. Phân tích giải thuật.
Giải thuật cho bài toán trên đợc xây dựng nh sau:
+ Cho trớc các thông số không đổi trong quá trình chuyển động nh vận tốc

của các tầu và quãng đờng di chuyển sau mỗi lần tiếp nhiên liệu.
+ Gieo một số ngẫu nhiên x=U(0,1):
Nếu x<hssst1 thì cả hai tàu xuất phát
Nếu x<hssst2 và x>hssst1 thì Tàu 2 xuất phát
+ Thực hiện so sánh quãng đờng di chuyển đợc của cả hai tàu so với độ
dài quãng đờng ấn định trớc nếu tầu nào về đích trớc thì kết thúc cuộc chơi.
5. Lu đồ thuật giải nh sau:
Nhập số liệu
hssst1,hssst2
Lấy x=U(0,1)
Nhập Y
x<=hssst1
hssst1<=x<=hssst2
Đ
Tau1cd,
Tàu 2cd
S
S
END
xt1=Y or
xt2=Y
xt1=xt1+x01 or
xt2=xt2+x02
Đ
Tàu 2cd
xt2=xt2+x02
xt2=Y
Đ
S
S

Đ
S
6.Điều kiện đầu, điều kiện dừng mô phỏng.
Điều kiện đầu ở đây là khoảng các của hành trình và các khoảng các
chuyển động của Tàu 1 và Tàu 2.
Điều kiện dừng chuyển động là khi một trong 2 tàu về đích.
7.Chơng trình mô phỏng.
Chơng trình sau đợc xây dựng trên ngôn ngữ C với mã nguồn đợc cho dới
đây.
#include"mou.h"
#include"iostream.h"
#include"math.h"
#include"conio.h"
#include"stdio.h"
#include"graphics.h"
#define Enter 13
void Write_text(int x1, int y1,double x,char*s,int color);
void write_value(int*x,int*y,char*s);
void read_value(int*x,int*y,char*s);
int Nhap_kc_tau1(void);
int Nhap_kc_tau2(void);
double Nhap_hsss_tau1(void);
double Nhap_hsss_tau2(void);
void Data_in(void);
void tau1_cd(void);
void tau2_cd(void);
void ve_tau1(void);
void ve_tau2(void);
void Win(void);
void Run(void);

void button(int x1,int y1,int x2,int y2,int color,char s[],int
color_s,int size);
void button1(int x1,int y1,int x2,int y2,int color);
void _line(int x1, int y1, int x2, int y2, int color);
void _line1(int x1,int y1,int x2,int y2);
void main_h(void);
void help(void);
int xt,yt,xr,yr,nt,nr,i,j=0,k=0,t=0,kct1,kct2,a=0;
double hst1,hst2;
char*pr,*pt;
void write_value(int*x,int*y,char*s)
{
outtextxy(*x,*y,s);
*x+=textwidth(s);
}
void read_value(int*x,int*y,char*s)
{
int i=0;char ch[2];
ch[1]=0;
while(1)
{
ch[0]=getch();
if(ch[0]==Enter) break;
write_value(x,y,ch);
s[i]=ch[0];
++i;
}
s[i]=0;
}
int Nhap_kc_tau1()

{
int xt1=220,yt1=380;
char*st1;
read_value(&xt1,&yt1,st1);
return atoi(st1);
}
int Nhap_kc_tau2()
{
int xt2=220,yt2=435;
char*st2;
read_value(&xt2,&yt2,st2);
return atoi(st2);
}
double Nhap_hsss_tau1()
{
int xst1=470,yst1=380;
char*sst1;
read_value(&xst1,&yst1,sst1);
return atof(sst1);
}
double Nhap_hsss_tau2()
{
int xst2=470,yst2=435;
char*sst2;
read_value(&xst2,&yst2,sst2);
return atof(sst2);
}
void Data_in(void)
{
char ch;

setcolor(1);
kct1=Nhap_kc_tau1();
kct2=Nhap_kc_tau2();
hst1=Nhap_hsss_tau1();
hst2=Nhap_hsss_tau2();
//ch=getch();
// if(ch==Enter) Run();
// Run();
}
void ve_tau1(void)
{
setcolor(RED);
ellipse(100,50,0,360,20,8);
ellipse(100,46,190,357,20,6);
line(107,44,110,38);
circle(110,38,2);
line(93,44,90,38);
circle(90,38,2);
setfillstyle(SOLID_FILL,BLUE);
floodfill(101,54,RED);
setfillstyle(SOLID_FILL,MAGENTA);
floodfill(94,45,RED);
}
void ve_tau2(void)
{
setcolor(7);
ellipse(100,50,0,360,20,8);
ellipse(100,46,190,357,20,6);
line(107,44,110,38);
circle(110,38,2);

line(93,44,90,38);
circle(90,38,2);
setfillstyle(SOLID_FILL,GREEN);
floodfill(101,54,7);
setfillstyle(SOLID_FILL,13);
floodfill(94,45,7);
}
void tau2_cd(void)
{
if(j<1)
putimage(50,170,pr,XOR_PUT);
j=j+1;
if(j>1)
putimage(xr,yr,pr,XOR_PUT);
xr=50+kct2*(j-1);yr=170;
while(xr<kct2*j+50)
{
putimage(xr,yr,pr,XOR_PUT);
delay(500);
putimage(xr,yr,pr,XOR_PUT);
xr=xr+7;
if(xr>=540)
{
setcolor(15);
outtextxy(250,250,"Tau 2 chien thang");
putimage(xr,yr,pr,XOR_PUT);
Win();
putimage(xr,yr,pr,XOR_PUT);
}
if(xr>=kct2*j+50)

putimage(xr,yr,pr,XOR_PUT);
}
}
void tau1_cd()
{
if(k<1)
putimage(50,250,pt,XOR_PUT);
k=k+1;
if(k>1)
putimage(xt,yt,pt,XOR_PUT);
xt=50+kct1*(k-1);yt=250;
while(xt<kct1*k+50)