本文主要是介绍matlab中带参数定积分的求解,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
1、首先,如果要对quad函数传递附加参数(Extra Parameters),应该先搞清楚quad函数的调用格式:
1 | 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得到的结果就挺平滑)。
这篇关于matlab中带参数定积分的求解的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!