You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
20 lines
646 B
20 lines
646 B
1 year ago
|
function val = dC_eval_dt(C,t0)
|
||
|
|
% Given a bernstein polynomial, evaluates it i at t0;
|
||
|
|
% Input parameters
|
||
|
|
% C: This is a bernstein polynomial's control points (aka coefficients)
|
||
|
|
% t0: This is the point at which we wish to split the polynomial
|
||
|
|
% Output parameters
|
||
|
|
% val: row vector of values of the bernstein polynomial defined by C at t0
|
||
|
|
[d,n] = size(C);
|
||
|
|
val=zeros(length(t0),d);
|
||
|
|
for j=1:d
|
||
|
|
% tempdivmat=repmat(C(j,:),1,length(t0));
|
||
|
|
cdiv_mat=C(j,:);
|
||
|
|
for i=2:(n-1)
|
||
|
|
cdiv_mat=cdiv_mat.*(1-t0') + circshift(cdiv_mat,-1,2).*t0';
|
||
|
|
end
|
||
|
|
% cdiv_mat=triu(fliplr(cdiv_mat));
|
||
|
|
|
||
|
|
val(:,j) = (n-1)*(cdiv_mat(:,2)-cdiv_mat(:,1));
|
||
|
|
end
|