Coefficents of piecewise polynomial in matlab

Question

Mukesh Singh Bisht 2023-11-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2050622-coefficents-of-piecewise-polynomial-in-matlab

编辑： Bruno Luong 2023-11-24

在 MATLAB Online 中打开

Hi I have a simple example of piecewise polynomial from MATLAB using the function 'cscvn'

cscvn( [1 0 -1 0 1;0 1 0 -1 0] )

ans : form: 'pp'

breaks: [0 1.1892 2.3784 3.5676 4.7568]

coefs: [8×4 double]

pieces: 4

order: 4

dim: 2

My question is : How do I read the coeffiecnt matrix as shown above. I can understand from the above coefs matrix that it has 4 columns corresponding to 4 coefficents of a 4th order polynomial. But why this matrix has 8 columns? I only have 4 pieces of polynomial, then why double the number of rows in coefs matrix?

Seocnd question is : How can I use the above coefs matrix to get the equation of piecewise poilynomial

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

请先登录，再进行评论。

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

Answer 1

John D'Errico 2023-11-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2050622-coefficents-of-piecewise-polynomial-in-matlab#answer_1357902

编辑：John D'Errico 2023-11-22

在 MATLAB Online 中打开

Do you understand the result will not be some simple polynomial like you are used to seeing?

plot([1 0 -1 0 1],[0 1 0 -1 0],'-o')

You want a curve that follows around this diamond shaped thing. So you CANNOT have a function of the form y(x). PERIOD. You cannot do so.

Instead, cscvn generates a curve as a function of a parameter, I'll call it t, where t is the arclength of the piecewise linear curve that connects those dots. So cscvn creates the pair of functions x(t), and y(t).

READ THE DOCS FOR CSCVN: cscvn

In there this scheme is explained fairly clearly.

Again, there is no simple equation of the form y(x). There cannot be, as if you look at the plot, the result cannot be a single valued function.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Torsten 2023-11-23

编辑：Torsten 2023-11-23

在 MATLAB Online 中打开

Does this help ?

c=cscvn( [1 0 -1 0 1;0 1 0 -1 0] );

c.coefs

ans = 8×4

0.2973 -1.0607 0.0000 1.0000 -0.2973 0.0000 1.2613 0 0.2973 0.0000 -1.2613 0 0.2973 -1.0607 -0.0000 1.0000 -0.2973 1.0607 -0.0000 -1.0000 0.2973 0 -1.2613 0 -0.2973 0.0000 1.2613 0 -0.2973 1.0607 0.0000 -1.0000

c.breaks

ans = 1×5

0 1.1892 2.3784 3.5676 4.7568

x = @(t)1-1.0607*t.^2+0.2973*t.^3;

y = @(t)1.2613*t-0.2973*t.^3;

t = linspace(0,1.1892,20);

hold on

plot(x(t),y(t))

x = @(t)-1.2613*(t-1.1892)+0.2973*(t-1.1892).^3;

y = @(t)1-1.0607*(t-1.1892).^2+0.2973*(t-1.1892).^3;

t = linspace(1.1892,2.3784,20);

plot(x(t),y(t))

x = @(t)-1+1.0607*(t-2.3784).^2-0.2973*(t-2.3784).^3;

y = @(t)-1.2613*(t-2.3784)+0.2973*(t-2.3784).^3;

t = linspace(2.3784,3.5676,20);

plot(x(t),y(t))

x = @(t)1.2613*(t-3.5676)-0.2973*(t-3.5676).^3;

y = @(t)-1+1.0607*(t-3.5676).^2-0.2973*(t-3.5676).^3;

t = linspace(3.5676,4.7568,20);

plot(x(t),y(t))

hold off

Walter Roberson 2023-11-23

(The rows appear to alternate, one row of x, then one row of y, with the pairs repeated according to the number of pieces. So x coordinates would be (1:2:end,:) and y coordinates would be (2:2:end,:)

Bruno Luong 2023-11-24

@John D'Errico "Instead, cscvn generates a curve as a function of a parameter, I'll call it t, where t is the arclength of the piecewise linear curve that connects those dots"

Actually it's the sum of the square-root of the linear arclength, this thing is called Eugene Lee centripetal parametrization according to cscvn doc

请先登录，再进行评论。