clear all
close all
clc
r=2;
for i=0:300;
t=0:0.1:10;
x0=10*sin(t);
theta11=pi/4;
iy=(x0+2*randn(size(t))).*(x0*sin(theta11)+2*randn(size(t)));
ix=(x0+2*randn(size(t))).*(x0*cos(theta11)+2*randn(size(t)));
iy1=sum(iy)/(101);
ix1=sum(ix)/(101);
theta10=atan2(iy1,ix1);
x1=cos(pi/4)*i;
y1=x1;
theta21=atan2(y1+4,x1);
iy0=(x0+2*randn(size(t))).*(x0.*sin(theta21)+2*randn(size(t)));
ix0=(x0+2*randn(size(t))).*(x0.*cos(theta21)+2*randn(size(t)));
iy10=sum(iy0)/(101);
ix10=sum(ix0)/(101);
theta20=atan2(iy10,ix10);
wucha2(i+1)=theta21-theta20;
wucha1(i+1)=theta11-theta10;
close all
clc
r=2;
for i=0:300;
t=0:0.1:10;
x0=10*sin(t);
theta11=pi/4;
iy=(x0+2*randn(size(t))).*(x0*sin(theta11)+2*randn(size(t)));
ix=(x0+2*randn(size(t))).*(x0*cos(theta11)+2*randn(size(t)));
iy1=sum(iy)/(101);
ix1=sum(ix)/(101);
theta10=atan2(iy1,ix1);
x1=cos(pi/4)*i;
y1=x1;
theta21=atan2(y1+4,x1);
iy0=(x0+2*randn(size(t))).*(x0.*sin(theta21)+2*randn(size(t)));
ix0=(x0+2*randn(size(t))).*(x0.*cos(theta21)+2*randn(size(t)));
iy10=sum(iy0)/(101);
ix10=sum(ix0)/(101);
theta20=atan2(iy10,ix10);
wucha2(i+1)=theta21-theta20;
wucha1(i+1)=theta11-theta10;