ch1example3prg2.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

% 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                        % 显示耗时