www.pudn.com > EDFA-980nm.zip > afs.m, change:2014-08-28,size:1915b


function [y]=afs(x) 
%输入x单位m,输出单位m^-2 
%1.46e-6m到 1.63e-6m 
fs= ((3921849949831840191*cos(32060000*x))/217780714829400616616559748756331655331840 +... 
    (26813651224398882309*cos(64120000*x))/4355614296588012332331194975126633106636800 +... 
    (19953487362310585047*cos(96180000*x))/43556142965880123323311949751266331066368000 ... 
    - (3430081931044148631*cos(128240000*x))/21778071482940061661655974875633165533184000 + ... 
    (39378324104424401667*sin(32060000*x))/4355614296588012332331194975126633106636800 +... 
    (37173564153526251173*sin(64120000*x))/4355614296588012332331194975126633106636800 + ... 
    (561435154782614531*sin(96180000*x))/174224571863520493293247799005065324265472 + ... 
    (18781440250866586829*sin(128240000*x))/43556142965880123323311949751266331066368000 + ... 
    5454526941720145553/435561429658801233233119497512663310663680).*((x>=1.460e-6) & (x<1.51e-6))... 
+(-1.862e-024 - 3.238e-024*cos(x*6.283e+007) + 2.736e-024*sin(x*6.283e+007)... 
               - 4.881e-025*cos(2*x*6.283e+007) + 3.272e-024*sin(2*x*6.283e+007)... 
               + 1.368e-024*cos(3*x*6.283e+007) + 1.806e-024*sin(3*x*6.283e+007)... 
               + 1.232e-024*cos(4*x*6.283e+007) + 1.801e-025*sin(4*x*6.283e+007)... 
               + 3.577e-025*cos(5*x*6.283e+007) -3.149e-025*sin(5*x*6.283e+007)... 
               - 6.656e-028*cos(6*x*6.283e+007) -1.285e-025*sin(6*x*6.283e+007)).*((x>=1.51e-6) & (x<1.56e-6))... 
  + (2.218e-025 + 1.875e-025*cos(x*4.488e+007) + 2.355e-026*sin(x*4.488e+007)... 
              + 4.929e-026*cos(2*x*4.488e+007) - 2.639e-026*sin(2*x*4.488e+007)... 
              - 2.444e-026*cos(3*x*4.488e+007) -2.147e-026*sin(3*x*4.488e+007)... 
              - 2.426e-026*cos(4*x*4.488e+007) -5.078e-027*sin(4*x*4.488e+007)... 
              - 6.194e-027*cos(5*x*4.488e+007) + 6.284e-028*sin(5*x*4.488e+007)).*((x>=1.56e-6) & (x<=1.63e-6)); 
y=fs./1.0048;