function?M?=?hunhebodao（M0hHH1HH2HH3slopeElvated1slopeElvated2）
%M0=330;
%HH1=40;HH2=100;HH3=150;
c1=slopeElvated1*1e-3;c2=slopeElvated2*1e-3;
z0=1.5e-4;
for?i=1:length（h）
?????if?h（i）<=HH1
?????M（i）=?M0+0.125*h（i）-0.125*（HH1*log（（z0+h（i））/z0））;
????elseif?（h（i）>HH1?&&?h（i）<=HH2）
?????M（i）=M0+0.125*HH1-0.125*（HH1*log（（z0+HH1）/z0））+c1*（h（i）-HH1）;
????elseif?（h（i）>HH2?&&?h（i）<=HH3）
?????M（i）=M0+0.125*HH1-0.125*（HH1*log（（z0+HH1）/z0））+c1*（HH2-HH1）+c2*（h（i）-HH2）;
????else
?????M（i）=M0+0.125*HH1-0.125*（HH1*log（（z0+HH1）/z0））+c1*（HH2-HH1）+c2*（HH3-HH2）+0.157*（h（i）-HH3）;?????
????end
end
M=（M*10^（-6））+1;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????625??2017-04-14?11:06??射線追蹤法\hunhebodao.m

?????文件???????7851??2017-05-19?21:49??射線追蹤法\Main_ray.m

?????文件????????170??2016-04-27?11:25??射線追蹤法\Sub_evaporation.m

?????文件????????560??2017-01-10?19:48??射線追蹤法\Sub_FourParaRefPro.m

?????文件?????????84??2016-10-27?20:20??射線追蹤法\Sub_StanRefPro.m

?????目錄??????????0??2017-05-28?10:32??射線追蹤法

-----------?---------??----------?-----??----

?????????????????9290????????????????????6