mmeasured.m

针对卡尔曼滤波的原创改进算法

% 功能          获得测量曲线
% 作者          刘宗尧
% 日期          2006-11-23
% 最后修改日期   2006-11-13
% 版本          v1.0


function y=mmeasured(y1,R)

%  y1  真实值
%  R   测量误差



f=R*randn(size(y1));

% a=1.2
% for i=1:length(x)
%     y(i)=a^(-2i-2)*randn;
% end

%-----------------------------------
%---随着时间的增加，噪声减小  CH(i)=1/sqrt(i)
%  CH=zeros(size(y1));
%  for i=1:length(y1)
%  CH(i)=i.^0.1;
%  end
%  f=f.*CH;
%  --------------------------
y=f+y1;