X
تبلیغات
رایتل
9 دی 1387 ساعت 01:20 ب.ظ

clc

clear

calib=input('enter calib matrix')

meas=input('enter meas matrix')

n=size(calib,2)/2

l=calib'

syms a b cx cy

for k=1:2:2*n

    x=a*meas(1,k)+b*meas(1,k+1)+cx

    y=-b*meas(1,k)+a*meas(1,k+1)+cy

   A(k,:)=jacobian(x,[a b cx cy])

    A(k+1,:)=jacobian(y,[a b cx cy])

end

A=eval(A)

xcap=inv(A'*A)*(A'*l)

r=l-(A*xcap) 

del.icio.us  digg  newsvine  furl  Y!  smarking  segnalo