Octave中高效递归计算数组值（无需for循环）

针对你给出的递推公式：

A(n) = A(n-1) * e^(-t(n)) + B(n) * (1 - e^(-t(n)))
初始条件：A(1) = 0

可以通过以下两种高效方法避免for循环：

方法一：利用累积乘积（cumprod）展开递推式 #

递推式展开后可转化为累加形式，借助Octave/Matlab的cumprod函数计算累积乘积，实现向量化运算：

t = 1:1000;
B = (1:1000).^2;
A = zeros(1, 1000);
A(1) = 0;

% 定义系数k(n) = e^(-t(n))
k = exp(-t);

% 构造乘积项：对应每个B(m)*(1-k(m))后续的k(p)累积乘积
if length(k) >= 3
    % 从k(3)到k(end)反向计算累积乘积，再反转回来
    prod_part = cumprod(k(end:-1:3));
    prod_part = [prod_part, 1];  % 最后添加1，对应最后一项无后续乘积
    prod_part = prod_part(end:-1:1);
else
    prod_part = [1];  % 处理t长度为2的边界情况
end

% 计算每个累加项并求和
terms = B(2:end) .* (1 - k(2:end)) .* prod_part;
A(2:end) = cumsum(terms);

plot(t, A)

方法二：使用filter函数处理线性递推 #

Octave和Matlab的filter函数原生支持线性递推计算，完全适配这类一阶递推场景：

t = 1:1000;
B = (1:1000).^2;
A = zeros(1, 1000);
A(1) = 0;

k = exp(-t);
% 递推式中的常数项C(n) = B(n)*(1 - k(n))
C = B(2:end) .* (1 - k(2:end));

% 构造filter的系数：对应递推式 A(n) - k(n)*A(n-1) = C(n)
a = [1, -k(2:end)];
b = [1];

% 计算递推结果，初始状态zi设为A(1)=0
A(2:end) = filter(b, a, C, 0);

plot(t, A)

两种方法都能得到和for循环完全一致的结果，且向量化运算的效率远高于循环，尤其在数组长度较大时优势明显。

内容的提问来源于stack exchange，提问作者Zoltan K.