Defining equations with interpolant cfit model (ex. linearinterp, pchipinterp)

2016 2 20
1 个回答
回答已采纳
更新时间：2017 11 12
8 次查看（30 天）

0 个投票

Hi,
does anyone know how to define equations to be used in solve or vpasolve using interpolant cfit model?
I have pchipinterp (pchip_int_fn) and linearinterp (lin_int_fn) cfit models from my raw data and I want to define an equation like the following:
syms x y z
eqn = x*pchip_int_fn(x)*lin_int_fn(x) + y*pchip_int_fn(y)*lin_int_fn(y) + x*pchip_int_fn(z)*lin_int_fn(z) == 0
solve(eqn)
Matlab gives me an error message saying
"Error using cfit/subsref>iParenthesesReference (line 46)
Cannot evaluate CFIT model for some reason.
Error in cfit/subsref (line 16)
out = iParenthesesReference( obj, currsubs );"
Does anyone know how to resolve this issue?
Thank you so much for your help!

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

 采纳的回答

Anish Mitra
Anish Mitra 2016-2-24

1 个投票

cfit objects do not work with inputs as symbolic variables. To work around this, you could re-write the cfit model equation using symbolic variables and then use it in 'eqn', instead of calling the cfit object.
The other option you could explore to find results for x,y,z is using the Optimization Toolbox, and minimizing the square of the value of 'eqn'. This would be a numerical optimization, and symbolic variables would not be used.

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Ho-Jun Suk
Ho-Jun Suk 2016-2-24
Thank you so much Anish! I could not figure out how to find the equations of interpolant cfit models (when I look at the information on the fit model, it returns ans(x) = piecewise polynomial computed from p Coefficients: p = coefficient structure), so I could not use your first option, but the second option works well! Thank you again!
Anish Mitra
Anish Mitra 2016-2-25
Yes, the first option would not be possible for piecewise polynomials. It would work for continuous functions (regular polynomials or exponentials).
Benedikt Trolp
Benedikt Trolp 2017-11-12
Hello Anish, I'm facing a similar Problem to the original poster, but having a hard time executing your second Option. Could you maybe explain it in Detail ? Thanks alot in advance

请先登录，再进行评论。

更多回答（0 个）

类别

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

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by