求解微分方程组的轨道以绘制相图

嘿，我来帮你把这个微分方程组的轨道求解和相图分析理清楚～首先你对**轨道（orbit）**的定义完全没问题，它就是解曲线$\gamma: t \mapsto (x(t), y(t))$在相平面上的点集$\gamma[I]$，也就是忽略时间参数后的轨迹，这正是相图分析里的核心概念。

你找的方向完全正确——利用$\frac{dy}{dx} = \frac{\dot{y}}{\dot{x}}$把原方程组转化为关于$x$和$y$的一阶微分方程，避开时间参数直接求轨线：
$$\frac{dy}{dx} = \frac{-xy}{y - 1}$$
这是个可分离变量的方程，我们一步步解它：

• 先分离变量，把含$y$的项放左边，含$x$的项放右边：
  $$\frac{y - 1}{y} dy = -x dx$$
• 拆分左边的分式，简化积分难度：
  $$\left(1 - \frac{1}{y}\right) dy = -x dx$$
• 两边同时积分：
  $$\int \left(1 - \frac{1}{y}\right) dy = -\int x dx$$
  计算后得到轨线的隐式方程：
  $$y - \ln|y| = -\frac{1}{2}x^2 + C$$
  这里的$C$是积分常数，由初值$(x_0, y_0)$确定：$C = y_0 - \ln|y_0| + \frac{1}{2}x_0^2$

除了上面的通解，还要注意几个特殊情况，这些是相图里的关键元素：

• 平衡点：令原方程组的导数都为0，即$\dot{x}=y-1=0$且$\dot{y}=-xy=0$，解得唯一平衡点$(0,1)$。我们可以用雅可比矩阵判断它的类型：
  原方程组的雅可比矩阵为：
  $$J(x,y)=\begin{pmatrix}0 & 1 \ -y & -x\end{pmatrix}$$
  代入平衡点$(0,1)$后得到：
  $$J(0,1)=\begin{pmatrix}0 & 1 \ -1 & 0\end{pmatrix}$$
  它的特征值是$\pm i$，说明这是一个中心型平衡点，周围的轨线是闭合的曲线，和我们得到的隐式方程趋势一致。
• $y=1$的特殊情况：当$y=1$时，$\dot{x}=0$意味着$x$保持常数，而$\dot{y}=-x$，所以：
  • 若$x=0$，就是刚才的平衡点$(0,1)$；
  • 若$x=a\neq0$，轨线是竖直线$x=a$，$y(t)=1 - a(t-t_0)$，也就是沿着竖直线上下移动。
• $y=0$的特殊情况：当$y=0$时，$\dot{x}=-1$，$\dot{y}=0$，轨线是水平线$y=0$，$x(t)=x_0 - (t-t_0)$，也就是沿着$y=0$向左移动的射线。

根据得到的隐式轨线方程$y - \ln|y| + \frac{1}{2}x^2 = C$，你可以给不同的$C$值绘制曲线：

• 当$C$接近$1 - \ln1 + 0 = 1$时，曲线靠近平衡点$(0,1)$，是闭合的小曲线；
• 当$C$增大时，曲线会向$y>0$的外侧扩展，$x$的取值范围也会变大；
• 对于$y<0$的区域，$\ln|y|=\ln(-y)$，轨线方程变为$y - \ln(-y) + \frac{1}{2}x^2 = C$，这里$y$为负，$y - \ln(-y)$是单调递增的，所以每个$C$对应一条横跨$x$轴左右的曲线。

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