操作方法
以f(x)=e^2x+sin(x+π/3) ,积分下限:a=0,积分上限:b=π/4 为例。 首先,建立被积函数M文件。 点击New,选择Function.
在编辑器中输入指令: function f=f(x) f=exp(2*x).*sin(x+pi/3); (注意,用 .* )
保存函数文件到自己的工作路径,比如G:\matlab\work。 输入指令: cd G:\matlab\work 是该路径成为当前路径。
matlab中计算积分的两种指令: 1. F=quad('fname',a,b,tol,trace) Simpson数值积分法 2. F=quad8('fname',a,b,tol,trace) Newton-Cotes数值积分法 其中: fname是被积函数表达式或函数名,a,b分别是上下限,tol可以控制积分精度,省略则取0.001;trace=1则用图形表示积分过程,trace=0,没有图形。 两者 比较,quad8精度更高。
调用积分函数squad进行计算。输入指令: F=quad('f(x)',0,pi/4) 如图,回车后可得到计算结果1.8612。 其他函数,只要修改函数文件中的表达式即可。