第 1 周作业:一维流的精确解

Views: --

源作业为教材练习 2.1.4,要求从隐式积分解出 x˙=sinx\dot x=\sin x,并处理指定初值与一般初值。

题目

已知

x˙=sinx\dot x=\sin x

的隐式解可写为

t=lncscx0+cotx0cscx+cotx,t=\ln\left| \frac{\csc x_0+\cot x_0}{\csc x+\cot x} \right|,

其中 x(0)=x0x(0)=x_0

  1. x0=π/4x_0=\pi/4 时,证明
x(t)=2arctanet1+2,x(t)=2\arctan\frac{e^t}{1+\sqrt2},

并推出 tt\to\inftyx(t)πx(t)\to\pi; 2. 对一般初值 x0x_0 写出解。

解答

利用恒等式

cscx+cotx=1+cosxsinx=cotx2,\csc x+\cot x =\frac{1+\cos x}{\sin x} =\cot\frac x2,

隐式解化为

cotx2=etcotx02.\cot\frac x2 =e^{-t}\cot\frac{x_0}{2}.

等价地

tanx2=ettanx02.\tan\frac x2 =e^t\tan\frac{x_0}{2}.

因此在不跨越相邻不动点的连续分支上,

x(t)=2arctan(ettanx02).x(t)=2\arctan\left(e^t\tan\frac{x_0}{2}\right).

x0=π/4x_0=\pi/4

tanπ8=21=11+2,\tan\frac\pi8=\sqrt2-1=\frac1{1+\sqrt2},

于是得到题目中的表达式。tt\to\infty 时反正切的自变量趋于 ++\infty,所以 x(t)πx(t)\to\pi

分支与特殊初值

x=kπx=k\pi 都是不动点,不能套入除以 sinx\sin x 后的公式;若 x0=kπx_0=k\pi,解就是常数 x(t)=kπx(t)=k\pi

一般初值应先判断它位于哪个区间 (kπ,(k+1)π)(k\pi,(k+1)\pi),再选择反正切分支。相线表明:

  • (2kπ,(2k+1)π)(2k\pi,(2k+1)\pi)sinx>0\sin x>0,轨道向右趋于 (2k+1)π(2k+1)\pi
  • ((2k1)π,2kπ)((2k-1)\pi,2k\pi)sinx<0\sin x<0,轨道向左趋于 (2k1)π(2k-1)\pi

这也说明偶数倍 π\pi 不稳定,奇数倍 π\pi 稳定。

评论