一、前言说明
今天看到微信群里一道六年级数学题,如下图,求阴影部分面积
看起来似乎并不是很难,可是博主添加各种辅助线,写各种方法都没出来,不得已而改用写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.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.下面有三种求解方法,第三种表示比大学高数还难看懂,你们呢?
总结
以上就是这篇文章的全部内容了,希望本文的内容对大家的学习或者工作具有一定的参考学习价值,如果有疑问大家可以留言交流,谢谢大家对的支持。
《魔兽世界》大逃杀!60人新游玩模式《强袭风暴》3月21日上线
暴雪近日发布了《魔兽世界》10.2.6 更新内容,新游玩模式《强袭风暴》即将于3月21 日在亚服上线,届时玩家将前往阿拉希高地展开一场 60 人大逃杀对战。
艾泽拉斯的冒险者已经征服了艾泽拉斯的大地及遥远的彼岸。他们在对抗世界上最致命的敌人时展现出过人的手腕,并且成功阻止终结宇宙等级的威胁。当他们在为即将于《魔兽世界》资料片《地心之战》中来袭的萨拉塔斯势力做战斗准备时,他们还需要在熟悉的阿拉希高地面对一个全新的敌人──那就是彼此。在《巨龙崛起》10.2.6 更新的《强袭风暴》中,玩家将会进入一个全新的海盗主题大逃杀式限时活动,其中包含极高的风险和史诗级的奖励。
《强袭风暴》不是普通的战场,作为一个独立于主游戏之外的活动,玩家可以用大逃杀的风格来体验《魔兽世界》,不分职业、不分装备(除了你在赛局中捡到的),光是技巧和战略的强弱之分就能决定出谁才是能坚持到最后的赢家。本次活动将会开放单人和双人模式,玩家在加入海盗主题的预赛大厅区域前,可以从强袭风暴角色画面新增好友。游玩游戏将可以累计名望轨迹,《巨龙崛起》和《魔兽世界:巫妖王之怒 经典版》的玩家都可以获得奖励。
更新动态
- DzintraErliha-Dreamscape(2024)[24Bit-96kHz]FLAC
- 群星.1999-上华世纪国语总冠军4CD【上华】【WAV+CUE】
- 叶丽仪.1992-仙鹤神针电视剧原声音乐专辑【星传真】【WAV+CUE】
- 陈柏宇.2011-QUINQUENNIUM【SONY】【WAV+CUE】
- 《51首热门当红歌曲 十大顶级发烧女声 3CD》[WAV+CUE][1.6GB]
- 《赵传 深受摇滚迷青睐的台湾音乐人 约定》[WAV+CUE][320MB]
- 《庄心妍 实力唱将 发烧音乐精选集》[WAV+CUE][380MB]
- [音展纪念碟]2006年《TAA台湾第十六届国际HI-ENDHI-FI音响展》2cd[WAV]
- 中国音乐地图之听见西藏山歌牧歌酒歌2020[WAV分轨]
- 王宏伟2017-永远的红飘带[太平洋影音][WAV+CUE]
- 费翔.1997-NOW现在【第一中国】【WAV+CUE】
- 新格合唱团.1981-我唱你和2辑(滚石新格民歌系列)【新格】【WAV+CUE】
- 陈柏宇.2011-QUINQUENNIUM【SONY】【WAV+CUE】
- 雨果唱片-《金唛系列-蔷薇依韵琼瑶》WAV
- 陈百强《凝望(2012华纳LPCD45限量版)[DMI]》[WAV整轨]