link管理
链接快照平台
- 输入网页链接，自动生成快照
- 标签化管理网页链接
相关文章推荐
|
飘逸的四季豆 · js--数组的filter()过滤方法的使用 ...· 4 月前 · |
|
|
豪情万千的热水瓶 · Can't Create new ...· 4 月前 · |
|
|
幸福的手术刀 · 在r语言中，如何使用支持向量机对多分类结局的 ...· 5 月前 · |
|
|
温文尔雅的生姜 · 【Python】解决pandas读取exce ...· 5 月前 · |
|
|
失望的篮球 · 后徐明时代实德背后关键人现身：新华人寿前董事 ...· 6 月前 · |
二维高斯与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);
推荐文章