Train.m

%% ========================================================================
% 模块名称：分布式时间同步算法训练主程序
% 功能描述：
%     1. 构建多跳网络拓扑（网格/随机2D/随机3D）。
%     2. 初始化物理层时钟模型（Skew/Offset）与逻辑层状态估计器。
%     3. 执行同步算法（ATS/MACTS/MRPCTS）并记录时域演进轨迹。
%% ========================================================================
%% 环境初始化
clc; clear; close all;

%% 节点结构体模板定义
WirelessNode_Str = struct( ...
    ...
    ... % =================================================================
    ... % 节点角色
    ... % =================================================================
    ...
    'NodeRole',                     'Follower', ... % 节点属性，默认为 跟随者 属性
    ...
    ... % =================================================================
    ... % 物理层属性
    ... % =================================================================
    ...
    'Hard_Clock_Rate',              1,  ...     % 硬件时钟速率
    'Hard_Clock_Rate_Record',       [],  ...    % 记录硬件时钟速率变化
    ...
    ... % =================================================================
    ... % 软件层属性
    ... % =================================================================
    ...
    'Hard_Clock_Rate_Comp',         1,   ...    % 硬件时钟速率补偿
    'Hard_Clock_Rate_Comp_Record',  [],   ...   % 记录硬件时钟速率补偿变化
    ...
    'Soft_Clock_Rate',              1,  ...     % 软件时钟速率
    'Soft_Clock_Rate_Record',       [],  ...    % 记录软件时钟速率变化
    ...
    'Soft_Offset',                  0,   ...    % 软件时间相位
    'Soft_Offset_Record',           [],  ...    % 记录软件时钟相位变化
    ...
    ... % =================================================================
    ... % 时间漂移和时间相位的耦合权重
    ... % =================================================================
    ...
    'rho_v',                        0.5,  ...   % 对邻居建议的时间漂移补偿的置信度
    'rho_o',                        0.5,  ...   % 对邻居建议的时间偏移补偿的置信度
    'leader_rho_v',                 0.5,  ...   % 对领导者建议的时间漂移补偿的置信度
    'leader_rho_o',                 0.5,  ...   % 对领导者建议的时间偏移补偿的置信度
    'follower_rho_v',               0.5,  ...   % 对跟随者建议的时间漂移补偿的置信度
    'follower_rho_o',               0.5,  ...   % 对跟随者建议的时间偏移补偿的置信度
    ...
    ... % =================================================================
    ... % 网络开销统计
    ... % =================================================================
    ...
    'BroadcastNumCNT',              0,  ...     % 记录一次同步周期内广播的总次数
    'BroadcastNumCNT_Record',       [], ...     % 记录每次同步周期内广播次数总次数的变化
    ...
    ... % =================================================================
    ... % 跳数控制
    ... % =================================================================
    ...
    'MaxMultiHop',                  1,  ...     % 最大通信跳数硬约束
    'CurrentHop',                   1,  ...     % 当前允许的动态广播跳数，用于跳数控制
    'Hop_Record',                   [], ...     % 记录跳数变化
    ...
    ... % =================================================================
    ... % 收敛性判断
    ... % =================================================================
    ...
    'ConvergentFlag',               0,  ...     % 收敛状态指示 (0: 瞬态/未收敛, 1: 稳态/已收敛)
    'Relative_Offset_TMP',          [], ...     % 存储与 1 跳邻居节点的相对时间相位差，用于局部收敛判断及动态多跳跳数控制
    ...
    ... %% ================================================================
    ... %% 1 跳邻居列表
    ... %% ================================================================
    ...
    'NeighborList',                 [] ...      % 1 跳邻居节点的 ID
);

%% 算法与实验环境配置
AlgList = {'ATS', 'MACTS', 'MRPCTS'};
WirelessNetwork.Alg = 'MRPCTS';

%% 开始仿真
fprintf(['**************** Simulation Start: ' WirelessNetwork.Alg ' ****************', '\n'])

%% 单轮同步次数和仿真轮数
WirelessNetwork.Sync_CNT_NUM     = 500;             % 单轮同步次数
WirelessNetwork.round_CNT_NUM    = 1;               % 仿真轮数

%% 网络生成
NetTypeList = {'GridNWK', 'RandomNWK2D', 'RandomNWK3D', 'DynamicRandomNWK3D'};
NWK.NetType = 'GridNWK';

switch NWK.NetType
    case 'GridNWK'
        % 配置栅格网络参数
        NWK.NodeNUM_Row = 1;        % 网格拓扑行数
        NWK.NodeNUM_Col = 7;        % 网格拓扑列数
        NWK.NodeNUM     = NWK.NodeNUM_Row * NWK.NodeNUM_Col;  % 节点数量计算
        
        % 生成栅格拓扑邻接表
        WirelessNetwork.NodeNeighborList = GenerateGridNWK(NWK.NodeNUM_Row, NWK.NodeNUM_Col);
    
    case 'RandomNWK2D'
        % 配置2维随机网络参数 (无 NWK.High 参数)
        NWK.Length         = 200;     % 场地长度 (m) 
        NWK.Width          = 200;     % 场地宽度 (m)
        NWK.NodeNUM        = 50;      % 节点总数 (n)
        NWK.nodeTXdistance = 50;      % 通信距离 (m)
        
        isGenerate = true;
    
        if isGenerate
            % 重新调用我们编写的 2D 网络生成函数
            WirelessNetwork.NodeNeighborList = GenerateRandomNWK2D(NWK);
        else 
            % 加载预存的邻接表
            matPath = fullfile('NWK2D', 'NodeNeighborList2D.mat');
        
            if exist(matPath, 'file')
                S = load(matPath, 'NodeNeighborList2D');
                WirelessNetwork.NodeNeighborList = S.NodeNeighborList2D;
                % 可视化处理
                if exist('NWK2D\NWK_2D.fig', 'file')
                    openfig('NWK2D\NWK_2D.fig', 'reuse'); 
                    drawnow;
                end
            else
                error('未找到网络数据文件: %s。请先将 isGenerate 设为 true 生成一次网络。', matPath);
            end
        end

    case 'RandomNWK3D'
        % 配置3维随机网络参数
        NWK.Length         = 200;     % 场地长度 (m) 
        NWK.Width          = 200;     % 场地宽度 (m)
        NWK.High           = 200;     % 场地高度 (m)
        NWK.NodeNUM        = 50;      % 节点总数 (n)
        NWK.nodeTXdistance = 50;      % 通信距离 (m)
        
        isGenerate = true;

        if isGenerate
            % 重新生成网络
            % 这里会调用 GenerateRandomNWK3D，按你现有函数逻辑会自动显示三维图
            WirelessNetwork.NodeNeighborList = GenerateRandomNWK3D(NWK);
        
        else 
            % 加载预存的邻接表
            matPath = fullfile('NWK3D', 'NodeNeighborList3D.mat');
        
            if exist(matPath, 'file')
                S = load(matPath, 'NodeNeighborList3D');
                WirelessNetwork.NodeNeighborList = S.NodeNeighborList3D;
                % 可视化处理
                if exist('NWK3D\NWK_3D.fig', 'file')
                    openfig('NWK3D\NWK_3D.fig', 'reuse'); 
                    drawnow;
                end
            else
                error('未找到网络数据文件: %s。请先将 isGenerate 设为 true 生成一次网络。', matPath);
            end
        end

    case 'DynamicRandomNWK3D'
        % 配置动态3维随机网络参数
        NWK.Length         = 200;     
        NWK.Width          = 200;     
        NWK.High           = 200;     
        NWK.NodeNUM        = 50;      
        NWK.nodeTXdistance = 50;      
        % 动态帧数与仿真同步总次数保持一致
        NWK.SnapshotNum    = WirelessNetwork.Sync_CNT_NUM; 
        
        isGenerate = false;

        if isGenerate
            WirelessNetwork.DynamicTopoSeq = GenerateDynamicNWK3D(NWK, NWK.SnapshotNum);
        else 
            matPath = fullfile('NWK3D_Dynamic', 'DynamicTopoSeq.mat');
            if exist(matPath, 'file')
                S = load(matPath, 'DynamicTopoSeq');
                WirelessNetwork.DynamicTopoSeq = S.DynamicTopoSeq;
            else
                error('未找到动态网络数据，请将 isGenerate 设为 true。');
            end
        end
        
        % 使用第 1 帧作为初始拓扑进行物理层初始化
        WirelessNetwork.NodeNeighborList = WirelessNetwork.DynamicTopoSeq{1};

    otherwise
        error('未知的网络类型配置: %s', NWK.NetType);
end

% --- 公共数据封装 ---
WirelessNetwork.NodeNUM  = NWK.NodeNUM;
WirelessNetwork.NetType  = NWK.NetType;

% 输出当前配置状态
fprintf('网络生成完毕。类型: %s, 无线节点数: %d\n', WirelessNetwork.NetType, WirelessNetwork.NodeNUM);

%% 网络全局特征度量-网络直径
WirelessNetwork.NetDiameter = CalculateNetDiameter(WirelessNetwork.NodeNUM, WirelessNetwork.NodeNeighborList);

%% 物理信道传输特性
% 通信延迟建模：采用高斯分布 N(mu, sigma^2) 叠加固定时延
WirelessNetwork.Delay_mu         = 3.33;               % [μs] 链路随机延迟期望
WirelessNetwork.Delay_sigma      = 0.0054;             % [μs] 链路随机延迟标准差
WirelessNetwork.Delay_fixed      = 3;                  % [μs] 链路固定处理时延

%% 时钟噪声与误差模型
% 时钟频率稳定性配置
WirelessNetwork.ppm              = 40;                 % [ppm] 最大频率偏移常数
% WirelessNetwork.MaxEstErro       = 1e8;                % [量化] 离散化估计精度极限
WirelessNetwork.MaxEstErroEnable = 1;                  % [布尔] 估计误差模型使能开关

WirelessNetwork.MaxMultiHop      = 4;                  % 算法多跳传播深度限制
WirelessNetwork.ClockOffsetMax   = 5e-6;               % [s] 同步收敛判据阈值

WirelessNetwork.SyncCycle        = 30;                 % [s] 同步周期
WirelessNetwork.CurrentSynNum    = 0;                  % 同步次数

%% 无线节点实例化
WirelessNetwork.WirelessNode = repmat(WirelessNode_Str, WirelessNetwork.NodeNUM, 1);

% 实例化跟随者节点的物理层参数：频率偏移 (Skew) 与 相位偏移 (Soft_Offset)
% Skew 遵循均匀分布：Skew ~ U[1-ppm, 1+ppm]
% Soft_Offset 遵循均匀分布：Offset ~ U[0, 500s]
SkewNoise_PPM       = randi([-WirelessNetwork.ppm, WirelessNetwork.ppm], WirelessNetwork.NodeNUM, 1); 
SkewHW              = 1 + SkewNoise_PPM * 1e-6;

InitPhaseOffset     = rand(WirelessNetwork.NodeNUM, 1) * 5e2; 

% 配置节点的相关信息
for i = 1 : WirelessNetwork.NodeNUM
    WirelessNetwork.WirelessNode(i).NodeRole                       = 'Follower';

    WirelessNetwork.WirelessNode(i).Hard_Clock_Rate                = SkewHW(i);
    WirelessNetwork.WirelessNode(i).Hard_Clock_Record              = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
    WirelessNetwork.WirelessNode(i).Hard_Clock_Record(1)           = WirelessNetwork.WirelessNode(i).Hard_Clock_Rate;

    WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp           = 1;
    WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp_Record    = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
    WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp_Record(1) = WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp;

    WirelessNetwork.WirelessNode(i).Soft_Clock_Rate                = WirelessNetwork.WirelessNode(i).Hard_Clock_Rate * WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp;
    WirelessNetwork.WirelessNode(i).Soft_Clock_Rate_Record         = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
    WirelessNetwork.WirelessNode(i).Soft_Clock_Rate_Record(1)      = WirelessNetwork.WirelessNode(i).Soft_Clock_Rate;

    WirelessNetwork.WirelessNode(i).Soft_Offset                    = InitPhaseOffset(i);
    WirelessNetwork.WirelessNode(i).Soft_Offset_Record             = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
    WirelessNetwork.WirelessNode(i).Soft_Offset_Record(1)          = WirelessNetwork.WirelessNode(i).Soft_Offset;

    WirelessNetwork.WirelessNode(i).rho_v                          = 0.5;
    WirelessNetwork.WirelessNode(i).rho_o                          = 0.5;

    WirelessNetwork.WirelessNode(i).leader_rho_v                   = 0.5;
    WirelessNetwork.WirelessNode(i).leader_rho_o                   = 0.5;

    WirelessNetwork.WirelessNode(i).follower_rho_v                 = 0.5;
    WirelessNetwork.WirelessNode(i).follower_rho_o                 = 0.5;

    WirelessNetwork.WirelessNode(i).BroadcastNumCNT                = 0;
    WirelessNetwork.WirelessNode(i).BroadcastNumCNT_Record         = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
    WirelessNetwork.WirelessNode(i).BroadcastNumCNT_Record(1)      = WirelessNetwork.WirelessNode(i).BroadcastNumCNT;

    WirelessNetwork.WirelessNode(i).MaxMultiHop                    = WirelessNetwork.MaxMultiHop;
    WirelessNetwork.WirelessNode(i).CurrentHop                     = WirelessNetwork.MaxMultiHop;
    WirelessNetwork.WirelessNode(i).Hop_Record                     = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
    WirelessNetwork.WirelessNode(i).Hop_Record(1)                  = WirelessNetwork.WirelessNode(i).CurrentHop;

    WirelessNetwork.WirelessNode(i).ConvergentFlag                 = 0;
    WirelessNetwork.WirelessNode(i).Relative_Offset_TMP            = zeros(length(WirelessNetwork.WirelessNode(i).NeighborList), 1);
    
    WirelessNetwork.WirelessNode(i).NeighborList                   = WirelessNetwork.NodeNeighborList(i).Neighbor;
end

if isfield(WirelessNetwork, 'Alg') && strcmp(WirelessNetwork.Alg, 'MRPCTS')
    %% 配置领导者节点
    % LeaderID = [1, 7, 18, 29, 35];   % 指定 1 号节点为领导者，可改成向量，例如 [1 5 9]
    % LeaderID = [1, 4, 7];
    LeaderID = 4;

    for k = 1:length(LeaderID)
        id = LeaderID(k);
    
        WirelessNetwork.WirelessNode(id).NodeRole = 'Leader';
    
        WirelessNetwork.WirelessNode(id).Hard_Clock_Rate                = 1;
        WirelessNetwork.WirelessNode(id).Hard_Clock_Record(1)           = WirelessNetwork.WirelessNode(id).Hard_Clock_Rate;
    
        WirelessNetwork.WirelessNode(id).Hard_Clock_Rate_Comp           = 1;
        WirelessNetwork.WirelessNode(id).Hard_Clock_Rate_Comp_Record(1) = WirelessNetwork.WirelessNode(id).Hard_Clock_Rate_Comp;
    
        WirelessNetwork.WirelessNode(id).Soft_Clock_Rate                = WirelessNetwork.WirelessNode(id).Hard_Clock_Rate * WirelessNetwork.WirelessNode(id).Hard_Clock_Rate_Comp;
        WirelessNetwork.WirelessNode(id).Soft_Clock_Rate_Record(1)      = WirelessNetwork.WirelessNode(id).Soft_Clock_Rate;
    
        WirelessNetwork.WirelessNode(id).Soft_Offset                    = 0;
        WirelessNetwork.WirelessNode(id).Soft_Offset_Record(1)          = WirelessNetwork.WirelessNode(id).Soft_Offset;
    
        % 领导者完全不相信周围邻居的建议（不论是领导者还是跟随者）
        WirelessNetwork.WirelessNode(id).leader_rho_v                   = 0;
        WirelessNetwork.WirelessNode(id).leader_rho_o                   = 0;
    
        WirelessNetwork.WirelessNode(id).follower_rho_v                 = 0;
        WirelessNetwork.WirelessNode(id).follower_rho_o                 = 0;
    end
    
    fprintf('领导者节点配置完成。LeaderID = %s\n', mat2str(LeaderID));
end

%% 存储算法基本信息
% 定义需提取的配置字段列表
configFields = {'Alg', 'NetType', 'NodeNUM', 'NetDiameter', ...
                'Delay_mu', 'Delay_sigma', 'Delay_fixed', ...
                'ppm', 'MaxEstErroEnable', ...
                'round_CNT_NUM', 'Sync_CNT_NUM', 'SyncCycle', ...
                'MaxMultiHop'};

% 提取存在的字段到 BasicParams
BasicParams = struct();
for i = 1:length(configFields)
    fieldName = configFields{i};
    if isfield(WirelessNetwork, fieldName)
        BasicParams.(fieldName) = WirelessNetwork.(fieldName);
    end
end

% 构建仿真标识符 (文件夹名)
% 格式：算法_网络类型_节点描述_跳数
% 使用 switch-case 结构提高可读性和扩展性
switch WirelessNetwork.NetType
    case 'GridNWK'
        % 确保 NWK 结构体存在，防止报错
        if exist('NWK', 'var')
            nodeDesc = sprintf('%dx%dNodes', NWK.NodeNUM_Row, NWK.NodeNUM_Col);
        else
            error('变量 "NWK" 未定义，无法获取网格行列信息。');
        end
        
    case 'RandomNWK2D'
        % 2D 随机网络的文件夹命名规则
        nodeDesc = sprintf('%dNodes', WirelessNetwork.NodeNUM);

    case 'RandomNWK3D'
        nodeDesc = sprintf('%dNodes', WirelessNetwork.NodeNUM);

    case 'DynamicRandomNWK3D'
        nodeDesc = sprintf('Dyn%dNodes', WirelessNetwork.NodeNUM);

    otherwise
        % 处理未知网络类型的默认命名
        nodeDesc = 'UnknownNodes';
        warning('未知的 NetType: %s，使用默认命名规则。', WirelessNetwork.NetType);
end

% 拼接最终文件夹名称
simNameStr = sprintf('%s_%s_%s_%dHops_%dRound_%dSyncNum', ...
                     WirelessNetwork.Alg, ...
                     WirelessNetwork.NetType, ...
                     nodeDesc, ...
                     WirelessNetwork.MaxMultiHop, ...
                     WirelessNetwork.round_CNT_NUM,...
                     WirelessNetwork.Sync_CNT_NUM);

% 创建目录结构
topDir = 'SimData';                 
fullDir = fullfile(topDir, simNameStr); 

% 仅当文件夹不存在时创建 (包含父级目录)
if ~exist(fullDir, 'dir')
    [status, msg] = mkdir(fullDir);
    if ~status
        error('创建文件夹失败: %s', msg);
    end
end

% 保存数据
fileName = 'BasicParams.mat';
fullFilePath = fullfile(fullDir, fileName);

try
    save(fullFilePath, 'BasicParams');
    fprintf('[System] 基本参数已保存至:\n  -> %s\n', fullFilePath);
catch ME
    fprintf('[Error] 保存失败: %s\n', ME.message);
end

%% 进行 WirelessNetwork.round_CNT_NUM 轮实验
for round_idx = 1 : WirelessNetwork.round_CNT_NUM
    WirelessNetwork.CurrentRound = round_idx;
    fprintf('>> Processing Trial: %d / %d \n', round_idx, WirelessNetwork.round_CNT_NUM);

    % ---------------------------------------------------------------------
    % 单轮仿真：执行 WirelessNetwork.Sync_CNT_NUM 次同步迭代
    % ---------------------------------------------------------------------
    for sync_idx = 1 : WirelessNetwork.Sync_CNT_NUM
        WirelessNetwork.CurrentSynNum = sync_idx;

        % 动态网络拓扑刷新
        if strcmp(WirelessNetwork.NetType, 'DynamicRandomNWK3D')
            % 提取当前时刻 (sync_idx) 的网络快照
            currentTopo = WirelessNetwork.DynamicTopoSeq{sync_idx};
            % 将新邻居关系覆盖至当前状态估计器中
            for i = 1 : WirelessNetwork.NodeNUM
                WirelessNetwork.WirelessNode(i).NeighborList = currentTopo(i).Neighbor;
            end
        end

        % 时间同步：ATS
        WirelessNetwork = ATS(WirelessNetwork);

        % 时间同步：MACTS
        WirelessNetwork = MACTS(WirelessNetwork);

        % 时间同步：MRPCTS
        WirelessNetwork = MRPCTS(WirelessNetwork);

        % 自适应网络跳数调节
        WirelessNetwork = HopController(WirelessNetwork);

        % 模拟物理时钟按 WirelessNetwork.SyncCycle 秒步进
        WirelessNetwork = ClockShift(WirelessNetwork);
        
        % randn 会生成均值为 0，方差为 1 的标准正态分布随机数
        % 某节点相位变化对收敛目标的影响
        % 示例 1：在第 400 次同步时，让 5 号节点的相位突然增加 100 秒
        randomPhaseDist = randn * 5e+2;
        WirelessNetwork = InjectDisturbance(WirelessNetwork, 50, 2, 'Phase', randomPhaseDist);
        WirelessNetwork = InjectDisturbance(WirelessNetwork, 250, 6, 'Phase', randomPhaseDist);
        WirelessNetwork = InjectDisturbance(WirelessNetwork, 350, 3, 'Phase', randomPhaseDist);

        % 示例 2：在第 600 次同步时，让 2 号和 3 号节点的硬件时钟速率发生 10 ppm 的漂移突变
        randomRateDist = randn * 10e-6;
        WirelessNetwork = InjectDisturbance(WirelessNetwork, 100, [2, 3], 'Rate', randomRateDist);
        WirelessNetwork = InjectDisturbance(WirelessNetwork, 200, [5, 6], 'Rate', randomRateDist);
        WirelessNetwork = InjectDisturbance(WirelessNetwork, 300, 3, 'Rate', randomRateDist);
    end

    % ---------------------------------------------------------------------
    % 保存单轮仿真的训练信息
    % ---------------------------------------------------------------------
    trainData = struct();
    
    trainData.NodeRole = {WirelessNetwork.WirelessNode.NodeRole};
    
    trainData.IsLeader = strcmp(trainData.NodeRole, 'Leader'); 

    trainData.Hard_Clock_Record             = reshape([WirelessNetwork.WirelessNode.Hard_Clock_Record], ...
                                                WirelessNetwork.Sync_CNT_NUM + 1, WirelessNetwork.NodeNUM)';

    trainData.Hard_Clock_Rate_Comp_Record   = reshape([WirelessNetwork.WirelessNode.Hard_Clock_Rate_Comp_Record], ...
                                                WirelessNetwork.Sync_CNT_NUM + 1, WirelessNetwork.NodeNUM)';
    
    trainData.Soft_Clock_Rate_Record        = reshape([WirelessNetwork.WirelessNode.Soft_Clock_Rate_Record], ...
                                                WirelessNetwork.Sync_CNT_NUM + 1, WirelessNetwork.NodeNUM)';
    
    trainData.Soft_Offset_Record            = reshape([WirelessNetwork.WirelessNode.Soft_Offset_Record], ...
                                                WirelessNetwork.Sync_CNT_NUM + 1, WirelessNetwork.NodeNUM)';
                                    
    trainData.BroadcastNumCNT_Record        = reshape([WirelessNetwork.WirelessNode.BroadcastNumCNT_Record], ...
                                                WirelessNetwork.Sync_CNT_NUM + 1, WirelessNetwork.NodeNUM)';
    
    trainData.Hop_Record                    = reshape([WirelessNetwork.WirelessNode.Hop_Record], ...
                                                WirelessNetwork.Sync_CNT_NUM + 1, WirelessNetwork.NodeNUM)';

    if ~isempty(fieldnames(trainData))
        fileName = sprintf('trainData_%d_%d.mat', WirelessNetwork.round_CNT_NUM, round_idx);
        fullFilePath = fullfile(fullDir, fileName);

        try
            save(fullFilePath, 'trainData', '-v7.3');
            fprintf('数据已保存: %s\n', fileName); 
        catch ME
            warning('文件保存失败: %s. 错误: %s', fileName, ME.message);
        end
    else
        warning('本轮 (%d) 未提取到数据，跳过保存。', round_idx);
    end

    % ---------------------------------------------------------------------
    % 信息重置
    % ---------------------------------------------------------------------
    if round_idx < WirelessNetwork.round_CNT_NUM
        for i = 1 : WirelessNetwork.NodeNUM
            WirelessNetwork.WirelessNode(i).Hard_Clock_Rate                = SkewHW(i);
            WirelessNetwork.WirelessNode(i).Hard_Clock_Record              = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
            WirelessNetwork.WirelessNode(i).Hard_Clock_Record(1)           = WirelessNetwork.WirelessNode(i).Hard_Clock_Rate;
        
            WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp           = 1;
            WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp_Record    = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
            WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp_Record(1) = WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp;
        
            WirelessNetwork.WirelessNode(i).Soft_Clock_Rate                = WirelessNetwork.WirelessNode(i).Hard_Clock_Rate * WirelessNetwork.WirelessNode(i).Hard_Clock_Rate_Comp;
            WirelessNetwork.WirelessNode(i).Soft_Clock_Rate_Record         = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
            WirelessNetwork.WirelessNode(i).Soft_Clock_Rate_Record(1)      = WirelessNetwork.WirelessNode(i).Soft_Clock_Rate;
        
            WirelessNetwork.WirelessNode(i).Soft_Offset                    = InitPhaseOffset(i);
            WirelessNetwork.WirelessNode(i).Soft_Offset_Record             = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
            WirelessNetwork.WirelessNode(i).Soft_Offset_Record(1)          = WirelessNetwork.WirelessNode(i).Soft_Offset;
        
            WirelessNetwork.WirelessNode(i).rho_v                          = 0.5;
            WirelessNetwork.WirelessNode(i).rho_o                          = 0.5;
        
            WirelessNetwork.WirelessNode(i).leader_rho_v                   = 0.5;
            WirelessNetwork.WirelessNode(i).leader_rho_o                   = 0.5;
        
            WirelessNetwork.WirelessNode(i).follower_rho_v                 = 0.5;
            WirelessNetwork.WirelessNode(i).follower_rho_o                 = 0.5;
        
            WirelessNetwork.WirelessNode(i).BroadcastNumCNT                = 0;
            WirelessNetwork.WirelessNode(i).BroadcastNumCNT_Record         = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
            WirelessNetwork.WirelessNode(i).BroadcastNumCNT_Record(1)      = WirelessNetwork.WirelessNode(i).BroadcastNumCNT;
        
            WirelessNetwork.WirelessNode(i).MaxMultiHop                    = WirelessNetwork.MaxMultiHop;
            WirelessNetwork.WirelessNode(i).CurrentHop                     = WirelessNetwork.MaxMultiHop;
            WirelessNetwork.WirelessNode(i).Hop_Record                     = zeros(WirelessNetwork.Sync_CNT_NUM + 1, 1);
            WirelessNetwork.WirelessNode(i).Hop_Record(1)                  = WirelessNetwork.WirelessNode(i).CurrentHop;
        
            WirelessNetwork.WirelessNode(i).ConvergentFlag                 = 0;
            WirelessNetwork.WirelessNode(i).Relative_Offset_TMP            = zeros(length(WirelessNetwork.WirelessNode(i).NeighborList), 1);
            
            WirelessNetwork.WirelessNode(i).NeighborList                   = WirelessNetwork.NodeNeighborList(i).Neighbor;
        end

        if isfield(WirelessNetwork, 'Alg') && strcmp(WirelessNetwork.Alg, 'MRPCTS')
            for k = 1:length(LeaderID)
                id = LeaderID(k);
    
                WirelessNetwork.WirelessNode(id).Hard_Clock_Rate               = 1;
                WirelessNetwork.WirelessNode(id).Hard_Clock_Record(1)          = WirelessNetwork.WirelessNode(id).Hard_Clock_Rate;
    
                WirelessNetwork.WirelessNode(id).Hard_Clock_Rate_Comp          = 1;
                WirelessNetwork.WirelessNode(id).Hard_Clock_Rate_Comp_Record(1)= WirelessNetwork.WirelessNode(id).Hard_Clock_Rate_Comp;
    
                WirelessNetwork.WirelessNode(id).Soft_Clock_Rate               = WirelessNetwork.WirelessNode(id).Hard_Clock_Rate * WirelessNetwork.WirelessNode(id).Hard_Clock_Rate_Comp;
                WirelessNetwork.WirelessNode(id).Soft_Clock_Rate_Record(1)     = WirelessNetwork.WirelessNode(id).Soft_Clock_Rate;
    
                WirelessNetwork.WirelessNode(id).Soft_Offset                   = 0;
                WirelessNetwork.WirelessNode(id).Soft_Offset_Record(1)         = WirelessNetwork.WirelessNode(id).Soft_Offset;
    
                WirelessNetwork.WirelessNode(id).leader_rho_v                  = 0;
                WirelessNetwork.WirelessNode(id).leader_rho_o                  = 0;
    
                WirelessNetwork.WirelessNode(id).follower_rho_v                = 0;
                WirelessNetwork.WirelessNode(id).follower_rho_o                = 0;
            end
        end
    end
end