第3回レポート問題

2016 年度・前期・数理解析・計算機数学3
1
● 第3回レポート問題
問題 C-1 は必ず解答すること. それ以外の問題から1問を選び解答すること.
★ 問題
すべての問題で, 適切な例に対するプログラムを書くこと. すべての計算(議論)は倍精度浮動
小数点数で計算すること. また, 行列 N × N 行列 DN , 行列 TN , ベクトル bN ∈ RN は以下で定
義されたものとし, プログラムでは, N = 10 および N = 20 で正しい結果が得られる必要がある.
さらに, DN は全ての固有値が相異なることは仮定してよい.
C-1 (20 点) 連立一次方程式 DN bN = x を, 以下の方法で解くプログラムをそれぞれ書きなさい.
方法: Gauss の消去法, LU 分解を行った上で解く方法, Jacobi の反復法, Gauss-Seidel 法,
SOR 法. ただし, SOR 法では ω = 1.75 と設定してよい.
C-2 (10 点) 行列 TN の核の基底を一組求めるプログラムを書きなさい.
C-3 (10 点) 行列 DN の最大固有値および最小固有値と, それに対する固有ベクトルを求めるプロ
グラムを書きなさい.
C-4 (20 点) 行列 DN の全ての固有値と固有ベクトルを求めるプログラムを書きなさい. ただし,
固有ベクトルは, 第 N 成分が 1 となるものを出力すること.
C-5 (25 点) 正定値エルミート行列 A に関する連立一次方程式を解く SOR 法は, ω ∈ (0, 2) なら
ば収束することを証明しなさい.






sin( Nπ+1 )
−2 1
−2 1
1






..


 1 −2 1

 1 −2 1

.












.
.
.
.
.
.
kπ

..
..
..
..
..
..
DN = 
 TN = 
 , bN =  sin( N +1 ) 







..






1 −2 1 
1 −2 1 


.


1 −2
1
1 −2
sin( N π )
(N +1
★ 締め切りと提出方法
締め切りは2016年7月27日(水曜日)17時とする.
• メールの送付先は [email protected] し, Subject には, 文字列「レポート提出」
を含むこと. また, Subject には問題番号を書くこと.
• 電子メールの発信者は, “nagoya-u.ac.jp” のアドレスであること.
• プログラム, PDF ファイル, 図版の電子データはメールで提出すること.
• 図版および数学的な内容を記述したレポートは, PDF ファイルにまとめて提出することが望
ましいが, それができない場合には紙に書いて(図版はプリンタで出力して)提出してもか
まわない.
Jul. 13, 2016, Version: 1.0
[email protected]
2016 年度・前期・数理解析・計算機数学3
2
• プログラム中には「問題番号, 学籍番号, 氏名」をコメントとして記入すること. PDF ファ
イルまたは紙で提出するレポートには, 「問題番号, 学籍番号, 氏名」を記入すること.
• PDF ファイルの用紙サイズ, および, 紙で提出する用紙サイズは A4 とする. (特に理由が
ない限り, A4 縦おきであること.)紙で提出する場合には, 片面であること.
• 図版を単独の電子データで送付する際のファイルフォーマットは, PDF, EPS, SVG, JPEG,
PNG のいずれかとする.
★ 採点基準
当初に説明した通り, 以下の基準で採点するが, 中間的な点をつけることもある. また, 非常に良
いレポートにはボーナス点を与える. (A:100%, B:60%, C:0% と考えて良い.)なお, 今回の満点
は 30 点とし, それを越える点はボーナス点とする.
A 数値計算が正しいプログラムによって行われ, その考察が適切であるもの.
B 数値計算は正しいプログラムによって行われているが, その考察に少々問題があるもの.
または, 数学的背景は正しく記述されているが, プログラムに小さな誤りがあるもの.
C それ以外.
★ 注意事項
以下の注意事項に沿わないレポートは採点から除外する.
• 提出方法に従わないレポートは採点の対象とはしない. 特に, 電子ファイルのファイルフォー
マットに注意すること.
• 提出するレポートおよびそのプログラムは, 自分で作成したものに限る. 特に, 他人書いても
らったプログラムや, 「インターネット」で拾ってきたプログラムをレポートとして提出し
てはならない.
このような場合には, 採点から除外するだけでなく, 不正行為として直ちに不可と判定する.
• プログラムは, 浮動小数点演算が IEEE の規定に沿う言語で記述すること. (C なら問題な
い.)また, プログラム中には, 何をしているプログラムか, 各部分が何を実行しているのかを
コメントとして記入すること. コメントが全くないプログラムは採点対象としない. コメン
トが不適切であったり, コメントが無いことにより何をしているのかが読み取れないプログ
ラムは減点対象とする.
★ 一言
レポートは, できたからと言って, 「よしできたぁ!」とすぐに出すのはやめましょう. 数日後に
もう一度自分の書いたレポート(プログラム)を見直して, 自分で何を書いているかがわかるかを
考えてください.
以上
Jul. 13, 2016, Version: 1.0
[email protected]