ch1example3prg1.m

单摆运动过程的建模和仿真

% ch1example3prg1.m
sita=0:0.01:2*pi;
x=sin(sita);
y=cos(sita);% 计算半径为1的圆周上的点，以便作出圆周观察
m=0;        % 在圆内在落点计数器
x1=2*rand(1000,1)-1;% 产生均匀分布于[-1, +1]直接的两个独立随机数x1,y1
y1=2*rand(1000,1)-1;
N=1000;     % 设置试验次数
for n=1:N   % 循环进行重复试验并统计
    p1=x1(1:n);
    q1=y1(1:n);
    if (x1(n)*x1(n)+y1(n)*y1(n))<1  % 计算落点到坐标原点的距离,判别落点是否在圆内
        m=m+1;                      % 如果落入圆中,计数器加1
    end
    plot(p1,q1,'.',x,y,'-k',[-1 -1 1 1 -1],[-1 1 1 -1 -1],'-k');
    axis equal;         % 坐标纵横比例相同
    axis([-2 2 -2 2]);  % 固定坐标范围
    text(-1,-1.2,['试验总次数 n=',num2str(n)]);% 显示试验结果
    text(-1,-1.4,['落入圆中数 m=',num2str(m)]);
    text(-1,-1.6,['近似圆面积 S_c=',num2str(m/n*4)]);
    set(gcf,'DoubleBuffer','on');      % 双缓冲避免作图闪烁
    drawnow;                           % 显示结果
end