matlab - 关于在 MATLAB 中使用四元组

Question

我在使用 quad 函数时遇到了一个奇怪的问题。 我正在使用 quad 计算简单积分，它工作了 10 到 20 次，然后 Matlab 发出以下错误:

Error using quad (line 75)
The integrand function must return an output vector of the same length as the input vector. 
yteor(k) = quad(@(q)(exp(-(q).^2).*q.^2/(k.^2+1)), 0, 1);

这里 q 和 k 是标量。 我不知道出了什么问题，也不知道为什么它在几个小时前就起作用了。

编辑

这是我的代码

for k=1:100,
    xteor(k)=step*k;
    yteor(k)=quad(@(q)(exp(-(q).^2).*q.^2/((step.*k+1).^2)),0,1);
end plot(xteor,yteor,'r');

Answer 1

以下代码片段适用于 Octave(Matlab GNU 克隆)

step = 1;
xteor = zeros(100,1);
yteor = zeros(100,1);

for k=1:100,
  xteor(k)=step*k;
  yteor(k)=quad(@(q)(exp(-(q).^2).*q.^2/((step.*k+1).^2)),0,1);
end
plot(xteor,yteor,'r');
pause

我的假设是您的错误是代码中较早发生的其他事情的结果(可能与 step 不是标量有关？)。而不是专注于出现错误的这一行。尝试在错误出现之前搜索您所做的更改。