function?dy=light_amplification（zy）
%?該函數的定義是為了調用ode函數計算微分方程
%%?光纖放大器的參數
Vp=3e8/980e-9;%泵浦光波長lambda=980nm??c/lambda?
sigma_p=2.0E-25;%波長為980nm時吸收截面面積
Vs=3e8/1531e-9;%信號光頻率??
sigma_as=7.2E-25;sigma_es=6.6E-25;%信號光的吸收和發射截面面積
Tp=0.6;Ts=0.6;%重疊因子
N0=2.0E24;%摻雜粒子數密度
Ac=1.26E-11;%光纖截面面積
t=1e-3;????%壽命
global?flag
if?flag==1
????Vp=3e8/1480e-9;%泵浦光波長lambda=1480nm
????sigma_p=7E-25;?%波長為1480nm時吸收截面面積
end
h=?6.626E-34;?%普朗克常數
%%?建方程
dy=zeros（21）;
N1=（Ac*h+t*Ts*sigma_es*y（2）/Vs）/（Ac*h+t*Tp*sigma_p*y（1）/Vp+t*Ts*sigma_as*y（2）/Vs+t*Ts*sigma_es*y（2）/Vs）*N0;?%下能級粒子數密度
%N1?=?N0*（1+Es*t*Ts*y（2）/（h*Vs））/（1+t*Tp*Ap*y（1）/（h*Vp）+t*Ts*As*y（2）/（h*Vs）+t*Ts*Es*y（2）/（h*Vs））
N2=N0-N1;%上能級粒子數密度
dy（1）=-Tp*sigma_p*N1*y（1）;%泵浦光功率微分方程
dy（2）=Ts*（sigma_es*N2-sigma_as*N1）*y（2）;%信號光功率微分方程

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????177074??2020-05-05?08:43??簡單光纖放大matlab模擬\980nm和1_480nm泵浦L波段EDFA的優化設計_關衛華.pdf
?????文件?????????926??2020-05-05?12:00??簡單光纖放大matlab模擬\light_amplification.m
?????文件?????????348??2020-05-05?12:02??簡單光纖放大matlab模擬\outP_information.m
?????文件?????????319??2020-05-05?12:02??簡單光纖放大matlab模擬\outP_pumping.m
?????文件????????2623??2020-05-05?12:29??簡單光纖放大matlab模擬\testmain.m
?????文件??????558843??2020-05-05?09:01??簡單光纖放大matlab模擬\摻鉺玻璃波導放大器中抽運光信號光重疊因子分析_顧浩然.pdf
?????文件??????421131??2020-05-05?08:52??簡單光纖放大matlab模擬\摻鉺光纖放大器的最佳光纖長度和增益特性_孫軍強.pdf
?????文件??????404089??2020-05-05?08:55??簡單光纖放大matlab模擬\抽運光與信號光的光強重疊因子和摻鉺玻璃波導放大器的增益特性_邵公望.pdf
?????文件??????504384??2020-05-05?08:48??簡單光纖放大matlab模擬\基于摻鉺光纖的高增益低噪聲光放大器實驗研究_張珂衛.pdf