嗨,杰夫,
很高兴看到我的贡献几乎从四年前,仅作为众多人对这个网站的贡献并得到引用。
除了提供很多查询的答案,网站似乎已经成为一个可搜索的信息来源。
之所以第一conv示例工作和第二并不是因为不同的x数组用于创建pdf文档。
假设fft变换y = y (x) z = z (ilam)
(ilam逆波长,它似乎没有常见的象征)。
fft将z (ilam = 0),直流值,产生的z数组中的第一个元素。
同样的fft
假设
y (x = 0)是在y数组第一个元素。
如果x = 0是别的地方,那么z ilam域中的拿起一个相移,会导致不正确的结果。
第一个例子使用了正确的x数组和方法。
第二个例子从一开始适应高斯- x的值,导致问题,正如你所看到的。
最直接的方法,允许负的x是开始一个数组在中点与x = 0(元素n / 2 + 1 n)。
定义y (x),然后使用ifftshift把直流y值y (1)。
我倾向于认为第一个数组作为策划组和移位的数组作为真正的数组。
经过fft算法,然后回来在x域的传输线,fftshift到把y (x = 0)的中心。
例如,
N = 10000;
x0 = 5;
delx = 2 * x0 / N;
delilam = 1 / (N * delx);
x = (- N / 2: N / 2 - 1) * delx;
μ= 1;
sig = 2;
m = 6;%的旋转
mp1 = m + 1;
y = (1 / (sqrt(π)* sig)) * exp (- (xμ)。^ 2 /团体^ 2);
%的结果
yconv =(1 /(√(π* mp1) * sig)) * exp (- (x-mp1 *μ)。^ 2 / (mp1 *团体^ 2));
%由fft
z = fft (ifftshift (y)) ^ mp1 * delx ^ mp1;
yconvfft = fftshift(传输线(z)) * N * delilam;
%以下情节原高斯缩减
图(1);网格在
情节(x, y / 2 x, yconvfft x(一60:结束),yconv(一60:结束),“o”)
I1 = trapz (x, y)
I2 = trapz (x, yconv)
I3 = trapz (x, yconvfft)
正常化,因为您正在使用fft(不仅资金)近似积分,有一个因素的阵列间距del_x每个变换。
回来有一个传输线有del_ilam的单一因素。
除了传输线自动划分的数量分N,所以你需要乘以N来弥补。
我有一个问题。在第二个例子中你有2 ^ 13 = 8192分,而不是说,10000。你能让我知道背后的基本原理用2 ^ n点吗?
