Download document (10.35 MB) - Escuela Superior de Informática ...

3.2. Técnicas de representación 3D 45anterior, obteniendo:⎛⎞ ⎛⎞⎛⎞⎜⎝1 w 01 w 11 w 2...1 w m−21 w m−11⎟ ⎜⎠ ⎝D 0D 1D 2...D m−2D m−1D m⎟⎠=⎜⎝δ 1δ 2δ 3...δ m−2δ m−1δ m⎟⎠así, la matriz que contiene las ω i únicamente tiene dos diagonales que no son ceros (acelerandolos cálculos un 33 %). Los valores de ω y δ se calculan mediante las expresiones:ω 0 = 1 21ω i =(4 − ω i−1 )1ω m =2 − ω m−1 )∀i ∈ [1, m − 1]δ 0 = 3(y 1 − y 0 )ω 0δ i = [3(y i+1 − y i−1 )δ i−1 ]ω i ∀i ∈ [1, m − 1]δ m = [3(y m − y m−1 )δ m−1 ]ω mLos casos particulares de cálculo son el primer y último punto de control. El resto, serán calculadosde forma general por la expresión de término i.Resolviendo la ecuación matricial, se aprecia que la derivada primera en el último puntode control se puede resolver directamente calculando el valor de δ. Será lo que se utilice en laimplementación del algoritmo; primero se resolverá el término m, y volviendo atrás hacia elprimero, se aplicará la ecuación general de término i. Se obtendrá:D m = δ mD m = δ i − ω i D i+1 ∀i ∈ [m − 1, 1]