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,’.-‘)

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

%d bloggers like this: