clc;
close?all;
clear?all;
hold?off;
r=2;n=20;%r為指數(shù)分布參數(shù)；
temp=0;?
Un=rand（n1）;%共產(chǎn)生n*1個隨機(jī)數(shù)
for?i=1:1:n
????S（i）=-log（Un（i）/r）/r;?%Y?=?log（X）：對參量X的每個元素計算自然對數(shù)（e為底數(shù)的對數(shù)函數(shù)），使X服從指數(shù)分布；
end
S=subs（S）;%S?=?subs（S）：替換表達(dá)式S中所有的符號變量，其值通過回調(diào)函數(shù)或MATLAB工作空間獲取。
for?i=1:1:n
????for?j=1:1:i
????????temp=temp+S（j）;????
????en