polynomial function 3D

Question

Christoph Thorwartl 2021-1-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/730203-polynomial-function-3d

评论： Christoph Thorwartl 2021-1-29

采纳的回答： Bjorn Gustavsson

Hello,

I would like to put a second order polynomial function through 9 points. The points are given in the form P(x,y,z).

Goal:

(1) Calculate polynomial function

(2) Calculate the mean curvature (= second derivative) of the function

Is there something similar to "p = polyfit(x,y,z,n)"?

Remark:

I can put a circle through three points in space and calculate the corresponding curvature of the circle. For this I use "circlefit3d(p1,p2,p3)". But I don't want to use only 3 of the 9 points for the curvature calculation.

I would be very pleased about a feedback.

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

Christoph Thorwartl 2021-1-29

3D markers were fixed on the surface of an object (in a row). A force is applied to the object which causes it to bend. The bending line does not have a completely homogeneous curvature.

Christoph Thorwartl 2021-1-29

Bending.pdf

Picture see attached.

请先登录，再进行评论。

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

Answer 1

Bjorn Gustavsson 2021-1-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/730203-polynomial-function-3d#answer_609263

You have the file exchange contribution polyfitn that should do what you want.

HTH

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Walter Roberson 2021-1-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/730203-polynomial-function-3d#answer_609268

在 MATLAB Online 中打开

Is there something similar to "p = polyfit(x,y,z,n)"?

No. However, polyfit() constructs the Vandermode companion matrix and then uses \ for linear least squared computation of the coefficients. You can do the same thing;

X = x(:); Y = Y(:); Z = Z(:);
A = [X.^2, X.*Y, Y.^2, X, Y, ones(size(X))];
b = Z;
A\b

Assuming a*x^2 + b*x*y + c*y^2 + d*x + e*y + f = z .

In particular, assuming that there is no z^2 or x*z or y*z .

If there are higher order z, then at the moment I am not sure how to set up the problem without ending up with a companion matrix and all zeros on the right hand side, which unfortunately leads you to an all-zero solution instead of something useful.