作者CRAZYAWIND (怒火烧不尽)
看板Grad-ProbAsk
标题Re: [理工] [工数]-98成大工科PDE
时间Sun Jan 31 23:46:45 2010
※ 引述《hihaka2001 (hihaka)》之铭言：
: Utt=Uxx+F(x,t)
: U(0,t)=U(L,t)=0
: U(x,0)=p(x)
: Ut(x,0)=0
吃光光 表演要舞台
令U(x t) = ψ(x t) + φ(x t)
代回pde中
ψ (x t) + φ (x t) = ψ (x t) +φ (x t) +F(x t)
tt tt xx xx
令ψ (x t) = ψ (x t)
tt xx
φ (x t) = φ (x t) +F(x t)
tt xx
U(0 t) = ψ(0 t) + φ(0 t) = 0 φ(0 t) = ψ(0 t)=0
U(L,t) = ψ(L t) + φ(L t) = 0 φ(L t)= ψ(L t)=0
Ut(x,0)= ψt(x 0) + φt(x 0) = 0 ψt(x 0) = φt(x 0)=0
U(x,0)=ψ(x 0) + φ(x 0) = p(x) 令ψ(x 0)=p(x) φ(x 0)=0
(1)
nπx ∞
先解φ(x t) 特徵函数{sin── }n=1
L
∞ nπx ∞ nπx
φ(x t)= Σ an(t)sin ── F(x t)=Σ qn(t)sin──
n=1 L n=1 L
代回PDE中
∞ nπ 2 nπx
Σ [an(t)" +(──) an(t) -qn(t)]sin── =0
n=1 L L
由an(0) = an'(0)
t L nπ
an(t) = ∫ ── sin──(t-τ)qn(τ) dτ
0 nπ L
(2) nπx ∞
解ψ(x t) 特徵涵数 {sin── }n=1
L
∞ nπx
ψ(x t)=Σ bn(t)sin ──
n=1 L
代回pde中
∞ nπ 2 nπx
Σ [bn(t)" +(──) bn(t) ]sin── =0
n=1 L L
nπ nπ
bn(t)= Ancos──t + Bnsin──t
L L
∞ nπ nπ nπx
ψ(x t)=Σ [Ancos──t + Bnsin──t]sin ──
n=1 L L L
ψt(x 0) = 0 =Bn
∞ nπx
ψ(x 0) = p(x) = Σ Ansin ──
n=1 L
2 L nπx
An = ──∫ p(x)sin──dx
L 0 L
∞ 2 L nπx nπ nπx
ψ(x t) = Σ [ ──∫ p(x)sin──dx] cos──t sin ──
n=1 L 0 L L L
∞ t L nπ nπx
φ(x t)= Σ [∫ ── sin──(t-τ)qn(τ) dτ]sin ──
n=1 0 nπ L L
U(x t) = ψ(x t) + φ(x t)
献丑了= =
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.105.159.190
1F：→ lalala419:我也是这样解 板上好像都没人补高成= = 01/31 23:47
2F：推 kagato:我也会这样解= = 修边界比较好算XD 01/31 23:51
3F：推 boy210637:耶!! 我算对了 我好感动 T_T 01/31 23:57
4F：→ CRAZYAWIND:听说这题会算的话就赢过1/2考成大工科的了= = 01/31 23:58
5F：→ boy210637:不过我没报工科ˊˋ 01/31 23:59
6F：→ doom8199:跟上一篇一样， g_n(t) 是自己假设的变数 02/01 00:00
7F：→ doom8199:应该要由 F(x,t) 表示才对。 02/01 00:01
8F：→ doom8199:这样最後表示出来的解比较完整 02/01 00:01
9F：→ CRAZYAWIND:应该没办法用F(x,t)来表示吧 @@~ F(x,t)都展成系数了 02/01 00:06
10F：→ CRAZYAWIND:而且都经过Dehamel的积分 而且他本身都是未知的函数了 02/01 00:07
11F：推 doom8199:可以阿，左右同乘 sin(nπt/L) 对t积分 02/01 00:07
12F：→ doom8199:就跟你在求 a_n(t) 的过程完全一样 02/01 00:08
13F：推 boy210637:恩恩 是这样没错 02/01 00:10
14F：→ doom8199:打错，是对x微分 02/01 00:10
15F：→ boy210637:2/L∫ F(x,t)sin(nπx/L)dx 是这样没错吧 02/01 00:12
16F：→ doom8199:嗯嗯，把它带进去 g_n(t) 才算完工 02/01 00:14
※ 编辑: CRAZYAWIND 来自: 59.105.159.190 (02/01 01:22)