zls.m

来自「最小二乘法中的递推最小二乘法」· M 代码 · 共 21 行

function[theta]=ZLS(N,n,In,Out)
%N为计算的数据点数
%n为模型的阶次
%In为输入 的数据向量
%Out为输入 的数据向量
x=rand(N,2*n+1);
for i=1:N
    for j=1:n
        if(i-j>0)
        x(i,j)=-Out(i-j);
        else x(i,j)=0;
        end
    end
    for k=(n+1):(2*n+1)
        if(i-k+n+1>0)
        x(i,k)=In(i-k+n+1);
        else x(i,k)=0;
        end
    end
end
theta=inv(x'*x)*x'*Out';