clear all;

close all;

clc;

p=ginput();

plot(p(:,),p(:,),'b-o');

N=length(p);    %确定贝塞尔阶数（控制点个数-1）

t=zeros(N,N);

M=;  %确定贝塞尔曲线点的个数

%计算杨辉三角

for i=:N

    t(i,) = ;

    t(i,i) = ;

end

if N>=

    for i=:N

        for j=:i-

            t(i,j) = t(i-,j-)+t(i-,j);

        end

    end

end

%根据公式计算贝塞尔曲线

re=zeros(M,);

for i=:M

    step = i/M;

    for k=:N-

        re(i,) = re(i,) + (-step)^(N-k-)*p(k+,)*step^k*t(N,k+); %t替换为nchoosek(N-,k)，不用计算杨辉三角了;

        re(i,) = re(i,) + (-step)^(N-k-)*p(k+,)*step^k*t(N,k+); %t替换为nchoosek(N-,k)，不用计算杨辉三角了;

    end

end

hold on;

plot(re(:,),re(:,),'r');