Python中黎曼积分字典键截断及DataFrame列序混乱问题求助

嘿，我来帮你搞定这个问题！你遇到的两个核心问题：一是小浮点数作为字典键时被自动转为科学计数法字符串，二是转DataFrame后列序混乱，其实都可以通过统一键的格式化方式和强制指定列顺序来解决，咱们一步步来：

问题根源拆解 #

1. 浮点数的字符串表示歧义：Python对极小的浮点数（比如0.00001）默认会用科学计数法（1e-05）显示，不管用str()还是repr()都没法改变这个默认行为，这就导致你的字典键出现了不符合预期的格式。
2. 列序混乱：即使字典在Python 3.7+是插入有序的，当你转DataFrame时，如果键是字符串，Pandas会按字符串的字典序排序（比如"0.00001"会和"0.001"的顺序不符合数值逻辑），这就需要手动指定列顺序来匹配你定义的h_range。

解决方案代码 #

方法一：统一格式化键为清晰的小数字符串 #

我们可以用format()函数强制把浮点数转为固定格式的小数字符串，避免科学计数法，同时在转DataFrame时指定列顺序：

import numpy as np
import pandas as pd

def myfunc(x, mu, sigma):
    px = np.exp(-(x-mu)**2/(2*sigma**2))
    return px

def get_area(h,mu,sigma):
    x = np.arange(-100,100+h,h)
    return sum([myfunc(xi,mu,sigma)*h for xi in x])

# 定义生成统一h键的函数
def get_h_key(h):
    # 格式化浮点数为小数形式，去掉末尾多余的0和小数点
    return format(h, '.10f').rstrip('0').rstrip('.')

sigma_range = [0.25,0.5,1,2]
h_range = [2,1,0.1,0.001,0.00001]
result_dict = {}

for sigma in sigma_range:
    sigma_dict = {}
    for h in h_range:
        h_key = get_h_key(h)
        # 直接赋值，每个h对应一个面积值（无需append，因为每个h只循环一次）
        sigma_dict[h_key] = get_area(h, 1, sigma)
    result_dict[str(sigma)] = sigma_dict

# 转DataFrame时指定列顺序，和h_range完全一致
columns = [get_h_key(h) for h in h_range]
df = pd.DataFrame.from_dict(result_dict, orient='index', columns=columns)
print(df)

方法二：用数值作为字典键，最后格式化列名 #

这种方法更简洁，直接保留h的数值类型作为字典键，避免字符串表示的问题，转DataFrame后再把列名改成友好的格式：

import numpy as np
import pandas as pd

def myfunc(x, mu, sigma):
    px = np.exp(-(x-mu)**2/(2*sigma**2))
    return px

def get_area(h,mu,sigma):
    x = np.arange(-100,100+h,h)
    return sum([myfunc(xi,mu,sigma)*h for xi in x])

sigma_range = [0.25,0.5,1,2]
h_range = [2,1,0.1,0.001,0.00001]
result_dict = {}

for sigma in sigma_range:
    sigma_dict = {}
    for h in h_range:
        # 直接用h的数值作为键
        sigma_dict[h] = get_area(h, 1, sigma)
    result_dict[sigma] = sigma_dict

# 转DataFrame时强制指定列顺序为h_range的原始顺序
df = pd.DataFrame.from_dict(result_dict, orient='index', columns=h_range)
# 把列名格式化为清晰的小数字符串
df.columns = [format(h, '.10f').rstrip('0').rstrip('.') for h in df.columns]
print(df)

额外小提示 #

你原来的代码里有个小细节：sigma_dict[str(repr(h))].append(get_area(h,1,sigma))其实不需要用append()，因为每个h在sigma的循环里只会出现一次，直接赋值就能得到对应的值，append反而会导致重复添加（如果不小心重复循环h的话）。

另外，关于你的黎曼积分实现：高斯函数的积分其实有解析解sigma*sqrt(2π)，你可以用这个结果来验证数值积分的准确性哦！

内容的提问来源于stack exchange，提问作者Michal B.