bsplinecurve.m

任意给定点坐标

% B样条曲线生成程序
% 说明：给定8个控制顶点{(3 5),(2 4),(3 2),(6 1),(5 8),(10 6),(8 1),(6 0)}
load c;x=800:25:1800;
a(1,:)=x;a(2,:)=c;
% a=[10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90;20 16 26 26 26 26 24 24 24 22 22 16 28 26 26 26 24];    % 控制顶点数组；
%a=[4 3 2 3 6 5 10 8 6 4;6 5 4 2 1 8 6 1 0 -1];%增加两个控制顶点（4 6）（4 －1）从而使曲线通过原来的首末点；
plot(a(1,:),a(2,:),':');                 % 绘制控制多边形；
hold on;                   % 默认为hold off，此命令用来保留控制多边形的图形；                               
for i=1:38;          % 由于4个顶点可以确定一条B样条，共有8个控制顶点，可确定5条B样条曲线；
%for i=1:7;
for u=0:0.001:1;    % 参数u；   
    
    b0=1.0./6.*(1-u).^3;                     % 基函数b0；
    b1=1.0./6.*(3.*u.^3-6.*u.^2+4);          % 基函数b1；
    b2=1.0./6.*(-3.*u.^3+3.*u.^2+3.*u+1);    % 基函数b2；
    b3=1.0./6.*u.^3;                         % 基函数b3；
    
    x=b0.*a(1,i)+b1.*a(1,i+1)+b2.*a(1,i+2)+b3.*a(1,i+3);   % 确定曲线的横坐标x；
    y=b0.*a(2,i)+b1.*a(2,i+1)+b2.*a(2,i+2)+b3.*a(2,i+3);   % 确定曲线的纵坐标y；
    
    line(x,y);       % 绘制曲线坐标点；
end
end
hold off;