MATLAB 常用函数简介
一、定积分
数值求解定积分基本方法有:欧拉法(矩形法)和梯形法; 1、矩形法数值积分公式为:1
1n i i
i s y x -==?∑
Matlab 命令:
s=sum(y1(1:n-1))*dx;
2、梯形法的公式为:1
1
1
1
2
2
2n n i i i n
i i i i i y y y y q x y x --+==++??
=?=+
? ??
?
∑
∑
Matlab 命令:
q=trapz(y1)*dx;
3、采用MA TLAB 中的quad 函数
quad 函数是自适应法(也称辛普生法),即采用二次曲线逼近被积分函数的算法。 它们的调用格式如下: 1) q=quad (fun,a,b,tol) 2) q=quadl(fun,a,b,tol)
其中,fun 是表示被积函数的字符串,a 是积为下限,b 是积分上限,tol 是规定计算的容差,其默认值为1e-6,通常不必输入。quadl 函数采用了更高阶的逼近方法。
二、常微分方程的数值求解
MA TLAB 中提供了专门求解常微分方程的功能函数—