【正文】 %有效長度Pth=16*ap/gr./Leff。 % 公式plot(L,Pth) % 輸出xlabel(39。長度 L/m39。)ylabel(39。泵浦功率 /w39。)附錄2L=(1:50)*1000。afa=。ap=50*10^(12)。gr=*10^(14)。Leff=(1exp(afa*L))/afa。Pth=20*ap/gr./Leff。plot(L,Pth)xlabel(39。長度 L/m39。)ylabel(39。泵浦功率 /w39。)附錄3Pp0=1。 % 輸入泵浦光功率，單位: WPs0=。 % 輸入信號光功率，單位: WP0=[Pp0。Ps0]。 %P0指向Pp0,Ps0[t,x]=ode45(39。cormp39。, [0,20000],P0)。 %調(diào)用函數(shù)cormp ，t的長度為20公里 x=10*log10(x*1000)。 %變量x的單位由公里轉(zhuǎn)換為dBplot(t, x) %輸出x波形，有兩條，一條是對信號進行放大的泵浦光功率，一條是被放大的信號光功率。t是長度。a=。 %光纖面積ahold onplot(t,10*log10(Pp0*exp(a*t)*1000)) %畫出只損耗而不對信號光進行放大的泵浦光功率波形title(39。同向拉曼39。)xlabel(39。長度 L m39。)ylabel(39。功率 p dB39。) % 標(biāo)注“同向拉曼”，以及橫，縱坐標(biāo)減小泵浦光功率程序Pp0=。 % 輸入泵浦光功率，單位: WPs0=。 % 輸入信號光功率，單位: WP0=[Pp0。Ps0]。 %P0指向Pp0,Ps0[t,x]=ode45(39。cormp39。, [0,20000],P0)。 %調(diào)用函數(shù)cormp ，t的長度為20公里 x=10*log10(x*1000)。 %變量x的單位由公里轉(zhuǎn)換為dBplot(t, x) %輸出x波形，有兩條，一條是對信號進行放大的泵浦光功率，一條是被放大的信號光功率。t是長度。a=。 %光纖面積ahold onplot(t,10*log10(Pp0*exp(a*t)*1000)) %畫出只損耗而不對信號光進行放大的泵浦光功率波形title(39。同向拉曼39。)xlabel(39。長度 L m39。)ylabel(39。功率 p dB39。) % 標(biāo)注“同向拉曼”，以及橫，縱坐標(biāo)附錄4Pp0=1。 % 泵浦光輸入功率，單位: WPs0=。 %信號光輸入功率，單位： Wa=。 % 泵浦光面積，單位mL=30000 % 定義常量a,LPp0=Pp0*exp(a*L)。 %定義Pp0P0=[Pp0。Ps0]。 % p0的范圍是在pp0到pso之間[t,x]=ode45(39。antirmp39。,[0,L],P0)。 %函數(shù)antirmp，長度為Lx=10*log10(x*1000)。 % 把x的單位由W轉(zhuǎn)變?yōu)閐Bplot(t,x,39。r39。) % 畫出t,x 的圖形hold onplot(t,10*log10(Pp0*exp(a*t)*1000)) %畫出只有損耗不放大的泵浦 功率曲線title(39。反 向拉 曼 39。) %標(biāo)注“反向拉 曼 ”xlabel(39。長度 L m39。) %橫坐標(biāo)ylabel(39。功率 p dB39。) %縱坐標(biāo)Cormp函數(shù)function xdot=cormp(t, x) %建立常系數(shù)微分函數(shù)cormpgs=*10^(14)。 % 拉曼增益系數(shù) m/WA=42*10^(12)。 % 光纖面積，單位 m^2a=。 % 泵浦光橫截面a,單位m^(1)g=gs/A。 % 定義g=gs/Axdot=[g*x(1)*x(2)a*x(1)。g*x(1)*x(2)a*x(2)]。 % ,對照公式() Antirmp函數(shù)function xdot=rmp(t,x) % 定義常系數(shù)微分函數(shù)antirmpgs=*10^(14)。 % 定義常量，單位:m/WA=42*10^(12)。 %光纖橫截面，單位: m^2a=。 %泵浦光面積單位: m^(1)g=gs/A。 %定義常量gxdot=[a*x(1)。g*x(1)*x(2)a*x(2)]。附錄5非線性光纖 Pp0=1。 % unit: WPs0=。 % unit: WP0=[Pp0。Ps0]。[t, x]=ode45(39。cormp839。,[0,5000],P0)。x=10*log10(x*1000)。plot(t, x)a=。hold onplot(t,10*log10(Pp0*exp(a*t)*1000))title(39。同向非線性光纖39。)xlabel(39。長度 L m39。)ylabel(39。功率 p dB39。) function xdot=cormp(t, x)gs=*10^(14)。 % unit: m/WA=*10^(12)。 % unit: m^2a=。 % unit: m^(1)g=gs/A。xdot=[g*x(1)*x(2)a*x(1)。g*x(1)*x(2)a*x(2)]。附錄6Pp0=1。 % unit: WPs0=。 % unit: WP0=[Pp0。Ps0]。[t, x]=ode45(39。cormp439。,[0,80000],P0)。x=10*log10(x*1000)。plot(t, x)a=。hold onplot(t,10*log10(Pp0*exp(a*t)*1000))title(39。同向拉曼39。)xlabel(39。長度 L m39。)ylabel(39。功率 p dB39。) function xdot=cormp(t, x)gs=*10^(14)。 % unit: m/WA=50*10^(12)。 % unit: m^2a=。 % unit: m^(1)g=gs/A。xdot=[g*x(1)*x(2)a*x(1)。g*x(1)*x(2)a*x(2)]。附錄7Pp0=1。 % unit: WPs0=。 % unit: WP0=[Pp0。Ps0]。[t, x]=ode45(39。cormp539。,[0,80000],P0)。x=10*log10(x*1000)。plot(t, x)a=。hold onplot(t,10*log10(Pp0*exp(a*t)*1000))title(39。同向拉曼39。)xlabel(39。長度 L m39。)ylabel(39。功率 p dB39。) function xdot=cormp(t, x)gs=*10^(14)。 % unit: m/WA=50*10^(12)。 % unit: m^2a=。 % unit: m^(1)g=gs/A。xdot=[g*x(1)*x(2)a*x(1)。g*x(1)*x(2)a*x(2)]。附錄8色散補償光纖Pp0=1。 % unit: WPs0=。 % unit: WP0=[Pp0。Ps0]。[t, x]=ode45(39。cormp739。,[0,8000],P0)。x=10*log10(x*1000)。plot(t, x)a=。hold onplot(t,10*log10(Pp0*exp(a*t)*1000))title(39。同向色散補償光纖39。)xlabel(39。長度 L m39。)ylabel(39。功率 p dB39。) function xdot=cormp(t, x)gs=*10^(14)。 % unit: m/WA=*10^(12)。 % unit: m^2a=。 % unit: m^(1)g=gs/A。xdot=[g*x(1)*x(2)a*x(1)。g*x(1)*x(2)a*x(2)]