- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在尝试构建二维框函数的 FFT 和二维高斯函数的 FFT 的乘积。之后,我找到了逆 FFT,结果我期待这两个函数的卷积。但是,我得到了一个奇怪的单方面结果,如下所示。结果显示在 subplot
的右下角。
我为重现上述 subplot
而编写的 Octave 代码以及我为构建卷积而执行的计算如下所示。谁能告诉我我做错了什么?
clear all;
clc;
close all;
% domain on each side is 0-9
L = 10;
% num subdivisions
N = 32;
delta=L/N;
sigma = 0.5;
% get the domain ready
[x,y] = meshgrid((0:N-1)*delta);
% since domain ranges from 0-(N-1) on both sdes
% we need to take the average
xAvg = sum(x(1, :))/length(x(1,:));
yAvg = sum(y(:, 1))/length(x(:,1));
% gaussian
gssn = exp(- ((x - xAvg) .^ 2 + (y - yAvg) .^ 2) ./ (2*sigma^2));
function ret = boxImpulse(a,b)
n = 32;
L=10;
delta = L/n;
nL = ((n-1)/2-3)*delta;
nU = ((n-1)/2+3)*delta;
if ((a >= nL) && (a <= nU) && ( b >= nL) && (b <= nU) )
ret=1;
else
ret=0;
end
ret;
endfunction
boxResponse = arrayfun(@boxImpulse, x, y);
subplot(2,2,1);mesh(x,y,gssn); title("gaussian fun");
subplot(2,2,2);mesh(x,y, abs(fft2(gssn)) .^2); title("fft of gaussian");
subplot(2,2,3);mesh(x,y,boxResponse); title("box fun");
inv_of_product_of_ffts = abs(ifft2(fft2(boxResponse) * fft2(gssn))) .^2 ;
subplot(2,2,4);mesh(x,y,inv_of_product_of_ffts); title("inv of product of fft");
最佳答案
让我们首先解决您最明显的错误。这是您在频域中计算卷积的方式:
inv_of_product_of_ffts = abs(ifft2(fft2(boxResponse) * fft2(gssn))) .^2 ;
第一个问题是您使用的是*
,即矩阵乘法。在频域中,element-wise乘法相当于频域中的卷积,所以需要用.*
代替。第二个问题是您也不需要 .^2
项和 abs
操作。您正在执行卷积,然后找到绝对值并对每个项进行平方。这不是必需的。删除 abs
和 .^2
操作。
因此,你需要的是:
inv_of_product_of_ffts = ifft2(fft2(boxResponse).*fft2(gssn)));
但是,您将要得到的是这个结果。让我们把它放在一个新图中而不是 subplot
*:
figure;
mesh(x,y,ifft2(fft2(boxResponse).*fft2(gssn)));
title('Convolution... not right though');
您可以看到这是正确的结果……但它没有居中……这是为什么?这实际上是在频域中计算卷积时最常见的问题之一。事实上,即使是最有经验的人也会遇到这个问题,因为他们不了解 FFT 工作原理的内部结构。
这是 MATLAB 和 Octave 如何计算 2D FFT 的结果。具体来说,MATLAB 和 Octave 定义了一个 meshgrid
坐标,水平坐标从 0,1,...M-1
到 0,1,.. .N-1
表示垂直,给我们一个 M x N
结果。这意味着原点/直流分量位于矩阵的左上角,而不是我们通常定义的中心。具体来说,传统的 2D FFT 定义坐标从 -(M-1)/2, ..., (M-1)/2
用于水平和 -(N-1)/2, ..., (N-1)/2
垂直。
引用上述内容,您将信号定义为中心,假设它是原点而不是左上角。为了弥补这一点,您需要添加一个 fftshift
(MATLAB doc,Octave doc),以便 FFT 的输出现在以原点为中心,而不是左上角。回到原来的样子,因此你真的需要:
figure;
mesh(x,y,fftshift(ifft2(fft2(boxResponse).*fft2(gssn))));
title('Convolution... now it is right');
如果您想仔细检查我们的结果是否正确,您可以在空间域中执行直接卷积,我们可以在频域中比较两种方法的结果。
这就是您直接在空间域中计算结果的方式:
figure;
mesh(x,y,conv2(gssn,boxResponse,'same'));
title('Convolution... spatial domain');
conv2
( MATLAB doc , Octave doc ) 在两个信号之间执行二维卷积,'same'
标志确保输出大小是最大的两个信号,可以是高斯滤波器或箱式滤波器。您会看到它是同一条曲线,为简洁起见,我不会在此处显示。
但是,我们可以比较这两个结果,看看它们在元素方面是否相同。执行此操作的一种方法是确定减去结果中的每个元素是否小于某个阈值……比如说……1e-10
:
>> out1 = conv2(boxResponse, gssn, 'same');
>> out2 = fftshift(ifft2(fft2(boxResponse).*fft2(gssn)));
>> all(abs(out1(:)-out2(:)) < 1e-10)
ans =
1
这意味着这两者确实是一样的。
希望这是有道理的!请记住,由于您是在原点位于中心的位置定义信号,一旦找到逆向 FFT,就必须移动,以便原点现在位于中心,而不是顶部 -左角。
*:小注 - 此答案中生成的所有图都是使用 MATLAB R2015a 生成的。但是,该代码在 Octave 中完全有效 - 使用 Octave 4.0.0 进行了测试。
关于matlab - 通过 FFT 执行 2D 卷积时的奇怪行为,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/34081933/
SQLite、Content provider 和 Shared Preference 之间的所有已知区别。 但我想知道什么时候需要根据情况使用 SQLite 或 Content Provider 或
警告:我正在使用一个我无法完全控制的后端,所以我正在努力解决 Backbone 中的一些注意事项,这些注意事项可能在其他地方更好地解决......不幸的是,我别无选择,只能在这里处理它们! 所以,我的
我一整天都在挣扎。我的预输入搜索表达式与远程 json 数据完美配合。但是当我尝试使用相同的 json 数据作为预取数据时,建议为空。点击第一个标志后,我收到预定义消息“无法找到任何内容...”,结果
我正在制作一个模拟 NHL 选秀彩票的程序,其中屏幕右侧应该有一个 JTextField,并且在左侧绘制弹跳的选秀球。我创建了一个名为 Ball 的类,它实现了 Runnable,并在我的主 Draf
这个问题已经有答案了: How can I calculate a time span in Java and format the output? (18 个回答) 已关闭 9 年前。 这是我的代码
我有一个 ASP.NET Web API 应用程序在我的本地 IIS 实例上运行。 Web 应用程序配置有 CORS。我调用的 Web API 方法类似于: [POST("/API/{foo}/{ba
我将用户输入的时间和日期作为: DatePicker dp = (DatePicker) findViewById(R.id.datePicker); TimePicker tp = (TimePic
放宽“邻居”的标准是否足够,或者是否有其他标准行动可以采取? 最佳答案 如果所有相邻解决方案都是 Tabu,则听起来您的 Tabu 列表的大小太长或您的释放策略太严格。一个好的 Tabu 列表长度是
我正在阅读来自 cppreference 的代码示例: #include #include #include #include template void print_queue(T& q)
我快疯了,我试图理解工具提示的行为,但没有成功。 1. 第一个问题是当我尝试通过插件(按钮 1)在点击事件中使用它时 -> 如果您转到 Fiddle,您会在“内容”内看到该函数' 每次点击都会调用该属
我在功能组件中有以下代码: const [ folder, setFolder ] = useState([]); const folderData = useContext(FolderContex
我在使用预签名网址和 AFNetworking 3.0 从 S3 获取图像时遇到问题。我可以使用 NSMutableURLRequest 和 NSURLSession 获取图像,但是当我使用 AFHT
我正在使用 Oracle ojdbc 12 和 Java 8 处理 Oracle UCP 管理器的问题。当 UCP 池启动失败时,我希望关闭它创建的连接。 当池初始化期间遇到 ORA-02391:超过
关闭。此题需要details or clarity 。目前不接受答案。 想要改进这个问题吗?通过 editing this post 添加详细信息并澄清问题. 已关闭 9 年前。 Improve
引用这个plunker: https://plnkr.co/edit/GWsbdDWVvBYNMqyxzlLY?p=preview 我在 styles.css 文件和 src/app.ts 文件中指定
为什么我的条形这么细?我尝试将宽度设置为 1,它们变得非常厚。我不知道还能尝试什么。默认厚度为 0.8,这是应该的样子吗? import matplotlib.pyplot as plt import
当我编写时,查询按预期执行: SELECT id, day2.count - day1.count AS diff FROM day1 NATURAL JOIN day2; 但我真正想要的是右连接。当
我有以下时间数据: 0 08/01/16 13:07:46,335437 1 18/02/16 08:40:40,565575 2 14/01/16 22:2
一些背景知识 -我的 NodeJS 服务器在端口 3001 上运行,我的 React 应用程序在端口 3000 上运行。我在 React 应用程序 package.json 中设置了一个代理来代理对端
我面临着一个愚蠢的问题。我试图在我的 Angular 应用程序中延迟加载我的图像,我已经尝试过这个2: 但是他们都设置了 src attr 而不是 data-src,我在这里遗漏了什么吗?保留 d
我是一名优秀的程序员,十分优秀!