Test of fits to conic curves

Different asphericities with constant radius (R=20000)

>> x=-3000:100:3000;
>> p=-50:10:10;
>> for c_p=1:length(p)
Z=ZConica_otraformula(x,20000,p(c_p));
conica=ajustaconica([x(:) Z(:)],[],1);
p_fit(c_p)=conica(2);
end
>> plot(p,p_fit,’.-‘)
>> xlabel(‘Real p’)
>> ylabel(‘Fitted p’)

Same, with noise:

>> for c_p=1:length(p)
Z=ZConica_otraformula(x,20000,p(c_p))+randn(1,length(x))*2;
conica=ajustaconica([x(:) Z(:)],[],1);
p_fit(c_p)=conica(2);
end
plot(p,p_fit,’.-‘)

Advertisement