How to fit a defined function?

13 次查看(过去 30 天)

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.

采纳的回答

Star Strider
Star Strider 2022-11-10
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 个评论
Sourabh Jain
Sourabh Jain 2022-11-10
Thank you Star Strider for your answer. It solved my problem.

请先登录,再进行评论。

更多回答(2 个)

VBBV
VBBV 2022-11-10
编辑:VBBV 2022-11-10
To define a and then call in function

Torsten
Torsten 2022-11-10
编辑:Torsten 2022-11-10
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 个评论
Sourabh Jain
Sourabh Jain 2022-11-10
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.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Mathematics 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by