B样条曲线的光顺设计 - 浙江大学本科生毕业论文

毕业论文

N2(i,i)=-6/(t12*t11)-6/(t21*t11); N2(i,i+1)=6 /(t21*t11); end

N2 = N2(2:n+1,:); t = t(3:n+2); DF = sparse(n - 2,n); for i = 1 : n - 2

DF(i, i) = 1 / (t(i+1) - t(i));

DF(i, i+1) = -(1 / (t(i+1) - t(i)) + 1 / (t(i+2) - t(i+1))) DF(i, i+2) = 1 / (t(i+2) - t(i+1)); end end

function [ matC] = getCurvatureMatrix( P,N1,N2 )

n = length(P); Px = P(1:n/2,:); Py = P(n/2+1:n,:); d1x = N1 * Px; d1y = N1 * Py; n=length(d1x); %% length(Px) - length(d1x)==2 M

=sparse(n,n+n); for i = 1 : n

d = [];

d(1) = d1x(i); d(2) = d1y(i);

M(i,i) = -d(2) / ((d(1)^2+d(2)^2)^1.5); M(i,i+n) = d(1) / ((d(1)^2+d(2)^2)^1.5);

end

N = [N2,sparse(n,n+2);sparse(n,n+2),N2]; matC = M * N; end

30

毕业论文

31