上海交通大学：《常微分方程》课程教学资源（讲座稿1）第七讲 一阶隐式微分方程（2/2）、高阶微分方程与Mathematica求解

第七讲、一阶隐式微分方程-2 高阶微分方程与Mathematica求解 张祥 xzhang@sjtu.edu.cn 答疑时间：周三晚上6：30-8：20点 答疑地点：老图书馆数学楼301 张样：上海交通大学数学系 第七讲、一阶隐式微分方程2、高阶微分方程与Mathematica求解

1‘˘!ò¤™á©êß-2 pá©êßÜMathematica¶) ‹ å xzhang@sjtu.edu.cn â¶ûmµ±n˛ 6:30–8:20 : â¶/:µP„÷,ÍÆ¢ 301 ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

本讲教学目的与目标 ·一阶隐式微分方程 x.y,dx =0 (1) 。不显含自变量x或y的方程的求解 。高阶微分方程的求解 。Mathematica求解 口0·4之·4生+2刀a0 张样：上海交通大学数学系 第七讲、一阶隐式微分方程-2、高阶教分方程与Mathematica求解

˘Æ8Ü8I ò¤™á©êß F  x, y, dy dx = 0. (1) ÿw¹gC˛ x ½ y ê߶) pá©ê߶) Mathematica ¶) ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

回顾：y可解出的方程的求解 设方程(1)可写成 y=f,p:p=衣 (2) 其中f(x,p)连续可微. 将p作为新的因变量，对方程(2)两边关于x求导，得到p关 于x的导数的显示方程 (f(x;p)-p)dx+p(x,p)dp=0. (3) 口164怎生12月风0 张样：上海交通大学数学系第七讲、一阶隐式微分方程2、高阶微分方程与Mathematica求解

£µy å)—ê߶) êß (1) å§ y = f(x,p), p = dy dx , (2) Ÿ• f(x,p) ÎYåá. Ú p äè#œC˛, Èêß (2) ¸>'u x ¶,  p ' u x Íw´êß (fx(x,p)−p)dx+fp(x,p)dp = 0. (3) ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

●若方程(3)有通解p=u(x,c),则方程(2)有通解 y=f(x,u(x,c),其中c是任意常数 对特解有类似的结论 ·若方程(3)有通解x=v(p,c),则方程(2)有含参数p的通解 x=v(p;c); 其中c是任意常数. y=f(v(p,c);p), 对特解有类似的结论。 4口6+4之·4生+2月a0四 张样：上海交通大学数学系 第七讲、一阶隐式微分方程-2、高阶教分方程与Mathematica求解

eêß (3) kœ) p = u(x, c), Kêß (2) kœ) y = f(x,u(x, c)), Ÿ• c ¥?ø~Í. ÈA)kaq(ÿ. eêß (3) kœ) x = v(p, c), Kêß (2) k¹ÎÍ p œ) x = v(p, c), y = f(v(p, c),p), Ÿ• c ¥?ø~Í. ÈA)kaq(ÿ. ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

例1：求解Clairaut方程 y=xp+f(p), 会ro0 p= 解：对Clairaut方程两边关于x求导得 d迎=0. (c+fp》 o由裴=0得到Clairaut方程的通解 y=cx+f(c), 其中c是任意常数.注意：Clairaut方程的通解由一族直线构 成. ●由x+f'(p)=0得到Clairaut方程的特解 x=-f(p),y=xp+f(p), 其中p是参数， 日回4他”“生12月风 张样：上海交通大学数学系 第七讲、一阶隐式微分方程2、高阶微分方程与Mathematica求解

~ 1: ¶) Clairaut êß y = xp+f(p), p = dy dx , f 00(p) 6= 0. ): È Clairaut ê߸>'u x ¶ (x+f 0 (p))dp dx = 0. d dp dx = 0  Clairaut êßœ) y = cx+f(c), Ÿ• c ¥?ø~Í. 5ø: Clairaut êßœ)dòxÜÇ §. d x+f 0 (p) = 0  Clairaut êßA) x = −f 0 (p), y = xp+f(p), Ÿ• p ¥ÎÍ. ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

进一步地，因f”(p)≠0，运用隐函数存在定理从 x=-f'(p) 解出 p=ω(x 则Clairaut方程的特解可写成 y=xo(x)+f(@(x)). (4) 思考：Clairaut方程在常微分方程基础理论中占有重要的地位。 ●Clairaut方程的特解与通解的关系如何？ 4口6+4之·4生+2风0 张样：上涛交通大学数学系 第七讲、一阶隐式微分方程-2、高阶教分方程与Mathematica求解

?ò⁄/, œ f 00(p) 6= 0, $^¤ºÍ3½nl x = −f 0 (p) )— p = ω(x). K Clairaut êßA)å§ y = xω(x) +f(ω(x)). (4) g: Clairaut êß3~á©ê߃:nÿ•”k­á/†" Clairaut êßA)Üœ)'XX¤º ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

Clairaut方程的特解与通解的关系如下： ●过Clairaut方程特解上任一点的切线是通解中的一条直线. 事实上，设(0，y(xo)是特解(4)上的任一点，则 有y(xo)=o(xo.所以特解过(x0,y(x0)的切线方程为 y-y(xo)=o(xo)(x-xo), 即 y=cox+f(co), c0=0(x0): 。特解(4)不能用通解表示，因为o(x)不是常数.事实上， 从x+f'(o(x)三0得 1 o国=o≠0. ●上述两点说明Clairaut方程特解对应的积分曲线上任一点都 有通解中的一条积分曲线通过，且两者在该点相切， 张样：上寿交通大学数学系第七讲、一阶隐式微分方程2、高阶微分方程与Mathematica求解

Clairaut êßA)Üœ)'XXeµ L Clairaut êßA)˛?ò:ÉÇ¥œ)•ò^ÜÇ. Ø¢˛,  (x0, y(x0)) ¥A) (4) ˛?ò:, K k y 0 (x0) = ω(x0). §±A)L (x0, y(x0)) ÉÇêßè y−y(x0) = ω(x0)(x−x0), = y = c0x+f(c0), c0 = ω(x0). A) (4) ÿU^œ)L´, œè ω(x) ÿ¥~Í. Ø¢˛, l x+f 0 (ω(x)) ≡ 0  ω 0 (x) = − 1 f 00(ω(x)) 6= 0. ˛„¸:`² Clairaut êßA)ÈA»©­Ç˛?ò:— kœ)•ò^»©­ÇœL, Ö¸ˆ3T:ÉÉ. ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

Clairaut方程是法国数学家Alexis Clairaut于1734年引入的. Alexis Clairaut(3 May 1713-17 May 1765)was a prominent French mathematician,astronomer,and geophysicist. o His paper procured his admission into the French Academy of Sciences in 1731,although he was below the legal age as he was only eighteen. o He was elected a Fellow of the Royal Society of London in 1737. 特别是在天文学的研究方面做出了杰出的贡献： o He obtained an ingenious approximate solution of the problem of the three bodies o In 1750 he gained the prize of theSt Petersburg Academy for his essay Theorie de la lune ●ln1759 he calculated the perihelion(近日点)of Halley's comet. 2月风0 张样：上涛交通大学数学系 第七讲、一阶隐式微分方程-2、高阶教分方程与Mathematica求解

Clairaut êߥ{IÍÆ[ Alexis Clairaut u 1734 c⁄\. Alexis Clairaut (3 May 1713–17 May 1765) was a prominent French mathematician, astronomer, and geophysicist. His paper procured his admission into the French Academy of Sciences in 1731, although he was below the legal age as he was only eighteen. He was elected a Fellow of the Royal Society of London in 1737. AO¥3U©ÆÔƒê°â— #—￾zµ He obtained an ingenious approximate solution of the problem of the three bodies In 1750 he gained the prize of theSt Petersburg Academy for his essay ThWorie de la lune In 1759 he calculated the perihelion(CF:) of Halley’s comet. and so on Clairaut seek to make his pleasures into his day to day work, he lost rest, health, and finally life at the age of fifty-two." ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

讨论+思考：在Clairaut方程中 通过f(x)的适当选取可以得到特解的哪些可能的几何形状？ 口1回1怎：“主12月双0 张样：上将交通大学数学系第七讲、一阶隐式微分方程2、高阶微分方程与Mathematica求解

?ÿ+g: 3 Clairaut êß• œL f(x) ·¿å±A)= åUA¤/G? ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)

例2.求解方程 =-p+式 p=y(x) 解：对上述方程两边关于x求导得 p-()=0 利用该方程的特解和通解得到原方程的特解和通解为 y=,y=+c+2 其中c是任意常数. 观察与思考： ·例2中的特解与通解的关系如何？ 口0·4之·4生+2刀a0 张样：上涛交通大学数学系 第七讲、一阶隐式微分方程-2、高阶微分方程与Mathematica求解

~ 2. ¶)êß y = p 2 −xp+ 1 2 x 2 , p = y 0 (x). ): È˛„ê߸>'u x ¶ (2p−x)  1− dp dx = 0. |^TêßA)⁄œ)êßA)⁄œ)è y = 1 4 x 2 , y = 1 2 x 2 +cx+c 2 , Ÿ• c ¥?ø~Í. * Üg: ~ 2 •A)Üœ)'XX¤º ‹å: ˛°œåÆÍÆX 1‘˘!ò¤™á©êß-2!pá©êßÜMathematica¶)