实值球谐函数表示结果不准确，调整系数数量无明显改善，问题出在哪？

错误原因及修复方案 #

你代码的问题主要出在以下三个地方：

1. 实球谐转换重复添加Condon-Shortley相位 #

Scipy自带的sph_harm函数默认已经包含了Condon-Shortley相位约定，你在转换实球谐时额外乘以(-1)**m相当于重复叠加了相位，直接导致基函数数值错误，系数计算精度失效。
修复方案：移除get_sph_harm_function方法里多余的相位项：

def get_sph_harm_function(self, l, m):
    def basis_function(theta, phi):
        Y = sph_harm(abs(m), l, phi, theta)
        
        if m < 0:
            Y = np.sqrt(2) * Y.imag
        elif m > 0:
            Y = np.sqrt(2) * Y.real
            
        return Y.real
    
    return basis_function

2. 测试用球坐标范围非法 #

球坐标的极角θ的取值范围是[0, π]，你在测试代码中将θ和方位角φ的范围都设为[0, 2π]，有一半的测试点不属于单位球面的合法坐标，属于无意义的取值区间，计算出来的误差自然会非常大。
修复方案：调整测试点生成逻辑，限制θ的范围：

theta_pts = np.linspace(0, np.pi, 50)
phi_pts = np.linspace(0, 2*np.pi, 50)
pts_2d = np.reshape(np.stack(np.meshgrid(theta_pts, phi_pts, indexing='ij'), axis=-1), (-1, 2))

3. 积分精度不足 #

默认的dblquad积分精度对于高阶球谐系数的计算可能不够，会引入额外的数值误差。
修复方案：调整积分精度参数：

def sph_harm_coeff(self, Y, f):
    def integrand(phi, theta):
        return f(theta, phi) * Y(theta, phi) * np.sin(theta)
    
    return integrate.dblquad(integrand, 0., np.pi, lambda x : 0., lambda x : 2*np.pi, epsabs=1e-8, epsrel=1e-8)[0]

验证效果 #

完成上述修改后，25个系数下的平均绝对误差会降到1e-3量级以下，调整系数数量可以进一步降低误差，符合球谐展开的收敛特性。

内容的提问来源于stack exchange，提问作者laura