python数据分析与展示(四)

实践:
课程的示例,引力波的绘制:
资源文件:http://python123.io/dv/grawave.html

import numpy as npimport matplotlib.pyplot as pltfrom scipy.io import wavfile#波形文件#速率 数据rate_h,hstrain= wavfile.read(r"H1_Strain.wav","rb")# print(rate_h,hstrain)#print(len(hstrain))#print(hstrain.shape)rate_l, lstrain = wavfile.read(r'L1_Strain.wav', 'rb')#读取出来 转置#理论模型reftime,ref_H1 = np.genfromtxt('wf_template.txt').transpose()#计算出时间间隔htime_interval = 1/rate_hltime_interval = 1 / rate_l#计算出单位时间的数据量htime_len = hstrain.shape[0]/rate_h# 时间坐标原点在中心 绘制时间与数据一一对应htime = np.arange(-htime_len/2,htime_len/2,htime_interval)# print(len(htime))ltime_len = lstrain.shape[0] / rate_lltime = np.arange(-ltime_len/2,ltime_len/2,ltime_interval)#绘制H1 Strainfig = plt.figure(figsize=(12,6))plth = fig.add_subplot(221)plth.plot(htime,hstrain,'y')plth.set_xlabel("时间(s)",fontproperties="Kaiti")plth.set_ylabel("H1 Strain")plth.set_title("探测器数据1",fontproperties="Kaiti",fontsize=20)#绘制 L1  Strainpltl = fig.add_subplot(222)pltl.plot(ltime,lstrain,'g')pltl.set_xlabel("时间(s)",fontproperties="Kaiti")pltl.set_ylabel("L1 Strain")pltl.set_title("探测数据2",fontproperties="Kaiti",fontsize=20)#绘制标准数据pltref = fig.add_subplot(212)pltref.plot(reftime,ref_H1)pltref.set_xlabel("时间(s)",fontproperties="Kaiti")pltref.set_ylabel("Template")pltref.set_title("标准数据",fontproperties="Kaiti",fontsize=20)fig.tight_layout()plt.show()

自己编写的例子

从 中国国家统计局获取分省经济增长数据绘制折线图网址:http://www.stats.gov.cn/tjsj/代码:

import numpy as npimport matplotlib.pyplot as pltimport matplotlibmatplotlib.rcParams['font.family']='SimHei'# 地区city_labels=['北京市','天津市','河北省','河北省','内蒙古自治区','辽宁省','吉林省']#获取数据data = np.loadtxt('分省月度数据.csv',delimiter=',',dtype=float)#获取数据形状print(data.shape)#抽取城市,因为城市过多,不利于分析,抽取7组city_data_style=data.shape[0]num=0for city in range(7):    plt.plot(data[city][::-1])plt.ylabel("较去年同期增长百分比(%)",fontproperties="Kaiti",fontsize=17)plt.xlabel("月份",fontproperties="Kaiti",fontsize=17)plt.legend(city_labels,loc=0,ncol=3,)plt.xticks([0,1,2,3,4,5,6,7,8,9,10,11,],           ['16年7月','16年8月','16年9月','16年10月','16年11月','16年12月','17年1月','17年2月','17年3月','17年4月','17年5月','17年6月',])plt.title("16年7月-17年6月,分省同比增长数据显示图",fontproperties="Kaiti",fontsize=20,)plt.show()