ECG-Kit 1.0
(582 bytes)
function f = objfun2(X,C,weight)
% error matrix to minimize: to find direction and point
%
% Rute Almeida 2.Dec.2004
% Last update:
%
% Designed for MATLAB Version R12; tested with MATLAB Version R13
%X=[a ;v]
if nargin<3
weight=ones(size(C));
end
a=NaN*ones(1,3);
v=NaN*ones(1,3);
a(1:size(X,2))=X(1,:);
v(1:size(X,2))=X(2,:);
%aux=[C-[C(:,1) v(2)/v(1)*(C(:,1)-a(1))+a(2) v(3)/v(1)*(C(:,1)-a(1))+a(3)]];
aux=[C(:,1) v(2)/v(1)*(C(:,1)-a(1))+a(2) v(3)/v(1)*(C(:,1)-a(1))+a(3)];
aux=weight.*(C-aux(:,1:size(C,2)));
f=aux(:);
%f =sqrt(sqrt(trace(aux'*aux)));