quadraticBetaSplines
Copyright © 1984 by Xerox Corporation. All rights reserved.
This file contains the derivation of the quadratic Beta-spline.
Last Edited by: Tso, July 18, 1984 10:37:34 am PDT
The Constraints:
0 = b1(0) 0 = b1'(0)
b1(1) = b0(0) bb1'(1) = b0'(0)
b0(1) = b-1(0) bb0'(1) = b-1'(0)
b-1(1) = 0 bb-1'(1) = 0
r = 1
E B[i+r](ui) = b1(0) + b0(0) + b-1(0)
r = -1
The Equations:
b1 = c0,1 + c1,1u + c2,1usq
b0 = c0,0 + c1,0u + c2,0usq
b-1 = c0,-1 + c1,-1u + c2,-1usq
b1' = c1,1 + 2c2,1u
b0' = c1,0 + 2c2,0u
b-1' = c1,-1 + 2c2,-1u
b1(0) = c0,1 b1(1) = c0,1 + c1,1 + c2,1
b0(0) = c0,0 b0(1) = c0,0 + c1,0 + c2,0
b-1(0) = c0,-1 b-1(1) = c0,-1 + c1,-1 + c2,-1
b1'(0) = c1,1 b1'(1) = c1,1 + 2c2,1
b0'(0) = c1,0 b0'(1) = c1,0 + 2c2,0
b-1'(0) = c1,-1 b-1'(1) = c1,-1 + 2c2,-1
b1(0) + b0(0) + b-1(0) = 1
The Derivation:
(1)
c0,1 = 0
c0,0 = c1,1 + c2,1
c0,-1 = c0,0 + c1,0 + c2,0
c1,1 = 0
c1,0 = 2bc2,1 9 equations in 9 unknowns
c1,-1 = bc1,0 + 2bc2,0
0 = c0,-1 + c1,-1 + c2,-1
0 = c1,-1 + 2c2,-1
1 = c0,1 + c0,0 + c0,-1
(2)
c0,0 = c2,1
c0,-1 = c0,0 + c1,0 + c2,0
c1,0 = 2bc2,1 7 equations in 7 unknowns, remove c0,1 and c1,1
c1,-1 = bc1,0 + 2bc2,0
0 = c0,-1 + c1,-1 + c2,-1
0 = c1,-1 + 2c2,-1
1 = c0,0 + c0,-1
(3)
c0,-1 = c0,0 + c1,0 + c2,0
c1,0 = 2bc0,0 6 equations in 6 unknowns, remove c2,1
c1,-1 = bc1,0 + 2bc2,0
0 = c0,-1 + c1,-1 + c2,-1
c1,-1 = -2c2,-1
1 = c0,0 + c0,-1
(4)
c0,-1 = c0,0 + c1,0 + c2,0
c1,0 = 2bc0,0 5 equations in 5 unknowns, remove c2,-1
c1,-1 = bc1,0 + 2bc2,0
0 = c0,-1 + (1/2)c1,-1
1 = c0,0 + c0,-1
(5)
c0,-1 = c0,0 + c1,0 + c2,0
c1,0 = 2bc0,0 4 equations in 4 unknowns, remove c1,-1
-2c0,-1 = bc1,0 + 2bc2,0
1 = c0,0 + c0,-1
(6)
1 = 2c0,0 + c1,0 + c2,0
c1,0 = 2bc0,0 3 equations in 3 unknowns, remove c0,-1
-2 = -2c0,0 + bc1,0 + 2bc2,0
(7)
1 = (2 + 2b)c0,0 + c2,0 2 equations in 2 unknowns, remove c1,0
-2 = (-2 + 2bsq)c0,0 + 2bc2,0
(8)
-2b = (-4b - 4bsq)c0,0 - 2bc2,0
-2 = (-2 + 2bsq)c0,0 + 2bc2,0
-----------------------------
-2(b + 1) = -2(1 + 2b + bsq)c0,0
(9)
c0,0 = (b + 1)/(1 + 2b + bsq) = 1/(b + 1), let d = 1 + b
(10)
c2,0 = -1
c1,0 = 2b/d
c0,-1 = b/d
c1,-1 = -2b/d
c2,-1 = b/d
c2,1 = 1/d
c0,1 = 0
c1,1 = 0
(11) +
b1 = 1/dusq
b0 = 1/d[1 + 2bu - usq]
b-1 = b/d[1 - 2u + usq]