python函数实验答案,实验二Python语言基础函数包练习

《实验二Python语言基础函数包练习》由会员分享，可在线阅读，更多相关《实验二Python语言基础函数包练习(13页珍藏版)》请在人人文库网上搜索。

1、实验二 Python语言基础函数包练习姓名： 季鹏 班级： 1208 学号：30实验目的1、熟练掌握Python语言基础函数包，如math、NumPy、SciPy和Matplotlib等函数包的使用实验内容 练习1 基本统计可视化：Step 1. 导入编程所需要函数包：Step 2. 柱状图可视化：可视化结果(将图粘贴到这个位置)：Step 3. 散点图可视化：可视化结果(将图粘贴到这个位置)：Step 4：box-plot可视化可视化结果(将图粘贴到这个位置)：Step 5: 三次样条插值可视化结果(将图粘贴到这个位置)：练习内容 1、要求读取某课程期末考试的成绩“scores.csv”，字。

2、段信息如下：ID: 学号；fscore：期末考试成绩groupe：分组组别class：班级score1：第一次平时成绩score2：第二次平时成绩score3：第三次平时成绩score：最终综合成绩要求：1) 分别画出四个成绩的柱状图，观察是否满足正态分布；2) 分别拟合第一次、第二次、第三次成绩和期末考试成绩之间的关系；3) 分别按照班级和组别做boxplot,观察每一组和每一班之间成绩的浮动区间。提示1： 读取csv文件import csvin_file = open(scores.csv, r)csv_reader = csv.reader(in_file, delimiter=,)# 。

3、read the header infoheader = csv_reader.next()print header代码#coding UTF-8import numpy as npimport scipy as spyimport matplotlib.mlab as mlabimport matplotlib.pyplot as pltfrom scipy.optimize import leastsq import pylab as plimport csvin_file =open(E:studtstudylearnpythongisExperiment2Experiment2scor。

4、es.csv,r)csv_reader=csv.reader(in_file,delimiter=,)# read the header infoheader= csv_reader.next()#print header#defineID=group=Class=fscore=score1=score2=score3=score=#read datafor line in csv_reader:#ingnore firstlineif csv_reader.line_num=0:continueID.append(line0)fscore.append(line1)group.append(。

5、float(line2)Class.append(float(line3)score1.append(float(line4)score2.append(float(line5)score3.append(float(line6)score.append(float(line7)#print ID#print group#print scoreplt.figure(1)#score1plt.subplot(411)n1,bins1,patches1=plt.hist(score1,25,normed=1,facecolor=green,alpha=0.8)y1=mlab.normpdf(bin。

6、s1,100,15)l=plt.plot(bins1,y1,r-,linewidth=1)plt.xlabel(ID)plt.ylabel(score1)plt.title(Histogram of score)#score2plt.subplot(412) n1,bins1,patches1=plt.hist(score2,25,normed=1,facecolor=red,alpha=0.8)y1=mlab.normpdf(bins1,100,15)l=plt.plot(bins1,y1,r-,linewidth=1)plt.xlabel(ID)plt.ylabel(score2)#sco。

7、re3plt.subplot(413)n1,bins1,patches1=plt.hist(score3,25,normed=1,facecolor=yellow,alpha=0.8)y1=mlab.normpdf(bins1,100,15)l=plt.plot(bins1,y1,r-,linewidth=1)plt.xlabel(ID)plt.ylabel(score3)#scoreax4 = plt.subplot(414)n1,bins1,patches1=plt.hist(score,25,normed=1,facecolor=blue,alpha=0.8)y1=mlab.normpd。

8、f(bins1,100,15)l=plt.plot(bins1,y1,r-,linewidth=1)plt.xlabel(ID)plt.ylabel(score)#plt.show()#222#2plt.figure(2)fig = plt.figure(2) def pic(score1,score,ax):import random order=9plt.plot(score,score1,color=m,linestyle=,marker=.) matA= for i in range(0,order+1): matA1= for j in range(0,order+1): tx=0.。

9、0 for k in range(0,len(score1):dx=1.0 for l in range(0,j+i): dx=dx*score1k tx+=dx matA1.append(tx) matA.append(matA1) #print(len(xa) #print(matA00) matA=np.array(matA) matB= for i in range(0,order+1): ty=0.0 for k in range(0,len(score1): dy=1.0 for l in range(0,i): dy=dy*score1k ty+=scorek*dy matB.a。

10、ppend(ty) matB=np.array(matB) matAA=np.linalg.solve(matA,matB) #print(matAA) yya=for i in range(0,len(score1): yy=0.0 for j in range(0,order+1): dy=1.0 for k in range(0,j): dy*=scorei dy*=matAAj yy+=dy yya.append(yy) ax.plot(score,yya,color=g,linestyle=-,marker=) ax.legend()return 1ax1 = fig.add_sub。

11、plot(311) pic(score1,score,ax1)ax2 = fig.add_subplot(312) pic(score2,score,ax2)ax3 = fig.add_subplot(313) pic(score3,score,ax3)#aa=for i in range(len(group):if groupi not in aa:temp=groupifor j in range(len(group):if groupj=temp:aa.append(temp)breakax=0 for col in range(5) for row in range(len(aa)m=。

12、0n=0for i in range(len(aa):for j in range(len(group):if groupj=aai:axmn=scorejn=n+1m=m+1n=0#print axi=0labels=list(aa)for i in range(len(aa):labels.append(str(aai)fs=10fig, axes=plt.subplots(nrows=1, ncols=len(aa), figsize=(40,10)for i in range(len(aa):axesi.boxplot(filter(lambda e:e!=0, axi),patch_。

13、artist=False,bootstrap=None)axesi.set_title(group+str(int(aai),fontsize=fs)i=0bb=list(set(Class)bx=0 for col in range(len(Class) for row in range(len(bb)m=0n=0for i in range(len(bb):for j in range(len(Class):if Classj=bbi:bxmn=scorejn=n+1m=m+1n=0#print bx#print the picturefig2, axes2=plt.subplots(nrows=1, ncols=len(bb), figsize=(10,10)i=0for i in range(len(bb):axes2i.boxplot(filter(lambda e:e!=0, bxi),patch_artist=False,bootstrap=None)axes2i.set_title(class+str(int(bbi),fontsize=fs)plt.show(。