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 各月睡眠时间条形图

1. #coding:utf8
3. import numpy as np
4. import pandas as pd
5. import matplotlib as mpl
6. import matplotlib.pyplot as plt
7. mpl.use("TkAgg")
8. import time
9. import datetime
12. """
13. 目标：统计分析出月入睡时间的数据
14. 大纲：
16. """
18. """
19. 读取excel
20. ------------------------------
21. 月份        入睡时间    起床时间
22. 2019年11月 1:44:00    9:11:00
23. 2019年12月 1:45:00    9:12:00
24. ------------------------------
25. """
26. df = pd.read_excel('./入睡时间月级.xls', sheet_name=0)
28. # 数据准备 按列获取数据作为x坐标轴
29. x = df['月份'].to_numpy()
30. y1 = np.array([])
31. for item in df['入睡时间'].to_numpy():
32.     y1 = np.append(y1,item.hour+(item.minute)/100)
34. y2 = np.array([])
35. for item in df['起床时间'].to_numpy():
36.     y2 = np.append(y2,item.hour+(item.minute)/100)
38. y3 = np.array([])
39. for item in (pd.to_timedelta(df['起床时间'].to_numpy().astype(str)) - pd.to_timedelta(df['入睡时间'].to_numpy().astype(str))):
40.     y3 = np.append(y3,(item.seconds)/3600)
43. # 画图1 月平均入睡时间频数图
44. 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]
45. segments = pd.cut(y2,bins,right=False)
46. counts = pd.value_counts(segments,sort=False)
47. b = plt.bar(counts.index.astype(str),counts)
48. plt.rcParams['font.sans-serif'] = ['SimHei']
49. plt.xticks(rotation=45)
50. plt.bar_label(b,counts)
51. plt.xlabel("月平均入睡时间")
52. plt.ylabel("次")
53. plt.title("入睡时间频数图")
54. plt.show()
57. # 画2 月平均起床时间频数图
58. # 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]
59. # segments = pd.cut(y2,bins,right=False)
60. # counts = pd.value_counts(segments,sort=False)
61. # b = plt.bar(counts.index.astype(str),counts)
62. # plt.rcParams['font.sans-serif'] = ['SimHei']
63. # plt.xticks(rotation=45)
64. # plt.bar_label(b,counts)
65. # plt.xlabel("月平均起床时间")
66. # plt.ylabel("次")
67. # plt.title("起床时间频数图")
68. # plt.show()
70. # 画图3 直方图
71. # plt.rcParams['font.sans-serif'] = ['SimHei']
72. # plt.bar(x,y2)
73. # plt.bar(x,y1)
74. # plt.plot(x,y3,'r')
75. # plt.xticks(rotation=60)
76. # plt.show()

复制代码

 

来源:https://www.cnblogs.com/Mozhiqin/p/17107455.html
免责声明：由于采集信息均来自互联网，如果侵犯了您的权益，请联系我们【E-Mail:cb@itdo.tech】 我们会及时删除侵权内容，谢谢合作！