optimalfun.asv

这是我用Delphi和Matlab写的一个程序

function Fun = OptimalFun(IniVector)
%  OPTIMALFUN  Return a vector of values: [f(x1), f(x2), ..., f(xn)]
%  IniVector: 七个参数(α、a、b、 eix、 eiy)的初始值组成的向量
%  [α, a, b, e0x, e0y, e1x, e1y];
%  Fun: [f(x1), f(x2), ..., f(xn)]
%
%  Ran Liu: [e-mail]
%  College of Computer Science, Chongqing University
%  Panovasic Technology Co.,Ltd

global I0P
global I1P

N = size(I0P, 2);  %  获得匹配点的对数
Fun = zeros(1, N);  %  将Fun变成1×N的零矩阵
%  计算Fun中的f(xi)
[H0, H1] = GetHiFrom7Params(IniVector);
Fun = H0 * I0P - H1 * I1P
for i = 1 : N
    Fun(i) = (H1(3, :) * I1P(:, i) * H0(2, :) - H1(2, :) * I1P(:, i) * H0(3, :)) * I0P(:, i);
end;