0
我對python非常陌生,並且正在使用它來分析仿真數據。我在Anaconda(python 2.7)中編寫的代碼工作了幾次,現在似乎變得無法響應。Spyder在Win 7上運行的Anaconda(py 2.7)中變得無響應64位
我試過重新啓動控制檯,甚至是卸載並重新安裝。但它不起作用。 有人可以指出問題嗎?
我對python非常陌生,並且正在使用它來分析仿真數據。我在Anaconda(python 2.7)中編寫的代碼工作了幾次,現在似乎變得無法響應。Spyder在Win 7上運行的Anaconda(py 2.7)中變得無響應64位
我試過重新啓動控制檯,甚至是卸載並重新安裝。但它不起作用。 有人可以指出問題嗎?
import os
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.backends.backend_pdf as bpdf
import pandas as pd
pdf = bpdf.PdfPages("Dictra_graphs.pdf")
% matplotlib inline
os.getcwd()
path= 'C:\\Users\\prajo25.AD\\Desktop\\D7X sim'
os.chdir(path)
a=os.listdir(os.curdir)
it=7 #number of iterations completed
Time=[1, 5, 10, 30, 50, 80, 100, 150, 300, 500, 800, 1000]
newlist_wp=[]
newlist_bph=[]
file_Al=[]
file_bph=[]
i=0
while i<5:
newpath=path+'\\'+a[i]
os.chdir(newpath)
for filename in os.listdir(os.curdir):
base_file, ext = os.path.splitext(filename)
if ext == '.exp':
os.rename(filename, base_file + '.txt')
List_wt=['wt1.txt','wt2.txt','wt3.txt','wt4.txt','wt5.txt','wt6.txt','wt7.txt']
List_bph=['npm1.txt','npm2.txt','npm3.txt','npm4.txt','npm5.txt','npm6.txt','npm7.txt']
I=range(it+1)
length=len(I)-1
j=0
while j < length:
file_Al.append(np.array(pd.read_table(List_wt[j],sep='\s',names=["Distance in microns","NAN","wt-% Al"], index_col=False, skiprows=17, skip_footer = 671,engine='python'))*1000000)
file_bph.append(np.array(pd.read_table(List_bph[j],sep='\s',names=["Distance in microns","NAN","wt-% Al"], index_col=False, skiprows=152, skip_footer = 68,engine='python'))*1000000)
j=j+1
newlist_wp.append(file_Al)
newlist_bph.append(file_bph)
file_Al=[]
file_bph=[]
i=i+1
i=0
k=0
r=7
while i<7:
list_f=[newlist_wp[0][i], newlist_wp[1][i],newlist_wp[2][i],newlist_wp[3][i],newlist_wp[4][i]]
list_bph=[newlist_bph[0][i], newlist_bph[1][i], newlist_bph[2][i], newlist_bph[3][i], newlist_bph[4][i]]
a0=np.array(list_f[0])
a1=np.array(list_f[1])
a2=np.array(list_f[2])
a3=np.array(list_f[3])
a4=np.array(list_f[4])
b0=np.array(list_bph[0])
b1=np.array(list_bph[1])
b2=np.array(list_bph[2])
b3=np.array(list_bph[3])
b4=np.array(list_bph[4])
plt.figure(k+1)
plt.grid()
plt.plot(a0[:,0],a0[:,2]/1000000,'r--', label='D71')
plt.title(str(Time[i])+'hours')
plt.plot(a1[:,0],a1[:,2]/1000000,'b--', label='D72')
plt.title(str(Time[i])+'hours')
plt.plot(a2[:,0],a2[:,2]/1000000,'m--', label = 'D73')
plt.title(str(Time[i])+'hours')
plt.plot(a3[:,0],a3[:,2]/1000000,'g--', label = 'D74')
plt.title(str(Time[i])+'hours')
plt.plot(a4[:,0],a4[:,2]/1000000,'k--', label = 'D75')
plt.title(str(Time[i])+'hours')
plt.legend(loc=2, ncol=1)
plt.ylabel('Wt.% of Al')
plt.xlabel('Distance in microns')
plt.figure(r+1)
plt.plot(b0[:,0],b0[:,2]/1000000,'r--',label='D71')
plt.title(str(Time[i])+'hours')
plt.plot(b1[:,0],b1[:,2]/1000000,'b--', label='D72')
plt.title(str(Time[i])+'hours')
plt.plot(b2[:,0],b2[:,2]/1000000,'m--', label='D73')
plt.title(str(Time[i])+'hours')
plt.plot(b3[:,0],b3[:,2]/1000000,'g--',label='D74')
plt.title(str(Time[i])+'hours')
plt.plot(b4[:,0],b4[:,2]/1000000,'k--', label='D75')
plt.title(str(Time[i])+'hours')
plt.legend(loc=2, ncol=1)
plt.ylim(ymax=1.2, ymin=0)
plt.ylabel('Amount of gamma phase')
plt.xlabel('Distance in microns')
k=k+1
r=r+1
plt.show()
i=i+1
q=0
for q in range(1,k+r-6): ## will open an empty extra figure :(
pdf.savefig(plt.figure(int(q)))
pdf.close()
歡迎來到Stack Overflow!您能否發佈不適合您的代碼文本? – gariepy
感謝您的回覆。文本超過100行。有沒有好的方法在這裏發佈? – prajo
100行是沒有問題的。將代碼粘貼到答案框中後,確保突出顯示該代碼,然後在格式工具欄中單擊花括號'{}'。這將使它以標準編碼格式顯示。 – gariepy