mmeasured.asv

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

% 功能          获得测量曲线
% 作者          刘宗尧
% 日期          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)=(R-1/(i+1);
 end
 f=f.*CH;
 %--------------------------
y=f+y1;