%E_G(100,:)
%Bestfi(100) CanshuDelta=[3 -3  -0.01 -0.0035 -2]-x_canshu;    
%load E_G E_G
%load Bestfi Bestfi   
%load    f_central    f_central

G=100
figure(1)
plot(E_G(:,1))
hold on
k=1:100;
plot(k,3,'r')
title('δH')
ylabel('單位（m）')
xlabel('代數(k)')
figure(2)
plot(E_G(:,2))
hold on

k=1:100;
plot(k,-3,'r')
title('δr1')
ylabel('單位（m）')
xlabel('代數(k)')
figure(3)
plot(E_G(:,3))
hold on
k=1:100;
plot(k,-0.01,'r')
title('δB')
ylabel('單位（m）')
xlabel('代數(k)')

figure(4)
plot(E_G(:,4))
hold on
k=1:100;plot(k,-0.0035,'r')
title('δθb ')
ylabel('單位（rad）')
xlabel('代數(k)')
figure(5)
plot(E_G(:,5))
hold on
k=1:100;
plot(k,-2,'r')
title('δ?( ')
ylabel('單位（rad）')
xlabel('代數(k)')
figure(6)
plot(-log(1-Bestfi)./2,'b')
hold on
plot(-log(1-f_central)./2,'r')
ylabel('單位（m）')
xlabel('代數(k)')
title('各定標點高程差的平方和')
Bestfi
E_G