关于变量替换x=2p后的常微分方程变换疑问
嗨,我来帮你把这个变换的细节拆解开——你困惑的点其实是没注意到高阶导数的链式法则,还有最后一步的方程化简,不是直接硬替换变量哦!
首先明确函数关系:原来$v$是$p$的函数,现在做替换$x=2p$,也就是$p=x/2$,所以$v$就变成了关于$x$的复合函数$v(p(x))$,求导的时候必须用链式法则,不能直接把导数符号里的变量换掉,这是关键!
我们一步步来算:
一阶导数的变换:
根据链式法则,$\frac{dv}{dp} = \frac{dv}{dx} \cdot \frac{dx}{dp}$,因为$x=2p$,所以$\frac{dx}{dp}=2$,代入后得到:
$$\frac{dv}{dp} = 2 \cdot \frac{dv}{dx}$$二阶导数的变换:
二阶导数是对一阶导数再求关于$p$的导数,也就是$\frac{d2v}{dp2} = \frac{d}{dp}\left(\frac{dv}{dp}\right)$,把刚才的一阶导数结果代入:
$$\frac{d2v}{dp2} = \frac{d}{dp}\left(2 \cdot \frac{dv}{dx}\right) = 2 \cdot \frac{d}{dp}\left(\frac{dv}{dx}\right)$$
这里又要用到链式法则:$\frac{d}{dp}\left(\frac{dv}{dx}\right) = \frac{d}{dx}\left(\frac{dv}{dx}\right) \cdot \frac{dx}{dp} = \frac{d2v}{dx2} \cdot 2$,所以最终二阶导数是:
$$\frac{d2v}{dp2} = 2 \cdot 2 \cdot \frac{d2v}{dx2} = 4 \cdot \frac{d2v}{dx2}$$代入原方程并化简:
现在把$p=x/2$、$\frac{dv}{dp}=2\frac{dv}{dx}$、$\frac{d2v}{dp2}=4\frac{d2v}{dx2}$全部代入原ODE:
$$\frac{x}{2} \cdot 4\frac{d2v}{dx2} + 2\left(l+1 - \frac{x}{2}\right) \cdot 2\frac{dv}{dx} + \left(k - 2(l+1)\right)v = 0$$
现在逐项计算:- 第一项:$\frac{x}{2} \times 4 = 2x$,所以第一项变为$2x \cdot \frac{d2v}{dx2}$
- 第二项:$2 \times \left(l+1 - \frac{x}{2}\right) \times 2 = 4(l+1) - 2x = 2\left(2(l+1) - x\right)$,所以第二项变为$2\left(2(l+1) - x\right) \cdot \frac{dv}{dx}$
- 第三项:保持不变,还是$\left(k - 2(l+1)\right)v$
现在方程变成:
$$2x \cdot \frac{d2v}{dx2} + 2\left(2(l+1) - x\right) \cdot \frac{dv}{dx} + \left(k - 2(l+1)\right)v = 0$$
最后一步,我们把方程两边同时除以2,就得到了你看到的变换后的方程:
$$x\frac{d2v}{dx2}+\big(2(l+1)-x\big)\frac{dv}{dx}+\left(\frac{1}{2}k-(l+1)\right)v = 0$$
这样你就明白为什么系数会有看起来“翻倍/减半”的变化了——一是高阶导数的链式法则带来的倍数,二是最后方程两边除以2的化简步骤,这两步是你之前没注意到的关键点!
备注:内容来源于stack exchange,提问作者AtomProgrammer

