Question

我想创建一个像这样的匿名函数：

n = 5;
x = linspace(-4,4,1000);

f = @(x,a,b,n) a(1)*exp(b(1)^2*x.^2) + a(2)*exp(b(2)^2*x.^2) + ... a(n)*exp(b(n)^2*x.^2);

我可以这样做，而不传递显式参数n：

f1  = @(x,a,b) a(1)*exp(-b(1)^2*x.^2);
for j = 2:n 
    f1  = @(x,a,b) f1(x,a,b) + a(j)*exp(b(j)^2*x.^2);
end

但似乎有点hacky。有人有更好的解决方案吗？我想知道其他人会如何对待这件事。

Answer 1

你的hacky解决方案肯定不是最好的，因为MATLAB中的递归函数调用效率不高，你可以快速进入最大递归深度（默认为500）。

您可以引入一个新维度，您可以在其中汇总数组a和b。假设x，a和b是行向量：

f = @(x,a,b,n) a(1:n)*exp((b(1:n).^2).'*x.^2)

这将使用第一个维度作为求和维度：(b(1:n).^2).'是一个列向量，当乘以x时会产生一个矩阵（确切地说，这是一个二元乘积）。生成的n * length(x)矩阵可以乘以a(1:n)，因为后者是大小为[1,n]的矩阵。此向量矩阵产品也将为我们执行求和。

迷你证明：

n = 5;
x = linspace(-4,4,1000);
a = rand(1,10);
b = rand(1,10);

y = 0;
for k=1:n
   y = y + a(k)*exp(b(k)^2*x.^2);
end

y2 = a(1:n)*exp((b(1:n).^2).'*x.^2);   %'

all(abs(y-y2))<1e-10

最后一个命令返回1，因此两者基本相同。

在输出向量的匿名函数中创建可变数量的项

1 个答案: