3

我需要对一组 2D 数组(如果需要,可以使用 3D 矩阵)执行两个参数(一组:t,另一个来自数组)的指数运算。f(t,x) = exp(tx) 然后我需要将第三维中每个值的结果相加。因为执行整个操作需要花费太多时间bsxfun,所以我正在考虑使用查找表。

我可以将表创建为矩阵LUT(由于两个参数,二维),然后我可以使用LUT(par1,par2). 但是使用循环访问第三维也很昂贵。

我的问题是:有没有办法实现这种机制(查找表)以具有预定义的值,然后只使用它们从矩阵元素(索引类型)访问而无需循环。或者,如何创建 MATLAB 自动处理的查找表以加快指数运算?

编辑:我实际上使用类似的方法来创建 LUT。现在,我的问题实际上是如何以有效的方式访问它。

可以说我有一个二维数组M。使用那些我想将函数f(t,M(i,j))应用于固定值的值t。我可以使用循环遍历 M 的所有值 (i,j)。但我想要一种更快的方法,因为我有一组 M,然后我需要将此过程应用于所有其他值.

我的函数比我给出的例子有点复杂:

pr = mean(exp(-bsxfun(@rdivide,bsxfun(@minus,color_vals,double(I)).^2,m)./2),3);

这是我的实际功能,正如您所看到的,它比我提供的示例更复杂。但想法是一样的。它在 M 集合的第三维中求两个数组之差的指数的平均值。

希望有帮助。

4

2 回答 2

5

我同意这个问题不是很清楚,并且显示一些代码会有所帮助。反正我会试试的。

为了让 LUT 完全有意义,t-x必须限制 by 获得的一组值,例如限制为整数。

假设指数可以是 -1000 到 1000 之间的任何整数,您可以像这样创建 LUT:

LUT = exp(-1000:1000);

然后创建索引(假设t是一维数组,x 是二维数组)

indexArray = bsxfun(@minus,reshape(t,[1,1,3]), x) + 1001; %# -1000 turns into 1

最后,你创建你的结果

output = LUT(indexArray);
%# sum along third dimension (i.e. sum over all `t`)
output = sum(output,3);
于 2010-10-15T20:18:16.343 回答
0

我不确定我是否理解你的问题,但我认为这就是答案。

x = 0:3
y = 0:2
z = 0:6

[X,Y,Z] = meshgrid(x,y,z)

LUT = (X+Y).^Z
于 2010-10-15T13:43:19.813 回答