You need to enable JavaScript to run this app.
最新活动
大模型
产品
解决方案
定价
生态与合作
支持与服务
开发者
了解我们

几何问题微分方程构建:求第一象限满足固定面积条件的曲线

求解第一象限内满足条件的曲线

别慌!这种看起来直白的几何微分方程问题,卡壳太正常了——我当初第一次做的时候也绕了好几个弯。咱们一步步拆解,把几何条件转化为可解的微分方程就行。

1. 先建立几何关系

设第一象限内的曲线为 y = y(x)(x>0,y>0),取曲线上任意一点T(x, y)。我们需要先写出T点处的切线方程,再找到它与x轴的交点,最后用三角形面积恒为p的条件列方程。

  • 切线方程:曲线在T点的导数为 y’ = dy/dx,用点斜式写出切线方程(用大写X、Y表示切线上的动点,避免和T点坐标混淆):
    Y - y = y'(X - x)
    
  • 切线与x轴的交点:令Y=0,解出对应的X值,得到交点A的坐标:
    X = x - y/y'
    
  • 三角形面积条件:原点O(0,0)、T(x,y)、A(X,0)围成的三角形面积恒为p。用三角形面积公式(底×高/2),这里底是OA的长度|X|,高是T点的纵坐标y,所以:
    (1/2) * |X| * y = p
    
    因为曲线在第一象限,且要保证三角形存在,切线与x轴的交点必然在正半轴(X>0),所以可以去掉绝对值,代入X的表达式:
    (1/2) * (x - y/y') * y = p
    

2. 转化为微分方程并求解

把上面的面积等式整理成标准微分方程形式:

(x y - y²/y') = 2p

两边乘以 y’ 后整理:

(x y - 2p) y' = y²

这个方程直接解 dy/dx 有点麻烦,我们换个思路,把它看成 关于x的一阶线性微分方程(即把x看作y的函数,求 dx/dy):

dx/dy = (x y - 2p)/y² = x/y - 2p/y²

整理成一阶线性非齐次方程的标准形式:

dx/dy - (1/y)x = -2p/y²

用积分因子法求解

对于一阶线性方程 dx/dy + P(y)x = Q(y),积分因子 μ(y) 的计算公式是:

μ(y) = e^(∫P(y) dy)

这里 P(y) = -1/y,所以:

μ(y) = e^(∫-1/y dy) = e^(-ln y) = 1/y

把方程两边乘以积分因子 1/y,左边会变成导数的形式:

d/dy (x/y) = -2p/y³

两边对y积分:

x/y = ∫-2p/y³ dy + C

计算右边的积分:

∫-2p/y³ dy = p/y² + C

(C为积分常数)

最后整理得到通解:

x y = p + C y²

或者写成更直观的形式:

y(x - C y) = p

3. 验证解的正确性

我们拿通解来验证是否满足原面积条件:
x y = p + C y² 可得 x = p/y + C y,求导得 dx/dy = -p/y² + C,因此 dy/dx = 1/(C - p/y²) = y²/(C y² - p)

切线与x轴的交点X:

X = x - y/(dy/dx) = (p/y + C y) - y*(C y² - p)/y² = 2p/y

三角形面积为:

(1/2) * X * y = (1/2)*(2p/y)*y = p

完全符合题目要求!

4. 第一象限内的曲线形式

因为限定在第一象限(x>0,y>0),我们可以分析通解的几种情况:

  • 当C=0时,方程为 x y = p,这是我们熟悉的反比例函数,是最典型的解;
  • 当C≠0时,曲线是双曲线的一支:
    • 若C>0,曲线位于直线 x = C y 的上方(保证x - C y > 0);
    • 若C<0,令C=-k(k>0),方程变为 y(x + k y) = p,曲线在第一象限内是向下凸的形态。

这些曲线都满足“切线、x轴、OT围成的三角形面积恒为p”的条件。

内容的提问来源于stack exchange,提问作者Nejc Kejzar

火山引擎 最新活动