loading毕业论文
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