作者yimao (MLP)
看板Grad-ProbAsk
标题Re: [理工] [工数]-一阶ODE问题~
时间Thu Jul 23 19:27:47 2009
※ 引述《mdpming (★ｐｉｇｍｉｎｇ★)》之铭言：
1.
2 2 2
(x + xy +y )dx-x dy=0
x
y= --------------
tan(ln c/x)
dy du
同除X^2 且令 u = y/x ---- = u + ----
dx dx
整理得
du
( 1 + u + u^2 ) -u -x---- = 0
dx
分离变数
dx du
----- = ------
x 1+u^2
-1
ln│x│ = tan u + c
u = tan(ln(x/c)) = y/x
移项就得到答案...但是X和C的位置颠倒...不知道是哪里算错(算是这样算)
2.
(x-2y+3)dx+(2x-4y-3)dy=0
答案是
3 9 4x 8y
(x-2y+ ---) e e =c
这种...一般用座标平移可转换为齐次
但这题是平行...故无法求交点
因此...
dy du
直接令 x-2y +3 = u .. 1- 2---- = ----
dx dx
带进去整理一下应该就会算的出来
另外我想问一个高阶ODE问题
3x 2
y''+2y'+3y=10e + 4x
这提到後面 答案要怎麽表示
用逆微分运算子
算阿算的
算到
10 3x 1 2 1 2 2
Y = ----- e + --- 4x + --- D 4x +.........
p 18 3 27
所以yp 的D 怎麽求出来@@
1 3x 10 3x
------------10e = ----e
D^2+2D+3 18
而X平方部分
1/3 -2/9D - 2/27 D^2
1 2 ------------------------
----------(4x ) 3 + 2D + D^2 │ 1
D^2+2D+3 1 + 2/3D + 2/3 D^2
------------------------
- 2/3D - 2/3 D^2
- 2/3D - 4/9 D^2 - 2/9D^3
-------------------------
- 2/9 D^2
- 2/9 D^2
2
= (1/3 -2/9D - 2/27 D^2) (4x ) ...变成用微的
D^3以後不考虑是因为微到4X^2为0
4 16 16
= ----- x^2 - ----- x - ---- ...Yp 应该长这样...假设没算错
3 9 27
最後Y=Yh+Yp
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.32.91.86
1F：→ doom8199:<2> 假设 u=x-2y <3>令 yp = Ae^(3x) + Bx^2 + Cx + D 07/23 18:31
2F：→ doom8199:带入解系数。 <1> 假设 u=y/x後 ，O.D.E.一定可以用 07/23 18:32
3F：→ doom8199:分离变数去解，原po可以在算一下 07/23 18:33
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 123.110.19.9
※ 编辑: yimao 来自: 123.110.19.9 (07/23 19:44)
4F：推 pimday:第一题算出来也跟原po一样 07/23 20:02
5F：推 doom8199:<1> 答案是一样的: tan(c1 + lnx) 07/23 20:33
6F：→ doom8199:= [tan(c1) + tan(lnx)]/[1 - tan(c1)*tan(lnx)] 07/23 20:34
7F：→ doom8199:= [1 + tan(lnx)/tan(c1)]/[1/tan(c1) - tan(lnx)] 07/23 20:35
8F：→ doom8199:= [1 + tan(lnx)*tan(c2)]/[tan(c2) - tan(lnx)] 07/23 20:36
9F：→ doom8199:= 1/tan(c2 - lnx) = 1/tan[ln(c3/x)] 07/23 20:37
10F：→ doom8199:where tan(c2) = 1/tan(c1) 、 ln(c3) = c2 07/23 20:38
11F：→ doom8199:或是直接对 O.D.E. 假设 u=x/y 也能得到那样的型态 07/23 20:39
12F：推 mdpming:太感谢了~~来研究一下!! 07/23 22:12
13F：推 mdpming:光那个长除法要打很久吧..=.= 07/23 22:21
14F：→ mdpming:D^3以後不考虑是因为微到4X^2为0 <--- 讲这句我就了解QQ 07/23 22:21