資源簡介
三次均勻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???????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
評論
共有 條評論