I’m in a context here that is very much less than clear.
Gauss = @(a,b,c,x) a*exp(-((x-b)/(2*c)).^2);
or in terms that the nlinfit function will understand:
Gfcn = @(p,x) p(1).*exp(-((x-p(2))./(2*p(3))).^2);
Use the ‘Gfcn’ function as your objective function with nlinfit and all should be well.
Your ‘mdl’ function makes absolutely no sense to me, so I won’t essay to correct it.
So just use ‘Gfcn’, and — with the correct initial parameter estimates vector — you should have no problems fitting it to your data, if they are truly described by a Gaussian function.
Note — I do not have your data to test this with, so I am listing it as UNTESTED CODE. It should work.
