بهينه سازي تقويت كنندههاي فيبر نوري رامان پمپ شده به وسيله پمپهاي گوسي
بهینه سازی تقویت کننده های فیبر نوری
چكيده
: استفاده از پمپهاي باند وسيع در طراحي تقويت كنندههاي فيبر نوري رامان اخيراً از اهميت ويژهاي برخوردار شده است. ازآنجا كه ليزرهاي نيمه هادي پر قدرت باند وسيع تقريباً داراي طيف گوسي هستند، در اين مقاله براي كمينه سازي اعوجاج بهره از
تركيب چند پمپ باند وسيع استفاده ميشود
. توان و طول موج مركزي هر كدام از پمپهاي گوسي براي حالت بهينه بر اساس روش50
و nm وردشي انتخاب ميشوند. محاسبات براي بهينه سازي اعوجاج بهره يك تقويت كننده رامان به طول 25 كيلومتر، پهناي باند0/25
بدست آمد. dB 4/1 انجام شد و اعوجاج بهره dB متوسط بهره -1 مقدمهتقويت كنندههاي رامان كه بر اساس اثر رامان القائي كار
ميكنند در سيستمهاي مخابرات نوري بسيار مورد توجه
هستند
[ 1]. مزيت اين تقويت كنندهها نسبت به تقويتكنندههاي نيمه رسانا و آلاييده به اربيوم اين است كه هر
سيگنال با هر طول موجي را ميتوان به وسيلهي اين تقويت
كنندهها تقويت كرده و از خود فيبر به عنوان يك محيط
بهره استفاده ميشود
. علاوه بر اين تقويت كنندههاي راماناز نسبت سيگنال به نوفه بالايي برخوردار هستند و با
استفاده از پمپهاي مناسب ميتوان پهناي باند طيف بهره
.[
آنها را گسترده كرد[ 2در گسترش پهناي باند بهره تقويت كنندهها مهمترين
مسئله تخت كردن و كاهش اعوجاجهاي طيف بهره است
.براي انجام اين امر معمولاً از تعداد زيادي پمپ استفاده مي
-20
پمپ). البته از هر پمپ با طول موج و - شود(تعداد 4توان دلخواه نميتوان استفاده كرد و بايد توان و طول موج
پمپها را به گونهاي انتخاب كرد كه اعوجاجهاي بهره كمينه
شود
[ 3]. به منظور كم كردن تعداد پمپها و كاهش هزينهيا پمپهايي با پهناي باند گسترده
WDM ميتوان از روشاستفاده كرد
[ 4]. راه ديگر براي تخت كردن تابع بهره راماناستفاده از پمپهايي است كه طول موج آنها بر روي يك
پهناي وسيع به سرعت قابل تنظيم باشد
[ 5]. با اين وجوداين روش هنوز به صورت واقعي و عملي مورد استفاده قرار
نگرفته است
.روشهاي استفاده از چند پمپ و پمپهايي با پهناي باند
گسترده در تقويت كنندههاي رامان به روش استفاده از يك
پمپ با طيف پيوسته پيشرفت كرد
. بهترين تخت كردن تابع.[
بهره را با استفاده از اين روش ميتوان بدست آورد[ 6] و[ 7اما شكل پمپ بدست آمده در مرجع
[ 6] بسيار پيچيده بودهو بدست آوردن آن به صورت عملي مشكل است
. بنابراينتلاش براي پيدا كردن شكل پمپ پيوستهاي كه اعوجاجهاي
طيف بهره را كمينه كند و از نظر عملي قابل دسترسي
باشد، جالب به نظر ميرسد
. معمولاً از سه شيوه براي بهينهسازي استفاده ميشود
: 1- ثابت نگهداشتن توان پمپها وتغيير طول موجها،
2- ثابت نگه داشتن طول موجها وتغيير توانها و
3- تغيير طول موجها و تغيير توانها و بهينهكردن آنها
[ 7]. از آنجا كه ليزرهاي نيمه هادي باند وسيع15-
داراي طيف گوسي با پهناي باند در نيمه بيشينه 10نانومتر هستند، در اينجا فرض ميشود كه توان اوليه پمپ
به صورت تركيبي از چند پمپ گوسي با پهناي باند
10نانومتر باشد
. براي شبيه سازي پمپ پيوسته از تعداد زياديپمپ تكفام
(گسسته) استفاده مي شود[ 8]. در كارهاي قبليانجام شده ابتدا شكل منحني پمپ پيوسته بدست ميآمد و
بعد با برازش گوسي روي منحني مورد نظر، شكل پمپي كه
2
از نظر عملي قابل دسترسي باشد، بدست ميآمد
[ 8] و [ 9] وبه
SMF- 10 ]. بر اساس گزارشات قبلي، براي يك فيبر 28 ]طول
100 كيلومتر، اعوجاج بهره بدست آمده براي دو شكل0/4
است [ 9]. در صورتي dB 0/2 و dB مختلف پمپ پيوسته، كه اگر از چهار پمپ گسسته در طول موجهاي
14201453
و 1480 نانومتر استفاده كنيم اعوجاج ،1435به طول
25 SMF- 0/8 است [ 9]. براي فيبر 28 dB بهرهكيلومتر اعوجاج بهره بدست آمده با تركيب چهار پمپ
0/3
با متوسط بهره بيشتر از صفر دسي بل dB گوسياست
[ 8]. اما در اين مقاله از ابتدا فرض ميشود كه تواناوليه پمپ به صورت تركيبي از چند پمپ گوسي با پهناي
باند
10 نانومتر باشد و بعد طول موجهاي مركزي و دامنهاين پمپهاي گوسي را به گونهاي انتخاب ميكنيم كه طيف
بهره داراي كمينه اعوجاج باشد
. براي اين منظور از روشوردشي و ضرايب نامعين لاگرانژ استفاده ميشود
.-2
معادلههاي حاكممعادلههاي مربوط به انتشار توان پمپ و سيگنال با در نظر
:[
گرفتن اثر رامان به صورت زير است[ 111
(1)2
1
( )
( ) ( ) ( , ) ( )
( , ) ( ) ; 1, 2,..,
s
m
N
i
i i i j j
j
N
i j j s
j
dS z
S z g S z
dz
g Q z i N
α ν ν ν
ν ν
=
=
⎛⎜
= − + +
⎜⎝
⎞⎟
=
⎟⎠
ΣΣ
1
(2)1
( )
( ) ( ) ( , ) ( )
( , ) ( ) ; 1,2,..,
s
m
N
i i i i jj
j
N
i j j m
j
dP z
P z g S z
dz
g Pz i N
α ν ν ν
ν ν
=
=
⎛⎜
= − −
⎜⎝
⎞⎟
=
⎟⎠
Σ
Σ
توان
Pi (z) ، ام i توان سيگنال Si (z) در معادلههاي بالانمونه
i پمپ ام ، iν
فركانس سيگنال i ، ام iν
فركانستعداد
Ns ، ضريب اتلاف در محيط فيبر α ، ام i نمونهثابت بهره رامان است كه
g تعداد نمونهها و Nm ، سيگنالها:[
به صورت زير تعريف ميشود[ 12(3)
2
( );
( , ) 0 ;
( )
. ;
i ref i j
i j
j ref eff
i j i j
j ref j i
i j
reff eff
g
A
g
g
A
ν ν ν
ν ν
ν ν
νν ν ν
ν ν ν
ν ν
ν
⎧ −
⎪− >
Γ ⎪⎪
= = ⎨⎪
⎪ − ⎪ < ⎩ Γ
j
ν و iν
فركانس موج هاي بر هم كنش كننده در اثر رامانبراي به حساب آوردن اثر كاتورهاي بودن
Γ هستند. عاملسطح مقطع موثر هسته فيبر و
Aeff . جهت قطبش استثابت رامان اندازه گيري شده در فركانس
gref (ξ −η )g
ref (Δν ) است. طيف ثابت بهره رامان ν reff پمپ مرجع.[
در شكل ( 1) نمايش داده شده است[ 120 5 10 15 20 25 30 35
0
1
2
3
4
5
6
x 10
-14Pump-signal frequency difference,THz
Raman gain cofficient,m / W
[12]200
THz شكل 1: طيف ثابت بهره رامان در فركانس پمپدر معادله
( 1) جمله اول تضعيف سيگنال در طول فيبر،جمله دوم برهم كنش سيگنال
- سيگنال و جمله سوم برهمكنش نمونههاي پمپ و سيگنال را نشان ميدهد
. در معادله2)
جمله اول تضعيف نمونهاي از پمپ در طول فيبر، جمله )دوم بر هم كنش سيگنال
– نمونهاي از پمپ و جمله سوم برهم كنش نمونهاي از پمپ با نمونهاي ديگر از پمپ را نشان
ميدهد
. همان گونه كه در معادله ( 2) ديده ميشود فرضبر اين است كه عمل پمپ كردن در جهت معكوس
1 انجامرا به صورت زير تعريف ميكنيم
: F ميشود. حال تابع هدف(4)
_
2
1
( )
N
sk
k
F G G
=
=
Σ −(5)
( )
(0)
k
k
k
S L
G
S
=
) 6 (
Σ==
N
sk
k
s
G
N
G
1
_
1نشان دهنده مجموع مربعات اختلاف بهره در هر فركانس
Fام،
k بهره مربوط به سيگنال Gk ، نسبت به متوسط بهرهطول فيبر است
. در واقع L متوسط بهره تقويت كننده و G1
backward3
با توجه به قيدهاي موجود است
. از F هدف ما كمينه سازيآنجا كه توان پمپها و نمونههايي از آن كه انتخاب ميكنيم
كميتي مثبت هستند براي جلوگيري از بدست آمدن جواب
را
Pk (z) = Qk2 (z) منفي براي نمونههاي پمپ، تغيير متغيريك تابع حقيقي است
. معادله- Qk (z) انجام ميدهيم كهباز نويسي ميشود
. در انتها Qk (z) هاي مربوطه بر حسب2
( ) را به عنوان نمونهاي از توان پمپ معرفي كرده و Qk zبراي اينكه از صحيح بودن جواب اطمينان حاصل شود آن را
در معادله
( 1) و ( 2) گذاشته و منحني بهره مربوطه رسممي شود
. پس معادلههاي( 1) و ( 2) را ميتوان به صورت زيربا زنويسي كرد
:1
(7)2
1
( )
( ) ( ) ( , ) ( )
( , ) ( ) ; 1, 2,..,
s
m
N
i
i i i j j
j
N
i j j s
j
dS z
S z g S z
dz
g Q z i N
α ν ν ν
ν ν
=
=
⎛⎜
= − + +
⎜⎝
⎞⎟
=
⎟⎠
Σ
Σ
1
(8)2
1
( ) ( ) ( ) ( , ) ( )
2
( , ) ( ) ; 1,2,..,
s
m
N
i i
i i j j
j
N
i j j m
j
dQ z Q z g S z
dz
g Q z i N
α ν ν ν
ν ν
=
=
⎛
= − − ⎜⎝
⎞
= ⎟⎠
Σ
Σ
S
i (0) = Si0 ;i = 1, 2,..., Ns (9)(10)
2 2 2
1
( ) exp( ( ) ); 1,2,...,
N
pi k i k m
k
Q L E
α λ λ i N=
=
Σ − − =رابطه
( 10 ) نشاندهنده اين است كه از ابتدا فرض شده طيفپمپ گوسي
Np توان پمپ ورودي به صورت تركيبي ازباشد در صورتي در مرجع
[ 10 ] از ابتدا چنين قيدي برمسئله اعمال نميشود
. براي معادلههاي ( 7) و ( 8) به ترتيبدر نظر گرفته و تابع هدف
βk ,γ k ضرايب نامعين لاگرانژتعميم يافته را به صورت حاصل جمع تابع هدف و
(
حاصلضرب ضرايب نامعين لاگرانژ در معادلههاي ( 7) و ( 8:[
تعريف ميكنيم[ 13(11)
(
0
1
2
1 1
0
1
2
1 1
( )
( ) ( ) ( )
( , ) ( ) ( , ) ( )
( ) 1 ( ) ( ) ( )
2
( , ) ( ) ( , ) ( )
s
s m
m
s m
N L i
i i
i i
N N
i j j i j j
j j
N L i
i i i
i
N N
i j j i j j
j j
dS z
J F z S z
dz
g S z g Q z dz
dQ z
z Qz
dz
g S z g Q z dz
γ αν
ν ν ν ν
β αν
ν ν ν ν
=
= =
=
= =
⎡
= + + − ⎢⎣
⎞⎤
− ⎟⎥
⎥ ⎟⎠
⎦
⎡ ⎛
+ ⎢ + ⎜⎜− ⎣ ⎝
⎞⎤
+ + ⎟⎥
⎥ ⎟⎠
⎦
Σ ∫
Σ Σ
Σ ∫
Σ Σ
نسبت به متغيرهاي موجود وردش گرفته و آن را
J حال ازاز
δ Qi (L) مساوي با صفر قرار ميدهيم. در اينجا به جايرابطه زير استفاده ميكنيم
:(12)
2
1
2
2
( ) exp( ( ) )
( )
( )exp( ( ));
( )
1, 2,...,
N
pk
i i k
k i
k
k i k i k
i
m
E
Q L
Q L
E
Q L
i N
δ
δ αλ λ
αδλ λ λ α λ λ
=
⎡
= − − + ⎢⎢⎣
⎤
− − − − ⎥⎥⎦
=
Σ
معادلههاي حاكم بر ضرايب
δ J با مساوي صفر قرار دادن.[
نامعين لاگرانژ و شرايط مرزي مربوطه بدست ميآيد[ 13(13)
1
2
1 1
1
( )
( ) ( ) ( , ) ( )
( , ) ( ) ( , ) ( ) ( )
1 ( , ) ( ) ( ) 0; 1,2,...,
2
s
m s
m
N
i
i i i j j
j
N N
i j j j i j j
j j
N
j i j j s
j
d z
z g S z
dz
g Q z g S z z
g Q z z i N
γ
γ αν ν ν
ν ν ν ν γ
ν ν β
=
= =
=
⎛⎜
− + − −
⎜⎝
⎞⎟
−
⎟⎠
+ = =
Σ
Σ Σ
Σ
(14)
2
1
1 1
1
( ) ( )
( ) ( , ) ( )
2
( , ) ( ) ( ) ( , ) ( ) ( )
2 ( , ) ( ) ( ) 0; 1,2,...,
m
s m
s
N
i i i i j j
j
N N
i j j i j i j j
j j
N
j i j j m
j
d z z
g Q z
dz
g S z Q z g zQ z
g S z z i N
β β
α ν ν ν
ν ν ν ν β
ν ν γ
=
= =
=
⎛⎜
− + − + +
⎜⎝
⎞ ⎛
⎟ + ⎜
⎟ ⎜
⎠ ⎝
⎞
− = = ⎟⎟⎠
Σ
Σ Σ
Σ
(15)
_
1
2 2 2 ( ) ( );
(0) (0) (0)
1, 2,...,
N
si
i j i
i s i j i
s
G
L G G G
S NS S
i N
γ
=
−
= + =− −
=
Σ
β
i (0) = 0;i = 1, 2,..., Nm (16)(17)
2
1
( )
exp( ( ) ) 0;
( )
1, 2,...,
N
mi
k i k
i i
p
L
E
Q L
k N
β
α λ λ
=
− − =
=
Σ
(18)
2 2
1
( )
( )exp( ( )) 0;
( )
1, 2,...,
N
mi
k i k i k
i i
p
L
E
Q L
k N
β
λ λ α λ λ
=
− − − =
=
Σ
4
از آنجا كه مقدار اوليه سيگنالها مشخص است
0)
يعني ( ) 0 i S δ = ( 0) هيچ شرط مرزي خاصي براي ) iγ
بدست نميآيد و مقدار آن با حل معادلههاي مربوطه بدست
ميآيد
. تعداد معادلههاي ديفرانسيل بدست آمده برابر با2
است (Nm + Ns ) 2 است. تعداد شرايط مرزي (Nm + Ns )ولي همانگونه كه مي دانيم در رابطه شرايط مرزي مربوط به
2
مجهول داريم و تعداد معادلههاي Np ها به تعداد Qi (L)2
است بنابراين اين دستگاه Np جبري ( 17 ) و ( 18 ) برابر بامعادله قابل حل است
.-3
محاسبات عدديبه طول
25 SMF − براي انجام محاسبات از يك فيبر 28كيلومتر استفاده ميشود كه منحني سطح مقطع و ضريب
.[
اتلاف آن در شكل 2 نشان داده شده است[ 81400 1450 1500 1550 1600
0
0.05
0.1
Wavelength,nm
Attenuation cofficint,1/km
7080
90
Effective area,
μm2effective area
attenuation cofficient
بر
SMF − شكل 2: منحني سطح مقطع و ضريب تضعيف براي 28[
حسب طول موج[ 8در انجام محاسبات فرض مي شود كه همه سيگنالها داراي
0.1
است mW توان اوليه مساوي هستند و مقدار آن برابر با.
Sk (0) = 0.1mW = −10dBm : يعنيمعادلههايي كه با آنها سر و كار داريم يك دستگاه معادله
-هاي ديفرانسيل و جبري غير خطي به هم جفت شده
هستند
. براي حل اين معادلهها به صورت زير عمل ميكنيم:در مرحله اول براي طول موج مركزي و دامنه پمپهاي
(
گوسي يك حدس ميزنيم بنابراين با توجه به معادله ( 10خواهيم داشت
. Qi (L) ما شرط مرزي را براي مقادير نهاييمجموع تعداد معادلههاي
( 9) و ( 10 ) و ( 15 ) و ( 16 ) كه در2(
Ns + Nm) واقع شرايط مرزي مسئله هستند برابر بااست
. تعداد معادلههاي ديفرانسيل موجود هم برابر با2
است. بنابراين تعداد معادلههاي ديفرانسيل و (Ns + Nm)شرايط مرزي برابر بوده و معادلههاي ديفرانسيل مربوطه
قابل حل است
. از آنجا كه بعضي از معادلهها داراي مقداراوليه، بعضي داراي مقدار نهايي و بعضي ديگر هم داراي
شرط مرزي مرتبط با شرايط مرزي معادلههاي ديگر هستند
براي حل اين دستگاه معادله ديفرانسيل از روش تفاضل
محدود استفاده ميشود
.2
است. اين Np تعداد معادلههاي ( 17 ) و ( 18 ) برابر بامعادلهها به صورت يك دستگاه معادلات غير خطي جبري
هستند
. بعد از حل معادلههاي ديفرانسيل از الگوريتم نيوتنرافسون استفاده ميشود و با توجه به مقادير بدست آمده
ميتوان مقادير طول موجهاي مركزي و دامنه ،
βi (L) برايپمپهاي گوسي را براي گام دوم بدست آورد
. با اين مقاديربدست آمده دستگاه معادله ديفرانسيل را دوباره حل مي-
كنيم
. آن قدر اين كار را ادامه ميدهيم تا( )
i L βها ، k E ، ها kλ
هاي بدست آمده از حل Qi (L) ها ومعادلههاي ديفرانسيل در روابط
( 17 ) و ( 18 ) صدق كنند.در انجام محاسبات تعداد سيگنالها را
40 ، تعداد نمونههايتعداد پمپهاي گوسي را
4 پمپ ، (Nm = 60) پمپ را 60-1575
nm و پهناي باند سيگنالها را از ( Np = 4 )1525
انتخاب ميكنيم. مقدار پهناي پمپهاي گوسي nm10
است. nm در نيمه بيشينه با هم مساوي بوده و برابر بادر شكل
3 طيف بهره بدست آمده بعد از بهينه سازي، نشانداده شده است در اين شكل متوسط بهره بدست
0/25
است. dB 4/1 و اعوجاج بهره برابر با dB آمده1520 1530 1540 1550 1560 1570 1580 1590
2.5
3
3.5
4
4.5
signal wavelength,nm
signal gain,dB
شكل
3 : طيف بهره سيگنال بدست آمده بر حسب طول موج5
در شكل
4 منحني تغيير توان سيگنالها در طول فيبر نشانداده شده است
. در اين شكل تقويت سيگنالها و در جهتعكس بودن پمپها به وضوح قابل مشاهده است
.0
10
20
30
1520
1540
1560
1580
0.05
0.1
0.15
0.2
0.25
0.3
signal wavelength,nm fiber length,km
signal power,mW
شكل
4 : تغيير توان سيگنالها بر حسب طول فيبردر شكل
5 طيف توان ورودي پمپ بهينه شده نشان دادهشده است
. طول موج مركزي پمپهاي گوسي برابر با1484
نانومتر است. / 1460/2 و 3 ، 1438/2 ،1421/61400 1420 1440 1460 1480 1500
0
0.2
0.4
0.6
0.8
1
pump power,nm
pump power,arb.units
شكل
5 : طيف توان پمپ ورودي656
است. روش mW در شكل 5 توان كل پمپ ورودي برابروردشي در كارهاي قبلي براي بهينه سازي تقويت كننده-
هاي رامان پمپ شده به وسيله پمپهاي تكفام
[ 13 ] وپمپهاي پيوسته با برازش گوسي مورد استفاده قرار گرفته
است
[ 10 ]. همان گونه كه در اين مقاله نشان داده شد برايبهينه سازي تقويت كنندههاي رامان پمپ شده به وسيله
پمپهاي گوسي نيز ميتوان از اين روش استفاده كرد ولي
بر خلاف روشهاي قبلي، از ابتدا فرض ميشود كه شكل
پمپ گوسي باشد ولي فركانسهاي مركزي و توان پمپها به
منظور كمينه سازي اعوجاج بهره بر اساس روش وردشي
تعيين ميشود و در اينجا از هيچ گونه برازش گوسي
استفاده نميشود
. اشكال برازش گوسي اين است كه ممكناست پهنا در نيمه بيشينه پمپهاي گوسي بدست آمده از
.[
مقدار آن در ليزرهاي نيمه هادي معمولي كمتر باشد[ 10همچنين در اين روش ميتوان قيدهاي مختلفي روي توان
پمپها، مجموع توان پمپها و بهره متوسط تقويت كننده
اعمال كرد
.-4
نتيجه گيريوردش روشي دقيق و كارامد در بهينه سازي تقويت كننده-
هاي رامان با هرگونه قيد است
. در اين مقاله نشان داده شدكه از اين روش ميتوان براي پيدا كردن شكل و مشخصهي
پمپهاي گوسي لازم براي كمينه كردن اعوجاج بهره
استفاده كرد
. نه تنها براي پيدا كردن شكل پمپهاي گوسيبلكه براي پيدا كردن پمپهايي با هرگونه شكل دلخواه مي-
توان از اين روش استفاده نمود
. در اين روش هيچ گونهمحدوديتي براي اعمال هر گونه قيد دلخواه بر روي مجموع
توان پمپ ها و متوسط بهره تقويت كننده وجود ندارد
. .