function e = creat_e(k,t)

%*********************************
% ** 鐘差模型 待查
% **********************************

b0=30;%單位是m
b1=0.3;%單位是m/s
b2=(1.0e-006)*3;%單位是m/s^2
e1=normrnd(0,0.3);
e=(b0+b1*t+b2*t^2+e1)*ones(k,1);

% e = zeros(k,1);