clear all
load seamount
N = 6;
v=[1 0.5;
0.78023 0.44678;
0.081126 0.30635;
0.6 0.2;
0.48679 0.81763;
0.2586 0.61483]
tri= delaunay(v(:,1),v(:,2))
z=zeros(N,1);
% Delaunay triangulation
figure;
subplot(2,1,1);triplot(tri,v(:,1),v(:,2));hold on;
A=[1,-1/3,0,-1/3,-1/3,0;
-1/4,1,-1/4,-1/4,0,-1/4;
0,-1/2,1,0,0,-1/2;
-1/4,-1/4,0,1,-1/4,-1/4;
-1/3,0,0,-1/3,1,-1/3;
0,-1/4,-1/4,-1/4,-1/4,1;
1,0,0,0,0,0;
0,0,0,1,0,0 ];
b=[0,0,0,0,0,0,1,0.6;
0,0,0,0,0,0,0.5, 0.2]’
newv=inv(A’*A)*A’*b
subplot(2,1,2);triplot(tri,newv(:,1),newv(:,2))
运行结果:
背景资料:变形后