c9_estimatepi.m

来自「％File: c9_estimatepi.m ％有5个pi的估计」· M 代码 · 共 28 行

%File: c9_estimatepi.m
%有5个pi的估计，每一个都是基于500次重复随机试验，所得的pi的五个估计值用以下向量表示
%pi的估计值=[3.0960 3.0720 2.9920 3.1600 3.0480]
%如果对5个结果进平均，则pi的估计值=3.0736,这样的结果等价于2500次的试验结果。
%用来产生该结果的MATLAB程序如下
m=input('Enter M,number of experiments>');   %pi的估计值个数
n=input('Enter N,number of trials/experiment>'); %pi的每个估计值重复试验次数
z=zeros(1,m);         %initialize array
data=zeros(n,m);      %initialize array
for j=1:m
    x=rand(1,n);
    y=rand(1,n);
    k=0;
    for i=1:n
        if x(n)^2+y(i)^2<=1  %Fall in pie slice?
            k=k+1;
        end
        data(i,j)=4*(k/i);  %jth estimate of pi  
    end
    z(j)=data(n,j);         %Store data
end
plot(data,'k')              %Plot curves
xlabel('Number of Trials')
ylabel('estimate of pi')
%End of script file.