PDEs
PDEs
ﺗﻌﺮﻳﻒ ۱.۱ﻣﻌﺎﺩﻟﻪ ﺍﻱ ﻛﻪ ﺩﺭ ﺁﻥ ﻣﺸﺘﻖ ﺍﺯ ﺗﺎﺑﻌﻲ ﺑﺎﺷﺪ ﻛﻪ ﺁﻥ ﺗﺎﺑﻊ ﺩﺍﺭﺍﻱ ﺩﻭ ﻳﺎ ﭼﻨﺪ ﻣﺘﻐﻴﺮ ﻣﺴﺘﻘﻞ ) (x, y, z, tﺑﺎﺷﺪ.
)∂U (x, t ∂U
→ )ex : U (x, t =4 =0
∂x ∂t
ﺗﻌﺮﻳﻒ P DE ۲.۱ﺭﺍ ﺧﻄﻲ ﻣﻲ ﻧﺎﻣﻴﻢ ﺍﮔﺮ P DEﻧﺴﺒﺖ ﺑﻪ ﺗﺎﺑﻊ ﻣﺠﻬﻮﻝ ﻭ ﻣﺸﺘﻘﺎﺕ ﺁﻥ ﺧﻄﻲ ﺑﺎﺷﺪ .ﻋﺎﻣﻞ ﻫﺎﻱ ﻏﻴﺮ
ﺧﻄﻲ ﺗﻮﺍﻥﻫﺎﻳﻲ ﺍﺯ Uﻭ ﻣﺸﺘﻘﺎﺕ -ﺩﺍﺧﻞ ﻛﻤﺎﻥ ﺗﻮﺍﺑﻊ ﻏﻴﺮ ﺧﻄﻲ Uﻭ ﻣﺸﺘﻘﺎﺕ ﺁﻥ ﺑﺎﺷﺪ -ﺣﺎﺻﻞ ﺿﺮﺏ ﺗﺎﺑﻊ ﻣﺠﻬﻮﻝ Uﺩﺭ
ﻣﺸﺘﻘﺎﺕ
∂U ∂2U ∂U ∂U
ex : = 2
ﺧﻄﻲ = ﺧﻄﻲ + ex
∂t ∂x ∂t ∂x
∂U ∂2U ∂2U ∂2U ∂U
ex : = 2
+ U 2
ﻏﻴﺮﺧﻄﻲ 2
= 2
= × U, ﻏﻴﺮﺧﻄﻲ
∂t ∂x ∂x ∂y ∂x
۱
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﺗﻌﺮﻳﻒ P DE ۳.۱ﺭﺍ ﻫﻤﮕﻦ ﮔﻮﻳﻴﻢ ﺍﮔﺮ ﻓﻘﻂ ﺷﺎﻣﻞ Uﻭ ﻣﺸﺘﻘﺎﺕ ﺁﻥ ﺑﺎﺷﺪ .ﻣﺜﻼ
∂2U ∂2U
= + sin(U ).
∂tx2 ∂x2
ﺗﻌﺮﻳﻒ P DE ۴.۱ﺭﺍ ﻧﺎﻫﻤﮕﻦ ﮔﻮﻳﻴﻢ ﺍﮔﺮ ﻋﻼﻭﻩ ﺑﺮ Uﻭ ﻣﺸﺘﻘﺎﺕ ﺁﻥ ﺩﺍﺭﺍﻱ ﺗﺎﺑﻌﻲ ﺍﺯ ﻣﺘﻐﻴﺮﻫﺎﻱ ﻣﺴﺘﻘﻞ ﺑﺎﺷﺪ
∂U ∂U
= + ex .
∂t ∂x
ﺗﻌﺮﻳﻒ ۵.۱ﻣﺮﺗﺒﻪ ﻱ P DEﺑﺰﺭﮔﺘﺮﻳﻦ ﻣﺮﺗﺒﻪ ﻱ ﻣﺸﺘﻖ ﻣﻮﺟﻮﺩ ﺩﺭ ﻣﻌﺎﺩﻟﻪ ﺍﺳﺖ.
ﺗﻌﺮﻳﻒ U ۶.۱ﺭﺍ ﭘﺎﺳﺦ ﻳﻚ P DEﺩﺭ ﻓﻀﺎﻱ Rﮔﻮﻳﻨﺪ ﺑﻪ ﺷﺮﻁ ﺁﻧﻜﻪ ﺍﻭﻻ ﺩﺭ ﺁﻥ ﻓﻀﺎ ،ﻣﺸﺘﻘﺎﺕ ﺗﻌﺮﻳﻒ ﺷﺪﻩ ﺩﺭ ﻣﻌﺎﺩﻟﻪ ﺭﺍ
ﺩﺍﺷﺘﻪ ﺑﺎﺷﺪ ﺛﺎﻧﻴﺎ ﺩﺭ ﻣﻌﺎﺩﻟﻪ ﺻﺪﻕ ﻛﻨﺪ.
∂U 2
)= −π 2 e−π t sin(πx
∂t
) ∂U = πe−π2 t cos(πx ∂U ∂2U
∂x → =
∂t ∂x2
∂2U
2
)= −π 2 e−π t sin(πx
∂x 2
ﺗﻌﺮﻳﻒ ۷.۱ﻫﻤﺎﻳﺶ ﺟﻮﺍﺏﻫﺎ :ﺍﮔﺮ U1 , U2 , ..., Unﺟﻮﺍﺏ P DEﺑﺎﺷﻨﺪ ﻭ ﺍﮔﺮ P DEﺗﻌﺮﻳﻒ ﺷﺪﻩ ﺧﻄﻲ ﻭ ﻫﻤﮕﻦ
ﺑﺎﺷﺪ ﺩﺭ ﺍﻳﻦ ﺻﻮﺭﺕ U = c1 U1 + c2 U2 + ... + cn Unﻧﻴﺰ ﺟﻮﺍﺏ ﻣﻌﺎﺩﻟﻪ ﻣﻲ ﺑﺎﺷﻨﺪ.
۲
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﻣﺜﺎﻝ ۱
ﺍﻟﻒ( ﺍﮔﺮ ← ∆(x, y) > 0ﻣﻌﺎﺩﻟﻪ ﻱ gﺩﺭ ﻧﻘﻄﻪ ﻱ ) (x, yﻫﺬﻟﻮﻟﻮﻱ ﺍﺳﺖ ) .ﻣﻌﺎﺩﻟﻪ ﻱ ﻣﻮﺝ(
ﺏ( ﺍﮔﺮ ← ∆(x, y) = 0ﻣﻌﺎﺩﻟﻪ ﻱ gﺩﺭ ﻧﻘﻄﻪ ﻱ ) (x, yﺳﻬﻤﻮﻱ ﺍﺳﺖ ) .ﻣﻌﺎﺩﻟﻪ ﻱ ﮔﺮﻣﺎ(
ﺝ( ﺍﮔﺮ ← ∆(x, y) < 0ﻣﻌﺎﺩﻟﻪ ﻱ gﺩﺭ ﻧﻘﻄﻪ ﻱ ) (x, yﺑﻴﻀﻮﻱ ﺍﺳﺖ.
→ Ut = Uxx
Uxx − Ut = 0
a=1 b=0 c=0 ﺳﻬﻤﻮﻱ → ∆ = 0
۳
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﺑﻪ ﺗﻌﺪﺍﺩ ﺑﺎﻻﺗﺮﻳﻦ ﻣﺮﺗﺒﻪ ﻱ ﻣﺸﺘﻖ ﺷﺮﻁ ﻻﺯﻡ ﺍﺳﺖ ﺗﺎ ﺟﻮﺍﺏ ﺧﺼﻮﺻﻲ ﺭﺍ ﺑﺪﺳﺖ ﺁﻭﺭﻳﻢ.
∂U ∂2U
= U (x, t) =? a ≤ x ≤ b ,t ≥ 0
∂t ∂x2
) Ux (a, t) = p(t
→ ﺷﺮﻁ ﻣﺮﺯﻱ ﻧﻴﻮﻣﺎﻥ
) Ux (b, t) = q(t
∂2U ∂2U
= a≤x≤b ,0 ≤ t ≤ T
∂t2 ∂x2
۴
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
) U (a, t) = p(t
→ ﺷﺮﻁ ﻣﺮﺯﻱ ﺩﻳﺮﻳﻜﻠﻪ
) U (b, t) = q(t
) Ux (a, t) = p(t
→ ﺷﺮﻁ ﻣﺮﺯﻱ ﻧﻴﻮﻣﺎﻥ
) Ux (b, t) = q(t
) U (x, 0) = f (t
→ ﺷﺮﺍﻳﻂ ﺍﻭﻟﻴﻪ
) Ut (x, 0) = g(t
ﻳﺎﺩﺁﻭﺭﻱ ۷.۱
ﻣﻌﺎﺩﻻﺕ ﻣﺮﺗﺒﻪ ﻱ ﺍﻭﻝ ﺧﻄﻲ
′
)y (x) + p(x)y(x) = g(x
)ﻫﻤﮕﻦ(
∫
− p(x)dx
g(x) = 0 → y(x) = ce
∫
− ∫ p(x)dx ∫
g(x) ̸= 0 → y(x) = −e ][ (g(x)e p(x)dx )dx + c
۵
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﻳﺎ
y(x) = c1 cosh r1 x + c2 sinh r2 x
-ﻣﺮﺣﻠﻪ ﻱ ﺍﻭﻝ :ﻓﺮﺽ ﻛﻨﻴﺪ ﺟﻮﺍﺏ ﻣﻌﺎﺩﻟﻪ ﻱ ﮔﺮﻣﺎ ﺑﻪ ﺻﻮﺭﺕ ﺭﻭﺑﺮﻭ ﺑﺎﺷﺪ:
۶
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
۷
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﭘﺲ
c2 = −c1 = 0 → F (x) = 0
: F (x) = eαx (c1 cos βx + c2 sin βx) = c1 cos µx + c2 sin µx
-ﻣﺮﺣﻠﻪ ﻱ ﺷﺸﻢ :ﻣﻌﺎﺩﻟﻪ ﻱ ﺯﻣﺎﻧﻲ ﺭﺍ ﺩﺭ ﺣﺎﻟﺘﻲ ﻛﻪ ﻣﻌﺎﺩﻟﻪ ﻱ ﻣﻜﺎﻧﻲ ﺟﻮﺍﺏ ﺩﺍﺩ ﺣﻞ ﺧﻮﺍﻫﻴﻢ ﻛﺮﺩ.
nπ 2 nπ
G′ (t) − kG(t) = 0 (k = − ) ) = −λ2n (λ2n = ( )2
L L
→ G′n (t) + λn Gn (t) = 0
Gn (t) = c1 e−λn t
۸
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
-ﻣﺮﺣﻠﻪ ﻱ ﻫﻔﺘﻢ :ﺗﻮﺍﺑﻊ ) Gn (t), Fn (tﺭﺍ ﺩﺭ ﻫﻢ ﺿﺮﺏ ﻛﺮﺩﻩ ﻭ ) U (x, tﺑﺪﺳﺖ ﻣﻲ ﺁﻳﺪ.
nπ 2 nπ 2
(U (x, t) = F (x)G(t) → Un (x, t) = c2n sin x) × c1n e−λn t → Un (x, t) = bn sin( x)e−λn t
L L
-ﻣﺮﺣﻠﻪ ﻱ ﻫﺸﺘﻢ :ﻫﻤﺎﻳﺶ ﺟﻮﺍﺏ ﻫﺎ
∑
∞ ∑
∞
2 nπ
= )U (x, t = )Un (x, t (bn e−λn t sin )x
n=1 n=1
L
∫ L
2 nπ
= bn (f (x) sin x)dx
L 0 L
۹
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
۱۰
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
: r2 − µ2 = 0 → r1 = µ, r2 = −µﻣﻌﺎﺩﻟﻪ ﻣﺸﺨﺼﻪ
: F (x) = c1 eµx + c2 e−µxﺟﻮﺍﺏ ﻋﻤﻮﻣﻲ
)ﺟﻮﺍﺏ ﺑﺪﻳﻬﻲ ﻣﻲ ﺑﺎﺷﺪ→ F (x) = 0 (.
: r2 + µ2 = 0 → r = ±iµﻣﻌﺎﺩﻟﻪ ﻣﺸﺨﺼﻪ
: F (x) = c1 cos µx + c2 sin µxﺟﻮﺍﺏ ﻋﻤﻮﻣﻲ
nπ
= → Fn (x) = c2n sin µn x, µn n = 1, 2, ...
L
-ﻣﺮﺣﻠﻪ ﻱ ﺷﺸﻢ :ﺣﺎﻻ ﻣﻌﺎﺩﻟﻪ ﻱ ﺯﻣﺎﻧﻲ ) (۱ﺭﺍ ﺣﻞ ﻣﻲ ﻛﻨﻴﻢ:
G′′ (x) − c2 kG(t) = 0 t ≥ 0
2 2
k = −µ2 = − n π
n
L2
n2 π 2 cnπ
→ G′′n (t) + c2 Gn (t) = 0 = → λn
L2 L
ﻓﺮﺽ ﻛﻨﻴﻢ G′′n (t) + λ2n Gn (t) = 0
ﺣﻞ ﻣﻌﺎﺩﻟﻪ ﻣﺮﺗﺒﻪ ﻱ ۲
-ﻣﺮﺣﻠﻪ ﻱ ﻫﻔﺘﻢ Gn (t), Fn (t) :ﺑﺪﺳﺖ ﺁﻣﺪﻩ ﺭﺍ ﺩﺭ ﻫﻢ ﺿﺮﺏ ﻣﻲ ﻛﻨﻴﻢ ﻭ ) U (x, tﺣﺎﺻﻞ ﻣﻲ ﺷﻮﺩ:
nπ
(U (x, t) = Fn (t)Gn (t) → Un (x, t) = crn sin ))x) × a1n cos(λn t) + a2n sin(λn t
L
-ﻣﺮﺣﻠﻪ ﻱ ﻫﺸﺘﻢ :ﻫﻤﺎﻳﺶ ﺟﻮﺍﺑﻬﺎ
∑
∞ ∑
∞
= )U (x, t = )Un (x, t
n=1 n=1
∑∞
nπ
= )U (x, t ((bn cos(λn t) + an sin(λn t)) sin )x
n=1
L
۱۱
ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ:ﻣﺪﺭﺱ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
۱۴۰۱-۱۴۰۲ ﺗﺮﻡ ﺩﻭﻡ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
: an , bn ﻣﺤﺎﺳﺒﻪ ﻱ
∑
∞
nπ
(۱) ﺷﺮﻁ ﺍﻭﻟﻴﻪf (x) = U (x, 0) = bn sin( x)
n=1
L
∫ L
2 nπ
ﺿﺮﺍﻳﺐ ﻓﻮﺭﻳﻪ ﺳﻴﻨﻮﺳﻲbn = f (x) sin(
x)dx
L 0 L
∑∞
nπ
(۲) ﺷﺮﻁ ﺍﻭﻟﻴﻪg(x) = Ut (x, 0) = λn an sin( x)
n=1
L
۱۲
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
1
])U (x, t) = [f (x − ct) + f (x + ct ﺍﮔﺮ ← g(x) = 0
2
∫ x+ct
1 1
U (x, t) = [f (x − ct) + f (x + ct)] + g(s)ds ﺍﮔﺮ ← g(x) ̸= 0
2 2c x−ct
ﻧﻜﺘﻪ ۱.۱ﺍﮔﺮ ﺷﺮﺍﻳﻂ ﻣﺮﺯﻱ ،ﺩﻳﺮﻳﻜﻠﻪ ﺑﺎﺷﺪ ،ﺑﻪ ﺑﺴﻂ sinﻣﻲ ﺭﺳﻴﻢ.
ﻧﻜﺘﻪ ۲.۱ﺍﮔﺮ ﺷﺮﺍﻳﻂ ﻣﺮﺯﻱ ،ﻧﻴﻮﻣﺎﻥ ﺑﺎﺷﺪ ،ﺑﻪ ﺑﺴﻂ cosﻣﻲ ﺭﺳﻴﻢ.
۱۳
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﺍﺯ ﺁﻧﺠﺎﻳﻲ ﻛﻪ ﻣﻲﺧﻮﺍﻫﻴﻢ ﻣﻌﺎﺩﻟﻪ ﻣﻮﺝ ﺟﺪﻳﺪ ،ﻫﻤﮕﻦ ﺑﺎﺷﺪ ،ﺷﺮﻁ ﺯﻳﺮ ﺭﺍ ﺍﻋﻤﺎﻝ ﺧﻮﺍﻫﻴﻢ ﻛﺮﺩ
−1
c2 w′′ (x) + H(x) = 0, → = )w′′ (x )H(x
c2
ﺍﻳﻦ ﻣﻌﺎﺩﻟﻪ ،ﻳﻚ ﻣﻌﺎﺩﻟﻪ ﺩﻳﻔﺮﺍﻧﺴﻴﻞ ﻣﺮﺗﺒﻪ ﺩﻭ ﻣﻲﺑﺎﺷﺪ ﻛﻪ ﺑﺮﺍﻱ ﺣﻞ ﺁﻥ ﺑﺎﻳﺪ ﺩﻭ ﺷﺮﻁ ﺩﺍﺷﺘﻪ ﺑﺎﺷﻴﻢ .ﺑﺮﺍﻱ ﻳﺎﻓﺘﻦ ﺍﻳﻦ ﺩﻭ
ﺷﺮﻁ ،ﺍﺯ ﺷﺮﻁ ﻣﺮﺯﻱ ﻣﻌﺎﺩﻟﻪ ﺍﺳﺘﻔﺎﺩﻩ ﺧﻮﺍﻫﻴﻢ ﻛﺮﺩ
ﺍﺯ ﻃﺮﻓﻲ ﺑﺮﺍﻱ ﻳﺎﻓﺘﻦ ﺷﺮﺍﻳﻂ ﺍﻭﻟﻴﻪ ﻣﻌﺎﺩﻟﻪ ﻣﻮﺝ ﺟﺪﻳﺪ ،ﺩﺍﺭﻳﻢ
۱۴
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﻧﺘﻴﺠﻪ ﮔﻴﺮﻱ :ﺑﺮﺍﻱ ﺣﻞ ﻣﻌﺎﺩﻟﻪ ﻣﻮﺝ ﻧﺎﻫﻤﮕﻦ ﺯﻳﺮ ،ﺑﺎﻳﺪ ﻣﺮﺍﺣﻞ ﺯﻳﺮ ﺭﺍ ﺍﻧﺠﺎﻡ ﺩﻫﻴﻢ:
Utt = c2 Uxx + H(x),
0 < x < L,
)U (x, 0) = f (x
)Ut (x, 0) = g(x
U (0, t) = A
U (L, t) = B
-ﻣﺮﺣﻠﻪ ﻱ ﺳﻮﻡ w :ﺑﺪﺳﺖ ﺁﻣﺪﻩ ﺍﺯ ﻣﺮﺣﻠﻪ ﻱ ﺍﻭﻝ ﻭ Vﺑﺪﺳﺖ ﺁﻣﺪﻩ ﺍﺯ ﻣﺮﺣﻠﻪ ﻱ ﺩﻭﻡ ﺭﺍ ﺑﺎ ﻫﻢ ﺟﻤﻊ ﻣﻲ ﻛﻨﻴﻢ ﻭ Uﺑﺪﺳﺖ
ﻣﻲ ﺁﻳﺪ.
)U (x, t) = V (x, t) + w(x
۱۵
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﭼﻮﻥ ﻣﻌﺎﺩﻟﻪ ﻧﺎﻫﻤﮕﻦ ﺍﺳﺖ ﺍﺯ ﺟﺪﺍﺳﺎﺯﻱ ﻣﺘﻐﻴﺮﻫﺎ ﻧﻤﻲ ﺗﻮﺍﻥ ﺍﺳﺘﻔﺎﺩﻩ ﻛﺮﺩ .ﺷﺮﺍﻳﻂ ﻣﺮﺯﻱ ﺩﻳﺮﻳﻜﻠﻪ ﺍﺳﺖ ﭘﺲ ﺟﻮﺍﺏ ﺑﺴﻂ
ﺳﻴﻨﻮﺳﻲ ﺍﺳﺖ ﻭ ﭼﻮﻥ Uﺑﺮ ﺣﺴﺐ x, tﺍﺳﺖ ﺑﺎﻳﺪ ﺿﺮﻳﺐ ﺑﺮ ﺣﺴﺐ tﺑﺎﺷﺪ.
∑
∞
nπ
= )U (x, t Un (t) sin x
n=0
L
ﻫﺪﻑ ﺑﺪﺳﺖ ﺁﻭﺭﺩﻥ ) Un (tﺍﺳﺖ .ﺑﻪ ﻋﻼﻭﻩ ﻓﺮﺽ ﻛﻨﻴﻢ ﺑﺴﻂ ﺳﻴﻨﻮﺳﻲ ﺗﺎﺑﻊ ) h(x, tﺑﻪ ﺻﻮﺭﺕ ﺯﻳﺮ ﺑﺎﺷﺪ.
∑
∞
nπ
= )h(x, t hn (t) sin x
n=1
L
∫ L
2 nπ
= )hn (t (h(x, t) sin x)dx
L 0 L
ﺭﻭﺍﺑﻂ ﺑﺪﺳﺖ ﺁﻣﺪﻩ ﺭﺍ ﺩﺭ ﻣﻌﺎﺩﻟﻪ ﺟﺎﻳﮕﺬﺍﺭﻱ ﻣﻲ ﻛﻨﻴﻢ
∑
∞
nπ ∑
∞
n2 π 2 nπ
= : Utt Un′′ (t) sin( x), Uxx = − (Un (t) sin )x
n=1
L n=1
L2 L
∑
∞
nπ ∑∞
c 2 n2 π 2 nπ ∑∞
nπ
Utt − c2 Uxx → )= h(x, t Un′′ (t) sin( x) + 2
U n )(t (sin )x = )hn (t) sin( x
n=1
L n=1
L L n=1
L
cnπ ∑∞
nπ ∑∞
nπ
= λn → = )(Un′′ (t) + λ2n Un (t)) sin( x )hn (t) sin( x
L n=1
L n=1
L
cnπ
= → Un′′ (t) + λ2n Un (t) = hn (t), λn
L
∫ t
1
→ Un (t) = c1n cos λn t + bn sin λn t + hn (s) sin λn (t − s)ds
λn 0
∑∞ ∑∞ ∫ t
nπ 1 nπ
= )U (x, t = )Un (t) sin( x [an cos λn t + bn sin λn t + )hn (s) sin λn (t − s)ds] sin( x
n=1
L n=1
λn 0 L
۱۶
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﻫﺪﻑ :ﺗﺸﻜﻴﻞ ﻳﻚ ﻣﻌﺎﺩﻟﻪ ﻱ ﻣﻮﺝ ﺟﺪﻳﺪ ﺑﻪ ﻓﺮﻡ ﺣﺎﻟﺖ ﺩﻭﻡ ﻧﺎﻫﻤﮕﻨﻲ ﺍﺳﺖ.
ﺑﺮﺍﻱ ﺍﻳﻦ ﻣﻨﻈﻮﺭ ﺗﻐﻴﻴﺮ ﻣﺘﻐﻴﺮ ﺯﻳﺮ ﺭﺍ ﺗﻌﺮﻳﻒ ﻣﻲ ﻛﻨﻴﻢ:
۱۷
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
)q(t) − p(t
w(x, t) = p(t) + x
L
ﺷﺮﺍﻳﻂ ﺍﻭﻟﻴﻪ ﻱ ﻣﻮﺝ ﺟﺪﻳﺪ : V
۱۸
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
۱۹
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
ﻣﺜﺎﻝ ۳ﻣﻌﺎﺩﻟﻪ ﺭﻭﺑﺮﻭ ﺭﺍ ﺑﺎ ﺭﻭﺵ ﺗﺒﺪﻳﻞ ﻓﻮﺭﻳﻪ ﺣﻞ ﻛﻨﻴﺪ .ﺗﺒﺪﻳﻞ ﻓﻮﺭﻳﻪ ﻣﺤﺪﻭﺩ ﺳﻴﻨﻮﺳﻲ
۲۰
ﻣﺪﺭﺱ :ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
ﺗﺮﻡ ﺩﻭﻡ ۱۴۰۱-۱۴۰۲ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
∑2
n̸=1
2 2 π
= )U (x, t Fs (U (x, t)) sin nx = Fs (U (x, t)) sin x = × cos t × sin nx
π n=0 π π 1
۲۱
ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ:ﻣﺪﺭﺱ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
۱۴۰۱-۱۴۰۲ ﺗﺮﻡ ﺩﻭﻡ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
U (x, 0) = e−x
x=0 − w
= Ux sin wx|x=∞ Ux cos wxdx = −w Ux cos wxdx
0
∫ ∞ 0
∫ ∞
= −w(U cos wx|x=∞
x=0 + w U sin wxdx = −w 2
U sin wxdx
∫ ∞ ∫ ∞ 0 0
w
te−x sin wxdx = t e−x sin wxdx = t
0 0 1 + w2
d w
→ Fs (U (x, t)) = −4w2 Fs (U (x, t)) + t
dt 1+w
d w
→ Fs (U (x, t)) + 4w2 Fs (U (x, t)) = t
dt w ∫
∫ ∫
→ ﻳﺎﺩﺁﻭﺭﻱy + p(x)y = q(x) → y(x) = e
′ − p(x)dx
[ e p(x)dx
qdx + c]
∫
2 2t w
: Fs (U (x, t)) = e−4w t [ e4w tdt + c]
1+w
2
2 w e4w t (4tw2 − 1)
e−4w t [ + c]
1 + w2 16w4
−w w w w 1
Fs (U (x, t)) = 2 2
+c→c= 2
+ 4 2
= 2
(1 + )
16w (1 + w ) 1+w 16w (1 + w ) 1+w 16w4
2 2
e−4w t w e4w t (4tw2 − 1) 1
Fs (U (x, t)) = 2
( 4
+1+ 4
)∗
∫
1 + w 16w 16w
2 ∞
U (x, t) = Fs (U (x, t)) sin wxdw
π 0
۲۲
ﻣﺼﻄﻔﻲ ﻋﺒﺎﺱﺯﺍﺩﻩ:ﻣﺪﺭﺱ ﻣﻌﺎﺩﻻﺕ ﺑﺎ ﻣﺸﺘﻘﺎﺕ ﺟﺰﺋﻲ
۱۴۰۱-۱۴۰۲ ﺗﺮﻡ ﺩﻭﻡ ﺭﻳﺎﺿﻲ ﻣﻬﻨﺪﺳﻲ
ﺩﺍﻧﺸﻜﺪﻩ ﺭﻳﺎﺿﻲ ﻭ ﻋﻠﻮﻡ ﻛﺎﻣﭙﻴﻮﺗﺮ
U (x, 0) = e−x
∫ ∞
F (w) = f (x)e−iwx dx
−∞
∫ ∞
1
f (x) = F (w)eiwx dw
2π −∞
۲۳