蒙特卡罗法计算椭圆c++

-lyslyslys

于 2022-08-24 10:30:33 发布

阅读量807

点赞数

分类专栏：考研文章标签： c++ 蓝桥杯算法

本文链接：https://blog.csdn.net/qq_43568078/article/details/126499543

版权

考研专栏收录该内容

7 篇文章 7 订阅

订阅专栏

题目：

3.已知一个椭圆的长轴长为10,，短轴长为8,请用蒙特卡罗法计算这个椭圆的面积，并编程实现。(20 分)（蒙特卡罗法）

思路：随机化算法
代码：

#include<bits/stdc++.h>
using namespace std;
int main() {
// y = 4*sqrt(1.0-nowx*nowx/25.0);根据 椭圆公式x*x/a*a + y*y/b*b == 1计算 还有y轴那种公式与这个同理 
	int sum = 0;
	for(int i = 0;i < 1000000;++i) {//测试1000000种情况 
		double nowx = double(rand()%500+1)/100.0;//忽略边 
		double nowy = double(rand()%400+1)/100.0;//忽略边 
		double y = 4.0*sqrt(1.0-nowx*nowx/25.0);
		if(nowy < y) {
			sum++;
		}
	}
//只计算1/4面积的正数部分更好计算 
	double ans = sum/50000.0;//根据在椭圆内部数/总数*实际要测试的面积 
	cout << ans << '\n';
	return 0;
}
//计算精度不够忽略了x轴y轴边的情况