Fit a curve based on its area?

2 次查看（过去 30 天）

Larissa Perez 2017-1-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/320696-fit-a-curve-based-on-its-area

评论： Larissa Perez 2017-1-19

Hi everyone! Basically, I need to estimate a curve based on its integral. I don't know the actual function.

I know my x values will be (0:1000). I know the area when x = 0:30 is 10% of the total area, x = 0:60 is 50% and x = 0:80 is 90%.

I have no idea where to start! I have tried using curve fitting tools but didn't succeed. Can someone help me with that?

Cheers!

13 个评论
显示 11更早的评论隐藏 11更早的评论

David Goodmanson 2017-1-18

Well, the problem now is that you have four conditions, assuming you are confident that the peak is right at x = 60. That means you have to fit an equation with four adjustable parameters. Walter has a solution technique, but so far it assumes a straight polynomial fit. The most obvious four parameter solution that has the shape you are talking about would be (ax + by^2 + cy^3)exp(-gy) [or, because of the nature of the conditions, (x + by^2 + cy^3 + dy^4))exp(-gy) might be needed, I haven't quite figured out which ].

There would be work involved to fit that, especially the coefficient g. If you knew g, which you don't, it would be a lot easier. Someone may have some tricks up his/her sleeve here. But it's a hard enough problem that it makes sense to ask, what is the context here? Are the conditions exact or pretty good but approximate? Is a not so simple many-parameter fit the way to go? Do you have any other information about where this function came from that might indicate what kind of function it is? Is the goal to do interpolation to find the area at other values of x?

David Goodmanson 2017-1-18

编辑：David Goodmanson 2017-1-18

在 MATLAB Online 中打开

Your best bet here may be the log normal distribution which is often used for particle sizes and a lot of other situations. You assume the log of particle size is normally distributed. For a good fit, log of the d50 size should be right in between log of the d10 size and log of the d90 size. So log(60) should be halfway between log(33) and log(95). Those three values are

ans = 3.4965 4.0943 4.5539

and not equally spaced. But fitting those three points with a straight line may be the best you can do. That fit is pretty decent, not super.

Larissa Perez 2017-1-19

I really appreciate your help! Thanks for that. :)

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Walter Roberson 2017-1-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/320696-fit-a-curve-based-on-its-area#answer_250865

编辑：Walter Roberson 2017-1-17

在 MATLAB Online 中打开

one of the solutions is

syms a b c d e F(x) f(x)
F(x) = a*x^4+b*x^3+c*x^2+d*x+e
sol = solve([F(30)-F(0)==10,F(60)-F(0)==50,F(80)-F(0)==90,F(1000)-F(0)==100])
f(x) = subs(diff(F(x),x),{a,b,c,d},{sol.a,sol.b,sol.c,sol.d})

For each additional known point you have, add another term to F(x) and add the point to the solve() expression.

This is based upon the fact that the indefinite integral is taken and becomes an actual function, and that the definite integral is determined by the subtraction of the value of the indefinite integral at the end points. F(x) is acting as the indefinite integral. You can then create simultaneous equations out of it and solve for the coefficients. Then substitute those into the derivative of F(x) to get the original f(x)