我的Python编程问题如下:
我想创建一个测量结果数组.每个结果可以描述为正态分布,其中平均值是测量结果本身,标准偏差是其不确定性.
伪代码可以是:
x1 = N(result1, unc1)
x2 = N(result2, unc2)
...
x = array(x1, x2, ..., xN)
比我想计算x 的FFT:
f = numpy.fft.fft(x)
我想要的是x中包含的测量的不确定性是通过FFT计算传播的,因此f是一个幅度阵列及其不确定性,如下所示:
f = (a +/- unc(a), b +/- unc(b), ...)
你能建议我这样做吗？