[論文レビュー] On the explicit representation of orthonormal Bernstein polynomials

Michael A. Bellucci|arXiv (Cornell University)|Apr 8, 2014

Advanced Mathematical Theories and Applications参考文献 27被引用数 14

ひとこと要約

本稿では、グラム・シュミット直交化を用いて導出された、標準的（非直交）ベルヌーイ多項式の線形結合としての、直交ベルヌーイ多項式の明示的解析的表現を初めて提示する。これらの直交多項式がストゥルム＝リウヴィル固有値問題を満たし、低次の基底集合を用いて高精度な曲線・曲面近似を可能にするとともに、ベジエ曲線・曲面モデリングにおける数値的安定性と制御点計算の改善を実現する。

ABSTRACT

In this work we present an explicit representation of the orthonormal Bernstein polynomials and demonstrate that they can be generated from a linear combination of non-orthonormal Bernstein polynomials. In addition, we report a set of $n$ Sturm-Liouville eigenvalue equations, where each of the $n$ eigenvalue equations have the orthonormal Bernstein polynomials of degree $n$ as their solution set. We also show that each of the $n$ Sturm-Liouville operators are naturally self-adjoint. While the orthonormal Bernstein polynomials can be used in a variety of different applications, we demonstrate the utility of these polynomials here by using them in a generalized Fourier series to approximate curves and surfaces. Using the orthonormal Bernstein polynomial basis, we show that highly accurate approximations to curves and surfaces can be obtained by using small sized basis sets. Finally, we demonstrate how the orthonormal Bernstein polynomials can be used to find the set of control points of Bezier curves or Bezier surfaces that best approximate a function.

研究の動機と目的

• 先行文献にない形式を欠いている現状を克服し、直交ベルヌーイ多項式の明示的・閉形式表現を確立すること。
• 直交ベルヌーイ多項式とストゥルム＝リウヴィル固有値問題との関係を確立し、自己随伴作用素構造を示すこと。
• 一般化フーリエ級数を用いた関数近似において、直交ベルヌーイ多項式の有用性を実証すること。
• 直交基底を活用することで、ベジエ曲線・曲面の最適制御点計算を簡素化すること。
• 直交化されたベルヌーイ基底に用いられる悪条件な変換行列の代替として、数値的に安定な手法を提供すること。

提案手法

• 次数を段階的に増加させた基底集合に対して、グラム・シュミット直交化を用いて、非直交ベルヌーイ多項式の線形結合として直交ベルヌーイ多項式を明示的に導出する。
• 直交ベルヌーイ多項式が満たすストゥルム＝リウヴィル微分方程式系を特定し、係数として p(x) = x(1−x)²、q(x) = n(n+2)(1−x)、w(x) = 1、固有値として λ = (n−j+1)(j−n) を得る。
• 区間 [0,1] 及び [0,1]² 上の関数近似に、直交ベルヌーイ多項式を用いた一般化フーリエ級数を定式化する。
• 直交基底を用いて、近似誤差を最小化するようにベジエ曲線・曲面の制御点を計算する。
• ベンチマーク関数（リサージュ曲線、sinc曲面、ランゲルマン曲面）を用いて性能を検証し、異なる次数における誤差を比較する。
• 最適近似精度が低次の基底（例：リサージュ曲線では n=20、曲面では n=12–13）で達成されることを示し、誤差は 5.2×10⁻⁶ 未満である。

実験結果

リサーチクエスチョン

• RQ1 反復的直交化に依存しない、直交ベルヌーイ多項式の明示的解析的表現を導出可能か？
• RQ2 直交ベルヌーイ多項式はストゥルム＝リウヴィル固有値問題を満たすか？ もしそうであれば、関連する微分作用素と固有値は何か？
• RQ3 曲線・曲面フィッティングにおいて、直交ベルヌーイ多項式の近似精度は非直交基底と比べてどの程度優れるか？
• RQ4 直交基底を用いることで、ベジエ曲線・曲面の最適制御点計算が簡素化可能か？
• RQ5 リサージュ曲線やsinc曲面のような複雑な関数を高精度に近似するための最小基底次数は何か？

主な発見

• 直交ベルヌーイ多項式は、非直交ベルヌーイ多項式の線形結合として明示的に表現され、文献上初の閉形式表現を提供する。
• 次数 n の各直交ベルヌーイ多項式は、自己随伴微分作用素を備えたストゥルム＝リウヴィル固有値問題の解であり、スペクトル的性質が裏付けられる。
• リサージュ曲線では、次数 n=20 の直交基底を用いた近似誤差が 2.1×10⁻⁷ にまで低下し、同じ次数の非直交基底と同等の精度を達成した。
• sinc曲面では、12×12 の直交基底を用いて誤差が 5.1×10⁻⁶ にまで低下したが、ランゲルマン曲面では n=13 が必要で誤差は 2.0×10⁻⁶ であった。
• 従来手法と比較して、著しく小さな基底集合で高精度な近似が可能となり、悪条件な変換行列に起因する数値的不安定性が低減される。
• 直交性のおかげで、最小二乗フィッティングにおける線形方程式系の条件数が改善され、ベジエ曲線・曲面の制御点最適化が簡素化される。