تابع گرین
تابع گرین به انگلیسی ( Green's function ) مفهومی است که اولین بار در دهه ۱۸۳۰ توسط جرج گرین ریاضی دان انگلیسی مطرح شد. بهطور کلی تابع گرین هسته انتگرال است که برای حل معادلات دیفرانسیل از جمله معادلات دیفرانسیل معمولی با شرایط اولیه یا مرزی و هم چنین مسائل دشوارتر مانند معادلات دیفرانسیل با مشتقات جزئی غیرهمگن با شرایط مرزی مورد استفاده قرار میگیرد.
بهطور دقیق فرض کنید عملگر دیفرانسیلی خطیL=L(x) که روی مجموعه ای از توزیع ها روی
D
⊆
R
n
{\displaystyle D\subseteq {R^{n}}}
عمل میکند داده شده است.یک تابع گرین در نقطه s ∈ Dمتناظر با عملگر L ، هر جواب معادله زیر است:
L
G
(
x
,
s
)
=
δ
(
x
−
s
)
{\displaystyle LG(x,s)=\delta (x-s)}
که
δ
{\displaystyle \delta }
تابع دلتای دیراک است.در واقع G تابعی است که با اثر کردن عملگر دیفرانسیلی L روی آن تابع دلتای دیراک حاصل می شود.
با ضرب طرفین این رابطه در f(s) و انتگرال گیری خواهیم داشت:
∫
L
G
(
x
,
s
)
f
(
s
)
d
s
=
∫
δ
(
x
−
s
)
f
(
s
)
d
s
{\displaystyle \int LG(x,s)f(s)ds=\int \delta (x-s)f(s)ds}
با استفاده از خواص تابع دلتای دیراک میدانیم عبارت سمت راست رابطه فوق برابر f(x) است، از طرفی چون عملگر L خطی است و فقط روی x اثر میکند عبارت سمت چپ رابطه فوق را میتوانیم به صورت زیر بنویسیم:
L
=
(
∫
G
(
x
,
s
)
f
(
s
)
d
s
)
{\displaystyle L=(\int G(x,s)f(s)ds)}
این روند برای حل معادلات دیفرانسیل زیر که u=u(x) بسیار مفید است:
L
u
(
x
)
=
f
(
x
)
{\displaystyle Lu(x)=f(x)}
⇒
L
u
(
x
)
=
L
(
∫
G
(
x
,
s
)
f
(
s
)
d
x
)
{\displaystyle Lu(x)=L(\int G(x,s)f(s)dx)}
u
(
x
)
=
∫
G
(
x
,
s
)
f
(
s
)
d
s
{\displaystyle u(x)=\int G(x,s)f(s)ds}
خواص تابع گرین در بعد یک[ ویرایش ]
تابع گرین مسئله اشتورم- لیوویل هموار پیوسته است یعنی:
[
G
(
x
,
s
)
]
x
=
s
−
x
=
s
+
=
0
{\displaystyle [G(x,s)]_{x=s^{-}}^{x=s^{+}}=0}
⇒
G
(
x
,
s
−
)
=
G
(
x
,
s
+
)
{\displaystyle G(x,s^{-})=G(x,s^{+})}
[
d
d
x
G
(
x
,
s
)
]
x
=
s
−
x
=
s
+
=
1
p
(
s
)
{\displaystyle [{\frac {\mathrm {d} }{\mathrm {dx} }}G(x,s)]_{x=s^{-}}^{x=s^{+}}={\frac {\mathrm {1} }{\mathrm {p(s)} }}}
به کمک این خواص میتوان تابع گرین مربوط به یک مسئله اشتورم- لیوویل هموار را مشخص کرد.
تابع گرین دو بعدی و سه بعدی[ ویرایش ]
اگر در بعد دو و سه دامنه مورد نظر ناحیه D باشد آنگاه شرایط مرزی روی مرز D یعنی
∂
D
{\displaystyle \partial D}
تعریف میشود.
فرض میکنیم D کراندار باشد:
مشابه قسمت یک بعدی میتوانیم تابع گرین دو بعدی را نیز به صورت زیر تعریف کنیم:
∇
2
G
(
x
,
s
)
=
δ
(
x
−
s
)
{\displaystyle {\nabla }^{2}G(x,s)=\delta (x-s)}
با شرایط مرزی
a
G
(
x
,
s
)
+
b
n
.
∇
x
G
(
x
,
s
)
=
0
{\displaystyle aG(x,s)+bn.{\nabla }_{x}G(x,s)=0}
که
∇
{\displaystyle \nabla }
تابع گرادیان است.
برای تعیین تابع گرین ، گرادیان را حول نقاطی بررسی میکنیم که به s نزدیک اند:
x
=
s
+
r
{\displaystyle x=s+r}
r
=
{
r
(
c
o
s
θ
,
s
i
n
θ
)
,
i
n
2
D
,
r
(
c
o
s
θ
s
i
n
φ
,
s
i
n
θ
s
i
n
φ
,
c
o
s
φ
)
,
i
n
3
D
.
{\displaystyle r={\begin{cases}r(cos\theta ,sin\theta ),&in2D,\\r(cos\theta sin\varphi ,sin\theta sin\varphi ,cos\varphi ),&in3D.\end{cases}}}
که
0
≤
θ
≤
2
π
{\displaystyle 0\leq \theta \leq 2\pi }
و
0
≤
φ
≤
π
{\displaystyle 0\leq \varphi \leq \pi }
S
β
{\displaystyle S_{\beta }}
رادایره ای به شعاع
β
{\displaystyle \beta }
در نظر میگیریم:
S
β
=
{\displaystyle S_{\beta }=}
{
x
:
|
x
−
s
|
=
r
≤
β
{\displaystyle {x:|x-s|=r\leq \beta }}
}
x
:
|
x
−
s
|
=
r
≤
β
{\displaystyle {x:|x-s|=r\leq \beta }}
∇
2
G
(
x
,
s
)
=
δ
(
x
,
s
)
{\displaystyle {\nabla }^{2}G(x,s)=\delta (x,s)}
∫
S
β
∇
2
G
(
x
,
s
)
d
x
=
∫
S
β
δ
(
x
,
s
)
d
x
{\displaystyle \int _{S_{\beta }}{\nabla }^{2}G(x,s)dx=\int _{S_{\beta }}\delta (x,s)dx}
طبق خاصیت تابع دلتای دیراک عبارت سمت راست رابطه فوق برابر یک است،از طرفی با توجه به قضیه دیورژانس خواهیم داشت :
∫
∂
S
β
δ
(
x
,
s
)
.
n
d
l
=
1
{\displaystyle \int _{\partial S_{\beta }}\delta (x,s).ndl=1}
⇒
∫
∂
S
β
∂
G
∂
n
(
x
,
s
)
d
l
=
1
{\displaystyle \int _{\partial S_{\beta }}{\frac {\mathrm {\partial G} }{\mathrm {\partial n} }}(x,s)dl=1}
که
n
=
e
r
{\displaystyle n=e_{r}}
بردار یکه در راستای شعاع و
d
l
=
r
d
θ
{\displaystyle dl=rd\theta }
است .لذا داریم :
lim
ϵ
→
0
∫
0
2
π
∂
G
∂
r
r
|
r
=
ϵ
d
θ
=
1
{\displaystyle \lim _{\epsilon \to 0}\int _{0}^{2\pi }{{\frac {\mathrm {\partial G} }{\mathrm {\partial r} }}r|_{r=\epsilon }d\theta =1}}
با محاسبه انتگرال فوق تابع گرین بدست خواهد آمد که به فرم زیر است :
G
(
x
,
s
)
=
1
2
π
l
n
|
x
−
s
|
{\displaystyle G(x,s)={\frac {\mathrm {1} }{\mathrm {2\pi } }}ln|x-s|}
با تکرار همین روند برای حالت سه بعدی داریم :
n
=
e
r
{\displaystyle n=e_{r}}
,
d
s
=
r
2
s
i
n
(
φ
)
d
θ
{\displaystyle ds=r^{2}sin(\varphi )d\theta }
:
lim
ϵ
→
0
∫
0
π
∫
0
2
π
∂
G
∂
r
r
2
|
r
=
ϵ
s
i
n
φ
d
θ
d
φ
=
1
{\displaystyle \lim _{\epsilon \to 0}\int _{0}^{\pi }\int _{0}^{2\pi }{{\frac {\mathrm {\partial G} }{\mathrm {\partial r} }}r^{2}|_{r=\epsilon }sin\varphi d\theta d\varphi =1}}
با محاسبه عبارت فوق تابع گرین سه بعدی در ناحیه بی کران به فرم زیر خواهد شد :
G
(
x
,
s
)
=
−
1
4
π
|
x
−
s
|
{\displaystyle G(x,s)={\frac {\mathrm {-1} }{\mathrm {4\pi |x-s|} }}}
روش تعیین پاسخ صریح برای تابع گرین[ ویرایش ]
1.تعیین دو جواب مستقل خطی از معادله همگن به نام
u
R
(
x
)
,
u
L
(
x
)
{\displaystyle u_{R}(x),u_{L}(x)}
بهطوری که
u
L
(
x
)
{\displaystyle u_{L}(x)}
در شرط مرزی مربوط به نقطه a و
u
R
(
x
)
{\displaystyle u_{R}(x)}
هم در شرط مرزی همگن مربوط به نقطه b صدق کند.
2.تابع گرین را به فرم زیر بنویسید ;
G
(
x
,
s
)
=
{
c
L
(
s
)
u
L
(
x
)
,
a
≤
x
≤
s
,
c
R
(
s
)
u
R
(
x
)
,
s
≤
x
≤
L
.
{\displaystyle G(x,s)={\begin{cases}c_{L}(s)u_{L}(x),&a\leq x\leq s,\\c_{R}(s)u_{R}(x),&s\leq x\leq L.\end{cases}}}
3.شرط پیوستگی
G(x,s)را در x=s اعمال کنید;
c
L
(
s
)
u
L
(
s
)
=
c
R
(
s
)
u
R
(
s
)
{\displaystyle c_{L}(s)u_{L}(s)=c_{R}(s)u_{R}(s)}
4.شرط پرش
d
d
x
G
(
x
,
s
)
{\displaystyle {\frac {\mathrm {d} }{\mathrm {dx} }}G(x,s)}
را در x=s اعمال کنید;
c
R
(
s
)
d
u
R
d
x
(
s
)
−
c
L
(
s
)
d
u
L
d
x
(
s
)
=
1
p
(
s
)
{\displaystyle c_{R}(s){\frac {\mathrm {du_{R}} }{\mathrm {dx} }}(s)-c_{L}(s){\frac {\mathrm {du_{L}} }{\mathrm {dx} }}(s)={\frac {\mathrm {1} }{\mathrm {p(s)} }}}
5.با حل روابط در مرحله 3و 4 مقادیر
c
L
(
s
)
,
c
R
(
s
)
{\displaystyle c_{L}(s),c_{R}(s)}
بدست خواهند آمد که به صورت زیر اند;
c
L
(
s
)
=
u
R
(
s
)
p
(
s
)
W
(
s
)
{\displaystyle c_{L}(s)={\frac {\mathrm {u_{R}(s)} }{\mathrm {p(s)W(s)} }}}
,
c
R
(
s
)
=
u
L
(
s
)
p
(
s
)
W
(
s
)
{\displaystyle c_{R}(s)={\frac {\mathrm {u_{L}(s)} }{\mathrm {p(s)W(s)} }}}
که W(s) رونسکین جواب های معادله همگن در نقطه s است.
معادله مستقل از زمان گرما را در نظر بگیرید.
d
2
u
(
x
)
d
x
2
=
f
(
x
)
{\displaystyle {\frac {\mathrm {d^{2}u(x)} }{\mathrm {d} x^{2}}}=f(x)}
با شرایط مرزی
u
(
0
)
=
u
(
L
)
=
0
{\displaystyle u(0)=u(L)=0}
.ما جواب این مسئله را میتوانیم به فرم زیر بنویسیم:
∫
0
L
f
(
s
)
G
(
x
,
s
)
d
s
{\displaystyle \int _{0}^{L}f(s)G(x,s)ds}
که تابع گرین در رابطه زیر صدق میکند:
d
2
G
(
x
,
s
)
d
x
2
=
δ
(
x
−
s
)
{\displaystyle {\frac {\mathrm {d^{2}G(x,s)} }{\mathrm {d} x^{2}}}=\delta (x-s)}
با شرایط مرزی
G
(
0
,
s
)
=
0
=
G
(
L
,
s
)
{\displaystyle G(0,s)=0=G(L,s)}
. حال با توجه به تعریف تابع دلتای دیراک در
x
=
s
{\displaystyle x=s}
داریم:
G
(
x
,
s
)
=
{
a
+
b
x
,
x
<
s
,
c
+
d
x
,
s
<
x
.
{\displaystyle G(x,s)={\begin{cases}a+bx,&x<s,\\c+dx,&s<x.\end{cases}}}
با اعمال شرایط مرزی مسئله خواهییم داشت
a
=
c
=
0
{\displaystyle a=c=0}
و تابع گرین به فرم زیر خواهد شد:
G
(
x
,
s
)
=
{
b
x
,
x
<
s
,
d
(
x
−
L
)
,
s
<
x
.
{\displaystyle G(x,s)={\begin{cases}bx,&x<s,\\d(x-L),&s<x.\end{cases}}}
با توجه به خاصیت پیوستگی تابع G در نقطه
x
=
s
{\displaystyle x=s}
داریم:
b
s
=
d
(
s
−
L
)
{\displaystyle bs=d(s-L)}
همچنین با شرط پرش روی مشتق در نقطه ی
x
=
s
{\displaystyle x=s}
داریم:
d
−
b
=
1
{\displaystyle d-b=1}
با حل دستگاه مقادیر b,d بدست خواهند آمد;
b
=
(
s
−
L
)
L
{\displaystyle b={\frac {\mathrm {(s-L)} }{\mathrm {L} }}}
و
d
=
s
L
{\displaystyle d={\frac {\mathrm {s} }{\mathrm {L} }}}
با جایگذاری مقادیر
b
{\displaystyle b}
و
d
{\displaystyle d}
در معادله داریم :
G
(
x
,
s
)
=
{
x
L
(
s
−
L
)
,
0
≤
x
≤
s
,
s
L
(
x
−
L
)
,
s
≤
x
≤
L
.
{\displaystyle G(x,s)={\begin{cases}{\frac {\mathrm {x} }{\mathrm {L} }}(s-L),&0\leq x\leq s,\\{\frac {\mathrm {s} }{\mathrm {L} }}(x-L),&s\leq x\leq L.\end{cases}}}