簡體   English   中英

制表球諧函數

[英]Tabulate Spherical Harmonic Functions

嘿,我正在嘗試以緊湊的方式表達球諧函數。 我不知道如何在表格中表達這一點,我在列表中得到了 l = 0,1,2 和 m = -l, ...,+l 的球面諧波。 但想要一個表格形式:

 m -2 -1 0 1 2
l        
0        x
1      x x x
2    x x x x x

有誰知道如何做到這一點,因為我只是想用普通標題創建更簡單的表格。 我有列表 sph_harmonics[] 中的函數。

from sympy import Ynm
from tabulate import tabulate

sph_harmonics = []
for l in range (0, 3):
            for m in range(-l, l+1):
                print('l =', l)
                print('m =',m)
                print(sp.simplify(Ynm(l, m, theta, phi).expand(func=True)))
                sph_harmonics.append(sp.simplify(Ynm(l, m, theta, phi).expand(func=True)))

Mathjax 本身沒有table / tabular環境,但仍然可以為此使用array環境( link1 ; link2 ; link3 )。

因此,基於此,我編寫了名為fmt_table的以下格式 function ,它像您提到的那樣格式化公式:

from sympy import Ynm, simplify, symbols, latex
theta, phi = symbols('\Theta \phi')
from IPython.display import HTML, display, Math, Latex

def fmt_table(data,center_data=False,add_row_nums=False):
    from math import ceil
    buf='''
\\newcommand\\T{\\Rule{0pt}{1em}{.3em}}
\\begin{array}{%s}    
'''
    max_cols = max(len(r) for r in data)
    column_spec = '|' + '|'.join(['c']*max_cols) + '|'
    buf = buf % column_spec
    row_idx = 0
    for row_data in data:
        row = ''
        if add_row_nums and row_idx > 0:
            row += str(row_idx) + " & "
        if center_data and row_idx > 0:
            to_add = ceil( (max_cols - len(row_data))/2 )
            row += ' & '.join([''] * to_add)
        row += ' & '.join(row_data)
        if row_idx == 0:
            row = '''\\hline''' + row + '''\\\\\hline'''
        else:
            row += '''\\\\\hline'''
        row += "\n"
        buf +=row
        row_idx += 1
    buf += '''\\end{array}'''
    # print(buf)
    return buf
    

tbl = []
tbl.append(['\\texttt{m/l}','-2','-1','0','1','2'])
for l in range (0, 3):
    new_row = []    
    for m in range(-l, l+1):
        h = simplify(Ynm(l, m, theta, phi).expand(func=True))
        new_row.append(latex(h))
    tbl.append(new_row)
    

display(Math(fmt_table(tbl,add_row_nums=True,center_data=True)))

OUTPUT:

在此處輸入圖像描述

此答案中的所有代碼也可在此處獲得

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM