请使用matlab写一段包含过载和欠载的负载容量模型代码，节点失效后将其负载分配给邻居节点，以失效节点数变化曲线为输出结果

好的，以下是包含过载和欠载的负载容量模型代码，节点失效后将其负载分配给邻居节点，并以失效节点数变化曲线为输出结果：

clear all;
clc;

% 初始化参数
numNodes = 100; % 节点数量
p = 0.2; % 节点失效概率
c = 0.8; % 节点负载容量系数
maxSteps = 100; % 最大模拟步数
failedNodes = zeros(1, maxSteps); % 记录每一步的失效节点数

% 初始化节点状态和负载
state = ones(1, numNodes); % 1表示节点正常，0表示节点失效
load = randi([5, 15], 1, numNodes); % 节点负载

% 模拟级联失效
for i = 1:maxSteps
    % 计算每个节点的失效概率
    prob = p + (1 - p) .* (1 - exp(-c .* load)); 
    
    % 随机失效节点
    for j = 1:numNodes
        if state(j) == 1 && rand() < prob(j)
            state(j) = 0;
            % 将失效节点的负载分配给邻居节点
            neighbors = find(load > 0 & state == 1 & j ~= 1:numNodes);
            if ~isempty(neighbors)
                % 判断失效节点的负载是否过载或欠载
                if load(j) > mean(load(neighbors))
                    load(neighbors) = load(neighbors) + (load(j) - mean(load(neighbors))) / length(neighbors);
                else
                    load(neighbors) = load(neighbors) + (load(j) - mean(load(neighbors))) * length(neighbors);
                end
            end
            load(j) = 0;
        end
    end
    
    % 统计失效节点数
    failedNodes(i) = numNodes - sum(state);
    
    % 如果所有节点都已失效，停止模拟
    if failedNodes(i) == numNodes
        break;
    end
    
    % 更新负载
    load = load + randi([-5, 5], 1, numNodes);
    % 将负载限制在[1, 20]范围内
    load(load < 1) = 1;
    load(load > 20) = 20;
end

% 绘制失效节点数变化曲线图
plot(1:i, failedNodes(1:i), '-o');
xlabel('模拟步数');
ylabel('失效节点数');
title('包含过载和欠载的负载容量模型的级联失效模拟（负载分配）');

代码中在节点失效后，将其负载分配给邻居节点。具体来说，找到所有负载大于0且状态为正常的邻居节点，根据失效节点的负载与邻居节点平均负载的关系，将失效节点的负载分配给这些邻居节点。如果失效节点的负载大于邻居节点平均负载，则将其超过平均负载的部分平均分配给邻居节点；如果失效节点的负载小于邻居节点平均负载，则将其不足平均负载的部分全部分配给邻居节点。最后，统计每一步的失效节点数，将其绘制为失效节点数变化曲线。