rls.rar

  • 李sira
    了解作者
  • matlab
    开发工具
  • 7KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 10 积分
    下载积分
  • 0
    下载次数
  • 2020-11-30 10:46
    上传日期
做回归用的代码,非常好用,亲测可用,非常好使,欢迎下载
rls.rar
  • rls
  • Data1.xlsx
    8.7KB
  • demo.m
    427B
  • rls.m
    835B
  • kalman.m
    460B
内容介绍
function [thetahat,xhat]=rls(x,y,N,lambda) % [thetahat,xhat]=rls(x,y,N,lambda) % % x - 完整序列 % y - 延迟序列 % N - 参数维度 % lambda - 遗忘因子 % thetahat - theta参数 % xhat - 估计值输出 % rls 递归最小二乘法 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% M = length(y); % 初始化 p = 10000*eye(N+1); thetahat = zeros(M+1,N+1); xhat = zeros(1,M+1); for n=1:M % 生成Y(n),设置Y中空值位为0; Y = zeros(N+1,1); for k = 1:N+1 if (n-k+1>0) Y(k) = y(n-k+1); else Y(k) = 0; end end % 对 x值进行估计 xhat(n+1) = Y'*thetahat(n,:)'; % 更新 K k = p*Y/(lambda+Y'*p*Y); % 更新 P p = (p-k*Y'*p)/lambda; % 更新thetahat thetahat(n+1,:)= thetahat(n,:) + k'*(x(n)-xhat(n+1)); end % 输出结果 thetahat=thetahat(2:M+1,:); xhat = xhat(2:end);
评论
    相关推荐