MATLAB solve/vpasolve将temp视为变量，求解结果异常问题排查

问题现象 #

通过for循环构建对应N=[3,5,7]的方程组（变量为a1至a7），使用MATLAB的solve或vpasolve求解时出现以下异常：

• 求解结果中缺失a3、a5、a7等目标变量，反而出现temp变量；
• 使用vpasolve指定求解Acou（如a1-a3、a1-a5）后，结果仍依赖temp，无法得到独立的变量解。

错误原因 #

1. temp未在循环内清零：temp作为全局符号变量，每次迭代构建方程时未重置为0，导致后续方程累加了之前的计算值，被MATLAB误判为待求解变量。
2. 方程构建逻辑冗余：偶数j的分支中，anslist(k)设为0但仍被累加到temp，导致方程出现无意义的累加，破坏了方程组的独立性，求解器被迫引入temp作为自由变量。

修正方案 #

1. 循环内初始化temp：每次进入j循环前，将temp重置为0，确保每个方程的累加都是独立的。
2. 优化偶数j的方程构建：针对偶数j的情况，直接跳过无效累加，保持temp为0后构建对应方程。
3. 明确目标变量范围：求解时仅传入目标变量组Acou，避免temp被纳入求解范围。

完整修正代码 #

cb = 1/5.25; % 展弦比倒数
cla = 2*pi; % 升力线斜率
aa = 5*pi/180; % 5度转换为弧度
N = [3 5 7]; % 方程数量列表
syms a1 a2 a3 a4 a5 a6 a7
A = [a1 a2 a3 a4 a5 a6 a7]; % 目标变量组

for i = 1:3
    counter = N(i);
    Acou = A(1:counter);
    eqncou = sym(zeros(1, counter)); % 初始化方程组
    spandiv = (pi/2)/counter;
    
    for j = 1:counter
        temp = 0; % 每次构建方程前重置temp
        for k = 1:counter
            if rem(j,2) ~= 0
                % 奇数j时计算累加项
                term = A(k)*sin(k*spandiv*j) + 0.25*cb*cla*sin(k*spandiv*j)/sin(spandiv*j);
                temp = temp + term;
            end
            % 偶数j时无累加项，temp保持0
        end
        % 构建第j个方程
        eqncou(j) = temp == 0.25*cb*2*pi*aa;
    end
    
    % 求解方程组，指定目标变量为Acou
    Asol = vpasolve(eqncou, Acou);
    disp(['N=', num2str(counter), ' 的求解结果：']);
    disp(Asol);
end

结果验证 #

修正后，求解结果将不再包含temp变量，每个目标变量（a1-a3、a1-a5、a1-a7）都会得到独立的数值解，例如N=3时的结果会输出a1、a2、a3的具体值，无依赖项。

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