python计算面积代码_利用Python求阴影部分的面积实例代码
利用Python求阴影部分的面积实例代码
来源:中文源码网 浏览: 次 日期:2019年11月5日
【下载文档: 利用Python求阴影部分的面积实例代码.txt 】
(友情提示:右键点上行txt文档名->目标另存为)
利用Python求阴影部分的面积实例代码一、前言说明
今天看到微信群里一道六年级数学题,如下图,求阴影部分面积看起来似乎并不是很难,可是博主添加各种辅助线,写各种方法都没出来,不得已而改用写Python代码来求面积了
二、思路介绍
1.用Python将上图画在坐标轴上,主要是斜线函数和半圆函数2.均匀的在长方形上面洒满豆子(假设是豆子),求阴影部分豆子占比*总面积三、源码设计
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.py2.计算源码,其中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计算身体质量指数BMI来判断体型
身体质量指数BMI:对身体质量的刻画(Body Mass Index) 国际上常用的衡量人体肥胖和健康程度的重要标准,主要用于统计分析. 定义: 有上图考虑BMI的值的结果因国内和国外的标准不同而不同 ...
- python正方形阴影面积计算_利用Python求阴影部分的面积实例代码
一.前言说明 今天看到微信群里一道六年级数学题,如下图,求阴影部分面积 看起来似乎并不是很难,可是博主添加各种辅助线,写各种方法都没出来,不得已而改用写Python代码来求面积了 二.思路介绍 1.用 ...
- python画pr曲线代码_利用Python中的numpy包实现PR曲线和ROC曲线的计算
闲来无事,边理解PR曲线和ROC曲线,边写了一下计算两个指标的代码.在 python 环境下,sklearn里有现成的函数计算ROC曲线坐标点,这里为了深入理解这两个指标,写代码的时候只用到numpy ...
- python解决物理问题_利用Python科学计算处理物理问题(和物理告个别)
背景: 2019 年初由于尚未学习量子力学相关知识,所以处于自学阶段.浅显的学习了曾谨言的量子力学一卷和格里菲斯编写的量子力学教材.注重将量子力学的一些基本概念了解并理解.同时老师向我们推荐了 Qua ...
- python计算密集型提速_利用Cython加速计算密集型python任务
何为计算密集型任务 下面贴上网上找到的描述计算密集型任务的特点是要进行大量的计算,消耗CPU资源,比如计算圆周率.对视频进行高清解码等等,全靠CPU的运算能力.这种计算密集型任务虽然也可以用多任务完成 ...
- python贪吃蛇最简单代码_利用python实现简易版的贪吃蛇游戏(面向python小白)
引言 作为python 小白,总是觉得自己要做好百分之二百的准备,才能开始写程序.以至于常常整天在那看各种语法教程,学了几个月还是只会print('hello world'). 这样做效率太低,正确的 ...
- python做游戏代码_利用Python基础代码语句,实现2G时代文字小游戏,世界如此简单!...
相信许多80,90后都玩过2G时代的文字小游戏,它是来自QQ家园的专属回忆.偷菜,美味小镇,大乐斗,还有精武堂等等,虽然只是文字的输出,但是留给我们这一代的人的印象却是最深刻的.曾经流量很少,响应很快 ...
- python制作图片墙_利用python生成照片墙的示例代码
PIL(Python Image Library)是python的第三方图像处理库,但是由于其强大的功能与众多的使用人数,几乎已经被认为是python官方图像处理库了.其官方主页为:PIL. PIL历 ...
- python 打卡记录代码_利用Python实现对考勤打卡数据处理的总结
利用Python实现对考勤打卡数据处理的总结 一.背景交代 二.说明 三. 8种方法 1. 查看文件是否存在 2. 导入excel文件,并把数据保存为dataframe格式 3. 计算程序运行时间 4 ...
最新文章
- 查询字符串中字母出现的个数
- 如何正确运用计算机,如何正确使用电脑
- 【备用】关于BOM替代物料与CK11N取数逻辑
- C++ STL : SGI-STL空间配置器源码剖析
- SQLServer2008安装教程
- named-config with name ‘c3p0-config.xml‘ does not exist. Using default-config
- 怎样设置计算机硬盘的空间,如何调整电脑硬盘的分区容量
- cmake 添加asan功能
- MYSQL导入数据出现ERROR 1290 (HY000): The MySQL server is running with the --secure-file-priv option so it
- Python-Pillow给图片加文字
- 数据挖掘:概念与技术 第五章-数据立方体技术
- selenuim爬虫实战(日lofter.com)
- unity中单位是米还是厘米_cm在单位里是厘米还是毫米
- 远程报:这可能是由于credssp加密oracle修正
- Java类型转换笔记
- 基于Pyhton的二维离散正弦变换(DST)及其反变换(IDST)
- 广和通L610+阿里云平台+钉钉机器人
- NJR日清纺微理光R5445系列 单节锂电池保护芯片,内置驱动器高位Nch FET开关温度保护
- 私有云的优缺点_什么是公有云、私有云、混合云 ?各自的优缺点有哪些?
- 贝叶斯方法---分段线性插值函数画图
热门文章
- ROS实战(一) 语音交互系统的学习:初步了解语音交互流程
- ListView-电影列表
- 查询IP地址的个性网站
- Python turtle库(螺旋正方形)
- 【正点原子FPGA连载】第四十一章OV5640摄像头RGB-LCD显示实验 -摘自【正点原子】新起点之FPGA开发指南_V2.1
- cornell数据集文件整理
- https://freeminecraftgiftcode.net/
- 平均年薪60.8万,拿下这个证书,数据挖掘就能起飞!
- 西瓜书习题 - 6.支持向量机
- docker swarm + nginx服务离线部署