更新：
以下模型已经上云实现，体验地址：
复杂网络网络计算平台

上一篇介绍了随机攻击网络节点与蓄意攻击节点的基本方法。

其中随机攻击的部分还有一些瑕疵，就是在实际的研究中，需要对网络进行多次（数十次甚至上百次）攻击后取指标变化平均值，这样的实验数据才具有一定的说服力。

其实这个问题乍一听起来，原理也比较简单：就是让一个程序运行指定的次数然后，累加程序中某一个变量后取平均值。

实现起来也并不费劲，首先就是把这个程序定义为一个函数，然后把所需要累加的变量作为函数的返回值。然后在另一个文件中创建循环，在循环中调用该函数，用一个变量接收该函数的返回值，以达到累加的效果，最后在循环外部取一个平均值即可。

代码如下所示：

定义函数文件名称（该名称需要与函数名称相同）

函数的参数介绍：

输入值： str：意为数据文件的路径；numDelete:删除节点的个数（这里的命名只是为了方便并不是必选项）
返回值：Eglob，即网络效率值的数组

function Eglob = ATestAver(str,numDelete)
%输入 ： 
%        str：意为数据文件的路径
%        numDelete:删除节点的个数

%返回值：Eglob，即网络效率值的数组

%加载数据文件
load(str);

Name_Struct = Node_Key_Sort_Descend;       %  Name_Struct 数据集名称，更换网络数据集时，需要更改此处结构体名称
A_Init =  Name_Struct.Adjacent_Matrix;      %% 网络邻接矩阵
N_Init = size(A_Init,1);                     %% 节点个数

NetEff_Init = zeros(1,numDelete);
Struct_Init = struct('Deg',NetEff_Init);

% 初始网络性能
%生成随机数，以此进行随机攻
Name_Struct.Node_Key_Degree = randperm(440);

%% 
%  按照 Degree 算法排序，删除节点
A = A_Init;          %% 网络邻接矩阵 A
for i = 1:numDelete
  % 按照 Degree 算法排序，删除节点
end

定义测试文件：

文件中需要定义随机攻击的次数和随机攻击节点的个数，具体参数设置应视具体网络而定。
在调用函数时，要传入文件路径，和删除节点的个数。

numDelete = 22;

netSum = zeros(1,numDelete);

numRandom = 50;
for i=1:numRandom
    
    netI = ATestAver('Data\test.mat',numDelete);
    netSum = netSum + netI;
end

%求出平均值
netAver = netSum/numRandom;

我的测试数据具有小世界特性与无标度特性。并得到了以下的结果:

大家可以看出，经过多次随机攻击取指标平均值之后，曲线近似于一条直线，下降速率较小且符合实际的网络情况。验证了无标度网络对于随机攻击拥有较好的鲁棒性。

复杂网络函数包分享
内含有各种复杂网络函数可直接使用。

希望本文对大家有所帮助，有任何问题或者是建议，欢迎大家与我交流。

上一篇文章：matlab实现随机攻击网络节点+蓄意攻击网络节点（1）

我们正在开发一个线上可视化复杂网络仿真平台，支持：
一键上传邻接矩阵（包含数据直接转换为数字邻接矩阵）

计算核心网络指标

模拟鲁棒性攻击

可视化网络结构变化

📌 点击加入内测 / 提出您的建议：

👉 平台入口（内测中）
如果需要本地化的解决方案，也可以联系我们。
并且我准备建立一个复杂网络交流群，想进群可以联系我