利用Python求阴影部分的面积实例代码


Posted in Python onDecember 05, 2018

一、前言说明

今天看到微信群里一道六年级数学题,如下图,求阴影部分面积

利用Python求阴影部分的面积实例代码

看起来似乎并不是很难,可是博主添加各种辅助线,写各种方法都没出来,不得已而改用写Python代码来求面积了

二、思路介绍

1.用Python将上图画在坐标轴上,主要是斜线函数和半圆函数

利用Python求阴影部分的面积实例代码

2.均匀的在长方形上面洒满豆子(假设是豆子),求阴影部分豆子占比*总面积

利用Python求阴影部分的面积实例代码

三、源码设计

1.做图源码

import matplotlib.pyplot as plt
import numpy as np


def init():
 plt.xlabel('X')
 plt.ylabel('Y')

 fig = plt.gcf()
 fig.set_facecolor('lightyellow')
 fig.set_edgecolor("black")

 ax = plt.gca()
 ax.patch.set_facecolor("lightgray") # 设置ax区域背景颜色    
 ax.patch.set_alpha(0.1) # 设置ax区域背景颜色透明度 
 ax.spines['right'].set_color('none')
 ax.spines['top'].set_color('none')
 ax.xaxis.set_ticks_position('bottom')
 ax.yaxis.set_ticks_position('left')
 ax.spines['bottom'].set_position(('data', 0))
 ax.spines['left'].set_position(('data', 0))


# 原下半函数
def f1(px, r, a, b):
 return b - np.sqrt(r**2 - (px - a)**2)


# 斜线函数
def f2(px, m, n):
 return px*n/m


# 斜线函数2
def f3(px, m, n):
 return n-1*px*n/m


if __name__ == '__main__':
 r = 4 # 圆半径
 m = 8 # 宽
 n = 4 # 高
 a, b = (4, 4) # 圆心坐标
 init()

 x = np.linspace(0, m, 100*m)
 y = np.linspace(0, n, 100*n)

 # 半圆形
 y1 = f1(x, r, a, b)
 plt.plot(x, y1)
 # 矩形横线
 plt.plot((x.min(), x.max()), (y.min(), y.min()), 'g')
 plt.plot((x.min(), x.max()), (y.max(), y.max()), 'g')
 plt.plot((x.max(), x.max()), (y.max()+2, y.max()+2), 'g') # 画点(8,6)避免图形变形
 # 矩形纵向
 plt.plot((x.min(), x.min()), (y.min(), y.max()), 'g')
 plt.plot((x.max(), x.max()), (y.min(), y.max()), 'g')
 # 斜线方法
 y2 = f2(x, m, n)
 plt.plot(x, y2, 'purple')

 # 阴影部分填充
 xf = x[np.where(x <= 0.5*x.max())]
 plt.fill_between(xf, y.min(), f1(xf, r, a, b), where=f1(xf, r, a, b) <= f2(xf, m, n),
      facecolor='y', interpolate=True)
 plt.fill_between(xf, y.min(), f2(xf, m, n), where=f1(xf, r, a, b) > f2(xf, m, n),
      facecolor='y', interpolate=True)
 # 半圆填充
 plt.fill_between(x, y1, y.max(), facecolor='r', alpha=0.25)
 plt.show()

Draw.py

2.计算源码,其中side是要不要计算图形边框上的点,理论上side只能为True;t设置越大运行时间越长也越精准

import numpy as np


def f1(px, r, a, b):
 return b - np.sqrt(r**2 - (px - a)**2)


def f2(px, m, n):
 return px*n/m


if __name__ == '__main__':
 r = 4 # 圆半径
 m = 8 # 宽
 n = 4 # 高
 a, b = (4, 4) # 圆心坐标
 t = 100 # 精度

 xs = np.linspace(0, m, 2*t*m)
 ys = np.linspace(0, n, t*n)

 # 半圆形
 y1 = f1(xs, r, a, b)
 # 斜线
 y2 = f2(xs, m, n)

 numin = 0
 numtotel = 0
 side = True # 是否计算边框
 for x in xs:
  for y in ys:
   if not side:
    if (x <= 0) | (x >= 8) | (y <= 0) | (y >= 4):
     continue
   numtotel += 1
   if x >= 4:
    continue
   y1 = f1(x, r, a, b)
   y2 = f2(x, m, n)
   if y1 - y2 >= 0:
    if y2 - y > 0:
     numin += 1
    if (y2 - y == 0) and side:
     numin += 1
   elif y2 - y1 > 0:
    if y1 - y > 0:
     numin += 1
    if (y2 - y == 0) and side:
     numin += 1

 print(32*numin/numtotel)

calc.py

四、最后小结

1.此种算法t为100时,阴影面积为1.268;t为1000时,阴影面积为1.253,已经非常接近正确答案(正确答案1.252)

2.举一反三,类似于这种不规则的面积,只要可以写出来函数,就可以求解面积.

2.下面有三种求解方法,第三种表示比大学高数还难看懂,你们呢?

利用Python求阴影部分的面积实例代码

利用Python求阴影部分的面积实例代码

利用Python求阴影部分的面积实例代码

总结

以上就是这篇文章的全部内容了,希望本文的内容对大家的学习或者工作具有一定的参考学习价值,如果有疑问大家可以留言交流,谢谢大家对三水点靠木的支持。

Python 相关文章推荐
Python类的多重继承问题深入分析
Nov 09 Python
用Python编写生成树状结构的文件目录的脚本的教程
May 04 Python
Python时间模块datetime、time、calendar的使用方法
Jan 13 Python
深入理解python中的select模块
Apr 23 Python
Python实现比较扑克牌大小程序代码示例
Dec 06 Python
Python实现生成随机数据插入mysql数据库的方法
Dec 25 Python
python Web开发你要理解的WSGI &amp; uwsgi详解
Aug 01 Python
python 多个参数不为空校验方法
Feb 14 Python
解决Python中定时任务线程无法自动退出的问题
Feb 18 Python
python循环定时中断执行某一段程序的实例
Jun 29 Python
python tkinter图形界面代码统计工具(更新)
Sep 18 Python
使用anaconda安装pytorch的实现步骤
Sep 03 Python
python之cv2与图像的载入、显示和保存实例
Dec 05 #Python
python存储16bit和32bit图像的实例
Dec 05 #Python
Python随机生成身份证号码及校验功能
Dec 04 #Python
python 从文件夹抽取图片另存的方法
Dec 04 #Python
使用Python实现微信提醒备忘录功能
Dec 04 #Python
flask-restful使用总结
Dec 04 #Python
Python读取YUV文件,并显示的方法
Dec 04 #Python
You might like
mysql5的sql文件导入到mysql4的方法
2008/10/19 PHP
驱动事件的addEvent.js代码
2007/03/27 Javascript
javascript 一个函数对同一元素的多个事件响应
2009/07/25 Javascript
JSON 数字排序多字段排序介绍
2013/09/18 Javascript
javascript四舍五入函数代码分享(保留后几位)
2013/12/10 Javascript
jquery获取URL中参数解决中文乱码问题的两种方法
2013/12/18 Javascript
JavaScript中的原型和继承详解(图文)
2014/07/18 Javascript
Jquery Mobile 自定义按钮图标
2015/11/18 Javascript
AngularJS单选框及多选框实现双向动态绑定
2016/01/13 Javascript
Bootstrap 实现查询的完美方法
2016/10/26 Javascript
JS比较两个数值的大小实例
2016/11/25 Javascript
jQuery Validate验证框架详解(推荐)
2016/12/17 Javascript
详解AngularJS通过ocLazyLoad实现动态(懒)加载模块和依赖
2017/03/01 Javascript
ES6解构赋值的功能与用途实例分析
2017/10/31 Javascript
10 种最常见的 Javascript 错误(频率最高)
2018/02/08 Javascript
如何在vue中使用百度地图添加自定义覆盖物(水波纹)
2020/11/03 Javascript
Python中用pycurl监控http响应时间脚本分享
2015/02/02 Python
python和bash统计CPU利用率的方法
2015/07/10 Python
python操作列表的函数使用代码详解
2017/12/28 Python
Python计算开方、立方、圆周率,精确到小数点后任意位的方法
2018/07/17 Python
python 对给定可迭代集合统计出现频率,并排序的方法
2018/10/18 Python
Django如何防止定时任务并发浅析
2019/05/14 Python
华为校园招聘上机笔试题 扑克牌大小(python)
2020/04/22 Python
100行Python代码实现每天不同时间段定时给女友发消息
2019/09/27 Python
使用IDLE的Python shell窗口实例详解
2019/11/19 Python
opencv+python实现均值滤波
2020/02/19 Python
如何在keras中添加自己的优化器(如adam等)
2020/06/19 Python
Html5新标签解释及用法
2012/02/17 HTML / CSS
Parfume Klik丹麦:香水网上商店
2018/07/10 全球购物
计算机毕业生求职信
2014/06/10 职场文书
学生会竞选演讲稿学习部
2014/08/25 职场文书
2014财务年终工作总结
2014/12/08 职场文书
幼儿园大班教师个人总结
2015/02/05 职场文书
酒店开业主持词
2015/07/02 职场文书
一次项目中Thinkphp绕过禁用函数的实战记录
2021/11/17 PHP
面试分析分布式架构Redis热点key大Value解决方案
2022/03/13 Redis