Demo This demo contains the calculation procedures to generate the results shown in Appendix A of the lab report. Degree Reduction and Raising Matrices for Degree M,4 = 3 4 4 4 3 4 3 4 4 4 3 4 4 M 4, = 3 3 4 8 3 9 M,3 = 3 3 3 7 37 4 3 7 3 7 37 4 3 7 3 3 9 8 3 3 3 3 M 3, = 3 3 3 3 3 3 8 7 8 9 3 8 7 8 9 4 7 M, = 9 7 4 3 9 8 7 8 3 8 3 7 3 8 3 M, = 3 3 3 3 3 8 M, = 4 4 8 M, = 4 3 3 4 M, = M, = ( )
QuadClip Applied to a Polynomial of th Degree with Two Roots X 3X 4 X 3 + 4X X + Called QuadClip with input polynomial on interval [, ]: p = X 3X 4 X 3 + 4X X +....3.4...7.8.9.. Recursion Branch for Input Interval [, ] Normalized monomial und Bézier representations and the Bézier polygon: p = X 3X 4 X 3 + 4X X + = B, (X) B, (X) B, (X) +.B 3, (X) + B 4, (X) + B, (X)...3.4...7.8.9 Best approximation polynomials of degree,,, 3 and 4: q =. =.B, q =.374X.487 =.487B, +.987B, q =.48X.7874X.748 =.748B,.448B, +.87B, q 3 =.X 3 +.47X.9X +.7349 =.7349B,3.7B,3 +.438B,3 +.7937B 3,3 q 4 = 7.X 4 7.X 3 +.8333X 7.97X +.99 =.99B,4 3.39484B,4 +.B,4 +.39484B 3,4 +.9794B 4,4
....3.4...7.8.9. q q q q 3 q 4 p Degree reduction and raising matrices:.849.487.743.349.874.743.49.487.487 M, =.487.487.743.874.349.743.487.849..3..74 M, = 8.43.4... 4.99.987 4..3. Degree reduction and raising: q =.48X.7874X.748 =.748B,.448B, +.87B, q =.877 X +.388 X 4.7937 X 3 +.48X.7874X.748 =.748B,.87B,.749B, +.B 3, +.874B 4, +.87B,...3.4...7.8.9 p q q The maximum difference of the Bézier coefficients is δ =.4. Bounding polynomials M and m: Root of M and m: Intersection intervals: M =.48X.7874X +.387 m =.48X.7874X.4743 N(M) = {} N(m) = {.93,.78} 3
...3.4...7.8.9 [, ] Longest intersection interval: = Bisection: first half [,.] und second half [., ]. Recursion Branch on the First Half [,.] Normalized monomial und Bézier representations and the Bézier polygon: p =.78X.87X 4.87X 3 + X 7.X + = B, (X).B, (X) B, (X).87B 3, (X).87B 4, (X) +.87B, (X).....3.4...7.8.9 Best approximation polynomials of degree,,, 3 and 4: q =.978 =.978B, q =.37944X.979 =.979B,.3844B, q = 4.839X.X +.84 =.84B,.837B, +.499B, q 3 = 4.798X 3 +.94X 7.997X +.4 =.4B,3.488B,3.4494B,3 +.8998B 3,3 q 4 =.3437X 4 3.X 3 +.X 7.93X +.3 =.3B,4.89B,4.83B,4.89B 3,4 +.B 4,4 4
....3.4...7.8.9. q q q q 3 q 4 p Degree reduction and raising matrices:.849.487.743.349.874.743.49.487.487 M, =.487.487.743.874.349.743.487.849..3..74 M, = 8.43.4... 4.99.987 4..3. Degree reduction and raising: q = 4.839X.X +.84 =.84B,.837B, +.499B, q =.7 X 3.37 X 4 +.77 X 3 + 4.839X.X +.84 =.84B,.399B,.7997B,.87B 3,.4734B 4, +.499B,...3.4...7.8.9 p q q The maximum difference of the Bézier coefficients is δ =.834. Bounding polynomials M and m: Root of M and m: Intersection intervals: M = 4.839X.X +.89 m = 4.839X.X +.3 N(M) = {.883,.79787} N(m) = {.84,.978}
.....3.4...7.8.9 [.84,.883], [.79787,.978] Longest intersection interval:.93 = Selective recursion: interval : [.9,.448], interval : [.398843,.48339],.3 Recursion Branch in Interval : [.9,.448] Normalized monomial und Bézier representations and the Bézier polygon: p =.78X.49X 4.338X 3 +.337X.933X +.83 =.83B, (X) +.799B, (X).88B, (X).89B 3, (X).89B 4, (X).373B, (X)....3.4...7.8.9. Best approximation polynomials of degree,,, 3 and 4: q =.849 =.849B, q =.387X +.43 =.43B,.4733B, q =.4X.894437X +.874 =.874B,.4473B,.374B, q 3 =.87X 3 +.49X.9397X +.8338 =.8338B,3.7787B,3.3783B,3.3733B 3,3 q 4 =.97X 4.3477X 3 +.34X.93X +.83 =.83B,4 +.74B,4.438B,4.8B 3,4.373B 4,4
....3.4...7.8.9..4 q q q q 3 q 4 p Degree reduction and raising matrices:.849.487.743.349.874.743.49.487.487 M, =.487.487.743.874.349.743.487.849..3..74 M, = 8.43.4... 4.99.987 4..3. Degree reduction and raising: q =.4X.894437X +.874 =.874B,.4473B,.374B, q =.499 3 X.88 X 4 + 8.4994 3 X 3 +.4X.894437X +.874 =.874B, +.388B,.947B,.87B 3,.884B 4,.374B,....3.4...7.8.9. p q q The maximum difference of the Bézier coefficients is δ =.8. Bounding polynomials M and m: Root of M and m: Intersection intervals: M =.4X.894437X +.8393 m =.4X.894437X +.887 N(M) = {.339, 3.38} N(m) = {.347349, 3.39} 7
....3.4...7.8.9. [.347349,.339] Longest intersection interval:.398 = Selective recursion: interval : [.88,.88],.4 Recursion Branch in Interval : [.88,.88] Found root in interval [.88,.88] at recursion depth 4!. Recursion Branch in Interval : [.398843,.48339] Normalized monomial und Bézier representations and the Bézier polygon: p =.78X +.7793X 4.878X 3 +.3977X +.347X.47 =.47B, (X).3B, (X).83B, (X) +.744B 3, (X) +.7978B 4, (X) +.3B, (X)....3.4...7.8.9.. Best approximation polynomials of degree,,, 3 and 4: q =.877 =.877B, q =.373X.44 =.44B, +.9B, q =.8X +.37X.8 =.8B,.39B, +.4B, q 3 =.94X 3 +.39X +.348X.48 =.48B,3.93B,3 +.387B,3 +.98B 3,3 q 4 =.743X 4.8994X 3 +.387X +.349X.4 =.4B,4.97B,4.9B,4 +.478B 3,4 +.B 4,4 8
....3.4...7.8.9.. q q q q 3 q 4 p Degree reduction and raising matrices:.849.487.743.349.874.743.49.487.487 M, =.487.487.743.874.349.743.487.849..3..74 M, = 8.43.4... 4.99.987 4..3. Degree reduction and raising: q =.8X +.37X.8 =.8B,.39B, +.4B, q =.734 4 X.8344 3 X 4 +.998 3 X 3 +.8X +.37X.8 =.8B,.33B,.47B, +.74878B 3, +.7889B 4, +.4B,....3.4...7.8.9.. p q q The maximum difference of the Bézier coefficients is δ =.773. Bounding polynomials M and m: Root of M and m: Intersection intervals: M =.8X +.37X.43 m =.8X +.37X.7 N(M) = { 3.79,.779} N(m) = { 3.737,.838} 9
....3.4...7.8.9.. [.779,.838] Longest intersection interval:.89 = Selective recursion: interval : [.44788,.4489],. Recursion Branch in Interval : [.44788,.4489] Found root in interval [.44788,.4489] at recursion depth 4!.7 Recursion Branch on the Second Half [., ] Normalized monomial und Bézier representations and the Bézier polygon: p =.78X +.787X 4.8X 3.937X +.3X +.87 =.87B, (X) +.B, (X) +.937B, (X) +.87B 3, (X) +.B 4, (X) + B, (X).8..4....3.4...7.8.9 Best approximation polynomials of degree,,, 3 and 4: q =.978 =.978B, q =.493X +.44494 =.44494B, +.9447B, q =.8473X +.37X +.39 =.39B, +.948B, +.8488B, q 3 =.794X 3.744X +.9873X +.939 =.939B,3 +.8B,3 +.498B,3 +.944444B 3,3 q 4 = 3.787X 4 4.48X 3.848X +.9384X +.8 =.8B,4 +.74B,4 +.433B,4 +.39B 3,4 +.999B 4,4
.8..4....3.4...7.8.9 q q q q 3 q 4 p Degree reduction and raising matrices:.849.487.743.349.874.743.49.487.487 M, =.487.487.743.874.349.743.487.849..3..74 M, = 8.43.4... 4.99.987 4..3. Degree reduction and raising: q =.8473X +.37X +.39 =.39B, +.948B, +.8488B, q = 3.79 X + 7.479 X 4.89 X 3.8473X +.37X +.39 =.39B, +.7B, +.778B, +.8893B 3, +.887B 4, +.8488B,.8..4....3.4...7.8.9 p q q The maximum difference of the Bézier coefficients is δ =.387. Bounding polynomials M and m: Root of M and m: Intersection intervals: M =.8473X +.37X +.7779 m =.8473X +.37X.937 N(M) = {.4774,.999} N(m) = {.49,.}
....3.4...7.8.9 [,.49] Longest intersection interval:.49 = Selective recursion: interval : [.,.333],.8 Recursion Branch in Interval : [.,.333] Normalized monomial und Bézier representations and the Bézier polygon: p =.794 7 X + 8.3 X 4.888X 3.49X +.947974X +.87 =.87B, (X) +.3779B, (X) +.4B, (X) +.74987B 3, (X) +.93B 4, (X) +.38B, (X).3.....3.4...7.8.9 Best approximation polynomials of degree,,, 3 and 4: q =.398 =.398B, q =.948X +.94 =.94B, +.3B, q =.44X +.9493X +.873 =.873B, +.7B, +.34B, q 3 =.998X 3.8X +.947998X +.87 =.87B,3 +.3B,3 +.8B,3 +.38B 3,3 q 4 = 8.89 X 4.89X 3.477X +.947974X +.87 =.87B,4 +.4449B,4 +.88B,4 +.887B 3,4 +.38B 4,4
.3.....3.4...7.8.9 q q q q 3 q 4 p Degree reduction and raising matrices:.849.487.743.349.874.743.49.487.487 M, =.487.487.743.874.349.743.487.849..3..74 M, = 8.43.4... 4.99.987 4..3. Degree reduction and raising: q =.44X +.9493X +.873 =.873B, +.7B, +.34B, q =.733 X +.7473 X 4.4 X 3.44X +.9493X +.873 =.873B, +.3779B, +.47B, +.7497B 3, +.933B 4, +.34B,.3.....3.4...7.8.9 p q q The maximum difference of the Bézier coefficients is δ =.94. Bounding polynomials M and m: Root of M and m: Intersection intervals: M =.44X +.9493X +.87 m =.44X +.9493X +.877 N(M) = {.8, 4.838} N(m) = {.8, 4.838} 3
.3.....3.4...7.8.9 No intersection intervals with the x axis. 4
.9 Result: Root Intervals Input Polynomial on Interval [, ] p = X 3X 4 X 3 + 4X X +....3.4...7.8.9. Result: Root Intervals [.88,.88], [.44788,.4489] with precision ε =..