Can I get the mathematical expression of an interpolation?

3 次查看（过去 30 天）
Jaime De La Mota Sanchis
评论： Steven Lord 2021-10-9
Hello everyone. I am working with interpolated data.
I want to obtain the mathematical expression of the interpollating polynomial. I have written
sfX1 = fit([positions(:,1), positions(:,2)],windX_in_the_node_vector,'poly34')
And the output returns the linear model poly which is what I am looking for, f(x,y)=a*x^3+b*y^3+...
However, I would be interested in using other kinds of interpolations. I have tried linear, lowess, linearinterp, cubicinterp and nearest and instead of poly34. Unfortunately, instead of a polynomial such as the one in poly34, the code returns
Linear interpolant:
sfX2(x,y) = piecewise linear surface computed from p
Coefficients:
p = coefficient structure
and so on. Is there a way for me to obtain the mathematical expressions as the ones in poly34?
Best regards.
Jaime.
  2 个评论
Walter Roberson
Walter Roberson 2021-10-8
For piecewise fits, the mathematical expressions would have to be written in terms of heaviside() or similar functions. Would that even be useful to you?
At the moment I am not aware of any function to convert piecewise fits to heaviside, but I believe it would be possible to write such a function... but is it worth doing?
Jaime De La Mota Sanchis
Jaime De La Mota Sanchis 2021-10-8
Hello. I think that would be very useful to me, as I feared, it seems like even if it can be done, a lot of work should be put in it.

请先登录，再进行评论。

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

采纳的回答

Steven Lord
Steven Lord 2021-10-8
编辑：Steven Lord 2021-10-8
Ran in:
Are you looking for a mathematical expression for the fitted model or do you just want to be able to evaluate the fitted model at new points? In the former case, how were you hoping to use such an expression? If the latter case, that's easy.
load census
f = fit(cdate, pop, 'linear', 'Normalize', 'on')
f =
Linear interpolant: f(x) = piecewise polynomial computed from p where x is normalized by mean 1890 and std 62.05 Coefficients: p = coefficient structure
plot(cdate, pop, 'r-+', 1977, f(1977), 'ko') % Evaluate f in 1977
Looks like a pretty good match.
  4 个评论
显示 2更早的评论
Jaime De La Mota Sanchis
Jaime De La Mota Sanchis 2021-10-8
I want to send it to another computer where pyomo, an optimizator made in python is used. I have sent the coefficients to said optimizator, but the result is numerically unfeasible.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Get Started with Curve Fitting Toolbox 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by