在MATLAB环境下插值拟合曲线

使用MATLAB，对图像中牙齿的轮廓曲线进行多种拟合，并且分析出误差大小，再评价出各种方法的优劣。

imshow（‘cast.bmp‘）
[xa]=ginput
n=numel（x）

% 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:n-3; % 由于4个顶点可以确定一条B样条，共有8个控制顶点，可确定5条B样条曲线；
%for i=1:7;
hold on
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；

s=b0.*x（i）+b1.*x（i+1）+b2.*x（i+2）+b3.*x（i+3）; % 确定曲线的横坐标x；
y=b0.*a（i）+b1.*a（i+1）+b2.*a（i+2）+b3.*a（i+3）; % 确定曲线的纵坐标y；

 plot（sy‘r‘）; % 绘制曲线坐标点；
end
end 
hold off;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1039  2009-10-10 00:15  牙齿\11912884mmdraw\B.m

     文件        485  2009-10-09 22:53  牙齿\11912884mmdraw\mmdraw.m

     文件       7159  2009-10-10 20:41  牙齿\abc.png

     文件       1037  2009-10-11 09:08  牙齿\area.m

     文件       1037  2009-10-11 09:08  牙齿\b样条.m

     文件     380406  2009-10-09 19:12  牙齿\cast.bmp

     文件        973  2009-10-10 23:43  牙齿\CubicNaturalSpline.asv

     文件        981  2009-10-10 23:46  牙齿\CubicNaturalSpline.m

     文件        794  2009-10-10 13:01  牙齿\FiveOrder.m

     文件       1235  2009-10-09 18:52  牙齿\PiecewiseFiveOrder.m

     文件        861  2009-10-09 18:52  牙齿\PiecewiseTwoOrder.m

     文件     221707  2009-10-10 15:42  牙齿\pp1.fig

     文件     251341  2009-10-10 23:30  牙齿\pp2.png

     文件       1652  2009-10-09 20:11  牙齿\sanci.m

     文件        293  2009-10-10 12:13  牙齿\try.asv

     文件        293  2009-10-10 12:13  牙齿\try.m

     文件        292  2009-10-09 18:37  牙齿\TwoOrder.m

     文件         31  2009-10-10 20:36  牙齿\untitled4.asv

     文件       1049  2009-10-09 20:59  牙齿\untitled4.fig

     文件        274  2009-10-11 00:25  牙齿\untitled4.m

     目录          0  2009-10-10 00:12  牙齿\11912884mmdraw

     目录          0  2009-10-11 09:08  牙齿

----------- ---------  ---------- -----  ----

               873157                    23