Question

matlab中带参数定积分的求解

程序和结果如下：
>> clear>> f=@(t,r) (142884+1.21*r.*r-831.6*r.*cos(2*t)).^(1/2);>> r=0:0.1:75;>> q=quad(f,0,90,r)Error using @(t,r)(142884+1.21*r.*r-831.6*r.*cos(2*t)).^(1/2)Not enough input arguments.
Error in quad (line 72)y = f(x, varargin{:});

求指导...万分感谢！

Answer 1

1、首先，如果要对quad函数传递附加参数（Extra Parameters），应该先搞清楚quad函数的调用格式：

q = quad(fun,a,b,tol,trace)

其中，第4、5个参数分别是tol和trace，不能想当然。你如果要传递附加参数，应该从第6个参数开始——尽管文档中并没有提及，但quad函数的确可以接受附加参数。类似地函数还有很多，例如fsolve、fmincon等函数的函数。

 

在传递附加参数时，需要把前面的参数补齐，即使你不知道也不关心那些参数到底是干什么用的。幸运的是，这类函数绝大多数（如果不是所有的话）都支持用空数组作为占位符，例如：

q = quad(fun,a,b,[],[],r)

这样，你也就可以不用管tol、trace到底该怎样取值了。

 

2、其次，quad函数并不支持一次性对多个附加参数求解。解决的办法一是用循环，另一是用arrayfun函数（其实后者本质上也是循环，只不过写法上简单些）。

 

参考代码

f=#(t,r) 142884+1.21*r.*r-831.6*r.*cos(2*t)).^(1/2);
r=0:0.1:75;
q=arrayfun(@(R)quad(f,0,90,[],[],R),r);
plot(r,q)

顺便提醒一句，求积分的计算量不小，r的间距没必要取那么小（取5得到的结果就挺平滑）。

 

PS: 上面代码中的#应为电子邮件的@，不知怎么回事贴不出来，严重鄙视百度知道开发人员这群猪。。。