clear
%二維的作勻速直線運動的水面目標的航跡示意圖：

%產生[0 1]間的正態分布的隨機數
a0=19;
C=51329;
x1(1)=47169; 
M=2^16;  
for i=1:5999
	  x1(i+1)=mod(x1(i)*a0+C,M);
end
for i=1:6000
	  fakerand(i)=x1(i)/M;
end
for i=1:500
	  normrnd1(i)=0;
	  for j=(12*(i-1)+1):(12*i)
		  normrnd1(i)=normrnd1(i)+fakerand(j);
      end
	  normrnd1(i)=normrnd1(i)-6;
end

%目標的參數，初始位置
x0=3000;
y0=4000;

v=-200;

Q1=120;
Q2=160;

b=pi/4;

vx=v*sin(b);
vy=v*cos(b);


t0=2;

t=[1:500]+t0;


for i=1:50
    x(i)=x0+vx*(t(i)-t0)+Q1*normrnd1(i);
    y(i)=y0+vy*(t(i)-t0)+Q2*normrnd1(i);
end
for i=1:50
    X(i)=x0+vx*(t(i)-t0);
    Y(i)=y0+vy*(t(i)-t0);
end
plot(x0,y0,'r^',[x0 x],[y0 y],'b.',[x0 X],[y0 Y],'r-')
hold on
legend('起始點','方差為Q1=120;Q2=160的觀測曲線','理論曲線');
title('二維勻速直線運動的連續模型目標軌跡示意圖');
gtext('起始點');
xlabel('X坐標軸');
ylabel('Y坐標軸');
grid on