原始数据如下：第一列为序号，第二列是主频值，第三列是振幅值。
import matplotlib.pyplot as pltimport pandas as pdimport numpy as npimport datetimex = []y = []z = []N = 13736def file_path():read_path = r'G:\testdatasets\1_weidatachuli\xuhao+KHz+zhenfu.txt'file_name = 'xuhao+KHz+zhenfu.txt'return read_path,file_namedef process_data():data_num = [line.strip().split('\t') for line in open(file_name)]for d1,d2,d3 in data_num:x.append(int(d1))y.append(float(d2))z.append(float(d3))# #进行排序# y.sort()#此时的y已经按照从下到大排好顺序# #print(y)#进行分类number_50 = 0number_50_150 =0number_150_250 = 0number_250_350 =0number_350 = 0for i in range (0,len(y)-1):if y[i] <= 50:number_50 += 1elif 50< y[i] <= 150:number_50_150 += 1if 150< y[i] <= 250:number_150_250 += 1elif 250< y[i] <= 350:number_250_350 += 1elif 350< y[i] :number_350 += 1print(('小于50KHz(低频)的个数为：{0}\n'+'50KHz--150KHz（中低频）的个数为：{1}\n'+'150KHz--250KHz（中频）的个数为：{2}\n'+'250KHz--350KHz（中高频）的个数为：{3}\n'+'大于350KHz（高频）的个数为：{4}\n').format(number_50 ,number_50_150 ,number_150_250 ,number_250_350 ,number_350))return x,y,z,number_50,number_50_150,number_150_250,number_250_350,number_350def main():x_label = ['0-50KHz(低频):0.88%','50-150KHz（中低频）:''58.94%','150-250KHz（中频）:14.97%','250-350KHz（中高频）:20.96%','350+KHz（高频）:4.25%']y_label = [number_50,number_50_150,number_150_250,number_250_350,number_350]plt.rcParams['font.sans-serif'] = ['SimHei'] # 显示中文plt.rcParams['axes.unicode_minus'] = Falseplt.figure(file_name+'莱州黄岗岩岩爆',figsize = (7,4)) # #plt.figure(figsize = (7,4))#设置图片大小plt.title('声发射主频折线图', fontsize=10)colors = ['#696969','#808080','#A9A9A9','#DCDCDC','#F5F5F5'] # 自定义颜色plt.subplot(211) #2行2列，第1个plt.title('声发射主频饼状图', fontsize=20)explode = [0,0,0,0,0] #扩大某一项plt.pie(x=y_label,labels=x_label,autopct='%3.2f%%',#图内显示百分数colors = colors,#设置饼图自定义颜色explode = explode,shadow = True,#shadow = True：阴影，有立体感textprops={'color':'black',#文本颜色'fontsize':'16',#文本大小#'fontfamily':'Microsoft JhengHe',#设置微软雅黑字体})plt.subplot(212)#2行2列，第2个plt.title('声发射中低频饼状分割图', fontsize=20)explode = [0,0.1,0,0,0] #扩大某一项plt.pie(x=y_label,labels=x_label,autopct='%3.2f%%',explode = explode,shadow = True, #shadow = True：阴影，有立体感textprops={'color': 'black', # 文本颜色'fontsize': '16', # 文本大小# 'fontfamily':'Microsoft JhengHe',#设置微软雅黑字体})if __name__ == '__main__':start_time = datetime.datetime.now()read_path, file_name = file_path()x,y,z,number_50,number_50_150,number_150_250,number_250_350,number_350 = process_data()main()end_time = datetime.datetime.now()run_time = end_time-start_timeprint(('程序运行时间{}').format(run_time))plt.show()
结果如图所示：
