非均匀三次B样条插值和插入节点算法.doc

非均匀三次B样条插值和插入节点算法  HYPERLINK /attachment.php?aid=1105k=6077916f21cc09646a3ff6ecb0f96832t=1239094195nothumb=yes \t _blank m.jpg (31.17 KB) 2009-2-18 15:33 复制内容到剪贴板 代码: %X:原始资料,d:控制顶点%n:数据条数,k:B样条的次数%%see also %X=load(data.txt);[n,numy]=size(X); %得数据维数;k=3;%弦长参数化u(k+1)=0;for i=1:(n-1)? ? u(k+i+1)=u(k+i)+abs(X(i+1,1)-X(i,1));end%====规范参数化=========================temp=u(n+k);for i=(k):(n+k) %(k+1):(n-1+k)足够了? ? u(i)=u(i)/temp;end%首末节点重复k+1个for i=1:(k+1)? ? u(i)=0;? ? u(n-1+k+i)=1;end%% ===============================%A:方程系数------采用递推算法-----------------------A=zeros(n+