如何解决Matlab中的显式IFFT2
我有一个2D矩阵,希望将其与一个明确的导数一起使用,因此具有矩阵的功能形式可能会有所帮助。我以为执行2D FFT,然后以傅立叶基础表示矩阵是个好主意,但我遇到了一些困难。我想如果给定x,y meshgrid的参数和在meshgrid上的矩阵f(是x和y的函数)的参数,则该函数将创建函数f2,使得f2(x(ii(j,jj)) y(ii,jj))= f(ii,jj)。
xc = -10;
yc = -10;
d = 0.1;
[x,y] = meshgrid(xc:d:10,yc:d:10);
f = sin(x + y);
n = size(f,2);
m = size(f,1);
coefs = fft2(f);
omega_x = exp((2*pi.*1i./n).*(0:(n-1)));
omega_y = exp((2*pi.*1i./m).*(0:(m-1)));
omega_y = omega_y';
func = @(x,y) (1/m).*(1/n).*sum(sum((omega_x.^((x-xc)./d).*(omega_y.^((y-yc)./d))).*coefs));
f2 = zeros(size(f));
for ii = 1:size(f,1)
for jj = 1:size(f,2)
f2(ii,jj) = func(x(ii,jj),y(ii,jj));
end
end
我们应该让矩阵f2和f具有可比性。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。