Matlab for循环代码问题排查：统计功率开启状态下变动超1W的次数

核心问题定位 #

• 索引越界错误：原循环范围设置为n=1:length(time)，当n遍历到数组最后一个元素时，访问power(n+1)会超出数组边界，直接导致程序运行中断。
• 状态判断逻辑缺失：原判断仅校验相邻采样对的前一个点状态为on，未校验后一个点的状态，会把on状态跳转off状态的跨段采样对误计入统计。即使修复了越界问题，按原逻辑会将3次on到off的功率跳变也算入，最终结果虚高到7，和预期的4次相差较大。
• 边界判断逻辑本身是正确的：原代码设置的power_change > 1符合“变动幅度超过1W”的要求（等于1W不计数），无需调整。

修正后的可运行代码 #

fid = fopen('power.txt');
% 跳过文件第一行的表头
fgetl(fid); 
data = textscan(fid, '%d %f %s');
fclose(fid);

time = data{1};
power = data{2};
status = data{3};

match_count = 0;
% 循环范围缩小到倒数第二个采样点，避免越界
for n = 1:length(time)-1
    delta_power = abs(power(n) - power(n+1));
    % 仅统计两个相邻点均处于on状态，且功率变动幅度超过1W的情况
    if status(n) == "on" && status(n+1) == "on" && delta_power > 1
        match_count = match_count + 1;
    end
end

disp(['统计得到的变动次数为：', num2str(match_count)]);

实现效果验证 #

逐段匹配采样数据后，符合要求的变动点共4处，和预期值完全一致：

• t=12（功率9W）到t=13（功率6W）：变动幅度3W，满足条件
• t=13（功率6W）到t=14（功率4W）：变动幅度2W，满足条件
• t=25（功率9W）到t=26（功率5W）：变动幅度4W，满足条件
• t=40（功率9W）到t=41（功率5W）：变动幅度4W，满足条件

可选的高效实现（向量化操作） #

Matlab支持向量化运算替代循环，代码更简洁且运行效率更高，最终输出结果和循环版本完全一致：

fid = fopen('power.txt');
fgetl(fid); 
data = textscan(fid, '%d %f %s');
fclose(fid);

power = data{2};
status = data{3};

% 计算所有相邻点的功率差
delta_power = abs(diff(power));
% 生成有效状态掩码：相邻两个点均为on
valid_on = strcmp(status(1:end-1), "on") & strcmp(status(2:end), "on");
% 统计满足条件的次数
match_count = sum(valid_on & (delta_power > 1));

内容的提问来源于stack exchange，提问作者skycharm