三次均匀B样条插值计算 matlab

三次均匀B样条插值计算 三次均匀B样条插值计算 三次均匀B样条插值计算 三次均匀B样条插值计算

function bs = uniformbspline（indexxsequence）
% 三次均匀B样条插值计算（CBI Cubic B-spline Interpolation）
% bs = uniformbspline（indexxsequence）
% 输入参数:  index     为代求点的横坐标向量（Y的坐标）
%               x     已知的点的坐标值
%        sequence     所有点的横坐标向量 sequence=1:n;
% 输出参数:
%              bs    三次均匀B样条插值结果
% eMail  clearsblue@163.com


%
% beta = 1;
% t=0.1;  % 可以更改
% p（t） = 1/6*[1 t t^2 t^3]*[0       6          0          0
%                           -3*beta 0          3*beta     0 
%                           6*beta  -18+3*beta 18-6*beta  -3*beta
%                           -3*beta 12-3*beta  -12+3*beta 3*beta  ]*[pL;pL1;pL2;pL3]
%
% 数据产生

if（nargin == 1）
     error（‘必须输入两个参数: bs = uniformbspline（indexx）‘）;
end

if length（index） ~=length（x）
   error（‘indexx向量长度不相等!‘）;
end;

% 均匀B样条插值计算
x=[x（1）x];
x=[x（1）x];
x=[x（1）x];
x=[xx（end）];
x=[xx（end）];
x=[xx（end）];

index = [index（1）index];
index = [index（1）index];
index = [index（1）index];
index = [indexindex（end）];
index = [indexindex（end）];
index = [indexindex（end）];

len = length（index）;
beta = 0.000001;
ind = 1;
for i = 1:len-3
    pL = x（i）;
    pL1 = x（i+1）;
    pL2 = x（i+2）;
    pL3 = x（i+3）;
    Step = abs（（index（i+2）-index（i+1）））;
    if Step == 0 
       Step = 1;
    end;
    k = 0:1/Step:1;
    for j = 1:length（k）
        t = k（j）;
        p（ind） = 1/6*[1 t t^2 t^3]*[0       6          0          0
                                   -3*beta  0          3*beta     0 
                                    6*beta  -18+3*beta 18-6*beta  -3*beta
                                   -3*beta   12-3*beta -12+3*beta  3*beta]*[pL;pL1;pL2;pL3];
        ind = ind + 1;
    end
    ind = ind  - 1;
end;
p = p（3:end-2）;
indexmin = min（index）;
if（nargin == 3）
   if indexmin < sequence（1）
      p = p（abs（indexmin）+2:end）;    
   end;
   bs = p（sequence）; 
else
   bs = p;
end;

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

     文件       9721  2008-12-16 19:13  使用帮助：新手必看.htm

     文件       2003  2009-01-18 13:11  uniformbspline.m

     文件        183  2008-09-03 11:09  Matlab中文论坛--助努力的人完成毕业设计.url

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

                11907                    3