待解方程：

% 輸入是 x_d 資料庫影象的 x 座標；y_d 查詢影象的 y 座標；a_d 查詢影象的視角

% 待求未知數為 x_query,初始值設為 [30 0 0];

[email protected](x_query)[(x_query(1)-x_d(1))*tan(x_query(3)-a_d(1))+(y_d(1)-x_query(2));
             (x_query(1)-x_d(2))*tan(x_query(3)-a_d(2))+(y_d(2)-x_query(2));
             (x_query(1)-x_d(3))*tan(x_query(3)-a_d(3))+(y_d(3)-x_query(2))];
x_query_0=[0.3 0 0];
options = optimoptions(@fsolve,'MaxFunctionEvaluations',3000,'MaxIterations',2000);

[Estimated_locazion]= fsolve(F,x_query_0,optimset('Display','iter'));


disp(Estimated_locazion);