From eba7c6428aebe1827d0858577ce2e52f06f7ba7b Mon Sep 17 00:00:00 2001 From: zhang Date: Sat, 16 Nov 2019 09:56:55 +0800 Subject: double 11 --- .idea/misc.xml | 4 ++++ .idea/modules.xml | 8 ++++++++ .idea/vlog.iml | 8 ++++++++ double_11.py | 19 +++++++++++++++++++ readme.md | 6 ++++++ 5 files changed, 45 insertions(+) create mode 100644 .idea/misc.xml create mode 100644 .idea/modules.xml create mode 100644 .idea/vlog.iml create mode 100644 double_11.py diff --git a/.idea/misc.xml b/.idea/misc.xml new file mode 100644 index 0000000..65531ca --- /dev/null +++ b/.idea/misc.xml @@ -0,0 +1,4 @@ + + + + \ No newline at end of file diff --git a/.idea/modules.xml b/.idea/modules.xml new file mode 100644 index 0000000..96c2fc7 --- /dev/null +++ b/.idea/modules.xml @@ -0,0 +1,8 @@ + + + + + + + + \ No newline at end of file diff --git a/.idea/vlog.iml b/.idea/vlog.iml new file mode 100644 index 0000000..f8446a8 --- /dev/null +++ b/.idea/vlog.iml @@ -0,0 +1,8 @@ + + + + + + + + \ No newline at end of file diff --git a/double_11.py b/double_11.py new file mode 100644 index 0000000..118cfb7 --- /dev/null +++ b/double_11.py @@ -0,0 +1,19 @@ + +import numpy as np +import matplotlib.pyplot as plt + +years = [y+2009 for y in range(11)] +gmvs = [0.5, 9.36, 52, 191, 350, 571, 912, 1207, 1682, 2135, 2684] + + +def polynomial(xs, ys, n): + f = np.polyfit(xs, ys, n) + def func(x): + return sum(f[i]*x**(n-i)for i in range(n+1)) + return func + +f = polynomial(years, gmvs, 3) +print(f(2020)) + +plt.plot(years, gmvs, marker='o') +plt.show() \ No newline at end of file diff --git a/readme.md b/readme.md index e69de29..8a990da 100644 --- a/readme.md +++ b/readme.md @@ -0,0 +1,6 @@ + +- monte_carlo_triangle.py, [【计算机科学】蒙特卡洛方法计算“一棍砍两刀”构成三角形的概率](https://www.bilibili.com/video/av75100858/) +- double_11.py + + + -- cgit v1.2.3