NumRecipe
目的
「数値計算」の基礎的な事項をPython codeを通して理解する. さらに,グループ評価を通じて,AM/PM的な学修スキルを身につける.
単位について
全般の注意
- フォーマンセルで取り組む.
座席指定. - フォーマンセルですが,基本は4人,4人以下ならOK, 5人以上は不可.
- 期末試験の結果をフォーマンセルで平均する.
- 授業中に与えた課題に解答あるいはまともな質問をしたフォーマンセルに加点(合計20点程度).
- 章末の演習課題を提出したフォーマンセルに加点(2点X9回程度),
- 現在の各チーム得点
- 過去年度の得点
- 相当量の印刷物を配布予定(100page以上).試験に備えてバインダーに保存しておくこと.
章末の演習課題について
- 次週の前日の夜中までにLUNAにpdfで提出.
- 全部やる必要はない.2,3問を指定します.
- 採点結果は私のURLへup.
- レポートっぽく色々試したり,追加してもらうと加点.
ペア(ツーマンセル)フォーマンセルで1つ出せばよい.- 採点の都合上,遅れたら,0点.
期末試験について(100点満点)
- NumRecipe_ExamTA_manual
- 7/10の講義時間中に,対面で実施.
- 講義室で実施
- 座席指定.時間中5分間だけ相談時間を設ける予定.
大学のコンピュータを使って実際に数値計算をする課題がある.- BYODでcodeを実行,LUNAへ提出.
- メール(含むLINE等のSNS)の使用は一切不可.見つかれば即退場.webの利用はOKだけど,「教えて」なんかはダメよ.
- 資料などの持ち込みは,プリント,ノート,コピーおよびレポート.
- USBメモリーの試験時間中の使用は不可.必要な資料や課題のファイルはあらかじめ自分のフォルダーに入れておくように
- 章末の演習課題と重複する問題が多い.よく消化しておくように.
講義の目的:数値計算をPythonで解説
数値計算
「数値計算」が対象とする分野は,関数の解,積分,微分方程式,固有値問題などで,数学の問題を解析的(analytical) ではなく,数値的 (numerical) に解く手法の集大成.数値計算には,料理と同じで,正しい調理法 (レシピ,recipe) がある.
基本的な戦術
数値計算が必要となる問題を解くときの戦術は,
- 手持ちのツールで解いてみる,
- 解析的な解を捜す,
- 既存の数値計算のサブルーチンを写して,利用.
- 使いやすいライブラリを捜す,
- 自分でサブルーチンを考える,
という順でおおむね進める. 数値計算のサブルーチンを実際に自分で考えるということはほとんどない. しかし,例え 3. や 4. の既存のサブルーチンに頼るときにも, ブラックボックスの中で何がおこなわれているかを大まかに理解しておかないと 大失敗をしでかすことがある.
といっても,数値計算手法の全てを理解しておく必要はない. 料理において基本となる,切る,煮る, 焼くなどを知っていれば,recipe を見ながら調理ができる. 数値計算においても同じで,基本となる誤差,精度,収束性,安定性,計算速度などの本質を理解しておけば, 新しい数値計算手法もだいたいの振る舞いを予測することが出来る.
本講義ではこのような視点にたって, いくつかの典型的な数値計算手法の基礎的な考え方と実際のプログラムを紹介する.
python
- jupyter notebookで資料を提供
- github の数値計算
- nbviewer notebookのviewer
- 自習するならanaconda
- https://www.anaconda.com でinstallするのが便利.
- 大学のpcにはサブセット(miniconda)が入っている.
- 情報準備室提供installer (Windows用)
- エディターはvscode(Microsoft謹製)が便利
- 環境設定はほぼ自動です.
- インストールしますかと尋ねられたのをyesしていく
スケジュール予定
初回環境整備と諸注意
講義内容
date | ipynb | viewer | マーク付きpdf | 解答例 | 動画 | |
day 0(4/10) | 授業進行説明 | テキスト | 講義手順 | |||
day 1(4/10) | 常微分方程式 | viewer | テキスト | テキスト | リンク | |
day 2(4/17) | 代数方程式 | viewer | テキスト | テキスト | テキスト | リンク |
day 3(4/24) | 誤差 | nbviewer | テキスト | テキスト | テキスト | リンク |
day 4(5/01) | 写像 | nbviewer | テキスト | テキスト | テキスト | リンク |
day 5(5/08) | 逆行列 | nbviewer | テキスト | テキスト | テキスト | リンク |
発展課題 | nbviewer | テキスト | ||||
day 6(5/15) | 固有値とライブラリ | nbviewer | テキスト | テキスト | テキスト | リンク |
計算速度競争 | nbviewer | テキスト | ||||
day 7(5/22) | 補間と数値積分 | nbviewer | テキスト | テキスト | テキスト | リンク |
day 8(5/29) | 線形最小二乗法 | nbviewer | テキスト | テキスト | テキスト | リンク |
day 9(6/05) | 非線形最小二乗法 | nbviewer | テキスト | テキスト | テキスト | link |
day 10(6/12) | FFT | nbviewer | テキスト | link | ||
day 11(6/19) | 常微分方程式 | nbviewer | テキスト | link | ||
MD and MC | ||||||
day 12(6/26) | 期末試験リハ | nbviewer | テキスト | 解答例 | ||
day 13(7/3) | 予備日 | |||||
day 14(7/10) | 期末試験本番 | nbviewer | テキスト | ipynb pdf | ||
- 第5回の発展課題(乳がん検出器)のdata train_A.data, train_b.data, validate_A.data, validate_b.data
- svdによるfitting
期末試験関連資料
年度 | 問題 | 内容 | 答案 |
23年度 | ipynb, pdf | 行列,差分商,GaussSeidel, PageRank | 答案, pdf |
22年度 | ipynb, pdf | 行列,漸化式の誤差,積分, Verlet | 答案, pdf |
21年度 | ipynb, pdf | 行列,差分商,積分, スムースな静⽌ | 答案, pdf |
20年度 | ipynb, pdf | fit, 数値解の収束性, 常微, fft | .ipynb, pdf |
19年度 | 行列, fft, 解の収束性, ODE | nbview, pdf | |
18年度 | 行列,解収束性,精度・誤差,最小二乗,常微 | nbview | |
17年度 | 行列,精度・誤差,積分,最小2乗法,常微 | nbview | |
16年度 | 収束性,Gauss-Seidel,FFT強度表示,tanニュートン差分商 | nbview | |
15年度 | 収束性,丸め誤差,Newton差分商,ページランク | nbview | |
14年度 | |||
11年度 |
来年度に向けて
- プリントアウトで不正が起こるかも.試験の提出はLUNAから.
- 学期の始めに,以下の項目を明言する
- コピペでいいこと.
- 何を削るべきか?
- どこを換えるべきか?
- どのような操作がJupyter Notebookで最低限修得しないといけないか?
- チームの結束.
- 課題の意味.
参考文献
- William H. Press 他著「ニューメリカルレシピ・イン・シー C言語による数値計算レシピ」(技術評論社,1993)
- 数値計算のバイブル.原著 Numerical recipeでは非常に広範な計算対象に対して,C, Fortran, C++, Pascal, Basic版が用意されており,数値計算プログラムをコーディングする際の洗練されたスタイルも提示している.記述は初学者には難しいが,ある程度経験を積んだプログラマには,手法を選ぶうえで非常に役に立つ情報である.
- 奥村晴彦著「C言語による最新アルゴリズム事典」(技術評論社,1991)
- いわゆる数値計算に限らず,いろいろな計算機問題とその解法が載っている.数値計算についても必要最小限の記述とプログラムがまとめられており,非常に便利.改訂版が2018に出た.
- 数値解析法,森正武著,1984, 朝倉書店朝倉現代物理学講座7
- 数値解析,一松信著,1982,朝倉書店新数学講座13
- 両書とも少し古いが,数値計算の基礎となる理論的な説明が明解.
- 河村哲也著「数値計算の初歩!」(山海堂,2002年)
- レベルは高くないが,重要な数値計算の初歩を丁寧に解説.
Keyword(s):
References:[LectureNotes]