clc,clear % 生成初始点 x=rand*100; y=rand*100; jiedian = [x,y]; % 初始化参数 p = 1e6; n = 0; % 计算初始覆盖率 area = mengtekaluo(jiedian); cover = area / 1e5; new_cover = []; new_cover = [new_cover;cover]; % 记录覆盖率的变化 cover_history = []; cover_history = [cover_history;cover]; % 随机生成点，进行循环判断覆盖率是否达到目标值 while new_cover < 0.95 flag1 = 1; x_new = rand * 100; y_new = rand * 100; jiedian = [jiedian; x_new, y_new]; for i = 1:p % 对于要投放的总共p个点 n = 0; flag = 0; px = rand*100; % 随机生成该点的横坐标 py = rand*100; % 随机生成该点的纵坐标 [a,~] = size(jiedian); for j = 1 : a if (px-jiedian(j,1))^2 + (py-jiedian(j,2))^2 < 100 % 检查点与已有节点的距离是否小于10 flag = flag + 1; end end if flag == a n = n + 1; end end area = (n/p)*10000; new_cover = area / 1e5; if (new_cover - cover) / cover < 0.1 flag1 = 0; % 用flag判断覆盖率变化 new_cover = cover; end % 记录覆盖率的变化 while flag1 cover_history = [cover_history; cover]; end end % 绘制覆盖率变化曲线 hold on; figure; plot(1:length(cover_history), cover_history); xlabel('生成节点数量'); ylabel('区域覆盖率'); title('区域覆盖率变化');

根据你提供的代码，我看到了问题所在。在while循环中，判断结束的条件应该是`new_cover >= 0.95`，而不是`new_cover < 0.95`。此外，在记录覆盖率变化时，需要在循环内部更新`cover`的值。

下面是修改后的代码段：

clc,clear

% 生成初始点
x=rand*100;
y=rand*100;
jiedian = [x,y];

% 初始化参数
p = 1e6;
n = 0;

% 计算初始覆盖率
area = mengtekaluo(jiedian);
cover = area / 1e5;
new_cover = [];
new_cover = [new_cover;cover];

% 记录覆盖率的变化
cover_history = [];
cover_history = [cover_history;cover];

% 随机生成点，进行循环判断覆盖率是否达到目标值
while new_cover >= 0.95
    flag1 = 1;
    x_new = rand * 100;
    y_new = rand * 100;               
    
    jiedian = [jiedian; x_new, y_new];
    for i = 1:p     % 对于要投放的总共p个点
        n = 0;
        flag = 0;
        px = rand*100;    % 随机生成该点的横坐标
        py = rand*100;    % 随机生成该点的纵坐标   
        [a,~] = size(jiedian);
        for j = 1 : a
            if (px-jiedian(j,1))^2 + (py-jiedian(j,2))^2 < 100      % 检查点与已有节点的距离是否小于10
                flag = flag + 1;
            end
        end
        if flag == a            
            n = n + 1;         
        end
    end
    area = (n/p)*10000;
    new_cover = area / 1e5;
    
    if (new_cover - cover) / cover < 0.1
        flag1 = 0;                % 用flag判断覆盖率变化
        new_cover = cover;
    end
    
    % 更新覆盖率和记录覆盖率的变化
    cover = new_cover;
    cover_history = [cover_history; cover];       
end

% 绘制覆盖率变化曲线
hold on;
figure;
plot(1:length(cover_history), cover_history);
xlabel('生成节点数量');
ylabel('区域覆盖率');
title('区域覆盖率变化');

这样修改后，代码应该能够正确计算和记录区域的覆盖率变化了。