作者CRAZYAWIND (怒火烧不尽)
看板Grad-ProbAsk
标题Re: [理工] [工数]-边界问题
时间Sun Jan 31 19:59:05 2010
※ 引述《hihaka2001 (hihaka)》之铭言:
: 请问大家
: Ut=4Uxx+xt
: U(0,t)=U(pi,t)=0
: U(x,0)=f(x)
: 可是这题是和t有关的
: 所以要修边界
: 但是边界没问题都为0
: 所以我要令
: U(x,t)=ψ(x,t)+ 0
: 好像这样带进去会有问题
: 可是如果直接用特徵函数展开
: 去解
: 也会遇到问题
: 请问这题该怎麽解决啊!!
: xt = sigma (1~inf)a(t)sin(nx)
: a(t)=(2 *t /pi) Int (0~pi) x sinnx dx
: 是这个意思吗
: 然後可以解出
: a(t)=2t/n * (1-(-1)^n)
: yp=1/ D^2 +4n^2 .....
: 这样解吗??
∞
令u(x t) = Σ an(t)sinnx
n=1
∞
x = Σ ψsinnx
n=1 n+1
2 π 2(-1)
ψ = ──∫ xsinnx dx = ────
π 0 n
因为ψ是常数我後面就用ψ带着跑了= =
代回PDE中
∞ ' 2
Σ [an(t) + 4n an(t) - ψt]sinnx = 0
n=1
' 2
解an(t) + 4n an(t) - ψt = 0
-4n^2t ψt ψ
an(t) = Ane + ─── - ───
4n^2 16n^2
∞ -4n^2t ψt ψ
u(x t) = Σ [Ane + ─── - ───]sinnx
n=1 4n^2 16n^2
∞ ψ
u(x 0) = f(x) = Σ [An - ───]sinnx
n=1 16n^2
ψ 2 π/4 40 π
An- ─── = ──∫ 20sinnxdx = ───(1-cos──)
16n^2 π 0 nπ 4
40 π ψ
An = ──── (1-cos──) + ────
nπ 4 16n^2
∞ 40 π ψ -4n^2t ψt ψ
u(x t) = Σ {[──(1-cos──) + ────]e + ─── - ───} sinnx
n=1 nπ 4 16n^2 4n^2 16n^2
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.105.159.190
1F:推 iyenn: m(_ _)m 01/31 20:02
2F:推 bjackal:感谢 刚刚卡在求an(t)那边 终於了解了 01/31 20:03
3F:→ bjackal:i大在偷拜拜 另一提一阶二次ODE 我不会! 01/31 20:04
※ 编辑: CRAZYAWIND 来自: 59.105.159.190 (01/31 20:04)
4F:→ mdpming:..=.= 01/31 20:05
5F:推 hihaka2001:ψ这个我算的是 2((-1)^n +1 )/n 01/31 20:06
6F:→ CRAZYAWIND:你应该算错了吧= = 下限是0 代入的话有个x在会是0 01/31 20:22
※ 编辑: CRAZYAWIND 来自: 59.105.159.190 (02/01 00:38)