有趣的问题.这是一个建议,这可能是相当困难的实现,取决于你有多少数学倾向. @H_404_2@请注意,我交换空间的速度,因为我建议可能是相当沉重的计算(但这是要针对真实数据进行测试).

@H_404_2@首先,使用功能性方法.概率分布是概率测度：

struct Distribution
{
    virtual ~Distribution() {};
    virtual double integrate(std::function<double(double)>) = 0;
};

@H_404_2@这样,您就可以从生成的样品中抽出,因为您不想存储它们.告诉自己,您可以使用“集成”方法做任何事情.

@H_404_2@当然,通过明确的样本,你会做一些类似的事情

struct SampledDistribution
{
    double integrate(std::function<double(double)> f)
    {
        double acc = 0;
        for (double x: samples) acc += f(samples);
        return acc / samples.size();
    }

    std::deque<double> samples;
};

@H_404_2@现在,存储部分：

@H_404_2@The usual representation — a histogram array — is undesirable mainly
because of the trade-off between quantization/resolution and space. I
imagine there must be a method of representation that adaptively
varies the bin size depending on local “complexity”.

@H_404_2@传统的方法是wavelets.你可以产生一个调用集成的系数,你可以串行化.如果它们产生的积分估计器的方差较大,则将系数除以.

@H_404_2@然后,为了反序列化,你生成一个Distribution对象,其集成方法执行与小波的集成.可以使用您最喜欢的正交方法进行集成.由于实际执行取决于您选择的小波系列(平滑,紧凑支持,正交或不等等),因此我在这里保持故意模糊.无论如何,您将需要潜入细节.

@H_404_2@这里的一点是,通常只需要很少的小波来表示具有很少特征的平滑功能(例如,几个峰,并且另外定期成形),不像更多的“常规”有限元(直方图是特定类型的有限元表示) .小波表示适应变形体的特征,无论其位置或大小.您还可以决定要保留多少系数,从而控制压缩比.

@H_404_2@另外,0.001的数字是相当高的数字：我怀疑你只需要少量的系数

@H_404_2@权衡就在于您使用的小波类：非常平滑的分布可能会很好地用平滑小波表示,但紧凑支持的小波可能更容易集成在一起等等.请注意,您不需要“小波变换”软件包：只有显式表示小波函数和正交例程(尝试Gauss-XXX程序进行重建,或者是高阶).

@H_404_2@我喜欢在傅立叶域中定义的小波(如Lemarie小波),因为它们在傅立叶空间中的值和导数在零是已知的,这允许您执行分布上的约束：概率测量必须集成到一个,您可能会事先知道预期的价值或差异.

@H_404_2@此外,您可能希望将变量更改为仅处理函数,例如.在[0,1]上.在间隔上有一个大小的小波.