1、问题
买了一副眼镜，花了 800 块钱，感觉很贵，然后我在心里突然思考，这个眼镜如果买了一天就不见了，那么相当于我这一天使用的成本是 800 块钱。如果使用了 4 天之后才不见，那么我这四天每天使用的成本就是 200 块钱。也就是说如果这一副眼镜我用的时间越多，每天所花费的使用成本就越低。那么我想要每天的使用成本低于 5 毛钱，至少要花费多少时间呢？
带着这个问题，我用 Python 写了一个简单的程序，使用 Matplotlib 来可视化，很直观地就看到要花多少时间才能够每天的平摊价格达到自己可接受的程度。
可以看到至少要使用 4 年半，才可以使得每天的使用成本低于 0.5 元。也就是说，我要好好地保养买的这个眼镜，至少要保证可以用到 4.5 年才行。
2、代码
下面是画出上面这个图的 Python 代码：
# 计算一个价格为 price 的物件, 在使用了 days 之后
# 平摊到每天所花的价格是多少.
# 例如一副眼镜是 800 元, 在第二天的时候平摊价格为 400 元/天
# 到第四天的时候平摊价格为 200 元/天. 有的时候我们想要知道
# 用到多少天的时候这个眼镜的平摊价格才能降到 0.5 元/天
import numpy as np
import matplotlib.pyplot as plt
price = 800 # 物品的价格
expected_prices = 0.5 # 预期可接受的平摊价格
max_day = int(price / expected_prices * 1.5) # 画出的最大时间为两年
amortized_prices = [price / day for day in range(1, max_day + 1)]
day_hit_expected = 0 # 达到预期平摊价格的天数
for i in range(max_day):
if amortized_prices[i] < expected_prices:
day_hit_expected = i + 1
break
# 每隔多少天设为一个区间单位, 用于画出多个虚线, 更加直观地看到自己关心经过饿了多少个时间区间
# 例如设定为 365 年的话可以看到大概要经过多少年才能达到自己预测可接受的平摊价格
bins_days = 365
bins_label = 'Year' # 相对地替换成自己所设定的天数代表的意思
xs = np.linspace(1, max_day, max_day)
plt.plot(xs, amortized_prices)
# axline 用于画一条无限长的线, 在这里用来标出一条代表达到 expected prices 的虚线
plt.axhline(expected_prices,
linestyle='--', linewidth=0.75, color='black')
for i in range(1, int(max_day / bins_days) + 1):
plt.axvline(i * bins_days,
linestyle='--', linewidth=0.75, color='orange')
plt.xlim(0, max_day)
plt.yscale('log')
# 改变一下坐标轴的显示, 以 bins_days 作为单位
old_xticks = [i * bins_days for i in range(int(max_day / bins_days) + 1)]
new_xticks = [i for i in range(int(max_day / bins_days) + 1)]
plt.xticks(old_xticks, new_xticks)
# 标出达到可接受平摊价格时的天数以及价格
annotate = "{:.2f}¥ on day {}".format(amortized_prices[day_hit_expected], day_hit_expected)
point_pos = (day_hit_expected, amortized_prices[day_hit_expected])
text_pos = (point_pos[0] * 1.05, point_pos[1] * 1.55)
plt.annotate(annotate,
xy=point_pos, xytext=text_pos,
arrowprops=dict(facecolor='k', width=0.5, headwidth=5, headlength=5))
plt.xlabel(bins_label)
plt.ylabel('Yuan')
plt.show()
对于 Matplotlib 有一些好用的函数可以讲一下：
plt.axline用于画直线（无限长的线）的函数：**kwargs是指定 2D 曲线的属性，比如线宽（linewidth）、颜色（color）、线型（linestyle）。plt.axline(xy1, xy2=None, *, slope=None, **kwargs)：其中xy1为直线通过的第一个点，xy2为第二个点，也可以不提供第二个点而是提供斜率slope；plt.axhline(y=0, xmin=0, xmax=1, **kwargs)：画一条值在y的平行线，xmin和xmax的值指定为 0 到 1 之间，其中 0 是图的最左边，1 是图的最右边；plt.axvline(x=0, ymin=0, ymax=1, **kwargs)：画一条值在x的平行线，ymin和ymax的值指定为 0 到 1 之间，其中 0 是图的最底端，1 是图的最上端。
plt.annotate(text, xy, xytext=None, ..., arrowprops=None)用于画出注释的函数：text是注释的内容，xy是注释指向的点的位置，xytext是注释文本的位置，arrowprops是指定箭头格式（arrowstyle）的一个字典（dict）；plt.xticks(ticks=None, labels=None, *, minor=False, **kwargs)改变坐标显示刻度：ticks为要替换的原刻度，labels为新的刻度，如果将ticks设为空则代表不显示坐标刻度。minor如果为False代表设定主坐标轴，为True代表设定次坐标轴。
在本文的代码中将原来表示天数的刻度 old_xticks 替代为代表年的 new_xticks。
