figure(2);
plot(Xn,Yn,'b.');
hold on;
plot(Xm,Ym,'r.');
hold on;
times
times = times + 1;
SATVs = SATV*ones(1,N);
Tpk = zeros(M,1); %代价函数
NEXT_ptr = 0;
NEXT_Set = ones(1,M);
while(NEXT_ptr<= M)
%所有活动访问点半径均被增加且所得为合适解
%计算代价函数
for j = 1:M
for i1 = 1:N
d = sqrt( (Xn(i1) - Xm(j))^2 + (Yn(i1) - Ym(j))^2 );
%判断是否在一定范围之内
if d <= Rs(j)
%进行资源分配
Tpk(j) = Tpk(j) + Requst(j,i1);
else
Tpk(j) = Tpk(j);
end
end
end
[A,I] = sort(Tpk);
if A > 0
%选择最小的一个
Tpk_min = A(1);
Tpk_ind = I(1);
NEXT_Set(Tpk_ind) = 0;
if feasible(A,rij) == 1
%没有被违反
Rs(Tpk_ind) = Rs(Tpk_ind) + Step;
if (NEXT_Set(Tpk_ind)) == 0
NEXT_ptr = NEXT_ptr;
else
NEXT_ptr = NEXT_ptr + 1;
end
else
%违反了,则直接退出进入下一个循环
NEXT_ptr = M+1;
end
else
%如果流量为0,则说明没有发生任何请求,其实半径自动递增
Tpk_min = A(1);
Tpk_ind = I(1);
Rs(Tpk_ind) = Rs(Tpk_ind) + Step;
end
end
%多个拥塞节点的重新设置影响半径
for j = 1:M
%表示该访问点处于第1阶段
if FLag(j) == 0
%计算每个节点到访问点的距离
for i1 = 1:N
d = sqrt( (Xn(i1) - Xm(j))^2 + (Yn(i1) - Ym(j))^2 );
%判断是否在一定范围之内
if d <= Rs(j)
%进行资源分配
SATVs(1,i1) = SATVs(1,i1) - Requst(j,i1);
else
SATVs(1,i1) = SATVs(1,i1);
end
%每次请求完之后,判断是否拥堵
if SATVs(1,i1) <= 0%表示拥堵
saturated_state{j,i1} = [1,Rs',Xm(j),Ym(j),Xn(i1),Yn(i1)];
FLag(j) = 1;
else
saturated_state{j,i1} = [0,zeros(1,M),0,0,0,0];
FLag(j) = FLag(j);
end
end
end
%********************************************************************
end
%绘制仿真结果
figure(3);
subplot(421);
plot(R(1,:),'b','linewidth',2);
xlabel('TIMES');
ylabel('Radius');
grid on;
title('资源点1半径请求变化');
subplot(422);
plot(R(2,:),'b','linewidth',2);
xlabel('TIMES');
ylabel('Radius');
grid on;
title('资源点2半径请求变化');
subplot(423);
plot(R(3,:),'b','linewidth',2);
xlabel('TIMES');
ylabel('Radius');
grid on;
title('资源点3半径请求变化');
subplot(424);
plot(R(4,:),'b','linewidth',2);
xlabel('TIMES');
ylabel('Radius');
grid on;
title('资源点4半径请求变化');
subplot(425);
plot(R(5,:),'b','linewidth',2);
xlabel('TIMES');
ylabel('Radius');
grid on;
title('资源点5半径请求变化');
subplot(426);
plot(R(6,:),'b','linewidth',2);
xlabel('TIMES');
ylabel('Radius');
grid on;
title('资源点6半径请求变化');
subplot(427);
plot(R(7,:),'b','linewidth',2);
xlabel('TIMES');
ylabel('Radius');
grid on;
title('资源点7半径请求变化');
subplot(428);
plot(R(8,:),'b','linewidth',2);
xlabel('TIMES');
ylabel('Radius');
grid on;
title('资源点8半径请求变化');
%绘制仿真结果
figure(4);
subplot(421);
plot(TPK(1,10:end),'b','linewidth',2);
xlabel('TIMES');
ylabel('代价函数');
grid on;
title('资源点1代价函数');
subplot(422);
plot(TPK(2,10:end),'b','linewidth',2);
xlabel('TIMES');
ylabel('代价函数');
grid on;
title('资源点2代价函数');
subplot(423);
plot(TPK(3,10:end),'b','linewidth',2);
xlabel('TIMES');
ylabel('代价函数');
grid on;
title('资源点3代价函数');
subplot(424);
plot(TPK(4,10:end),'b','linewidth',2);
xlabel('TIMES');
ylabel('代价函数');
grid on;
title('资源点4代价函数');
subplot(425);
plot(TPK(5,10:end),'b','linewidth',2);
xlabel('TIMES');
ylabel('代价函数');
grid on;
title('资源点5代价函数');
subplot(426);
plot(TPK(6,10:end),'b','linewidth',2);
xlabel('TIMES');
ylabel('代价函数');
grid on;
title('资源点6代价函数');
subplot(427);
plot(TPK(7,10:end),'b','linewidth',2);
xlabel('TIMES');
ylabel('代价函数');
grid on;
title('资源点7代价函数');
subplot(428);
plot(TPK(8,10:end),'b','linewidth',2);
xlabel('TIMES');
ylabel('代价函数');
grid on;
title('资源点8代价函数');
12_008m
- 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.
- 144.
- 145.
- 146.
- 147.
- 148.
- 149.
- 150.
- 151.
- 152.
- 153.
- 154.
- 155.
- 156.
- 157.
- 158.
- 159.
- 160.
- 161.
- 162.
- 163.
- 164.
- 165.
- 166.
- 167.
- 168.
- 169.
- 170.
- 171.
- 172.
- 173.
- 174.
- 175.
- 176.
- 177.
- 178.
- 179.
- 180.
- 181.
- 182.
- 183.
- 184.
- 185.
- 186.