How to fit a defined function?

I want to fit a custom function to my experimental data. For simplicity, I have some arbitary x & y values and a very simple linear function. I write the following code:

clear all;

x = [1 2 3 4 5]'; % x data
y = [.8 4 10 18 23]'; % y data

Y = lsqcurvefit(fun,1,x,y) % fitting function 'fun' defined below to find parameter 'a'

function y = fun(a,x)
y = a*x; % just a simple function for example, in actual problem, it is a long complicated function with various parameters
end

I get the following error: 'Not enough input arguements.'

I know this particular simple function can be defined as anonymous function and be fitted but I don't it that way.

 采纳的回答

The ‘fun’ function must be presented to lsqcurvefit as a function handle using the ‘@’ operator —
x = [1 2 3 4 5]'; % x data
y = [.8 4 10 18 23]'; % y data
Y = lsqcurvefit(@fun,1,x,y) % fitting function 'fun' defined below to find parameter 'a'
Local minimum found. Optimization completed because the size of the gradient is less than the value of the optimality tolerance.
Y = 4.1055
function y = fun(a,x)
y = a*x; % just a simple function for example, in actual problem, it is a long complicated function with various parameters
end
See What Is a Function Handle? for details.
.

2 个评论

Thank you Star Strider for your answer. It solved my problem.
As always, my pleasure!

请先登录,再进行评论。

更多回答(2 个)

VBBV
VBBV 2022-11-10
编辑:VBBV 2022-11-10

1 个投票

To define a and then call in function

2 个评论

Then use @
Y = lsqcurvefit(@fun,1,x,y)
Thank you for your respone, it helped me.

请先登录,再进行评论。

A simpler way for this problem, but I guess your "real" model is more complicated:
x = [1 2 3 4 5]'; % x data
y = [.8 4 10 18 23]'; % y data
a = x\y
a = 4.1055

1 个评论

Thank you Torsten for this intelligent answer, I really like your perspective.
Yes, my actual problem is not so simple, it is a transfer matrix model.

请先登录,再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息

产品

版本

R2021b

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by