来自其他函数的输入元素 1 by 1 w/o for 循环

Question

我有一个函数可以计算一个数字数组（randparam），我想将它逐个元素地输入到另一个进行模拟的函数中。

例如

function [randparam] = arraycode

code 

randparam = results of code 
% randparam is now a 1x1001 vector. 

end

接下来我想将 randparam 1 by 1 输入到我的模拟函数中

function simulation 

x = constant + constant * randparam + constant

return 
end

让我觉得困难的是因为模拟函数中的return命令，它只计算上面方程x的一步，将结果返回到另一个函数，调用它积分器，然后积分器函数会调用模拟再次使用函数来计算 x。

所以积分函数可能看起来像

function integrator (x)
y = simulation(x) * 5
u = y+10
yy = simulation(x) * 10 + u

end

如您所见，积分函数调用模拟函数两次。这给我带来了两个问题：

1. 如果我在模拟函数中创建一个 for 循环，在该循环中我使用类似的方式逐个元素地输入：

   对于 i = 1:100

   x = 常数 + 常数 * randparam(i) + 常数

   返回
   结束

然后每次我的积分器函数再次调用我的模拟函数时，我的 for 循环都会从 1 重新开始。

2.如果我以某种方式将 i 保留在我的基础工作区中，以便我的模拟函数中的 for 循环知道从 1 开始，那么我的 y 和 yy 函数将具有不同的 x 输入，因为只要它会第二次为 yy 打电话，那么由于 y 的电话，我现在是 i+1。

在这种情况下有没有办法避免 for 循环？第二个问题的一个潜在解决方案是复制脚本但使用不同的名称，并让我的 for 循环使用不同的变量，但这似乎效率很低。

希望我说清楚了。

谢谢。

Answer 1

首先，如果您通常希望对数组的每个元素应用相同的函数，并且还没有内置的矢量化方法，您可以 使用arrayfun（虽然通常一个简单的for 循环更快、更易读）：

%# randparam is a 1x1001 vector. 
%#next I want to input randparam 1 by 1 into my simulation function
function simulation 
    x = constant + constant * randparam + constant
    return 
end

（注意：问问自己这个函数可能在做什么，因为它没有返回值，而且 MATLAB 也没有通过引用传递。）这就是 arrayfun 的用途：将函数应用于一个数组（或向量，在这种情况下）。同样，您应该确保在您的情况下这样做是有意义的，而不是显式循环。

function simulation(input_val)
#% your stuff
end


sim_results = arrayfun( @simulation, randparam);

当然，你写的方式，行

x = constant + constant*randparam + constant;

可以（并且将）进行矢量化 - 如果你给它一个矢量或矩阵，一个矢量或矩阵将是结果。

第二您似乎对MATLAB中函数变量的“范围”不清楚。如果调用函数，则会创建一个干净的工作区。因此，来自一个函数的x 在您调用的另一个函数中不会自动可用。变量也会在函数末尾超出范围，因此在函数内使用 x 不会更改/覆盖该函数外部存在的变量 x。并且对一个函数的多次调用每个都有自己的工作区。

Answer 2

积分器级别的循环有什么问题？

function integrator (x)

for i=1:length(x)
    y = simulation(x(i)) * 5
    u = y+10
    yy = simulation(x(i)) * 10 + u
end

然后将您的整个 randparm 传递给积分器？从您的问题中不清楚您是否希望 simulation 在给定相同的输入时返回相同的值，或者您是否希望它使用相同的输入执行两次，或者您是否希望每次调用都有新的输入。也不清楚simulation 是否保持任何内部状态。根据您编写示例的方式，模拟仅取决于输入值，而不取决于任何先前的输入或输出，这使得矢量化变得微不足道。如果我们都错过了这条船，请使用更精致的示例代码编辑您的问题。