actual_fun.m

来自「简单的最小二乘逼近算法」· M 代码 · 共 28 行

function y_return = actual_fun(k)
%真实模型
%输入参数：
%输出参数：真实模型计算出的准确值.
%调用关系：被以下函数调用：递推最小二乘法系统辨识sys_identify

%初值
y(1) = 0;
y(2) = 1;

%参数
a1 = 2;
a2 = 1.3;
b0 = 0.4;
b1 = 0.88;
b2 = 2.2;

if k == 1
    y_return = y(1);
elseif k == 2
    y_return = [y(1),y(2)];
else
    for i = 3:k
        y(i) = -a1 * y(i-1) - a2 * y(i-2) + b0 * u(i) + b1 * u(i-1) + b2 * u(i-2) + xi;
    end
    y_return = y;
end