Python分析睡眠数据

在19年11月的时候买了一个运动手环，然后时不时会用它来记录睡眠数据；积累到现在已经有40个月了。现在想要调整作息，分析一下这些数据，来制定合理的作息计划。

图1 月平均入睡时间

从图1可以看出，我最经常的入睡时间是（02:00:00~02:10:00)之间；

现在我想要早睡，逐步调整，第一个目标值就是（00:50:00~01:00:00）之间，也就意味着我要在（00:20:00~00:30:00）停止刷手机【躺下到入睡需要30分钟】

 图2 月平均入睡时间
从图2可以看出，我经常在（08:00:00~08:10:00）以及（08:40:00~08:50:00）之间醒来。
为了保证睡眠充足，选（08:40:00~08:50:00）作为起床闹钟时间。

 图3 各月睡眠时间条形图

#coding:utf8

import numpy as np
import pandas as pd
import matplotlib as mpl
import matplotlib.pyplot as plt
mpl.use(\"TkAgg\")
import time
import datetime


\"\"\"
目标：统计分析出月入睡时间的数据
大纲：

\"\"\"

\"\"\"
读取excel 
------------------------------
月份        入睡时间    起床时间
2019年11月 1:44:00    9:11:00
2019年12月 1:45:00    9:12:00
------------------------------
\"\"\"
df = pd.read_excel(\'./入睡时间月级.xls\', sheet_name=0)

# 数据准备 按列获取数据作为x坐标轴
x = df[\'月份\'].to_numpy()
y1 = np.array([])
for item in df[\'入睡时间\'].to_numpy():
    y1 = np.append(y1,item.hour+(item.minute)/100)

y2 = np.array([])
for item in df[\'起床时间\'].to_numpy():
    y2 = np.append(y2,item.hour+(item.minute)/100)

y3 = np.array([])
for item in (pd.to_timedelta(df[\'起床时间\'].to_numpy().astype(str)) - pd.to_timedelta(df[\'入睡时间\'].to_numpy().astype(str))):
    y3 = np.append(y3,(item.seconds)/3600)


# 画图1 月平均入睡时间频数图
bins = [0.1,0.2,0.4,0.5,1.0,1.1,1.2,1.3,1.4,1.5,2.0,2.1,2.2,2.3,2.4,2.5,3.0,4.0,5.0,6.0]
segments = pd.cut(y2,bins,right=False)
counts = pd.value_counts(segments,sort=False)
b = plt.bar(counts.index.astype(str),counts)
plt.rcParams[\'font.sans-serif\'] = [\'SimHei\']
plt.xticks(rotation=45)
plt.bar_label(b,counts)
plt.xlabel(\"月平均入睡时间\")
plt.ylabel(\"次\")
plt.title(\"入睡时间频数图\")
plt.show()


# 画2 月平均起床时间频数图
# bins = [7.4,7.5,8.0,8.1,8.2,8.3,8.4,8.5,9.0,9.1,9.2,9.3,9.4,9.5,10.0,10.1,10.2,10.3,11.0]
# segments = pd.cut(y2,bins,right=False)
# counts = pd.value_counts(segments,sort=False)
# b = plt.bar(counts.index.astype(str),counts)
# plt.rcParams[\'font.sans-serif\'] = [\'SimHei\']
# plt.xticks(rotation=45)
# plt.bar_label(b,counts)
# plt.xlabel(\"月平均起床时间\")
# plt.ylabel(\"次\")
# plt.title(\"起床时间频数图\")
# plt.show()

# 画图3 直方图
# plt.rcParams[\'font.sans-serif\'] = [\'SimHei\']
# plt.bar(x,y2)
# plt.bar(x,y1)
# plt.plot(x,y3,\'r\')
# plt.xticks(rotation=60)
# plt.show()

来源：https://www.cnblogs.com/Mozhiqin/p/17107455.html
本站部分图文来源于网络，如有侵权请联系删除。