|
|
相关文章推荐
|
淡定的猕猴桃 · JavaScript字符串MD5 - ...· 8 月前 · |
|
|
乖乖的芹菜 · 用atob解码前端的JWT令牌会使我在“窗口 ...· 8 月前 · |
|
|
曾深爱过的手链 · PHP 将HTML转成PDF文件,HTML ...· 10 月前 · |
|
|
博学的圣诞树 · 使用 JavaScript ...· 10 月前 · |
|
|
豪气的煎鸡蛋 · 为什么在python中for-range比w ...· 1 年前 · |
二维高斯与matlab？
问
二维高斯与matlab？
EN
EN
Stack Overflow用户
提问于
2022-04-17 19:35:46
回答 1
查看 1.2K
关注 0
票数 1
我想用matlab绘制一个以矩阵中间为中心的二维高斯。现在我在和
代码语言：
javascript
复制
for i = 1:size2
for j = 1:size2
gauss(i,j)=A*exp(-1/(sigma^2)*((i-round(size2/2))^2+(j-round(size2/2))^2)); %gaussiana
end但它可能会非常慢。
在matlab中我没有找到任何函数，最快的方法是什么？mex文件是一种有价值的选择吗？
EN
回答 1
Stack Overflow用户
回答已采纳
发布于 2022-04-17 20:44:16
用 网格 对高斯计算进行矢量化比较简单。
在 维基百科 中甚至有一个八度实现。
执行工作是：
代码语言：
javascript
复制
[X, Y] = meshgrid(1:size2, 1:size2);
G = A*exp(-1/(sigma^2)*((Y-size2/2).^2 + (X-size2/2).^2));
注意:我用
round(size2/2)
替换了
size2/2
，因为舍入不是公式的一部分。
将一维高斯变换为二维高斯 的技巧
为了使计算速度稍微快一点，我们可以创建一维高斯，然后用它计算2D高斯：
代码语言：
javascript
复制
x = 1:size2;
G1 = sqrt(A)*exp(-1/(sigma^2)*(x-size2/2).^2); % Create 1D gaussian
G2 = G1'*G1; % Compute the 2D gaussian out of 1D gaussian. 测试：
代码语言：
javascript
复制
size2 = 101;
A = 10;
sigma = 50;
gauss = zeros(size2);
for y = 1:size2
for x = 1:size2
gauss(y,x)=A*exp(-1/(sigma^2)*((y-size2/2)^2+(x-size2/2)^2)); %gaussiana
[X, Y] = meshgrid(1:size2, 1:size2);
G = A*exp(-1/(sigma^2)*((Y-size2/2).^2 + (X-size2/2).^2));
disp(['sum abs diff (G, gauss) = ' num2str(sum(abs(G(:) - gauss(:))))])
x = 1:size2;
G1 = sqrt(A)*exp(-1/(sigma^2)*(x-size2/2).^2);
推荐文章