⚡

边界值问题Skill boundary-value-problems

这个技能提供解决常微分方程和偏微分方程中边界值问题的策略和方法，包括问题分类、射击方法、有限差分法、配置法以及特征值问题求解。适用于数学建模、科学计算和工程分析等领域，关键词包括边界值问题、ODE、PDE、数值方法、射击法、有限差分、特征值问题。

预测建模 0 次安装 0 次浏览 更新于 3/14/2026

name: 边界值问题 description: “用于常微分方程和偏微分方程中边界值问题的解决策略” allowed-tools: [Bash, Read]

边界值问题

何时使用

在常微分方程和偏微分方程中处理边界值问题时使用此技能。

决策树

  1. 问题分类

    • 两点边界值问题：条件在 x=a 和 x=b？
    • Sturm-Liouville 问题：特征值问题？
    • 混合条件：Dirichlet、Neumann、Robin？

  2. 射击方法

    • 将边界值问题转换为初值问题
    • 猜测缺失的初始条件
    • 迭代以满足边界条件
    • scipy.integrate.solve_ivp + 根查找

  3. 有限差分方法

    • 离散化域：x_i = a + i*h
    • 用差分替换导数：y’’ ~ (y_{i+1} - 2y_i + y_{i-1})/h^2
    • 求解得到的线性系统
    • sympy_compute.py linsolve "tridiagonal_matrix" "boundary_vector"

  4. 配置法/边界值问题求解器

    • scipy.integrate.solve_bvp(ode, bc, x, y_init)
    • 提供初始网格和猜测
    • 检查残差以获取准确性

  5. 特征值问题

    • Sturm-Liouville 形式：-(p(x)y’)’ + q(x)y = lambda*w(x)*y
    • 如果 p, w > 0，特征值为实数
    • 特征函数在权重 w 下正交
    • sympy_compute.py eigenvalues "sturm_liouville_matrix"

工具命令

Scipy_Solve_Bvp

uv run python -c "from scipy.integrate import solve_bvp; import numpy as np; ode = lambda x, y: [y[1], -y[0]]; bc = lambda ya, yb: [ya[0], yb[0]-1]; x = np.linspace(0, np.pi, 10); y = np.zeros((2, 10)); sol = solve_bvp(ode, bc, x, y); print('Solution at pi/2:', sol.sol(np.pi/2)[0])"

Sympy_Linsolve

uv run python -m runtime.harness scripts/sympy_compute.py linsolve "tridiagonal_matrix" "boundary_vector"

Z3_Sturm_Liouville

uv run python -m runtime.harness scripts/z3_solve.py prove "eigenvalue_real"

关键技术

来自索引教材：

  • [Elementary Differential Equations and… (Z-Library)] 边界值问题和偏微分方程 (第6版，波士顿：学术出版社，2012年8月7日 21:05 c10 页数 88 第 676 页 青色黑色 2012年8月7日 21:05 c11 页数 1 第 677 页 青色黑色 第十一章 边界值问题和Sturm–Liouville理论 作为在第十章偏微分方程中分离变量的结果，我们反复遇到微分方程 X + λX = 0, 0 < x < L 边界条件 X (0) = 0, X (L) = 0。这个边界值问题是一大类在应用数学中重要问题的原型。
  • [Elementary Differential Equations and… (Z-Library)] 非齐次边界值问题 本节讨论如何解决常微分方程和偏微分方程中的非齐次边界值问题。我们的大部分注意力集中在微分方程本身是非齐次而边界条件是齐次的问题上。我们假设解可以展开为相关齐次问题的特征函数级数，然后我们确定这个级数中的系数，以便满足非齐次问题。
  • [Elementary Differential Equations and… (Z-Library)] 考虑边界条件 y, y 当 x → −1, −1 m = n。 2012年8月7日 21:05 c11 页数 46 第 722 页 青色黑色 第十一章 边界值问题 一般微分方程或边界条件。
  • [An Introduction to Numerical Analysis… (Z-Library)] 现代常微分方程数值方法 Wiley, 纽约。 DVERK用户指南：用于求解非刚性常微分方程的子程序。 Keller (1966), 数值方法分析。
  • [Elementary Differential Equations and… (Z-Library)] 用几句话描述解如何随时间演变。一种非反应性示踪剂以浓度 c0 连续引入到长度为 L 的填充有均匀颗粒介质的柱的上游端。假设柱中示踪剂浓度初始为零，建模此过程的边界值问题是 0 < x < L, t > 0, t > 0, 0 < x < L, 其中 c(x, t), v, 和 D 如问题27中所述。

认知工具参考

查看 .claude/skills/math-mode/SKILL.md 获取完整工具文档。