bounded biharmonic weights deformation，alec jacobson的论文matlab函数biharmonic_bounded注解

V = [ 0.41089 3.2532;
2.8251 2.7639;
0.14673 0.96887;
2.0997 1.4961];
F=delaunay(V)
C=[1,2;1,1.5];
bc=[0,1;1,0]
n = size(V,1);
m = size(bc,2);
L = cotmatrix(V,F)
M = massmatrix(V,F,‘voronoi’)
low = 0;
up = 1;
type = ‘conic’;
b=[3,4]

param = [];
Qi = L*(M\L)

Q = sparse(mn,mn)

for ii = 1:m
d = (ii - 1)*n + 1;
Q(d:(d + n-1), d:(d + n-1)) = Qi;
end
Q
PA = repmat(speye(n,n),1,m)
Pb = ones(n,1)

BCAi = speye(n,n)
BCAi = BCAi(b,:)
BCA = sparse(msize(BCAi,1),msize(BCAi,2))

for ii = 1:m
di = (ii - 1)*size(BCAi,1) + 1;
dj = (ii - 1)*size(BCAi,2) + 1;
BCA(di:(di + size(BCAi,1)-1), dj:(dj + size(BCAi,2)-1)) = BCAi;
end

BCb=bc(😃

ux = up.ones(mn,1)
lx = low.ones(mn,1)

pabca=[PA;BCA]
pabca=[Pb;BCb]
W = quadprog(Q,zeros(n*m,1),[],[],[PA;BCA],[Pb;BCb],lx,ux,[],param)

F =

 2     1     4
 1     3     4

bc =

 0     1
 1     0

L =

(1,1) -0.78572
(2,1) 0.14676
(3,1) 0.27439
(4,1) 0.36457
(1,2) 0.14676
(2,2) -0.86944
(4,2) 0.72268
(1,3) 0.27439
(3,3) -0.68713
(4,3) 0.41274
(1,4) 0.36457
(2,4) 0.72268
(3,4) 0.41274
(4,4) -1.5
M =

(1,1) 1.1267
(2,2) 0.6081
(3,3) 0.78498
(4,4) 1.349
b =

 3     4

Qi =

(1,1) 0.77778
(2,1) -0.11689
(3,1) -0.32
(4,1) -0.3409
(1,2) -0.11689
(2,2) 1.6494
(3,2) 0.25685
(4,2) -1.7893
(1,3) -0.32
(2,3) 0.25685
(3,3) 0.79458
(4,3) -0.73143
(1,4) -0.3409
(2,4) -1.7893
(3,4) -0.73143
(4,4) 2.8617

Q =
全零稀疏矩阵: 8×8

Q =

(1,1) 0.77778
(2,1) -0.11689
(3,1) -0.32
(4,1) -0.3409
(1,2) -0.11689
(2,2) 1.6494
(3,2) 0.25685
(4,2) -1.7893
(1,3) -0.32
(2,3) 0.25685
(3,3) 0.79458
(4,3) -0.73143
(1,4) -0.3409
(2,4) -1.7893
(3,4) -0.73143
(4,4) 2.8617
(5,5) 0.77778
(6,5) -0.11689
(7,5) -0.32
(8,5) -0.3409
(5,6) -0.11689
(6,6) 1.6494
(7,6) 0.25685
(8,6) -1.7893
(5,7) -0.32
(6,7) 0.25685
(7,7) 0.79458
(8,7) -0.73143
(5,8) -0.3409
(6,8) -1.7893
(7,8) -0.73143
(8,8) 2.8617
PA =
(1,1) 1
(2,2) 1
(3,3) 1
(4,4) 1
(1,5) 1
(2,6) 1
(3,7) 1
(4,8) 1
Pb =

 1
 1
 1
 1

BCAi =

(1,1) 1
(2,2) 1
(3,3) 1
(4,4) 1
BCAi =

(1,3) 1
(2,4) 1
BCA =

全零稀疏矩阵: 4×8
BCb =

 0
 1
 1
 0

ux =

 1
 1
 1
 1
 1
 1
 1
 1

lx =

 0
 0
 0
 0
 0
 0
 0
 0

pabca =

(1,1) 1
(2,2) 1
(3,3) 1
(5,3) 1
(4,4) 1
(6,4) 1
(1,5) 1
(2,6) 1
(3,7) 1
(7,7) 1
(4,8) 1
(8,8) 1
pabca =

 1
 1
 1
 1
 0
 1
 1
 0

W =

  0.58858
        1
        0
        1
  0.41142
  7.6593e-11
        1
        0

help quadprog
quadprog - Quadratic programming

Solver for quadratic objective functions with linear constraints.

x = quadprog(H,f)
x = quadprog(H,f,A,b)
x = quadprog(H,f,A,b,Aeq,beq)
x = quadprog(H,f,A,b,Aeq,beq,lb,ub)
x = quadprog(H,f,A,b,Aeq,beq,lb,ub,x0)
x = quadprog(H,f,A,b,Aeq,beq,lb,ub,x0,options)
x = quadprog(problem)
[x,fval] = quadprog(___)
[x,fval,exitflag,output] = quadprog(___)
[x,fval,exitflag,output,lambda] = quadprog(___)