matlab練習程序（B樣條反算控制點）

如果要准確反求B樣條的控制點，有幾個參數還是要事先知道的：

1. 樣條的控制點個數。

2. B樣條曲線的所有點坐標和個數。

3. B樣條基函數。

一般條件2容易知道一些，1和3還是比較難事先知道的。

如果待求控制點為四個，B樣條曲線點個數為n個，並且已知基函數形式如下面代碼中[b0 b1 b2]。

則可以列出類似下面這樣一個方程組求解，其中px,py為帶求解控制點，x、y為B樣條曲線點：

matlab代碼如下：

clear all;close all;clc;

p =ginput();  
plot(p(:,1),p(:,2),'g-o');  

p
%二次均勻b樣條
re=[];
for i=1:length(p)-2  
    for t=0:0.1:1   
        
        b0 = 1/2*(1-t)^2;
        b1 = 1/2*(-2*t^2+2*t+1);
        b2 = 1/2*t^2;
        
        x=b0*p(i,1)+b1*p(i+1,1)+b2*p(i+2,1);
        y=b0*p(i,2)+b1*p(i+1,2)+b2*p(i+2,2);        
        
        re=[re;x y];
    end
end

hold on;  
plot(re(:,1),re(:,2),'r-o');

%構造方程求解
n = length(p);
step = length(re) / (n-2);
A = zeros(length(re)*2,n*2);
B = zeros(length(re)*2,1);
t = linspace(0,1,step); 
for i=1:n-2
    for j=1:length(t)
        b0 = 1/2*(1-t(j))^2;
        b1 = 1/2*(-2*t(j)^2+2*t(j)+1);
        b2 = 1/2*t(j)^2;  
        
        A((j-1)*2+1+(i-1)*length(t)*2,1+(i-1)*2) = b0;
        A((j-1)*2+1+(i-1)*length(t)*2,3+(i-1)*2) = b1;
        A((j-1)*2+1+(i-1)*length(t)*2,5+(i-1)*2) = b2;
        A((j-1)*2+2+(i-1)*length(t)*2,2+(i-1)*2) = b0;        
        A((j-1)*2+2+(i-1)*length(t)*2,4+(i-1)*2) = b1;
        A((j-1)*2+2+(i-1)*length(t)*2,6+(i-1)*2) = b2;        
        
        B((j-1)*2+1+(i-1)*length(t)*2) = re(j+(i-1)*length(t),1);
        B((j-1)*2+2+(i-1)*length(t)*2) = re(j+(i-1)*length(t),2); 
    end
end

X = inv(A'*A)*A'*B;

%畫出新曲線看看
figure;
p = reshape(X,[2 length(X)/2])'
plot(p(:,1),p(:,2),'g-*');  
re=[];
for i=1:length(p)-2  
    for t=0:0.1:1   
        
        b0 = 1/2*(1-t)^2;
        b1 = 1/2*(-2*t^2+2*t+1);
        b2 = 1/2*t^2;
        
        x=b0*p(i,1)+b1*p(i+1,1)+b2*p(i+2,1);
        y=b0*p(i,2)+b1*p(i+1,2)+b2*p(i+2,2);        
        
        re=[re;x y];
    end
end

hold on;  
plot(re(:,1),re(:,2),'r-*');

結果如下：

原始B樣條曲線：

反求解控制點后生成的B樣條曲線：