angular data spline interpolation

1 次查看（过去 30 天）

mukesh bisht 2021-8-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1438384-angular-data-spline-interpolation

评论： darova 2021-8-23

Hi, I have the following query using spline interpolation for angular data. The code attached below is from matlab documentation:

x = pi*[0:.5:2];

y = [0 1 0 -1 0 1 0;

1 0 1 0 -1 0 1];

pp = spline(x,y);

yy = ppval(pp, linspace(0,2*pi,101));

plot(yy(1,:),yy(2,:),'-b',y(1,2:5),y(2,2:5),'or')

axis equal

This code gives the closed piecewise polynomial for circle of radius 1unit. Now if I calculate the value of plynomial for a given query point, it doesnt give me the correct answer. why??

For ex: xq = 0; yq = ppval(pp,xq);

It gives yq = [1;0] but it should be yq = [1,-1];

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

Julius Muschaweck 2021-8-23

If you think my answer is correct, please make it the "accepted answer".Thanks

请先登录，再进行评论。

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

采纳的回答

Julius Muschaweck 2021-8-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1438384-angular-data-spline-interpolation#answer_772204

Your x has five values, your y has seven values. This is a mismatch.

I think you should use

y = [0 1 0 -1 0;

1 0 -1 0 1];

and when you do that,

xq = 0; yq = ppval(pp,xq);

gives you [0;1] as it should. Your "circle" looks more triangular, but that's ok since the sine is badly approximated by the spline with only 5 support points.