MATLAB 代码错误解决：索引无效问题

在使用MATLAB进行遗传算法优化方向图时，代码出现了“位置 3 的索引无效。数组索引必须为正整数或逻辑值。”的错误。

错误原因

该错误通常出现在对数组进行循环访问时，索引值超出了数组的范围。

本例中，错误发生在第15行的循环中，代码为 f(i,:)，其中i是从1到NP的整数。然而，在第11行初始化f的时候，第一行已经被赋值为fff了，因此在循环的时候i从2开始才是正确的。

改正方案

为了解决该错误，需要修改循环的起始值，从2开始遍历f数组。

正确代码

clear all;
close all;
clc

%输入基础数据

N = 33; % x方向阵元个数
M = 53; % y方向阵元个数
f = 12; % 信号频率 单位GHz
c = 2.99792 * 10^8; % 光速 单位m/s
lambda = c / (f * 10^9); % 波长，需要将GHz转换为Hz
lambdag = lambda * 0.8; % 波导中波长
dx = 7.5*0.001; % x方向单元间距
dy = 3.5*0.001; % y方向单元间距
phi = linspace(-90, 90, 181); % 方位角范围
theta = linspace(-90, 90, 181); % 俯仰角范围
theta0 = 30 * pi/180; % 目标俯仰角
phi0 = 90 * pi/180; % 预先设定的方向，目标方位角

%生成单元位置矩阵
x = -(N + 1) / 2 * dx + (1:N)' * dx; % x方向阵元位置
y = -(M + 1) / 2 * dy + (1:M)' * dy; % y方向阵元位置
[X, Y] = meshgrid(x, y);
Z = zeros(size(X)); % z方向阵元位置全为0
p = [X(:), Y(:), Z(:)]' ;% 阵元的位置信息，每一列是一个阵元的位置
q = sqrt(sum(p.^2, 1));

%%%%%第1步，初始化参数%%%%%
NP = 50;          %种群数量
Pc = 0.8;         %交叉率
Pm = 0.1;         %变异率
G = 100;          %最大遗传代数
Xs = 10;          %在区间[0,10]求函数的极值
Xx = 0;           %定义域下限
fff = m0(lambda,lambdag,theta0,phi0,p,M,N);  %随机产生第一代随机种群
ff = randi([0 1],NP-1,M*N);
f = cat(1 , fff ,ff);
trace=zeros(1,G);     %创建数组，用于记录每一代，适应度最好的函数值
for i=2:NP
    i
    C = f(i,:)
    %获取方向图
    B = B( M, N ,C, lambda , lambdag , theta , phi, q ,p, theta0, phi0);
    %获取主瓣与副瓣位置
    ppeaks = ppeaks( B, phi0 , theta0 );     %适应度，即目标函数的值
    fit = -ppeaks
end

通过修改循环起始值，代码可以正常运行，并避免了“索引无效”的错误。