Can I get the mathematical expression of an interpolation?

Question

Jaime De La Mota Sanchis 2021-10-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1469771-can-i-get-the-mathematical-expression-of-an-interpolation

评论： Steven Lord 2021-10-9

采纳的回答： Steven Lord

在 MATLAB Online 中打开

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 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 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.

请先登录，再进行评论。

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

Answer 1

Steven Lord 2021-10-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1469771-can-i-get-the-mathematical-expression-of-an-interpolation#answer_804481

编辑：Steven Lord 2021-10-8

在 MATLAB Online 中打开

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更早的评论隐藏 2更早的评论

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.