问题描述
我有一个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 (将#修改为@)