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

 题目:

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

思路:随机化算法
代码:

#include
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轴边的情况 

你可能感兴趣的:(考研,c++,蓝桥杯,算法)