lsqcurvefit

Would LSQCURVEFIT work with comparing two complex functions? It looks like the metric for success is min sum {(FUN(X,XDATA)-YDATA).^2}, where it should be min sum {abs(FUN(X,XDATA)-YDATA).^2)} for complex-valued functions, or maybe sum {(FUN(X,XDATA)-YDATA).conj((FUN(X,XDATA)-YDATA))} Thanks Matt

Answer 1

[Andrew Newell]

Since the distance between two complex numbers a+i*b and c+i*d is (a-c)^2 + (b-d)^2, you could treat each complex number as a pair of independent real numbers. Wrap fun in a function rfun that returns the real and complex parts of the output. For example, assuming rfun returns a column vector,

function y = rfun(x)
  y = rfun(x);
  y = [real(y); imag(y)];
end

Similarly, you could replace ydata by

rdata = [real(ydata); imag(ydata)].

Then you have the first two inputs for lsqcurvefit.