作者deepwoody (快回火星吧)
看板Math
标题[微积] 级数解微分方程
时间Sun Dec 26 22:46:43 2010
2
解 8x u''+ 6xu'+ (x-1)u = 0
我的解法: 用method of Frobenius
2
x u''+ (3x/4)u'+ (x-1)u/8 = 0
P_0(x) = 3/4 => b_0 = 3/4 otherwise b_n = 0
P_1(x) = (x-1)/8 => c_0 = -1/8, c_1 = 1/8 otherwise c_n = 0
r ∞ n
则解 u(x) = x Σ (a_n)x , a_n ≠ 0
n=0
indicial eq : F(r) = r(r-1) + (3/4)r - (1/8) = 0 => r = 1/2 or -1/4
n
解递回 F(n+r)a_n + Σ [G_n(k)][a_(n-k)] = 0, (n = 1,2,3,...) ………(*)
k=1
2
F(n+(1/2)) = (n+(1/2))(n-(1/2)) + (3/4)(n+(1/2)) - (1/8) = n + (3/4)n
G_n(k) = (n-k+r)b_k + c_k
=> G_n(1) = (n-1+(1/2))*b_1 + c_1 = c_1 = 1/8, otherwise G_n(k) = 0
(*) 2
===> (n +(3/4)n)a_n + (1/8)a_(n-1) = 0
n
-1 (-1) a_0
=> a_n = ────── a_(n-1) => a_n = ───────，令 a_0 = 1
8n^2 + 6n n 2
Π (8i + 6i)
i=1
n
(1/2) ∞ (-1) n (1/2)
所以一解为 u_1(x) = x Σ ─────── x + x
n=1 n 2
Π (8i + 6i)
i=1
因为(1/2)-(-1/4) = 3/4，非整数，所以同理
n
(-1/4) ∞ (-1) n (-1/4)
另外一解为 u_2(x) = x Σ ─────── x + x
n=1 n 2
Π (8i - 6i)
i=1
则 u(x) = Au_1 + Bu_2
------------------------------------------------------------------------------
我用mathematica解，他给我Bessel和Gama function的型态
我不知道该怎麽转换，还是我算错了呢?
--
╔《新版十二生肖》═════════════════════════════╗
║
◣◣ ˍ ║
●●╰
‧‧╯
◣ ◢ [ ] ιι . .◣
- - ◣ ˍ▁
◢﹎◣▁▂ ║
◢'' .. '〒' '. ' ' ' ξ ▌
. . ≡◢'◎@@ ●﹏●★' ◣
ˊ▄ˋ ★︰ ║
‵′/ ██╯
█/ ▲
╭≡╮ ▃ ￣█ㄟ @@@@ ◣◢ ◥█◤
ˋ▄ˊ ‥ ║
╚═══" " ════════════
"═
"══"═"══
╯══════
liszt1025╝
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.211.87