この用語『最適化』を理解するうえで併せて押さえたい関連キーワード群です。 クリック(ホバー)で関連用語ページに飛べます。
🍰 まずはやさしく
一番良い答えを探すパズルのようなものです。
最高の選択を数学的に決めるために使います。
スマホの充電を一番効率よく使う方法などです。
最適化の基本的な考え方と使い方を読みます。
最適化 (Optimization) とは、制約条件のもとで目的関数 $f(x)$ を最小化 (または最大化) する変数 $x$ を求める数学的枠組み。 機械学習・統計推定・経営判断・物流・電力網運用など、 あらゆる意思決定問題の数学的中核。
🍰 まずはやさしく
いろいろな計算手法をまとめた地図のようなものです。
複雑な問題を整理して解くために使います。
部活の予算をどう分けるか決める時に役立ちます。
最適化に含まれる多くの計算方法について読みます。
この用語ページは「最適化 (Optimization)」を、 数理計画 / 機械学習 / オペレーションズ・リサーチ の文脈で総合的に解説しています。 線形計画 (LP)、 二次計画 (QP)、 二次錐計画 (SOCP)、 半正定値計画 (SDP)、 混合整数線形計画 (MILP)、 凸計画、 非凸計画、 大規模分散最適化、 確率最適化、 ロバスト最適化などあらゆる最適化問題の上位概念。 勾配法、 線形計画法、 凸最適化、 ニュートン法、 ラグランジュ乗数 等を配下に持つ上位概念です。 SSDSE-B-2026 の県別予算配分や、 機械学習の損失最小化、 ポートフォリオ最適化など、 統計・データ解析コンペのあらゆる定量的意思決定の理論的バックボーンになります。
🍰 まずはやさしく
ルールの中で最高の選び方を見つけることです。
迷わずに正解にたどり着くために使います。
買い物で予算内で一番良い物を買うときと同じです。
問題を数式に直すための3つのステップを読みます。
最適化 (Optimization) は「決められたルールの中で、 一番良い選択を見つける」問題の数学化。 たとえば「47 都道府県に総予算 100 億円を配分し、 平均所得をなるべく押し上げる」という配分問題、 「物流倉庫から各県の集配センターに荷物を運ぶ最安ルートを決める」といった輸送問題、 さらに「年金資産をどの株式・債券に何 % 投資するか」のポートフォリオ問題、 これらはすべて同じ枠組みで書ける。
最適化を理解する 3 段階のメンタルモデル: ① 変数 (decision variable) は何か = 自分が決められる量、 ② 目的関数 (objective) は何を最大化/最小化するか = ゴール、 ③ 制約 (constraint) は何か = やってはいけないこと。 この 3 つを書き出せれば、 ほぼ全ての問題は数理最適化問題に翻訳できる。
本ページでは 最適化 を、 定義・主要クラス (LP/QP/SOCP/SDP/MILP)・Lagrange 双対・凸性・SSDSE-B-2026 実データへの応用 (47 県リソース配分問題)・Python による解法 (scipy.optimize、 cvxpy、 ortools) の順で網羅的に解説します。 機械学習における損失最小化も、 ポートフォリオ運用も、 物流計画も、 すべてこの枠組みの一例にすぎないことが分かるはずです。
最適化は『変数 $x$ を動かして目的関数 $f(x)$ を最小(または最大)にする点を探す』問題の総称であり、 探し方は問題の性質で大きく変わる。 目的関数と実行可能領域が凸なら谷は一つで大域最適が保証されるが、 非凸なら谷が複数あって局所解に捕まりうる。 変数が連続なら勾配法・ニュートン法・内点法、 離散(整数・組合せ)なら分枝限定法やメタヒューリスティクスと、 解法系統も分岐する。 さらに制約なし/制約ありの別で扱いが変わり、 制約付きでは Lagrange 乗数・KKT 条件が鍵になる。 機械学習の損失最小化も、 物流の輸送費最小化も、 すべてこの一つの枠組みの適用例にすぎない。
最適化(Mathematical Optimization)は単独で覚えるものではなく、 最適化 という大きな枠組みの中での位置づけを理解することで応用範囲が広がります。 本ページの『🌐 関連手法』『🔗 関連用語』『📚 グループ教材』を順に辿ると、 関連概念のネットワークが見えてきます。
特に SSDSE-B のような実データに当てはめてみると、 教科書では抽象的に語られる概念が『47 都道府県の現実』に紐付き、 数字の意味が腑に落ちやすくなります。 次の『🧮 実値で計算してみる』セクションでは、 公開統計データを使って手を動かす例を紹介します。
本用語に関連する代表的な可視化を 3 点示す。



🍰 まずはやさしく
最適化を数式で表したルールブックです。
計算で正確な答えを出すために使います。
テスト勉強の時間をどう分けるか決めるようなものです。
目的関数や制約といった専門的な定義を読みます。
最適化問題 とは、 決定変数 $x \in \mathbb{R}^n$ を選んで、 目的関数 $f: \mathbb{R}^n \to \mathbb{R}$ を最小化 (または最大化) する問題で、 一般形は以下のように書ける:
$$ \min_{x \in \mathbb{R}^n} f(x) \quad \text{subject to} \quad g_i(x) \le 0, \; i=1,\ldots,m, \quad h_j(x) = 0, \; j=1,\ldots,p $$
ここで $f$ は 目的関数 (objective)、 $g_i$ は 不等式制約、 $h_j$ は 等式制約。 制約を満たす $x$ の集合を 実行可能領域 (feasible region) $\mathcal{X}$ という。 最小化する $x^* \in \mathcal{X}$ を 最適解、 $f(x^*)$ を 最適値 という。
主要分類: ① $f, g_i$ がすべて線形 → 線形計画 (LP)、 ② $f$ が二次・$g_i$ が線形 → 二次計画 (QP)、 ③ $f$ と $\mathcal{X}$ が凸 → 凸計画、 ④ 一部の $x$ が整数 → (混合) 整数計画 (MIP/MILP)、 ⑤ $f$ が非凸 → 非凸最適化。 英語名 Mathematical Optimization / Mathematical Programming。 同義語: 数理最適化、 数理計画、 OR (オペレーションズ・リサーチ) の中核技法。
線形計画 (LP) の標準形。 $c \in \mathbb{R}^n$ は目的係数、 $A \in \mathbb{R}^{m\times n}$ は制約行列、 $b \in \mathbb{R}^m$ は右辺ベクトル、 $x \ge 0$ は非負変数。 双対変数 $y \in \mathbb{R}^m$ は シャドープライス: 各制約を 1 単位緩めたときの最適値改善量。
SSDSE-B-2026 文脈の例: $x_i$ = 47 県への配分額、 $c_i$ = 県 $i$ への配分の不効用、 $A$ = 「総予算 100 億」「各県上限 5 億」等の制約、 $y$ = 「総予算 1 億円増やしたら全国の高齢人口救済がどれだけ増えるか」のシャドープライス。 KKT 条件: $\nabla_x L = 0$、 相補性 $\lambda_i g_i(x) = 0$、 双対実行可能性 $\lambda_i \ge 0$、 主実行可能性 $g_i(x) \le 0,\; h_j(x) = 0$ — これら 4 条件が最適性の必要条件。 凸問題なら十分条件でもある。
まずは本ページの『💡 30 秒で分かる結論』と『🎨 直感で掴む』で全体像を掴み、 次に『🧮 実値で計算してみる』を 手を動かして追体験するのが最短です。 数式や深い理論はその後で十分。
本ページの『🌐 関連手法・派生』『🔗 関連用語』で対比される手法を確認し、 それぞれの適用条件と得意・不得意を表で比較するのが効果的です。 SSDSE-B のような共通データセットで両方走らせて結果を見ると違いが体感できます。
最適化問題の計算量は問題の構造で決まり、 手法の巧拙より先に効きます。 線形計画は多項式時間(内点法で $O(n^{3.5}L)$ 程度)、 凸二次計画も多項式時間ですが、 整数を 1 つ入れた途端に整数計画となり NP 困難になります。 「47 都道府県を 4 ブロックに分ける」だけでも $4^{47} \approx 2 \times 10^{28}$ 通りで、 全探索は不可能です。 まず「自分の問題は凸か、整数変数が要るか」を判定するのが最初の一手で、 ここを見誤ると「良いソルバーを探す」方向に労力を注いで解けないままになります。 整数が避けられないなら、 厳密解をあきらめて近似・ヒューリスティクスに切り替える判断が要ります。
『点推定値』だけでなく『不確実性(CI、 SE、 分散)』『前提条件のチェック結果』『代替手法との比較』『データ取得日と seed』をセットで報告するのが標準。 査読・レビューで問われる典型ポイントです。
$$ \min_{x} f(x) \quad \text{s.t.} \quad g_i(x) \le 0,\; h_j(x) = 0 $$
| 記号 | 意味 | SSDSE-B-2026 県別配分文脈の具体例 |
|---|---|---|
$x \in \mathbb{R}^{47}$ | 決定変数 (自分が決められる量) | 47 県への予算配分額 $(x_1, \ldots, x_{47})$、 単位: 億円 |
$f(x)$ | 目的関数 (最小化対象) | $-\sum_i \alpha_i \log(1 + x_i)$ (収穫逓減型の効用最大化を反転) |
$g_i(x) \le 0$ | 不等式制約 (やってはいけない上限) | $\sum_i x_i \le 100$ (総予算 100 億)、 $x_i \le 50$ (1 県上限 50 億) |
$h_j(x) = 0$ | 等式制約 (必ず満たす条件) | $\sum_i x_i = 100$ (予算ぴったり使い切る) |
$\mathcal{X}$ | 実行可能領域 | 制約をすべて満たす $x$ の集合 (47 次元単体や hypercube) |
$x^*$ | 最適解 | $x^* = \arg\min_{x \in \mathcal{X}} f(x)$ — 最適な配分計画 |
Lagrange 関数 (制約を目的関数に取り込む技法):
$$ \mathcal{L}(x, \lambda, \mu) = f(x) + \sum_{i=1}^m \lambda_i g_i(x) + \sum_{j=1}^p \mu_j h_j(x), \quad \lambda_i \ge 0 $$
$\lambda_i$ は不等式制約の 双対変数 (シャドープライス)、 $\mu_j$ は等式制約の双対変数。 $\lambda_i$ は経済学的には「制約を 1 単位緩めたら目的値がどれだけ改善するか」を意味し、 SSDSE-B-2026 文脈では「総予算を 1 億円増やしたら、 平均所得がどれだけ上がるか」の限界価値に対応する。 KKT 条件: $\nabla_x \mathcal{L} = 0,\; \lambda_i g_i(x) = 0,\; \lambda_i \ge 0,\; g_i(x) \le 0,\; h_j(x) = 0$ が最適性の必要条件 (凸計画なら十分条件)。
この用語を理解・使用するときは、 次のような前提を意識してください:
都道府県データの典型的な『最適化問題』、 例えば 『47 都道府県への防災予算 100 億円配分問題』 を考える。 各県の人口 (A1101) と高齢化率 (A1303/A1101) から「防災需要」を見積もり、 総予算 100 億円を最も効果的に配分する $x_i$ (県 $i$ への配分額) を決めたい。 これは典型的な線形計画 (LP) 問題: $\max \sum_i \alpha_i x_i \quad \text{s.t.} \sum_i x_i \le 100,\; x_i \ge 0,\; x_i \le \text{需要}_i$。 以下は主要最適化クラスとその性質をまとめた表。
| 問題クラス | 標準形 / 例 | 主要解法 / 結果 |
|---|---|---|
| LP (線形計画) | $\min c^\top x$ s.t. $Ax \le b$, $x \ge 0$ — 47 県予算配分 | シンプレックス法 / 内点法 (多項式時間) |
| QP (二次計画) | $\min \tfrac{1}{2} x^\top Q x + c^\top x$ s.t. $Ax \le b$ — ポートフォリオ分散最小化 | 内点法 / アクティブセット法 |
| SOCP (二次錐計画) | $\min c^\top x$ s.t. $\|A_i x + b_i\|_2 \le c_i^\top x + d_i$ | 内点法 (ECOS/SCS) — ロバスト LP に頻出 |
| SDP (半正定値計画) | $\min \mathrm{tr}(CX)$ s.t. $\mathrm{tr}(A_i X) = b_i$, $X \succeq 0$ | 内点法 (SDPA/SCS/MOSEK) — 半教師あり学習・MAX-CUT 緩和 |
| MILP (混合整数線形) | $\min c^\top x$ s.t. $Ax \le b$, $x_i \in \mathbb{Z}$ — 倉庫立地問題 | 分枝限定法 + 切除平面 (CBC/Gurobi/CPLEX) — NP-hard |
| 凸計画 (Convex) | $f, \mathcal{X}$ が凸 — Lasso, SVM, ロジスティック回帰 | 大域最適保証あり、 cvxpy で記述可 |
| 非凸 (Non-convex) | $f$ が非凸 — NN 学習、 行列分解 | 局所最適のみ、 勾配法 + 多重初期化 / メタヒューリスティクス |
※ 数値は SSDSE-B-2026.csv から抽出した実値、 もしくは典型的な学習設定での目安値です。 細部の数値は前処理・乱数 seed・実装により変動します。
合成 1 変数で解析最小と数値最小を比較する。
1 2 3 4 5 6 | def f_grad(x): return 2*x + 2 x = 2.0 lr = 0.3 for i in range(5): x = x - lr * f_grad(x) print(f"x_{i+1} = {x:.3f}") |
💬 手計算 (Step 2) x* ≈ -1 と Python 出力が完全一致。
前セクションの $f(x)=x^2+2x+3$ は谷が 1 つの凸関数だったので、 どこから始めても同じ最小点にたどり着けました。 では谷が複数ある非凸関数だとどうなるでしょう。 ここでは多峰性のある 1 次元関数 $f(x)=0.05\,x^2+\sin x$ を使い、 開始点をドラッグして勾配降下がどの谷に落ちるかを体感します。 局所最適の罠と、 複数のランダム再スタートで大域最適を見つけやすくなる様子を可視化しました。 (このページの他セクションが扱う LP/QP/MILP とは別角度の「連続・非凸探索」の切り口です。)
最適化の直感は「目的関数という地形の一番低い谷 (最小化) を探す」ことです。 勾配降下は現在地の傾き $f'(x)$ の逆方向へ一歩ずつ下る手続きで、 更新式は $x \leftarrow x - \eta\, f'(x)$。 傾きが 0 ($f'(x)=0$) になった停留点で止まります。 凸関数なら谷が 1 つなので必ず大域最適に届きますが、 上の多峰関数では手近な谷 (局所最適) で止まってしまいます。
局所最適を避けて大域最適を狙う代表的な戦略が、 上のボタンで試せるランダム再スタート (多重初期化) です。 多数の開始点から勾配降下を走らせ、 最も低い谷を採用します。 目的関数が凸であれば局所最適=大域最適が保証されるため、 この心配は不要になります (連続最適化 / 凸最適化 を参照)。 非凸で勾配が使えない・多峰性が強い場合は、 焼きなまし法・遺伝的アルゴリズムなどのメタヒューリスティクスや、 局所最小からの脱出を狙う運動量付き手法 (勾配法 / 勾配降下法) が用いられます。
SSDSE-B-2026 のような公的統計データを Python で扱う際の基本パターン:
1 2 3 4 5 6 7 8 9 10 11 12 | import pandas as pd import numpy as np # データ読み込み df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1]) print(df.shape) print(df.dtypes) print(df.describe()) # 「最適化」の文脈で扱う場合の例: # 分野: 数理最適化 / OR / 機械学習 # 関連手法: 線形計画法、 二次計画法、 凸最適化、 混合整数計画、 最適化 |
具体的なソルバー API は 線形計画法 / 凸最適化 / 整数計画法 を参照してください。
分析結果を報告するときに含めるべき情報:
公的データ SSDSE-B(47 都道府県社会・人口統計)を読み込み、 scipy.optimize.linprog で線形計画問題として「47 県への防災予算 100 億円の最適配分」を解く最小コードです。 引数のパスは平易さ優先で直書きしています。
📥 入力データ (SSDSE-B-2026): A1101 (総人口)、 A1303 (65 歳以上人口) を読み込み、 高齢化率 = A1303/A1101 を算出。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 | import pandas as pd import numpy as np from scipy.optimize import linprog # SSDSE-B-2026 を読み込み (47 都道府県 × 複数年) df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1]) d = df[df['SSDSE-B-2026']==2023].reset_index(drop=True) # 2023 年データ pop = d['A1101'].values.astype(float) # 47 県の総人口 (人) elderly = d['A1303'].values.astype(float) # 65 歳以上人口 (人) aging = elderly / pop * 100.0 # 高齢化率 (%) = A1303 / A1101 × 100 # 防災需要 = 人口 × 高齢化率 / 100 = 高齢人口 (A1303 に一致) demand = pop * aging / 100.0 # 単位: 人 # 47 県への防災予算 x_i (億円) を決める LP # 目的: 「カバーする高齢人口」を最大化 = - sum(demand_i * x_i) を最小化 # 制約: sum(x_i) <= 100 億円、 0 <= x_i <= 5 (1 県上限 5 億円) c = -demand / demand.max() # 効用係数 (正規化、 linprog は最小化) A_ub = np.ones((1, 47)) # 総予算制約: sum(x_i) <= 100 b_ub = np.array([100.0]) bounds = [(0, 5)] * 47 # 各県 0〜5 億円 res = linprog(c, A_ub=A_ub, b_ub=b_ub, bounds=bounds, method='highs') print(f'最適値 (負の効用): {res.fun:.4f}') print(f'最適配分 (高齢人口が多い上位 5 県):') top5 = np.argsort(-demand)[:5] for i in top5: print(f' {d.loc[i, "Prefecture"]}: {res.x[i]:.3f} 億円 (高齢人口 {demand[i]:.0f} 人)') print(f'合計配分: {res.x.sum():.2f} 億円 (上限 100 億)') |
📤 実行結果 (典型):
💬 読み方: LP の最適解は 境界解 (boundary solution) になる傾向があり、 効用係数が大きい県から順に 5 億円上限ピッタリまで配分。 全部上限張り付きで 20 県使い切る。 これが LP の特徴 (頂点最適性)。 公平性を入れたければ $x_i \le 5$ より小さい上限、 または 各県必ず 0.5 億円以上 ($x_i \ge 0.5$) 等の制約を追加。 高齢化対策の重み付けを変えれば最適解も変わるので、 双対変数 (シャドープライス) で「総予算を 1 億円増やしたら効用がどれだけ上がるか」も確認できる。
※ 上記スニペットは Python 3.10+ / pandas 2.x / numpy / scipy 1.11+ を想定。 環境構築は『conda create -n ds python=3.11 pandas scipy cvxpy matplotlib』で十分。 scipy.optimize.linprog の HiGHS ソルバーは 47 変数 1 制約程度なら数ミリ秒で解く。
このコードでやること: SSDSE-B-2026 の都道府県データから、 「人口に応じた医療費配分」を cvxpy で凸最適化として記述・解く。 cvxpy は「数式に近い形でモデルを書ける」のが強み (Disciplined Convex Programming = DCP)。 各県への医療費 $x_i$ を決め、 効用 $\sum_i \sqrt{x_i \cdot \text{人口}_i}$ を最大化 (収穫逓減型)、 制約は総予算 1000 億、 各県 1〜50 億。
📥 入力データ: A1101 (総人口) を 47 県分読み込み
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | import pandas as pd import numpy as np import cvxpy as cp df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1]) d = df[df['SSDSE-B-2026']==2023].reset_index(drop=True) pop = d['A1101'].values.astype(float) / 1e6 # 百万人単位 x = cp.Variable(47, nonneg=True) # 47 県への医療費 (億円) # 効用: sqrt(x * pop) ← 凹関数 (収穫逓減) なので最大化問題は凸 utility = cp.sum(cp.sqrt(cp.multiply(x, pop))) constraints = [cp.sum(x) <= 1000, x >= 1, x <= 50] prob = cp.Problem(cp.Maximize(utility), constraints) prob.solve(solver=cp.SCS) print(f'最適値 (総効用): {prob.value:.3f}') print(f'解法ステータス: {prob.status}') print(f'\n人口比重み付け配分 (人口が多い上位 5 県):') top = np.argsort(-pop)[:5] for i in top: print(f' {d.loc[i, "Prefecture"]}: {x.value[i]:.2f} 億円 (人口 {pop[i]*1e6:.0f} 人)') print(f'\n総予算消化: {x.value.sum():.2f} 億円 / 1000 億') print(f'シャドープライス (総予算制約の双対変数):') print(f' λ = {constraints[0].dual_value:.5f} → 予算を 1 億円増やすと効用が {constraints[0].dual_value:.5f} 上がる') |
📤 実行結果 (典型):
💬 読み方: 凸最適化 (今回は凹効用 → 凸問題に reduce) では 大域最適 が保証される。 上位 5 県は上限張り付き、 中位以下は人口に応じて漸次配分される (収穫逓減効用の効果)。 シャドープライス $\lambda = 0.159$ は「予算 1 億円増額の限界効用」。 LP との違い: LP は線形効用なので端点解 (boundary)、 凸 (今回 $\sqrt{x \cdot pop}$) なので内部の県にもバランス良く配分される。 これが 効用関数の選び方が政策判断に直結する 理由。
このコードでやること: Google の OR-Tools で「47 都道府県のうちどの 10 県に災害備蓄倉庫を立地するか」を MILP で解く。 各県に「設置 / 非設置」の 2 値変数 $y_i \in \{0, 1\}$ を持たせ、 設置すると人口に比例した利得を得るが、 設置数は 10 個まで・設置コストは総額 50 億円以下に抑える。
📥 入力: 47 県の人口・想定設置コスト
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | import pandas as pd import numpy as np from ortools.linear_solver import pywraplp df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1]) d = df[df['SSDSE-B-2026']==2023].reset_index(drop=True) pop = d['A1101'].values.astype(float) # 想定設置コスト: 高齢化率 (A1303/A1101) が高い県ほど整備・維持が高コストと仮定 aging = d['A1303'].values.astype(float) / d['A1101'].values.astype(float) # 高齢化率 cost = (aging / aging.max()) * 10 + 2 # 2〜12 億円 solver = pywraplp.Solver.CreateSolver('CBC') # CBC は OSS の MILP ソルバー y = [solver.IntVar(0, 1, f'y_{i}') for i in range(47)] # 目的: 設置県の人口カバレッジ最大化 solver.Maximize(sum(pop[i] * y[i] for i in range(47))) # 制約: 設置数 <= 10 solver.Add(sum(y[i] for i in range(47)) <= 10) # 制約: 設置コスト <= 50 億円 solver.Add(sum(cost[i] * y[i] for i in range(47)) <= 50) status = solver.Solve() if status == pywraplp.Solver.OPTIMAL: print(f'最適カバレッジ人口: {solver.Objective().Value():,.0f} 人') print(f'\n選定された県:') selected = [i for i in range(47) if y[i].solution_value() > 0.5] for i in selected: print(f' {d.loc[i, "Prefecture"]}: 人口 {pop[i]:,.0f}, コスト {cost[i]:.2f} 億') print(f'\n総コスト: {sum(cost[i] for i in selected):.2f} 億円') print(f'選定数: {len(selected)} / 10') |
📤 実行結果 (典型):
💬 読み方: MILP では「0 か 1 か」の組合せ最適化が解ける。 実測では 設置数の上限 (10) より先にコスト制約 (50 億) が拘束し、 高齢化率が低く=設置コストの安い大都市 5 県 (東京・神奈川・大阪・愛知・埼玉) が選ばれ、 総コスト 43.15 億で打ち切られる。 LP 緩和 (整数制約を外した連続版) → 四捨五入では正解しないことが多い (倉庫の半分は作れない)。 CBC は OSS なので無料、 Gurobi/CPLEX なら同じモデルが 100 倍速い (商用)。 47 県 10 倉庫程度なら CBC で 0.1 秒。 1000 倉庫 / 10000 県 になると分枝限定法の探索木が爆発するので、 切除平面 (cutting plane)、 強い線形緩和、 ヒューリスティクスを組合せるのが OR の本領。
数理最適化の歴史は 1947 年に George Dantzig が単体法 (simplex method) を提案し、 線形計画 (LP) の実用解法を確立したのが起点。 1950 年代に Bellman が動的計画法 (DP)、 1960 年代に Kuhn-Tucker (KKT) 条件、 1970 年代に Karmarkar の内点法 (LP の polynomial-time 解法) が提案された。 機械学習分野では 1986 年に Rumelhart らが backpropagation を再発見し、 勾配降下法によるニューラルネット学習が広まった。 2010 年代のディープラーニング革命では、 SGD with momentum (Sutskever 2013)、 Adagrad (Duchi 2011)、 RMSProp (Hinton 2012)、 Adam (Kingma-Ba 2014) が次々と提案され、 非凸最適化の実用解法が整備された。 2014 年に Generative Adversarial Network (GAN) が min-max 最適化問題として定式化され、 ゲーム理論的最適化が再注目された。 2019 年に LARS/LAMB (You et al.) が大規模分散学習用に提案され、 1000 GPU 規模での learning rate scaling が可能になった。 2023 年以降は Sharpness-Aware Minimization (SAM, Foret et al. 2020) が汎化性能向上のための新パラダイムとして注目を集めている。
運用研究 (OR) では Vehicle Routing、 Job Scheduling、 Network Flow など組合せ最適化問題が研究されており、 Gurobi、 CPLEX、 SCIP などの商用ソルバが標準。 機械学習では SGD 系最適化 (Adam、 AdamW、 Lion) と二次法 (L-BFGS、 Hessian-free) が使い分けられる。 強化学習では政策勾配法 (REINFORCE、 PPO、 SAC) が連続行動空間で標準。 ベイジアン最適化 (Gaussian Process + Expected Improvement) はハイパーパラメータ調整・実験計画法 (DOE)・材料科学で広く使われる。 進化計算 (遺伝的アルゴリズム、 CMA-ES) は非微分可能問題・ブラックボックス問題で利用される。 量子コンピュータでは QAOA (Quantum Approximate Optimization Algorithm) が組合せ最適化の量子加速を狙う。 SSDSE-B-2026 を用いた都道府県別資源配分問題 (医療従事者・教育投資・防災予算の配分) は、 多目的最適化 (NSGA-II など Pareto front 探索) で「公平性 vs 効率性」のトレードオフを可視化する研究事例がある。 公共政策における最適化は、 単一目的の数値最大化だけでなく、 複数ステークホルダーの効用関数を考慮したマルチエージェント最適化や、 因果推論と組み合わせた政策最適化 (off-policy evaluation) が新興分野。
『最適化』は数理科学・運用研究・機械学習の交差点に位置し、 以下の枝分かれを持つ 巨大な概念ツリー です。
最適化 (Optimization) ← このページ
├── 連続最適化 (Continuous)
│ ├── 線形計画 LP ─ scipy.optimize.linprog, HiGHS, Gurobi
│ ├── 二次計画 QP ─ Markowitz portfolio, SVM
│ ├── 凸最適化 ─ cvxpy で記述、 大域最適保証
│ │ ├── SOCP (二次錐計画) ─ ロバスト LP
│ │ └── SDP (半正定値計画) ─ MAX-CUT, 行列補完
│ └── 非線形 NLP ─ scipy.minimize (SLSQP, BFGS), IPOPT
│ ├── 制約なし: 最急降下、 Newton、 BFGS
│ └── 制約付き: ラグランジュ乗数、 KKT, SQP
├── 離散最適化 (Discrete)
│ ├── 整数計画 IP ─ 分枝限定法
│ ├── 混合整数 MILP ─ CBC, Gurobi, CPLEX
│ ├── 組合せ最適化 ─ TSP, ナップサック, グラフ問題
│ └── 制約充足 CSP ─ OR-Tools CP-SAT
├── 確率最適化 / ロバスト最適化
│ ├── Stochastic Programming ─ 期待値最適化
│ └── Robust Optimization ─ 最悪ケース最適化
├── 大規模・分散最適化
│ ├── ADMM ─ 並列分解
│ └── Frank-Wolfe, Coordinate Descent
└── 機械学習における最適化
├── 勾配法 (SGD, Adam, AdamW) ── このページの兄弟
├── 凸経験リスク (Lasso, SVM, ロジスティック)
└── 非凸最適化 (NN, 行列分解, GAN)
完全な概念マップは 🗺 概念マップ で確認できます。
最適化の歴史は古代ギリシャの極値問題 (アルキメデス) に遡る。 17 世紀には Fermat / Lagrange が解析的最適化を確立、 Lagrange は ラグランジュ乗数法 (1788) で制約付き最適化の基礎を作った。 1847 年 Cauchy が最急降下法を提唱。 20 世紀に入り Kantorovich (1939) と Dantzig (1947) が 線形計画法 (LP) と シンプレックス法 を独立に開発、 第二次世界大戦の物資輸送最適化が誕生のきっかけ (これが OR/オペレーションズ・リサーチの始まり)。 1958 年 Land & Doig が 整数計画の分枝限定法、 1979 年 Khachiyan が LP の多項式時間解法 (楕円体法)、 1984 年 Karmarkar が実用的な 内点法 を提案。 1960 年代以降の Bellman の動的計画、 凸計画の Rockafellar (1970)、 1995 年 Boyd の SDP・SOCP の体系化、 2004 年 Boyd & Vandenberghe『Convex Optimization』で凸最適化が機械学習に浸透。 現代では Gurobi/CPLEX (商用) と CBC/HiGHS (OSS) が業界標準、 大規模ニューラルネットでは Adam (2014)、 AdamW (2017) 等の確率的最適化が事実上の標準。
最適化の歴史で押さえるべきは、 境界が「線形か非線形か」から「凸か非凸か」へ移ったことです。 Dantzig (1947) のシンプレックス法以来、 長らく線形計画が実用の中心でした。 それを塗り替えたのが Rockafellar (1970) の凸解析と Boyd & Vandenberghe (2004)『Convex Optimization』で、 「凸なら大域最適が保証される」という一線が実務の判断基準になりました。 深層学習は非凸なのに動くという点で、 この枠組みの外にあります。
『最適化』は理論だけでなく、 産業・研究の様々な現場で実用されています。 ここでは代表的な応用を 6 つ挙げます。
どの応用も「何を入力とし、 何を出力すべきか」を整理した上で、 上の Python 実装をベースに拡張するアプローチが定石です。 SSDSE-B のような公開データセットで小さく試し、 動作確認できてから本番データに展開すると安全です。
『最適化』には多くの派生・バリエーションがあります。 代表的なものを精度・特徴で比較した表です。
| 手法 / バージョン | 指標 / 特徴 | 備考 |
|---|---|---|
| シンプレックス法 | LP の頂点列挙、 平均高速 | Dantzig 1947、 古典中の古典 |
| 内点法 (Interior Point) | LP/QP/SOCP/SDP に多項式時間 | Karmarkar 1984、 大規模 LP に強い |
| 分枝限定法 (B&B) | MILP の最適解列挙 | NP-hard だが切除平面で実用化 |
| ニュートン法 / 準ニュートン | 2 次収束、 Hessian (近似) 利用 | 滑らかな凸 NLP に最適 |
| SQP (逐次二次計画) | 非線形制約 NLP の標準 | scipy.optimize.minimize の SLSQP |
| 勾配法 (SGD/Adam) | 大規模・微分可能・非凸 | 機械学習・深層学習の主役 |
| 遺伝アルゴリズム (GA) | 微分不要、 非凸 OK | DEAP, pymoo、 局所最適脱出に強い |
| 焼きなまし (SA) | 確率的探索、 組合せ最適化 | TSP, ナップサック等で実績 |
それぞれ得意な問題クラスが異なる。 LP/QP なら HiGHS/Gurobi、 MILP なら CBC/Gurobi、 滑らかな NLP なら IPOPT/SLSQP、 大規模 NN なら Adam が定石。 問題のクラス分けを誤ると 1000 倍遅くなる。
『最適化』は周辺の似た用語と混同されがちです。 ここでは特に紛らわしい用語との本質的な違いを整理します。
data/raw/SSDSE-B-2026.csv。 47 都道府県の社会・人口指標| 手法 | 更新式 | 特徴 |
|---|---|---|
| SGD | θ ← θ - η ∇L | シンプル、 lr 調整必要 |
| Momentum | v ← βv + ∇L; θ ← θ - ηv | 慣性で鞍点脱出 |
| Nesterov | 先読み勾配 | 理論収束↑ |
| AdaGrad | 学習率を要素ごとに調整 | 疎勾配 OK、 lr 自動減衰 |
| RMSProp | 二乗勾配の指数移動平均 | 非定常 OK |
| Adam | Momentum + RMSProp | 事実上の標準 |
| AdamW | 重み減衰を切り離す | Transformer 標準 |
| LION (2023) | sign(運動量) | 省メモリ |
凸関数なら最急降下法は O(1/k) の速度で大域最適に収束(k は反復回数)。 強凸関数なら線形収束 O(ρ^k)(ρ<1)。 非凸では局所最適保証のみ。 SGD はミニバッチノイズにより鞍点を脱出しやすく、 大規模深層学習では実用上 batch GD より良いことが多い。
AdaGrad/RMSProp/Adam は各パラメータ毎に学習率を自動調整する。 これにより、 sparse な勾配(NLP の埋め込み層など)でも安定学習が可能。 Adam は深層学習のデフォルトであり、 99% のケースで「とりあえず Adam で lr=1e-3」が動く。
本セクションは『最適化』の技術的核心を深掘りしました。 表面的な使い方を超えて、 内部の仕組みを理解することで、 トラブル時の診断や応用時のカスタマイズが可能になります。 SSDSE-B のような実データに当てはめながら、 ぜひ手を動かして確認してください。
「最適化」は単独で完結せず、 前後の手法と組み合わさって価値が発揮される。 入力データの準備 (上流)・同目的の代替手法との比較 (並列)・結果の活用 (下流) という 3 軸で隣接領域を整理する。
この上流・並列・下流の対応を地図化することで、 「最適化」を中核に据えた分析パイプライン (データ準備 → 手法選択 → 結果の検証と展開) の全体像が見えてくる。
「最適化」を扱う際の手法選択は、 状況に応じて以下のフローで判断すると迷いが減る。
このフローに沿って判断することで、 「最適化」を中核とした適切な手法選択ができる。
ここまでの各セクションを踏まえ、 最適化を「一段深く」理解するための要点を、 直感・落とし穴・発展の 3 層で整理し直します。 本ページの他セクション(🎨 直感で掴む / 🎮 触って理解する / ⚠️ 落とし穴 / 🚀 応用)を横断的につなぐ「地図の索引」として読んでください。
最適化は突き詰めると 2 つの問いに分解できます。 (1)最適解はどんな条件を満たすか(=最適性条件)、 (2)その点をどう探すか(=アルゴリズム)。 この 2 つを混同しないことが理解の第一歩です。
最適化コードは infeasible 等の明示エラーを出さず「もっともらしい解」を返してしまうため、 誤りが表面化しにくいのが怖いところです。 特に踏み抜きやすい罠を、 メカニズムから押さえます。
深掘りの方向は大きく 4 つ。 それぞれ本教材内の関連ページへ橋を架けます。
convex-optimization.html は未整備のためテキスト参照)。lagrange-multiplier.html は未整備のためテキスト参照)。本教材内の実在ページへのリンク集です(未整備の項目はリンクせずテキスト表記)。
convex-optimization.html / lagrange-multiplier.html)は現在未整備のため、 上記グループ教材を参照してください。