含双函数的二阶微分方程组：求解可行性及方法问询

这个二阶耦合非线性微分方程组是可以求解的，但得看你要的是解析解还是数值解——两种情况的方法差异不小，我给你详细拆解：

首先：解的存在性肯定没问题 #

根据微分方程里的皮卡-林德洛夫定理，只要方程组右侧的函数满足利普希茨条件，给定合适的初值（比如指定某个点$x_0$处的$h(x_0)$、$h'(x_0)$、$f(x_0)$、$f'(x_0)$），局部范围内一定存在唯一解。而这里的$\sin(h(x))$和$\cos(h(x))$都是光滑、有界的函数，所以不仅局部有解，甚至可以把解延拓到整个实轴上，不存在“无解”的情况。

求解方法分两类 #

1. 解析解：仅特殊情况可行，通用解很难找 #

这是强耦合的非线性系统，一般情况下没有初等解析解——也就是没法用多项式、三角函数、指数函数这类常见初等函数组合出解的表达式。不过可以试试找一些特殊解：

• 假设$h(x)$是常数$h_0$，代入第一个方程得$0 = f(x)\sin(h_0)$：
  • 如果$\sin(h_0)≠0$，那$f(x)=0$，但代入第二个方程会得到$0=\cos(h_0)$，矛盾；
  • 如果$\sin(h_0)=0$（即$h_0=k\pi$，$k$为整数），第二个方程就变成$f''(x)=\cos(k\pi)=(-1)^{k$，积分两次可得$f(x)=(-1)}k·\frac{x^2}{2} + C_1x + C_2$（$C_1$、$C_2$是任意常数），这是一类简单的特殊解。
• 也可以尝试把二阶方程组转化为一阶方程组来寻找守恒量：
  令$u=h(x)$，$p=h'(x)$，$v=f(x)$，$q=f'(x)$，原方程组就变成4维一阶系统：

  u' = p
  p' = v·sin(u)
  v' = q
  q' = cos(u)
  

  不过尝试推导首次积分（守恒量）的过程会很繁琐，目前看来很难找到通用的解析表达式，除非给初值加特殊限制或者利用对称性。

2. 数值解：通用可行，方法成熟 #

如果不需要解析表达式，数值解是完全靠谱的选择，常用方法有：

• 龙格-库塔法（比如经典RK4）：把上面的4维一阶方程组作为输入，给定初值后，就能逐点迭代计算$h(x)$、$h'(x)$、$f(x)$、$f'(x)$的值，精度和效率都不错，是最常用的方法。
• 多步法（如亚当斯系列方法）：适合需要高精度或者长时间积分的场景，比如亚当斯-巴什福思（显式）、亚当斯-穆尔顿（隐式），能利用之前计算的点来提高精度。
• 有限差分法：把$x$轴离散成网格点，用差分近似二阶导数（比如$h''(x)≈\frac{h(x+\Delta x)-2h(x)+h(x-\Delta x)}{(\Delta x)^2}$），把原方程组转化为非线性代数方程组，再用牛顿法迭代求解，适合需要处理边界值问题的情况。

内容的提问来源于stack exchange，提问作者Mathy1002