%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%蒙特卡洛仿真
%N：仿真實(shí)驗(yàn)光子數(shù)
%L：仿真基線長(zhǎng)度
%thgma：消光系數(shù)
%theta_trans：發(fā)射機(jī)發(fā)散角
%theta_rthv：接收機(jī)視場(chǎng)角
%g：非對(duì)稱因子
%PP：散射光子數(shù)占接受光子數(shù)的百分比
%num_of_scat：限定最大散射次數(shù)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function?[PPT]=MonteCarlo2（NLthgmatheta_transtheta_rthvgdnum_of_scat）%


num_of_photon?=?0;%光子計(jì)數(shù)器初始化
num_of_scath?=?zeros（1num_of_scat+1）;%多次散射到達(dá)接收機(jī)的光子數(shù)

for?i=1:1:N
????
????P=zeros（13num_of_scat+2）;%存放每個(gè)散射點(diǎn)位置
????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）;%第一個(gè)散射點(diǎn)位置
????%test_2=acos（P（112）/sqrt（P（112）