Matlab数组计算结果错误

Question

我正在尝试从论文中复制结果，并给出链接以避免写下所有需要的数学公式：

基于博弈论的DoS和DDoS攻击防御机制建模与仿真

更具体地说，我遇到的问题是图3的图。 该图在z轴上给出了给定两个变量m和M的方程式3的结果。将需要的其他方程式是5,6,7，并且方程式6之前的段落中还有两个小方程式。看看Xi是什么，检查4.2部分。 所需的所有变量值都在绘图之前给出。

现在，我要尝试在matlab中创建完全相同的图，但是我失败了，并且我需要帮助，因为我的matlab技能​​不是很好。

我有一个脚本文件，其中包含以下内容：

w1  = 1000;
w2  = 1000;
w3  = 10;  
B   = 2000;
n   = 20;  
r_l = 60;  
s_l = 20;  
g   = 10;  
a_f = 5000;
b   = 20;  

vx = 0 : 1 : 500;
vy = 0 : 1 : 90;

[x,y] = meshgrid(vx,vy);

z = payoff(w1, w2, w3, y, r_l, n, g, B, b, x, s_l, a_f);

h = surfc(x,y,z);
set(h, 'edgecolor','none')
xlabel('Firewall Midpoint (M)')
ylabel('Number of zombies')
zlabel('Attackers payoff')
view(-41,11);

支付功能如下：

function out = payoff(w1, w2, w3, m, r_l, n, g, B, b, M, s_l, a_f)

r_a = a_f./ m;
r_a_dash = r_a.*(1-Fx(r_a, b, M, B));
r_l_dash = r_l.*(1-Fx(r_l, b, M, B));

v_b = ( m .* r_a_dash ) ./ ( n .* r_l_dash + m .* r_a_dash );
v_n = normcdf(( g .* ( n .* r_l_dash + m .* r_a_dash ) ./ B ), r_l, s_l);

out = w1 * v_b + w2 * v_n - w3 * m;

Fx再​​次执行以下功能：

function out = Fx(x,b,M,B)

out=1./(1+exp(-b.*(x-M)./B)); 

我不知道确切的错误在哪里，但是下面得到的图与本文中的图不同。

本文中的图沿Firewall Midpoint轴呈U形曲线，而矿山呈单调递增。

谁能发现我有任何错误？ 提前致谢。

Answer 1

我注意到的最大的事情是在您使用的代码中：

v_n = normcdf(( g .* ( n .* r_l_dash + m .* r_a_dash ) ./ B ), r_l, s_l);

什么时候应该使用（我认为）：

v_n = normcdf(( g .* ( n .* r_l_dash + m .* r_a_dash ) ./ B ), r_l_dash, s_l);

他们在论文中指出：

represents the expected rate of a legitimate flow. 回想一下， 代表合法流量的预期比率。 ′. 令通过防火墙的合法流量的平均速率为 '。

在normcdf函数中，第二个参数应该是平均值mu 。 这使我沿着防火墙中点形成了U形曲线，但是我看到它与图片不完全一样，并且我相信这是由于b的值所致，因为有人已经给出了该值。

希望这可以帮助。 由于我使用不同的b值进行游戏时可能仍然存在计算错误，并且仍然无法匹配纸张中的图像。