根据第二列的值计算表达式(方程式)的 Python 程序

Question

我的值(value)观是:

0.263 0

0.265 0

0.267 0

0.269 0.0001

0.271 0.0003

0.273 0.0006

0.275 0.0011

0.277 0.0021

0.279 0.0029

0.281 0.0046

0.283 0.0072

0.285 0.0113

我想计算表达式，即 PMF(W_r)= k_b T ln g(r)。其中 g(r) 是第二列的值。

这是我试过的代码:

import numpy as np
#import panda as pd
import scipy.constants as sc
#from astropy import constants as const
import matplotlib.pyplot as plt
import math
A=open('rdf_CaOw.dat','r')
B=open('pmf.dat','w')
for column in A:
    c=column.strip().split()
    B.write(column[6:11]+'\n')
B.close()
A.close()

C=open('pmf.dat', 'r')
D=open('pmf1.dat','w')
for line in C:
    W = (- float(sc.Boltzmann * 298 * float (math.log (C))))
    print (W)

1) 对此有何建议？

2) 我可以直接取第二列插入/放入方程来计算吗？

我该怎么做？因为在这里我将第二列写入一个输出文件，然后再次从该输出文件输入方程式。

非常感谢任何建议。

Answer 1

可以直接拿第二列

with open('rdf.dat') as f:
    for line in f:
        cols = [float(x) for x in line.split()]
        print (cols[1])

输出

0.0
0.0
0.0
0.0001
0.0003
0.0006
0.0011
0.0021
0.0029
0.0046
0.0072
0.0113

您现在可以将其用于进一步的数学运算。