著者

沼田宗敏 輿水 大和 秦野 やす世 神谷 和秀 野村 俊 二宮 市三

雑誌

情報処理学会論文誌 (ISSN:18827764)

巻号頁・発行日

vol.48, no.12, pp.4051-4058, 2007-12-15

---

従来,直線,円,楕円などの幾何学図形の生成法が多数報告され,傾斜楕円についても若干の研究業績がある.本論文では,1 点あたり2 回の乗算で傾斜楕円を生成する画期的な三項漸化算法を提案する.三角関数の間に成り立つ一対の漸化式が同型であることが,その秘訣である.初期設定や楕円弧描画は従来法と比較して簡単であり,それを従来の高速生成法と同等以上の計算速度で実現する.また,実験によって新算法の有効性を検証する.

著者

吉田 年雄 二宮 市三

雑誌

情報処理学会論文誌 (ISSN:18827764)

巻号頁・発行日

vol.21, no.3, pp.238-245, 1980-05-15

---

Xが小さい場合について 第2種変形ベッセル関数瓦K_ν(x)の新しい数値計算法を提案している.K_ν(x)の定義式K_ν(x)=(π/2)・(I_<-ν>(x)-I_ν(x))/sinνπに 第1種変形ベッセル関数I_ν(x)およびムI_ν(x)のベキ級数展開を代入し 適当に項をまとめ 桁落ちを生ずる部分を所要の精度を有する最良近似式にて計算することにより K_ν(x)を精度良く しかも能率的に計算している.例えば (0 0.5)では ν=0を中心とする展開として得られるK_ν(x)=(π/(2sinνπ))・㊥^^∞__<k=0>{A_k(ν x)+B_k(ν x)}により計算を行う.ただし A_k(ν x)の式中には 関数(1/Γ(k+1-ν)-1/Γ(k+1+v))/(k!)を含むが そのまま式どおりに計算を行うと桁落ちを生ずる.そこで これに対する最良近似式を新たに作成し それによりその関数値を計算する.また A_k(ν x)には Φ_1(ν x)=(x/2)^<ν-1>-1およびΦ_2(ν x)=1-(x/2)^νを含むが そのまま計算すると桁落ちを生ずるときには f(t)=e^t-1なる関数の最良近似式によりΦ_1(ν x)=f(-νln(x/2)) Φ_2(ν x)=-f(νln(x/2))として計算を行う.このようにして 桁落ちが無く しかも能率的にK_ν(x)の値を計算する方法を述べている.

著者

吉田 年雄 二宮 市三

出版者

一般社団法人情報処理学会

雑誌

情報処理学会論文誌 (ISSN:18827764)

巻号頁・発行日

vol.21, no.3, pp.238-245, 1980-05-15

---

Xが小さい場合について 第2種変形ベッセル関数瓦K_ν(x)の新しい数値計算法を提案している.K_ν(x)の定義式K_ν(x)=(π/2)・(I_<-ν>(x)-I_ν(x))/sinνπに 第1種変形ベッセル関数I_ν(x)およびムI_ν(x)のベキ級数展開を代入し 適当に項をまとめ 桁落ちを生ずる部分を所要の精度を有する最良近似式にて計算することにより K_ν(x)を精度良く しかも能率的に計算している.例えば (0 0.5)では ν=0を中心とする展開として得られるK_ν(x)=(π/(2sinνπ))・㊥^^∞__<k=0>{A_k(ν x)+B_k(ν x)}により計算を行う.ただし A_k(ν x)の式中には 関数(1/Γ(k+1-ν)-1/Γ(k+1+v))/(k!)を含むが そのまま式どおりに計算を行うと桁落ちを生ずる.そこで これに対する最良近似式を新たに作成し それによりその関数値を計算する.また A_k(ν x)には Φ_1(ν x)=(x/2)^<ν-1>-1およびΦ_2(ν x)=1-(x/2)^νを含むが そのまま計算すると桁落ちを生ずるときには f(t)=e^t-1なる関数の最良近似式によりΦ_1(ν x)=f(-νln(x/2)) Φ_2(ν x)=-f(νln(x/2))として計算を行う.このようにして 桁落ちが無く しかも能率的にK_ν(x)の値を計算する方法を述べている.