計算アルゴリズム

計算アルゴリズム 計算理工学専攻　張研究室山本有作

物理における数値計算 • 解析的に求められない式の計算 • 非線形方程式の解 • 積分で表示された式の値の計算 • データ解析 • 最小二乗法による実験データのフィッティング • フーリエ変換を用いた結晶解析 • シミュレーション • 分子動力学シミュレーション • 流体シミュレーション数値計算の必要性

数値計算とは • 数値計算とは • 方程式など数学的に定式化された問題の解を計算機を用いて求めること • 計算の対象 • 計算機のメモリ内に格納された数値 • 計算の手順 • 有限回の演算を用いて（近似的な）解を返すように設計された計算手順（アルゴリズム）

数値計算における有限性 • 表現できる数値の有限性 • 表現できる数値は実数の一部のみ • 丸め誤差の存在　→　今回の授業で解説 • 無限大・無限小が扱えない • 無限級数の和は有限級数で近似 • 微積分で現れる無限小は有限の微小量で近似 • 打ち切り誤差の存在　→　今回の授業で解説 • 計算能力・記憶領域の有限性 • 実用的な時間・記憶領域で計算するためには，効率的なアルゴリズムが必要 • 収束の次数の概念 • 計算量のオーダーの概念

授業で学ぶこと • 基本的な数値計算法 • アルゴリズム • 数学的原理 • アルゴリズムの振る舞いの解析 • 計算量 • 記憶領域 • 誤差 • 各アルゴリズムの適用範囲 • 各アルゴリズム使用の際の注意事項

授業の構成 • 数値計算における誤差 • 非線形方程式の解法 • 関数の補間 • 数値積分法 • 数値微分法と加速法 • 常微分方程式の解法 • 偏微分方程式の解法 • 連立一次方程式の解法 • 固有値の計算法

テーマ間の関連 偏微分方程式固有値の計算常微分方程式の解法数値積分法連立一次方程式非線形方程式関数の補間数値微分法と加速法数値計算における誤差

授業の進め方(1) • レポート • ２回出題 • 提出は授業後10日以内に5号館6階625号室まで • 中間テスト • 前半の講義内容について簡単な問題を出題 • 11/19 に実施予定

授業の進め方(2) • 期末テスト • 通常の筆記テスト（昨年度のHPを参照） • 出題範囲は本講義全体 • 出題募集（希望者のみ） • 期末試験には，受講者が出題した問題を含める。 • 出題者には，出題された問題の良さに応じて成績（A，B，Cなど）を保証 • 質問シート • 授業の終わりの５分間に，感想，質問，コメントなどを紙に書いて提出 • 全部で２、３回実施予定

授業の進め方(3) • 授業のホームページ • http://www.na.cse.nagoya-u.ac.jp/~yamamoto /lectures/algorithms2007/algorithms2007.html • 授業ノート，C/Javaプログラム，授業に関するお知らせ等は，すべてホームページに掲載の予定 • 質問など • 5号館6階625号室　山本まで • あるいはメールで yamamoto@na.cse.nagoya-u.ac.jp　まで

数値計算に関する参考書 • 水島二郎，柳瀬眞一郎：　「理工学のための数値計算法」，数理工学社，2002. • 杉浦洋：　「数値計算の基礎と応用」，サイエンス社，1997. • 伊理正夫，藤野和建：　「数値計算の常識」，共立出版，1985. • 山本哲朗：　「数値解析入門」，サイエンス社，1976. • 森正武：　「数値解析（第２版）」，共立出版，2002.

情報落ちの例 • S =Σi=1n(1/i2)の計算結果（単精度）・n = 4097で情報落ち。・nをそれ以上増やしても，　和の値は増加しない。 1.64493 0668 n = ∞の真値

情報落ちを防ぐ方法 (1) • S =Σi=n1(1/i2)の計算結果（単精度）・小さい値の項から足して　いくことで，情報落ちが　生じにくくなる。・nを増やすと，和の値も　順調に増加する。 1.64493 0668 n = ∞の真値

情報落ちを防ぐ方法 (2) • S =Σi=1n(1/i2)の計算結果（単精度と倍精度） 1.64493 0668 n = ∞の真値

計算アルゴリズム

計算アルゴリズム

Presentation Transcript