1 # -*- coding: utf-8 -*-
2
3 ###############################################################################
4 ####################### 正文代码 #######################
5 ###############################################################################
6
7 # 代码 3-5
8 import numpy as np
9 import matplotlib.pyplot as plt
10 plt.rcParams['font.sans-serif'] = 'SimHei' ## 设置中文显示
11 plt.rcParams['axes.unicode_minus'] = False
12 data = np.load('../data/国民经济核算季度数据.npz')
13 name = data['columns'] ## 提取其中的columns数组，视为数据的标签
14 values = data['values']## 提取其中的values数组，数据的存在位置
15
16 plt.figure(figsize=(8,7))## 设置画布
17 plt.scatter(values[:,0],values[:,2], marker='o')## 绘制散点图
18 plt.xlabel('年份')## 添加横轴标签
19 plt.ylabel('生产总值（亿元）')## 添加y轴名称
20 plt.xticks(range(0,70,4),values[range(0,70,4),1],rotation=45)
21 plt.title('2000-2017年季度生产总值散点图')## 添加图表标题
22 plt.savefig('../tmp/2000-2017年季度生产总值散点图.png')
23 plt.show()
24
25
26
27 # 代码 3-6
28 plt.figure(figsize=(8,7))## 设置画布
29 ## 绘制散点1
30 plt.scatter(values[:,0],values[:,3], marker='o',c='red')
31 ## 绘制散点2
32 plt.scatter(values[:,0],values[:,4], marker='D',c='blue')
33 ## 绘制散点3
34 plt.scatter(values[:,0],values[:,5], marker='v',c='yellow')
35 plt.xlabel('年份')## 添加横轴标签
36 plt.ylabel('生产总值（亿元）')## 添加纵轴标签
37 plt.xticks(range(0,70,4),values[range(0,70,4),1],rotation=45)
38 plt.title('2000-2017年各产业季度生产总值散点图')## 添加图表标题
39 plt.legend(['第一产业','第二产业','第三产业'])## 添加图例
40 plt.savefig('../tmp/2000-2017年各产业季度生产总值散点图.png')
41 plt.show()
42
43
44 # 代码 3-7
45 plt.figure(figsize=(8,7))## 设置画布
46 ## 绘制折线图
47 plt.plot(values[:,0],values[:,2],color = 'r',linestyle = '--')
48 plt.xlabel('年份')## 添加横轴标签
49 plt.ylabel('生产总值（亿元）')## 添加y轴名称
50 plt.xticks(range(0,70,4),values[range(0,70,4),1],rotation=45)
51 plt.title('2000-2017年季度生产总值折线图')## 添加图表标题
52 plt.savefig('../tmp/2000-2017年季度生产总值折线图.png')
53 plt.show()
54
55
56
57 # 代码 3-8
58 plt.figure(figsize=(8,7))## 设置画布
59 plt.plot(values[:,0],values[:,2],color = 'r',linestyle = '--',
60 marker = 'o')## 绘制折线图
61 plt.xlabel('年份')## 添加横轴标签
62 plt.ylabel('生产总值（亿元）')## 添加y轴名称
63 plt.xticks(range(0,70,4),values[range(0,70,4),1],rotation=45)
64 plt.title('2000-2017年季度生产总值点线图')## 添加图表标题
65 plt.savefig('../tmp/2000-2017年季度生产总值点线图.png')
66 plt.show()
67
68
69
70 # 代码 3-9
71 plt.figure(figsize=(8,7))## 设置画布
72 plt.plot(values[:,0],values[:,3],'bs-',
73 values[:,0],values[:,4],'ro-.',
74 values[:,0],values[:,5],'gH--')## 绘制折线图
75 plt.xlabel('年份')## 添加横轴标签
76 plt.ylabel('生产总值（亿元）')## 添加y轴名称
77 plt.xticks(range(0,70,4),values[range(0,70,4),1],rotation=45)
78 plt.title('2000-2017年各产业季度生产总值折线图')## 添加图表标题
79 plt.legend(['第一产业','第二产业','第三产业'])
80 plt.savefig('../tmp/2000-2017年季度各产业生产总值折线图.png')
81 plt.show()
82
83
84
85 ###############################################################################
86 ####################### 任务实现 #######################
87 ###############################################################################
88
89 # 代码 3-10
90 import numpy as np
91 import matplotlib.pyplot as plt
92 plt.rcParams['font.sans-serif'] = 'SimHei' ## 设置中文显示
93 plt.rcParams['axes.unicode_minus'] = False
94 data = np.load('../data/国民经济核算季度数据.npz')
95 name = data['columns']## 提取其中的columns数组，视为数据的标签
96 values = data['values']## 提取其中的values数组，数据的存在位置
97 p = plt.figure(figsize=(12,12)) ##设置画布
98
99 ## 子图1
100 ax1 = p.add_subplot(2,1,1)
101 plt.scatter(values[:,0],values[:,3], marker='o',c='r')## 绘制散点
102 plt.scatter(values[:,0],values[:,4], marker='D',c='b')## 绘制散点
103 plt.scatter(values[:,0],values[:,5], marker='v',c='y')## 绘制散点
104 plt.ylabel('生产总值（亿元）')## 添加纵轴标签
105 plt.title('2000-2017年各产业季度生产总值散点图')## 添加图表标题
106 plt.legend(['第一产业','第二产业','第三产业'])## 添加图例
107
108 ## 子图2
109 ax2 = p.add_subplot(2,1,2)
110 plt.scatter(values[:,0],values[:,6], marker='o',c='r')## 绘制散点
111 plt.scatter(values[:,0],values[:,7], marker='D',c='b')## 绘制散点
112 plt.scatter(values[:,0],values[:,8], marker='v',c='y')## 绘制散点
113 plt.scatter(values[:,0],values[:,9], marker='8',c='g')## 绘制散点
114 plt.scatter(values[:,0],values[:,10], marker='p',c='c')## 绘制散点
115 plt.scatter(values[:,0],values[:,11], marker='+',c='m')## 绘制散点
116 plt.scatter(values[:,0],values[:,12], marker='s',c='k')## 绘制散点
117 ## 绘制散点
118 plt.scatter(values[:,0],values[:,13], marker='*',c='purple')
119 ## 绘制散点
120 plt.scatter(values[:,0],values[:,14], marker='d',c='brown')
121 plt.legend(['农业','工业','建筑','批发','交通',
122 '餐饮','金融','房地产','其他'])
123 plt.xlabel('年份')## 添加横轴标签
124 plt.ylabel('生产总值（亿元）')## 添加纵轴标签
125 plt.xticks(range(0,70,4),values[range(0,70,4),1],rotation=45)
126 plt.savefig('../tmp/2000-2017年季度各行业生产总值散点子图.png')
127 plt.show()
128
129
130
131 # 代码 3-11
132 p1 = plt.figure(figsize=(8,7))## 设置画布
133 ## 子图1
134 ax3 = p1.add_subplot(2,1,1)
135 plt.plot(values[:,0],values[:,3],'b-',
136 values[:,0],values[:,4],'r-.',
137 values[:,0],values[:,5],'g--')## 绘制折线图
138 plt.ylabel('生产总值（亿元）')## 添加纵轴标签
139 plt.title('2000-2017年各产业季度生产总值折线图')## 添加图表标题
140 plt.legend(['第一产业','第二产业','第三产业'])## 添加图例
141 ## 子图2
142 ax4 = p1.add_subplot(2,1,2)
143 plt.plot(values[:,0],values[:,6], 'r-',## 绘制折线图
144 values[:,0],values[:,7], 'b-.',## 绘制折线图
145 values[:,0],values[:,8],'y--',## 绘制折线图
146 values[:,0],values[:,9], 'g:',## 绘制折线图
147 values[:,0],values[:,10], 'c-',## 绘制折线图
148 values[:,0],values[:,11], 'm-.',## 绘制折线图
149 values[:,0],values[:,12], 'k--',## 绘制折线图
150 values[:,0],values[:,13], 'r:',## 绘制折线图
151 values[:,0],values[:,14], 'b-')## 绘制折线图
152 plt.legend(['农业','工业','建筑','批发','交通',
153 '餐饮','金融','房地产','其他'])
154 plt.xlabel('年份')## 添加横轴标签
155 plt.ylabel('生产总值（亿元）')## 添加纵轴标签
156 plt.xticks(range(0,70,4),values[range(0,70,4),1],rotation=45)
157 plt.savefig('../tmp/2000-2017年季度各行业生产总值折线子图.png')
158 plt.show()
