Fig2_叶绿素浓度和散射吸收系数的关系图_源码

%本程序的计算需参考以下文献： %1、“Chlorophyll-based model of seawater optical properties” %2、“An algorithm to restore spectral signatures of all inherent optical properties of seawater using a value of one property at one wavelength” %3、“Absorption spectrum：380–700 nm of pure water II. Integrating cavity measurements” %4、“An optical classification of coastal and oceanic waters based on the specific spectral absorption curves of phytoplankton pigments, dissolved organic matter, and other particulate materials” clear all; clc C=[0.31,0.83]; %叶绿素浓度：mg/m^3 lamda=[400:5:690]; %波长：nm %aw和ac这两个因子，似乎都是根据文献4中的数据计算所得 aw=[0.018,0.018,0.017,0.017,0.016,0.016,0.015,0.015,0.015,0.015,0.015,0.016,0.016,0.016,0.016,0.017,0.018,0.019,0.020,0.023,... 0.026,0.031,0.036,0.042,0.048,0.050,0.051,0.054,0.056,0.060,0.064,0.068,0.071,0.076,0.080,0.094,0.108,0.133,0.157,0.201,... 0.245,0.268,0.290,0.300,0.310,0.315,0.320,0.325,0.330,0.340,0.350,0.380,0.410,0.420,0.430,0.440,0.450,0.475,0.500]; %纯水吸收系数，波长：[400:5:690] ac=[0.687,0.781,0.828,0.883,0.913,0.939,0.973,1.001,1.000,0.971,0.944,0.928,0.917,0.902,0.870,0.839,0.798,0.773,0.750,0.717,... 0.668,0.645,0.618,0.582,0.528,0.504,0.474,0.444,0.416,0.384,0.357,0.321,0.294,0.273,0.276,0.268,0.291,0.274,0.282,0.249,... 0.236,0.279,0.252,0.268,0.276,0.299,0.317,0.333,0.334,0.326,0.356,0.389,0.441,0.534,0.595,0.544,0.502,0.420,0.329]; %叶绿素特别吸收系数,波长：[400:5:690] af=35.959; %黄腐酸特别吸收系数：m^2/mg kf=0.0189; %nm^-1 ah=18.828; %腐殖酸特别吸收系数：m^2/mg kh=0.01105; %nm^-1 Cf=1.74098*C.*exp(0.12327*C); %黄腐酸浓度：mg/m^3 Ch=0.19334*C.*exp(0.12343*C); %腐殖酸浓度：mg/m^3 Cs=0.01739*C.*exp(0.11631*C); %小粒子浓度：g/m^3 Cl=0.76284*C.*exp(0.03092*C); %大粒子浓度：g/m^3 bw=0.005826*(400./lamda).^4.322; %纯水散射系数：m^-1 bs=1.151302*(400./lamda).^1.7; %小粒子物质特别散射系数：m^2/g bl=0.341074*(400./lamda).^0.3; %大粒子物质特别散射系数：m^2/g for ii=1:length(C) a(ii,:)=aw+ac*0.06*C(ii)^0.602+af*Cf(ii)*exp(-kf*lamda)+ah*Ch(ii)*exp(-kh*lamda); %吸收系数：文献1和2中计算吸收系数的公式第二项中少乘了一个0.06的因子，但文献2附录所给Fortran程序中有，故修正 b(ii,:)=bw+bs*Cs(ii)+bl*Cl(ii); %反射系数 c=a+b; %损耗系数 end plot(lamda,a(1,:),'bs--',lamda,b(1,:),'ro--',lamda,c(1,:),'kd--','MarkerSize',4);hold on; grid on;set(gca,'GridLineStyle',':','GridAlpha',1); plot(lamda,a(2,:),'bs-','MarkerFaceColor','b','MarkerSize',4); plot(lamda,b(2,:),'ro-','MarkerFaceColor','r','MarkerSize',4); plot(lamda,c(2,:),'kd-','MarkerFaceColor','k','MarkerSize',4); legend('a, C=0.31','b, C=0.31','c, C=0.31','a, C=0.83','b, C=0.83','c, C=0.83','Location','northwest'); xlabel('\lambda(nm)'); ylabel('Coef(m^{-1})');