-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathMSN_CALCULATE4.m
143 lines (129 loc) · 8.17 KB
/
MSN_CALCULATE4.m
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
141
142
143
%%%%%%%%%%%%%%%% MSN system %%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%% MSN_CALCULATE.m %%%%%%%%%%%%%%%%%%%%
% This script is used to calcualte the mobile of MNs%
function [MN_DATA_SOCIA] = MSN_CALCULATE4(input_settings,MN_DATA,s_data_day,ageing,protocol)
%myFun - Description
%
% Syntax: [MN_DATA] = myFun(input_settings,MN_DATA)
%
% Long description
%根据论文2，计算移动节点的社交属性
clear MN_DATA_SOCIA_temp;
%The Role Playing Mobility Model.
global MN_DATA_SOCIA_temp; %global MN_DATA Index.
MN_DATA_SOCIA_temp = MN_DATA;
%ROUTING_TABLE_temp = ROUTING_TABLE;
%创建节点相遇的字段
for MN_INDEX = 1 : input_settings.MN_N
%记录相遇
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX).MEET_ALLTIMES = 0; %用于记录节点相遇总次数
%MN_DATA_temp.VS_NODE(MN_INDEX).MEET_ALLID = [];%用于记录所有节点相遇的ID
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX).MEET_ID = []; %用于记录相遇过节点的ID
%每一个节点都创建50个记录与对应节点相遇的储存空间
for INTERMEET_INDEX = 1 : input_settings.MN_N
%节点1与节点INTERMEET_INDEX的相遇记录，社交参数
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX).SOCIAL_CONTACT(INTERMEET_INDEX).ID = INTERMEET_INDEX;
%用于记录节点1与INTERMMET_INDEX的相遇次数
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX).SOCIAL_CONTACT(INTERMEET_INDEX).MEET_TIMES = 0;
%用于记录节点1与INTERMMET_INDEX的相遇时间
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX).SOCIAL_CONTACT(INTERMEET_INDEX).MEETING_TIME = [];
%两节点之间的相遇概率 = 1节点遇到2节点次数/1节点遇到所有节点的次数
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX).SOCIAL_CONTACT(INTERMEET_INDEX).ENCOUNTER_PROBABILITY = 0;
%修正 ：1节点遇到2节点次数/最大相遇次数
%MN_DATA_temp.VS_NODE(MN_INDEX).SOCIAL_CONTACT(INTERMEET_INDEX).ENCOUNTER_PROBABILITY = 0;
%TODO:待补充
end
end
%count = 0
%对所有时间，计算移动节点相遇状况
%waitbar
wait_bar = waitbar(0, 'Meet Situation');
set(wait_bar, 'name','Meeting...');
wb = 1;
for time = 1 : 24 * 60 * s_data_day%考虑到抽样
%对于每一个移动节点
for MN_INDEX_1 = 1 : input_settings.MN_N
%计算所有时刻和移动节点的相遇情况
%创建节点相遇数组
for MN_INDEX_2 = 1 : input_settings.MN_N
%获取当前节点1的 x,y
temp_x_1 = MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).X_POSITION(time);
temp_y_1 = MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).Y_POSITION(time);
%获取当前节点2的 x,y
temp_x_2 = MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).X_POSITION(time);
temp_y_2 = MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).Y_POSITION(time);
%计算两个节点的距离
inter_distance = sqrt( (temp_y_2 - temp_y_1)^2 + (temp_x_2 - temp_x_1)^2 );
%如果节点之间的距离小于定义的距离，那么定义为两个节点相遇，节点1相遇次数加一
%然后检查他们是否已经相遇，如果没有，那么修改两个节点之间的相遇记录
if(inter_distance < input_settings.MN_R)
%如果节点1从未与节点2相遇，然后将其添加到相遇ID数组中，记录相遇时间
if( isempty( find( MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).MEET_ID == MN_INDEX_2, 1)))
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).MEET_ALLTIMES = ... %节点1总相遇次数+1
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).MEET_ALLTIMES + 1;
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).MEET_ID(end + 1) = MN_INDEX_2; %记录相遇对象的ID
%节点1 与 节点2 的社会联系
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).SOCIAL_CONTACT(MN_INDEX_2).MEET_TIMES = ...
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).SOCIAL_CONTACT(MN_INDEX_2).MEET_TIMES + 1;%节点1与节点2相遇次数+1
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).SOCIAL_CONTACT(MN_INDEX_2).MEETING_TIME(end + 1) = ...
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).V_TIME(time);%记录节点1与节点2的相遇时刻
else
%若节点1以及和节点2相遇过，为了防止相同的时间多次相遇（这不符合常理）
%若时间与最后一次相遇的时间不同
if( (time * 60 ~= MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).SOCIAL_CONTACT(MN_INDEX_2).MEETING_TIME(end) ))
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).MEET_ALLTIMES = ... %节点1总相遇次数+1
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).MEET_ALLTIMES + 1;
%此时无需记录节点ID，因为已经有过记录
%MN_DATA_temp.VS_NODE(MN_INDEX_1).MEET_ID(end + 1) = MN_INDEX_2; %记录相遇对象的ID
%节点1 与 节点2 的社会联系
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).SOCIAL_CONTACT(MN_INDEX_2).MEET_TIMES = ...
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).SOCIAL_CONTACT(MN_INDEX_2).MEET_TIMES + 1;%节点1与节点2相遇次数+1
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).SOCIAL_CONTACT(MN_INDEX_2).MEETING_TIME(end + 1) = ...
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_1).V_TIME(time);%记录节点1与节点2的相遇时刻
end
end
%对节点2做同样的操作
if( isempty( find( MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).MEET_ID == MN_INDEX_1, 1)))
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).MEET_ALLTIMES = ... %节点1总相遇次数+1
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).MEET_ALLTIMES + 1;
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).MEET_ID(end + 1) = MN_INDEX_1; %记录相遇对象的ID
%节点1 与 节点2 的社会联系
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).SOCIAL_CONTACT(MN_INDEX_1).MEET_TIMES = ...
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).SOCIAL_CONTACT(MN_INDEX_1).MEET_TIMES + 1;%节点1与节点2相遇次数+1
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).SOCIAL_CONTACT(MN_INDEX_1).MEETING_TIME(end + 1) = ...
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).V_TIME(time);%记录节点1与节点2的相遇时刻
else
%若节点1以及和节点2相遇过，为了防止相同的时间多次相遇（这不符合常理）
%若时间与最后一次相遇的时间不同
if( (time * 60 ~= MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).SOCIAL_CONTACT(MN_INDEX_1).MEETING_TIME(end) ))
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).MEET_ALLTIMES = ... %节点1总相遇次数+1
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).MEET_ALLTIMES + 1;
%此时无需记录节点ID，因为已经有过记录
%MN_DATA_temp.VS_NODE(MN_INDEX_1).MEET_ID(end + 1) = MN_INDEX_2; %记录相遇对象的ID
%节点1 与 节点2 的社会联系
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).SOCIAL_CONTACT(MN_INDEX_1).MEET_TIMES = ...
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).SOCIAL_CONTACT(MN_INDEX_1).MEET_TIMES + 1;%节点1与节点2相遇次数+1
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).SOCIAL_CONTACT(MN_INDEX_1).MEETING_TIME(end + 1) = ...
MN_DATA_SOCIA_temp.VS_NODE(MN_INDEX_2).V_TIME(time);%记录节点1与节点2的相遇时刻
end
end
end
end
end
str_bar = ['The Meet ' num2str(wb) ' Time'];
waitbar(wb/(24 * 60 * s_data_day), wait_bar, str_bar);
wb = wb + (24*60*s_data_day)/(24*60*s_data_day);
end
close(wait_bar);
if strcmp(protocol, 'SCPR')
protocol_SCPRclc;
elseif strcmp(protocol, 'SimBet')
protocol_SimBetUtil_calc;
elseif strcmp(protocol, 'Prophet')
protocol_prophetMetricSetting2;
elseif strcmp(protocol, 'NULL')
disp('NULL')
end
MN_DATA_SOCIA = MN_DATA_SOCIA_temp;
end
%计算社交活跃度