我是StackOverflow的新手,我对Python非常陌生.
我的问题是这个……我需要写一个双和,如下：
动机是这是对用于大地水准面的重力势的角度修正.
>我写这笔钱很困难.并且,在你说“去这样一个这样的资源”或者对我不耐烦之前,这是我第一次完成编码/编程/不管这是什么.
>这是一个使用“for”循环的好地方吗？
>我有两个索引(n,m)的数据以及.txt文件中的系数c_ {nm}和s_ {nm}的数据.每个项目都是一列.当我说usecols时,我会将它们编号为0到3,还是1到4？
(上面的等式)
\begin{equation}
V(r,\phi,\lambda) = \sum_{n=2}^{360}\left(\frac{a}{r}\right)^{n}\sum_{m=0}^{n}\left[c_{nm}*\cos{(m\lambda)} + s_{nm}*\sin{(m\lambda)}\right]*\sqrt{\frac{(n-m)!}{(n+m)!}(2n + 1)(2 - \delta_{m0})}P_{nm}(\sin{\lambda})
\end{equation}
最佳答案
(2)是的,“for”循环没问题.正如@jpmc26指出的那样,生成器表达式是“for”循环的一个很好的替代品. IMO,如果效率对你很重要,你会想要使用numpy.
(3)正如@askewchan所说,“usecols”是指genfromtxt的参数;如该文档中所指定,列索引从0开始,因此您将要使用0到3.
一个天真的实现可能是好的,因为较大的阶乘是分母,但如果遇到数字问题我不会感到惊讶.这是让你入门的东西.请注意,您需要定义P()和a.我不明白“0到3”是如何与c和s相关的,因为它们的索引范围更远.我将假设每个(和delta)都有自己的值文件.
import math
import numpy
c = numpy.getfromtxt("the_c_file.txt")
s = numpy.getfromtxt("the_s_file.txt")
delta = numpy.getfromtxt("the_delta_file.txt")
def V(r,phi,lam):
ret = 0
for n in xrange(2,361):
for m in xrange(0,n + 1):
inner = c[n,m]*math.cos(m*lam) + s[n,m]*math.sin(m*lam)
inner *= math.sqrt(math.factorial(n-m)/math.factorial(n+m)*(2*n+1)*(2-delta[m,0]))
inner *= P(n,m,math.sin(lam))
ret += math.pow(a/r,n) * inner
return ret
一定要写unittests来检查数学.请注意,“lambda”是保留字.