%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%蒙特卡洛仿真
%N：仿真實驗光子數
%L：仿真基線長度
%thgma：消光系數
%theta_trans：發射機發散角
%theta_rthv：接收機視場角
%g：非對稱因子
%PP：散射光子數占接受光子數的百分比
%num_of_scat：限定最大散射次數
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function?[PPT]=MonteCarlo2（NLthgmatheta_transtheta_rthvgdnum_of_scat）%


num_of_photon?=?0;%光子計數器初始化
num_of_scath?=?zeros（1num_of_scat+1）;%多次散射到達接收機的光子數

for?i=1:1:N
????
????P=zeros（13num_of_scat+2）;%存放每個散射點位置
????D=zeros（13num_of_scat+1）;
????P（::1）?=?[000];%初始位置
????theta_0?=?acos（1?-?rand*（1-cos（theta_trans）））;%初始條件抽樣
????phi_0?=?2*rand*pi;
????D（::1）?=?[sin（theta_0）*cos（phi_0）sin（theta_0）*sin（phi_0）cos（theta_0）];%初始方向

????l?=?-（1/thgma）*log（rand）;%自由路程抽樣
????P（::2）?=?P（::1）+l.*D（::1）;%第一個散射點位置
????%test_2=acos（P（112）/sqrt（P（112）