Home > bioelectromagnetism > elec_fit_ellipse_optim.m

elec_fit_ellipse_optim

PURPOSE ^

elec_fit_ellipse_optim - Optimization for elec_fit_ellipse.m

SYNOPSIS ^

function [f] = elec_fit_ellipse_optim(r, X, Y, Z, xo, yo, zo)

DESCRIPTION ^

 elec_fit_ellipse_optim - Optimization for elec_fit_ellipse.m

 Called from elec_fit_ellipse.m

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [f] = elec_fit_ellipse_optim(r, X, Y, Z, xo, yo, zo)
0002 
0003 % elec_fit_ellipse_optim - Optimization for elec_fit_ellipse.m
0004 %
0005 % Called from elec_fit_ellipse.m
0006 %
0007 
0008 % $Revision: 1.1 $ $Date: 2004/11/12 01:32:33 $
0009 
0010 % Licence:  GNU GPL, no implied or express warranties
0011 % History:  02/2002, Darren.Weber_at_radiology.ucsf.edu
0012 %
0013 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0014 
0015 % r is a 3x1 vector of radius values for each x,y,z axis component of ellipse
0016 %
0017 % equation of ellipsoid with center (xo,yo,zo) and radius for each axis (x,y,z) = (a,b,c):
0018 % (( x - xo )^2 / a^2) + (( y - yo )^2 / b^2) + (( z - zo )^2 / c^2) = 1
0019 %
0020 % This function below creates a scalar value to
0021 % return to the fminsearch function in elec_fit_sphere.
0022 
0023 E = ( (X-xo).^2 )/r(1).^2 + ((Y-yo).^2)/r(2).^2 + ((Z-zo).^2)/r(3).^2  - 1;
0024 
0025 f = sum( E .* E );  % sum of squares returned
0026 
0027 return

Generated on Mon 15-Aug-2005 15:36:19 by m2html © 2003