MATLAB 中的复制求和运算符

我正在尝试在 MATLAB 中对以下时间相关函数进行建模：

其中 gamma、beta 和 delta 是常数。 我不确定进行求和的最佳方法是什么。我可以制作 u(t) 和 y(t) 函数句柄，但 symsum 只能用于符号。我唯一的选择是使用 for 循环手动添加 y(t-i) 和 u(t-m) 吗？到目前为止，我为生成 y(t) 所做的工作如下所示：

N = 100;
y_out = zeros(N,1);
u = zeros(N,1);

for t = 1:1:N
u(t) = cos(beta*t);
end

y_out(1) = NARMA(beta,delta,gamma,true);
for t=2:N
    y_out(t) = NARMA(beta,delta,gamma);
end


function y_out = NARMA(beta,delta,gamma,first_call)
    % Set first_call to true only the first time you call the function
    persistent y
    if nargin>3 && first_call
        y = 0;
    end

    for i=1:1:m
    y = y + gamma*u + delta;
    end

    y_out = y;

end

这是一个不需要循环的可能解决方案（未经测试）：

t = 1:N;
u = @(t) cos(beta * t);
y = cumsum(gamma * u(t - m) .* u(t) + delta);
y(m+1:end) = y(m+1:end) - y(1:end-m);

EDIT:

事实证明，上面的解决方案是错误的。我认为这是正确的：

u = @(t) cos(beta*t);
y = zeros(N+m, 1);
for t = 1:N
    y(t+m) = sum(y(t:t+m-1)) + gamma * u(t-m) * u(t) + delta;
end
y = y (m+1:end);

EDIT2:

通过进一步简化，求和可以简化为简单的 for 循环：

u = @(t) cos(beta*t);
y = zeros(N+m, 1);
e = gamma * u((1:N)-m) .* u(1:N) + delta;
de = diff(e);
y(1+m) = e(1);
for t = 2:N
    y(t+m) = 2 * y(t+m-1) - y(t-1) + de(t-1);
end
y = y (m+1:end);