问题根源 #

你遇到的核心问题是Simulink的运行逻辑和你预期的不一样：Simulink是逐时间步迭代运行的，默认情况下Matlab Function模块在每个时间步只会接收到单个采样点（标量），而非整个输入信号数组。所以length(y)返回1是正常的——它只看到了当前时刻的单个值，不是你以为的完整信号。

而且RLS算法本身是递归迭代的，其实根本不需要提前知道信号总长度N，每一步用当前的u和y就能完成参数更新。

解决方案 #

1. 适配Simulink逐步运行的RLS实现（推荐）

修改代码，用persistent变量保存跨时间步的状态（比如协方差矩阵C、迭代步数），去掉对总N的依赖，让算法逐步更新：

function [P, current_step] = fcn(u,y)
    persistent C step_count
    % 初始化状态（仅第一次运行时执行）
    if isempty(C)
        sigma = 1;
        C = sigma*eye(2); % 初始协方差矩阵
        step_count = 0;
    end
    
    step_count = step_count + 1;
    current_step = step_count; % 如果你需要跟踪当前迭代步数，替代原来的N
    
    lamda = 1;
    z = [y; u];
    % 补全RLS核心更新公式（你原来的代码在这里截断了）
    K = (C*z)/(lamda + z'*C*z);
    C = (C - K*z'*C)/lamda;
    P = ...; % 根据你的RLS需求计算最终输出P
    
end

• persistent变量会在模块多次运行时保留值，实现递归状态的传递。
• 这种方式完全符合Simulink的实时运行逻辑，不管信号多长都能持续迭代。

2. 一次性输入完整信号（仅适合离线仿真）

如果你确实需要一次性处理整个信号，需要调整模型设置：

• 将输入信号改为一维数组（比如用Constant模块导入Workspace里的完整信号数组）。
• 打开Matlab Function模块的端口设置，把y的维度指定为:（即任意长度的一维数组），这样模块会一次性接收整个信号，此时length(y)就能返回正确的总长度。

但这种方式只适合离线批量处理，不支持实时或连续运行的仿真场景。

关键提醒 #

RLS算法的设计初衷就是逐步递归更新，提前知道总长度N并不是必要条件。适配Simulink的逐步运行逻辑，才是最合理的解决方式。

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