理学部理学科数学プログラム(数理物質科学専攻 数理科学コース) 田中研究室(最適化理論,オペレーションズ・リサーチ)

最適化数学A / B (251S1517 / 252S1527) (Optimization Mathematics A / B)
使用教科書:田中環著,「もっと知りたくなる 最適化数学の基礎」,培風館,ISBN:978-4-563-01171-0 (発売が少し遅れているかもしれません。それまではプリントを配布します。)
最適化数学A (251S1517) のまとめと筆記試験は 6月4日(水) の予定です。
最適化数学B (252S1527) のまとめと筆記試験は 7月30日(水) の予定です。
[理学部棟 B303]
小テストの解答や課題レポートについては,学務情報システムの 小テスト機能,レポート機能を利用して提出して下さい。
小テストの解答用紙やレポートをPDFで提出する仕方
  1. 表紙は不要。解答は手書き。 小テストの添削を提出するときは用意された解答用紙を使うこと。
  2. レポートの場合はノート(ルーズリーフなど)の一番上に 在籍番号,氏名を記載すること。
  3. スキャンは,プリンタかコピー機でできる。 スマホしかなくてもOfficeLensなどのアプリを使用すれば可能。 ただし,出力形式をPDFに指定すること。 なお,パソコン上で他の形式のファイルをPDFにする方法は色々あるが, 印刷の時にプリンターのところで「Microsoft Print to PDF」を選んべばよい。 なお,その際に用紙のサイズを指定すること。
  4. 出力サイズはB5用紙かA4用紙サイズ。 PDFを表示させた時,ズームで100%表示を選ぶと読み取った元の大きさが分かる。 プロパティからもページサイズを確認することが可能。
  5. 提出するレポートはなるべく1つのファイルに結合したり, ファイル容量を小さくする工夫を行うこと。 仮に,PDFの結合用のソフトが無くても, OfficeLensの場合は追加のボタンで複数のスキャンを1つのPDFとして出力できる。 また,OfficeLensを使わなくても,既に複数の画像ファイル(JPEG,GIF,PNGなど)があった 場合には,それらを複数選択して(Ctrlキーで選択), 右クリックメニューでポップアップメニューから「印刷」を選択すれば,上記に書いたように プリンターのところで「Microsoft Print to PDF」を選ぶと自動的に1つのPDFが作成できる。
  6. 提出するファイル名は, (半角英数の)学籍番号で始まる名前にすること。
  7. 提出締め切りは「日曜日まで」。 ただし,小テストやその添削については,学務情報システムに記載してある日時が優先される。 もちろん,間に合わなかったときはメールで提出しても構わない。

WordやExcel等のファイルをPDFに変換する仕方
  1. WordやExcelを使ってPDF変換:保存するときに,PDFとして保存すればよい。
  2. Windows10の仮想プリンター [Print to PDF] を使ってPDF変換: WordやExcelで印刷するときに,「Microsoft Print to PDF」を選択すればよい。

[1]配布資料と学習内容
  1. 講義案内と授業予定
  2. 4/9 は,高等学校で習う,内積の定義を見直し, スライド[1]に基づいて, 連立方程式(不等式)の解の集合を手際よく図示できるように考えるコツを学習します。 また,超平面と半空間の概念を学習します。
    なお,教科書がまだ入荷されていないかもしれませんので, 参考資料として,参考資料[1] も利用して下さい。
  3. 4/16 は,前回の復習を簡単に行って, すぐに 小テスト[A1](4/16) を行いました。 その後, 参考資料[2](教科書37--48頁の内容)に基づいて, 数理計画問題の紹介とその内,制約関数に制限をつけて線形計画問題を説明します。 あらかじめ,該当するところをよく読んで理解してきてください。 それと,スラック変数を入れて不等式を等式に変形する事も学びます。 次回は,基底解,実行可能解,実行可能基底解などの概念とその線形代数的な意味付けを学びます。 また,連立一次方程式や連立一次不等式の表す領域(実行可能領域)についての 基本的な考え方とその上での線形関数の最大・最小を学習します。
  4. 4/23 は, 参考資料[2]の残りを説明し線形計画問題の」標準形について学習しました。 教科書の写し(pp.41--48)を使用して, 線形計画問題の実行可能領域が有限個の頂点を持つ凸多面体(polytope)である時は,最適解はその頂点の いずれかとなり,最大値あるいは最小値は必ず存在することを説明しましたが, 証明は次回となりました。教科書の写しをよく読んでみてください。 また, 実行可能領域が凸多面体でない場合,つまり,有界な集合でなく,無限遠まで領域が 広がっている場合は,目的関数の法線ベクトル(勾配ベクトル)の向きによって, 最大値や最小値が存在しないこともあることについても次回となります。
    その後,具体例として 小テスト[A2]で 線形計画問題の実行可能領域が有限個の頂点を持つ凸多面体である時には,端点で最小値,最大値をとる ことを確認しました。
    その他,スラック変数を挿入して標準形に直す方法を確認します。 その他の用語(基底行列,基底変数,非基底行列,非基底変数など)については 各自で復習しておいて下さい。
    次回は,4月30日(水)です。 それまでに, シンプレックス法の手続きとその具体例について, 教科書の写し(pp.49--54) を参考にして理解して下さい。
  5. 4/30 は, 線形計画問題の標準形に直す方法を復習してから, 小テスト[A3] を実施します。
  6. 6月4日(水)まとめと筆記試験

[2]確認問題(小テスト)とレポート問題や練習問題
  1. レポート課題[A1](4/9)
    スライド[1]の課題をレポートとしています。 線形計画問題の不等式を満足する領域を内積の考え方を使って図示し(各不等式の法線ベクトルを記載する), 目的関数の等高線を表示して(法線ベクトルを記載する),最大値を与える解とその時の値を求めよう。
  2. 小テスト[A1](4/16)その解答
    内積を利用して不等式領域(半空間)を図示する問題。 法線ベクトル(勾配ベクトルともいう)と半空間の関係が理解できているかを見る問題。 半空間の境界(平面の時は直線)を超平面というが,それと法線ベクトルは直交している。
  3. 小テスト[A2](4/23)
    線形計画問題の実行可能領域が有限個の頂点を持つ凸多面体である時の最大・最小に関する問題。
  4. レポート課題[A2](授業中に途中まで作業します)
    簡単な線形計画問題に対してスラック変数を用いて標準形へ変形したり,実行可能基底解をすべて求めたり, グラフから最適解を直感的に求めたりする問題。 また,シンプレックス法で解いて最適解と最適値が他の方法で求めたものと同じくなることを確かめてみよう。
  5. 小テスト[A3](4/30)その解答例
    不等号の向きが反対だったり,非負条件のない一般の線形計画問題を 標準形(最小化問題,等式制約,すべての変数に非負条件あり)に直す問題。
  6. 小テスト[A4]と その解答例(answer example)
    初期の実行可能基底解を適当に見つけ,シンプレックス法で最適解とその時の最適値を求める問題。 実行可能基底解は3通りあるので,どれから出発しても良いので,一通りだけの解答でよい。
  7. 課題[A3][Report-A3]と その解答例
    線形計画問題をシンプレックス法で解いて最適解と最適値を求める問題。
  8. 小テスト[A5]と その解答

    人為変数を挿入して2段階法のシンプレックス法で解く線形計画問題。 最初に,2段階法で省エネ型として,既に与えられている。 2本の連立方程式なので,どの2つの列を組合せて基底行列になるように選べばよいかを考えればよい。
    第1段階の人工問題の最適値が0となる問題となっているのだが,ちゃんと0となることを シンプレックス法で確かめて,その後,その最後の表の元の変数に当たる部分の係数行列と右辺を 取り出し,目的関数を元の関数の係数に取り換えて,あらためて第2段階をシンプレックス法で 解くことになる。
  9. 練習問題と その最適解・最適値のみ
  10. 小テスト[B1]と その解答および その解説
    与えられた線形計画問題からその双対問題を作る問題。(2)についてはいくつも方法があるが, 式を整理すると同じものになる。昨年の授業の動画に詳しい解説があります。
  11. 小テスト[B2]と その解答
    与えられた線形計画問題を標準形に直し,その双対問題を作る問題。 また,双対問題が2変数の線形計画問題になるので,容易に最適解と最適値を求めることができる。 もちろん,シンプレックス法で解いても答えは同じになる。 授業では与えられた問題(LP)を標準形にして双対問題を求めたが, 一般形の不等号の場合に対する最大化の双対問題にしてもよい。この場合は, 見かけは異なって見えるが,同値な問題である。
  12. 課題[B2]
    具体的な行列$A$とベクトルを与えて, ファルカスの補題とゴルダンの定理について, 二者択一の定理の一方が成立して他方が成立しない例をそれぞれ図を書いて説明せよ。
  13. 小テスト[B3]と その解答例<
    複数の2変数実数値関数からなる制約条件の下で,2変数実数値関数の最適性条件を確認する問題。

[3]2022年のオンライン授業動画(履修登録者のみ視聴可能) YouTubeの32名までの制限により,学務情報システムの学生のアカウントを各回ごと2つのリンク先のどちらかに登録してあります。 両方には登録できないので,一方がだめなら他方を試してみてください。
  1. [A]第1回(令和4年4月13日)のYouTube動画
  2. [A]第1回(令和4年4月13日)のYouTube動画
  3. [A]第2回(令和4年4月20日)のYouTube動画
  4. [A]第2回(令和4年4月20日)のYouTube動画
  5. [A]第3回(令和4年4月27日)のYouTube動画
  6. [A]第3回(令和4年4月27日)のYouTube動画
  7. [A]第4回(令和4年5月11日)のYouTube動画
  8. [A]第5回(令和4年5月18日)のYouTube動画
  9. [A]第6回(令和4年5月25日)のYouTube動画
  10. [A]第7回(令和4年6月1日)のYouTube動画
  11. [B]第2回(令和4年6月22日)のYouTube動画
  12. [B]第3回(令和4年6月29日)のYouTube動画
  13. [B]第4回(令和4年7月6日)のYouTube動画
  14. [B]第5回(令和4年7月13日)のYouTube動画
  15. [B]第6回(令和4年7月20日)のYouTube動画
  16. [B]第7回(令和4年7月27日)のYouTube動画

[4]関連リンク
  1. 数理計画法パッケージ Nuorium Optimizer
    株式会社「NTT DATA 数理システム」で開発・販売されている数理計画ソフトウェアで, 「数式ベースの自然で簡潔な記述ができるモデリング言語」や「様々な実用的な求解アルゴリズム」 などを含んだ汎用数理最適化パッケージです。
  2. 数式処理・数式モデル設計環境 Maple (メイプル)
    Maplesoft社の数式処理ソフトMapleの概要と事例を見ることができます。
  3. (社)日本オペレーションズ・リサーチ学会
    日本OR学会のホームページでは, 日本でのオペレーションズ・リサーチ(OR)関係の活動情報が見れます。 2022年度の秋季研究発表会は新潟でハイブリッド開催されました。
  4. 第28回RAMPシンポジウム
    日本OR学会の数理計画研究部会 (RAMP: Research Association of Mathematical Programming) が 開催する,最適化や数理計画に関するシンポジウムです。2016年度は新潟大学で開催されました。

トップへ
トップへ


新潟大学 理学部数学プログラム/大学院自然科学研究科数理物質科学専攻 数理科学コース
E-mail: prtana@gs.niigata-u.ac.jp