MatLab递归错误（初学者）

Question

好的。 所以我在 MatLab 中有两个互相调用的函数。

Riemann.m

function I = Riemann(f, dx, a, b)
   x = a:dx:b;
   fx = f(x).*dx;
   I = sum(fx);

和 myfunc.m

function f = myfunc(x)
   f = sin(1./x);
   for n=1:100
        I = Riemann(@myfunc, 0.001, 1/n, 1);
   end
   plot(I)

问题是让它运行。 我如何调用 myfunc 来从中获取任何内容。 我尝试过的所有内容最终都会陷入无休止的递归调用堆栈（这是有道理的）。

Answer 1

您的问题在于函数的定义：为了能够使用递归定义，您必须能够计算两个函数中的至少一个而无需另一个，至少对于某些值而言。 您还必须确保每次计算最终都依赖于无需递归即可获得的结果。

对于您的具体问题，我感觉您想要积分函数 f(x)=sin(1./x)。 如果是这样，第二个函数的代码应为：

function f = myfunc(x)
   fct = @(x) sin(1./x);
   f = fct(x);
   for n=1:100
        I = Riemann(fct, 0.001, 1/n, 1);
   end
   plot(I)

Answer 2

函数 myfunc 没有在 f = sin(1./x); 之后结束。 在那里终止该函数并从其他地方（单独的文件）调用绘图代码。

来自手册：
您可以使用 end 语句终止任何函数，但在大多数情况下，这是可选的。 仅在使用一个或多个嵌套函数的 M 文件中才需要 end 语句。 在这样的 M 文件中，每个函数（包括主函数、嵌套函数、私有函数和子函数）都必须以 end 语句终止。 您可以使用 end 终止任何函数类型，但除非 M 文件包含嵌套函数，否则不需要这样做。

Answer 3

当 myFunc 停止调用黎曼时，您需要输入 x 的最终条件。 将实际函数（在本例中为 sin）发送给 Riemann 比调用 myFunc 更好。