ch1example3prg2.m

来自「一些有用的MATLAB通信信号处理的源代码程序」· M 代码 · 共 10 行

% ch1example3prg2.m
tic                             % 计时器启动
n=10000;                        % 每次随机落点10000个
for k=1:1000                    % 重复试验1000次
    x1=2*rand(n,1)-1;
    y1=2*rand(n,1)-1;           % 随机落点产生
    m(k)=sum((x1.*x1+y1.*y1)<1);% 求落入圆中的点数和
end
S_c=mean(m).*4./n               % 计算并显示结果
time=toc                        % 显示耗时