% 2floaters each one 2 pts
ls
load obs_float_test.dat
whos;
X=obs_float_test (:,:)
[m,n]=size(X)
x1=obs_float_test (1:2:m/2,4)
y1=obs_float_test (2:2:m/2,4)
x1=x1+1;
y1=y1+1;
plot(x1(1),y1(1),'+')
hold on
plot(x1(4),y1(4),'>')
hold on
plot(x1(1:4),y1(1:4),'-r')
hold on
x2=obs_float_test (m/2+1:2:m,4)
y2=obs_float_test (m/2+2:2:m,4)
x2=x2+1;
y2=y2+1;
plot(x2(1),y2(1),'+')
hold on
plot(x2(4),y2(4),'>')
hold on
plot(x2(1:4),y2(1:4),'-b')
hold on
load r_float_optim25float-filtre15-inc_16.dat;
whos;
r_float_optim25float_filtre15_inc_16(1:end,:)
x3=r_float_optim25float_filtre15_inc_16(1:4:16,4)
y3=r_float_optim25float_filtre15_inc_16(3:4:16,4)
x3=x3+1;
y3=y3+1;
plot(x3(1),y3(1),'+')
hold on
plot(x3(4),y3(4),'>')
hold on
plot(x3(1:4),y3(1:4),'-g')
hold on
x4=r_float_optim25float_filtre15_inc_16(2:4:16,4)
y4=r_float_optim25float_filtre15_inc_16(4:4:16,4)
x4=x4+1;
y4=y4+1;
plot(x4(1),y4(1),'+')
hold on
plot(x4(4),y4(4),'>')
hold on
plot(x4(1:4),y4(1:4),'-y')
hold on
ls;
load uoptim25float-filtre15-inc16.dat;
whos;
X=uoptim25float_filtre15_inc16(:,end-2);
Y=uoptim25float_filtre15_inc16(:,end-1);
U=uoptim25float_filtre15_inc16(:,end);
ls;
load voptim25float-filtre15-inc16.dat;
whos;
V=voptim25float_filtre15_inc16(:,end);
quiver(X,Y,U,V,'k');
