如何在matlab中写“圆上任意两点的中点”的程序代码？？求指点呀！！！

关键字：matlab 更新时间：2019-04-11 19:46:31 209次访问

知道怎么用蒙特卡洛仿真求圆的面积吗？下面有一个pudn上的代码，使用蒙特卡洛方法求圆心在原点，半径为1的圆的面积。你可以用类似的方法解决你的问题，只要把条件改成落点同时在三个圆里面就行了
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

相关标签：matlab