M17
一维加权 Gagliardo–Nirenberg 不等式最优常数的数值变分测定：以解析已知情形为标尺
1 · 研究问题
用数值变分法（梯度流/约束优化）求解一维 Gagliardo–Nirenberg（GN）不等式的最优常数与极值函数，在解析已知情形（Nagy 1941 一维族）校准后，对幂权 |x|^a 加权族，最优常数随权指数 a 的曲线是什么形状？极值函数的对称性在 a 变化中是否保持？
2 · 研究背景与空白
GN 插值不等式 ‖u‖_q ≤ C‖u'‖₂^θ‖u‖_p^(1−θ) 的最优常数在一维由 Nagy 1941 完整解决，极值函数为显式（如 sech 型）；高维族由 Del Pino–Dolbeault 2002 给出一支显式解。加权情形:Caffarelli–Kohn–Nirenberg（CKN）不等式族的最优常数只在部分参数区已知，对称性破缺边界（Felli–Schneider 曲线）是已发表的深刻结果；近期有最优输运方法处理带权 GN（2022，ResearchGate 所见《Weighted Gagliardo-Nirenberg inequalities via Optimal Transport》，权相等时常数显式）。
已有工作给出:无权情形显式常数、加权情形部分区域的显式常数与对称性理论。空白在参数维度：一维加权族中"常数曲线的完整数值形状与对称性行为"在已知解析点之间缺少公开数值表（检索未见系统表列；未找到不等于不存在）。变分问题一维、极值函数光滑衰减，数值上是笔记本级任务；解析锚点密集，校验极硬。
3 · 可检验假设
- H1 数值管线在 ≥ 3 个解析已知情形（Nagy 一维族取三组 (p,q)）复现最优常数至相对误差 ≤ 0.1%，极值函数与显式 sech 型解的 L² 距离 ≤ 10⁻³。
- H2 加权族中最优常数 C(a) 随 a 连续单调（方向可测），且在扫描区间内极值函数保持偶对称——或在某个 a* 处出现对称破缺（数值判据:非对称初值收敛到能量更低的非对称解，能量差 > 10 倍离散误差）。两种结果都构成结论。
4 · 量化验收标准
- 方法学校验（硬门槛）：H1 即硬门槛——三个解析情形复现 ≤ 0.1%；不过关则加权族全部结果无效。
- 离散化收敛检查:网格 N ∈ {10³, 4×10³, 1.6×10⁴} 三档 Richardson 外推，常数外推余项 ≤ 10⁻⁴；区间截断 [−L, L] 两档验证边界误差 ≤ 10⁻⁵。
- 每个 a 值从 ≥ 5 个独立初值（含刻意非对称者）启动优化，报告收敛到的极小值个数与能量差。
- 权指数扫描 ≥ 20 点，疑似破缺点二分加密至 Δa ≤ 0.01。
- 独立交叉校验:变分法结果与打靶法解 Euler–Lagrange ODE（scipy.integrate + 二分打靶）在 ≥ 5 个 a 点互验，常数差 ≤ 0.05%。
- 代码开源、一键重跑。
5 · 数据与工具
| 用途 | 来源 / 工具 |
|---|---|
| 变分优化 | 自建:归一化梯度流（离散 L² 梯度 + 投影）或 scipy.optimize.minimize(L-BFGS-B)；离散化用有限差分（numpy） |
| ODE 打靶交叉校验 | scipy.integrate.solve_ivp + 自建二分打靶 |
| 解析锚点 | Nagy 1941 一维显式常数（Frank 等综述有现代转写，人工转录双核对）；Del Pino–Dolbeault 2002，仅用于校验，不计入贡献 |
| 符号验证 | sympy 核验 Euler–Lagrange 方程推导与显式解代入 |
| 对称性理论对照 | Felli–Schneider 曲线文献（CKN 情形；一维加权对应关系需核实——若一维退化无破缺，H2 按单调性结论走） |
| 算力 | 单次优化 N=1.6×10⁴ 网格 ~ 秒–分钟级；全项目 CPU < 10 小时 |
6 · 方法路径
- 推导离散能量泛函与梯度（sympy 核验），实现梯度流。
- 无权解析情形三点校准（第 1 条验收）。
- 完成收敛性检查与截断误差控制，冻结数值参数。
- 加权族 a 扫描，多初值启动，记录常数曲线与极值函数形状。
- 疑似破缺区二分加密；对称/非对称解能量差判定（判据写死）。
- 打靶法独立互验 5 个点。
- 与已发表加权情形（权相等特例）的显式常数对照，成文。
7 · 新颖性边界
本课题不声称证明任何不等式或新常数的解析值；Nagy 1941、Del Pino–Dolbeault 2002、CKN 对称破缺（Felli–Schneider、Dolbeault–Esteban–Loss）及最优输运加权 GN（2022）均为已发表理论。已有数值工作多为验证性个例。本项目的主结论是一维加权族最优常数曲线与对称性行为的系统数值表，解析已知点作为标尺嵌入表中——属"参数维度扩展 + 可复现性贡献"。风险坦白:一维情形可能理论上排除对称破缺（需第 8 周核实文献），届时 H2 自动退为常数曲线的单调性与凸性测量，仍是有效主结论。历届丘奖不等式类获奖（2024 bronze 等周问题）为纯理论，无重合。
8 · 决策门槛（go / no-go）
- 第 4 周末:无权情形 0.1% 校准通过；不通过先查离散化（分析类课题此步失败率高，预留 2 周缓冲）。
- 第 8 周末:完成"一维加权是否存在破缺"的文献核实。若被理论排除，H2 改述（见上），框架不变。
- 第 16 周末:加权扫描中若梯度流在大 |a| 处不收敛（权奇性致刚性），降级路径 A：缩小 a 区间并改用打靶法为主、变分为辅——常数曲线主结论保留。
- 第 30 周末:若打靶与变分互验差 > 0.05% 无法消除，报告两法差异本身并以外推余项为误差estimate，主结论降精度保留。
- 适合对象:分析口味最重的一题，适合喜欢证明与推导、同时能写代码的学生；数学写作门槛高于 M11–M16。