MATLAB绘制B样条曲线[通俗易懂]

2 B样条曲线的绘制

2.1 节点矢量的确定

不同类型的B样条曲线区别主要在于节点矢量，对于具有$n+1$个控制顶点$(P_0,P_1,{...},P_n)$的 $k$ 次B样条曲线，无论是哪种类型都具有$n+k+2$个节点$([u_0,u_1{...}u_{n+k+1}])$。

根据图示，三种类型的B样条曲线对应的节点矢量分别为：

$$[0 \dfrac 17 \dfrac 27\dfrac 37\dfrac 47\dfrac 57\dfrac 67 1]$$

$$[0 \space 0\space 0 \dfrac 13\dfrac 23 1 \space 1\space 1]$$

$$[0 \space 0\space 0 \dfrac 12\dfrac 12 1 \space 1\space 1]$$ 需要注意的是分段Bezier曲线必须满足 $\dfrac nk=\text 正整数$。

这里给出准均匀B样条和分段Bezier曲线的生成节点矢量的代码，均匀B样条的很简单就不列出了。假设共n+1个控制顶点，k次B样条，输入参数为 n, k ，输出节点矢量NodeVector。

准均匀B样条曲线的节点矢量生成：

% U_quasi_uniform.m文件
function NodeVector = U_quasi_uniform(n, k)
% 准均匀B样条的节点向量计算，共n+1个控制顶点，k次B样条
NodeVector = zeros(1, n+k+2);
piecewise = n - k + 1;       % 曲线的段数
if piecewise == 1       % 只有一段曲线时，n = k
    for i = n+2 : n+k+2
        NodeVector(1, i) = 1;
    end
else
    flag = 1;       % 不止一段曲线时
    while flag ~= piecewise
        NodeVector(1, k+1+flag) = NodeVector(1, k + flag) + 1/piecewise;
        flag = flag + 1;
    end
    NodeVector(1, n+2 : n+k+2) = 1;
end

分段Bezier曲线的节点矢量生成：

% U_piecewise_Bezier.m文件
function NodeVector = U_piecewise_Bezier(n, k)
% 分段Bezier曲线的节点向量计算，共n+1个控制顶点，k次B样条
% 分段Bezier端节点重复度为k+1，内间节点重复度为k,且满足n/k为正整数

if ~mod(n, k) && (~mod(k, 1) && k>=1)   % 满足n是k的整数倍且k为正整数
    NodeVector = zeros(1, n+k+2);   % 节点矢量长度为n+k+2
    NodeVector(1, n+2 : n+k+2) = ones(1, k+1);  % 右端节点置1

    piecewise = n / k;      % 设定内节点的值
    Flg = 0;
    if piecewise > 1
        for i = 2 : piecewise
            for j = 1 : k
                NodeVector(1, k+1 + Flg*k+j) = (i-1)/piecewise;
            end
            Flg = Flg + 1;
        end
    end

else
    fprintf('error!\n');
end

2.2 B样条曲线的绘制

根据B样条曲线的定义公式，曲线上任一点坐标值是参数变量$u$的函数，用矩阵形式表示

$$p(u)=\pmatrix{d_0 & d_1 & \ldots & d_n} \pmatrix{N_{0,k}(u)\crN_{1,k}(u)\cr \vdots \crN_{n,k} (u)}$$ 只需要确定控制顶点$d_i$、曲线的次数$k$ 以及基函数$N_{i,k}(u)$，就完全确定了曲线。

B样条曲线的绘制函数：

% DrawSpline.m文件
function DrawSpline(n, k, P, NodeVector)
% B样条的绘图函数
% 已知n+1个控制顶点P(i), k次B样条，P是2*(n+1)矩阵存控制顶点坐标, 节点向量NodeVector
plot(P(1, 1:n+1), P(2, 1:n+1),...
                    'o','LineWidth',1,...
                    'MarkerEdgeColor','k',...
                    'MarkerFaceColor','g',...
                    'MarkerSize',6);
line(P(1, 1:n+1), P(2, 1:n+1));
Nik = zeros(n+1, 1);
for u = 0 : 0.005 : 1-0.005
    for i = 0 : 1 : n
        Nik(i+1, 1) = BaseFunction(i, k , u, NodeVector);
    end
    p_u = P * Nik;
    if u == 0
        tempx = p_u(1,1);
        tempy = p_u(2,1);
        line([tempx p_u(1,1)], [tempy p_u(2,1)],...
            'Marker','.','LineStyle','-', 'Color',[.3 .6 .9], 'LineWidth',3);
    else
        line([tempx p_u(1,1)], [tempy p_u(2,1)],...
            'Marker','.','LineStyle','-', 'Color',[.3 .6 .9], 'LineWidth',3);
        tempx = p_u(1,1);
        tempy = p_u(2,1);
    end
end

调用 DrawSpline(n, k, P, NodeVector) 函数就能绘制曲线，注意输入变量要正确。

下面给出绘制三种不同B样条曲线的命令流，可以参考比较每种类型之间的区别。

% 绘制三种类型的B样条曲线，需要前面所给的所有.m文件
clear all;
%控制顶点
P = [9.036145, 21.084337, 37.607573, 51.893287, 61.187608;
    51.779661, 70.084746, 50.254237, 69.745763, 49.576271];

n = 4; k = 2;

flag = 2;
% flag = 1，绘制均匀B样条曲线
% flag = 2, 绘制准均匀B样条曲线
% flag = 3, 绘制分段Bezier曲线

switch flag
    case 1
        NodeVector = linspace(0, 1, n+k+2); % 均匀B样条的节点矢量

        % 绘制样条曲线
        plot(P(1, 1:n+1), P(2, 1:n+1),...
                        'o','LineWidth',1,...
                        'MarkerEdgeColor','k',...
                        'MarkerFaceColor','g',...
                        'MarkerSize',6);
        line(P(1, 1:n+1), P(2, 1:n+1));
        Nik = zeros(n+1, 1);
        for u = k/(n+k+1) : 0.001 : (n+1)/(n+k+1)
            % for u = 0 : 0.005 : 1
            for i = 0 : 1 : n
                Nik(i+1, 1) = BaseFunction(i, k , u, NodeVector);
            end
        p_u = P * Nik;
        line(p_u(1,1), p_u(2,1), 'Marker','.','LineStyle','-', 'Color',[.3 .6 .9]);
        end
    case 2
        NodeVector = U_quasi_uniform(n, k); % 准均匀B样条的节点矢量
        DrawSpline(n, k, P, NodeVector);
    case 3
        NodeVector = U_piecewise_Bezier(n, k);  % 分段Bezier曲线的节点矢量
        DrawSpline(n, k, P, NodeVector);
    otherwise
        fprintf('error!\n');
end

三种类型的B样条曲线：
1. 均匀B样条曲线

2. 准均匀B样条曲线

3. 分段Bezier曲线

参考文献：

[1] 施法中. 计算机辅助几何设计与非均匀有理B样条(修订版)[M]. 北京: 高等教育出版社, 2013 : 217-248.