link之家
链接快照平台
- 输入网页链接，自动生成快照
- 标签化管理网页链接
相关文章推荐
|
想出国的饼干 · 【PEL夏季赛总决赛主kv】 ...· 7 月前 · |
|
|
暗恋学妹的佛珠 · 原神mana的画师的画从哪里看 - 抖音· 9 月前 · |
|
|
高大的包子 · 合肥轨道4号线南延线开始铺轨 _光明网· 10 月前 · |
|
|
大气的芒果 · 相守相成·见证者说（4）这里，将成中国自动驾 ...· 11 月前 · |
|
|
傻傻的毛豆 · ubuntu限制syslog文件大小_rsy ...· 12 月前 · |
二维高斯与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);