-
Notifications
You must be signed in to change notification settings - Fork 0
/
2019 China mathematical modeling Porblem C
140 lines (133 loc) · 4.18 KB
/
2019 China mathematical modeling Porblem C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
%% %导入爬取得数据
load data_time
load data_taxi
%% %初始化数据
[m,~]=size(data_time); %计算航班数
pl=randi([100,180],1,m); %随机生成一组数据，每个航班的乘客数
pl_flag=1;
tx_flag=1;
pl_sum=zeros(24,60); %每个时间段的人数变化
plane=data_time; %飞机的航班数据
num_pl=0; %初始化站台上剩余人数
pl_num=zeros(1,m); %记录每分钟站台上剩余人数
num_take=0 ;%初始化离开的车辆数
num_tpl=10;
num_wait=0;
taxi=data_taxi;
decih=randi([1,22],1,500);
decim=randi([1,59],1,500);
all_data=zeros(500,5);
income_total=zeros(500,1);
% k=1;
de_flag=0;
now_num_tx=nan;
now_num_pl=0;
now_num_take=0;
%begin_time=0;
for k=1:500
pl_flag=1;
tx_flag=1;
pl_sum=zeros(24,60); %每个时间段的人数变化
plane=data_time; %飞机的航班数据
num_pl=0; %初始化站台上剩余人数
pl_num=zeros(1,m); %记录每分钟站台上剩余人数
num_take=0 ;%初始化离开的车辆数
num_tpl=10;
num_wait=0;
de_flag=0;
now_num_tx=nan;
now_num_pl=0;
now_num_take=0;
for i=1:23
for j=1:60
%%模拟登机
if ismember(i*100+j,plane)
num_pl=num_pl+pl(1,pl_flag);
pl_num(1,pl_flag)=floor(pl(pl_flag)*0.08);% 假定在地铁工作时间内，乘坐公交车的人数有1%
num_tpl=num_tpl+ pl_num(1,pl_flag);
pl_flag=pl_flag+1;
end
%%
%模拟停车场剩余车辆
if ismember(i*100+j,taxi(:,1))
num_tx=taxi(tx_flag,2);
tx_flag=tx_flag+1;
end
%%
%模拟接客
if num_tpl>0
pl_take=randi([2,3]);% 设定每分钟有1名或者5名乘客被劫走
num_tpl=num_tpl-pl_take;
if num_tpl<0
num_tpl=0;
end
num_take=num_take+1;
else
num_wait=num_wait+1;
end
pl_sum(i,j)=num_tpl;
tx_sum(i,j)=num_tx;
%% %模拟司机决策
if i==decih(1,k)&&decim(1,k)==j %假定某一时刻司机到达
now_num_tx=num_tx;
now_num_pl=num_pl;
now_num_take=num_take;
begin_time=100*i+j;
begin_timeg=[i,j];
end_time=0;
de_flag=1;
end
if num_take-now_num_take==now_num_tx&& de_flag==1 %计算排序可能花费的时间
end_time=100*i+j;
end_timeg=[i,j];
cost_time=begin_time-end_time;
cost_timeg=end_timeg- begin_timeg;
de_flag=0;
end
end
end
%% %模拟在机场拉客产生的收益与成本
dis=normrnd(36.6,9) ; %随机生成正态分布的距离
%% %价格
st_pri=6;
pri=1.5;
back_pri=0.5;
oil_pri=0.5;
dis_line=2;
%% 速度与计算
speed=53; %根据百度地图导航生成的速度
if dis<12
income=(dis-dis_line)*pri+st_pri-oil_pri*dis;
else
income=((dis-dis_line)*pri+st_pri)*(1+back_pri)-oil_pri*dis;
end
time=dis/speed;
total_time=cost_timeg(1,1)+cost_timeg(1,2)/60+time;
per_income= income/ total_time;
income_total(k,1)= per_income;
%% &模拟回市区拉客产生的收益与成本
back_dis=36.6;
back_time=back_dis/speed;
back_cost=dis*oil_pri;
back_income=getincome(begin_timeg(1,1),begin_timeg(1,2),total_time-back_time)-back_cost;
%% %做决策
if back_income>income
disp(['现在时间是:',num2str(begin_timeg(1,1)),':',num2str(begin_timeg(1,2))]);
disp(['决策为','回去拉客',',','此时每小时收益为:',num2str(back_income/total_time),'元',' 耗费时间为：',num2str(total_time),'小时']);
all_data(k,1)=0;
all_data(k,2)=begin_timeg(1,1)+begin_timeg(1,2)/60;
all_data(k,3)=back_income/total_time;
all_data(k,4)=total_time;
all_data(k,5)=dis;
else
disp(['现在时间是:',num2str(begin_timeg(1,1)),':',num2str(begin_timeg(1,2))]);
disp(['决策为','过去排队',',',' 此时每小时收益为:',num2str(income/total_time),'元',' 耗费时间为：',num2str(total_time),'小时']);
all_data(k,1)=1;
all_data(k,2)=begin_timeg(1,1)+begin_timeg(1,2)/60;
all_data(k,3)=income/total_time;
all_data(k,4)=total_time;
all_data(k,5)=dis;
end
end
sum(income_total)/500
sum(pl_num);