第2章节矩阵数组符号运算幻灯片.pptVIP

数值积分 例8-2 求定积分。 (1) 被积函数文件fx.m。 function f=fx(x) f=x.*sin(x)./(1+cos(x).*cos(x)); (2) 调用函数quadl求定积分。 I=quadl(fx,0,pi) I = 2.4674 数值积分 例 分别用quad函数和quadl函数求定积分的近似值，并在相同的积分精度下，比较函数的调用次数。 调用函数quad求定积分： format long; fx=inline(exp(-x)); [I,n]=quad(fx,1,2.5,1e-10) I = 0.28579444254766 n = 65 数值积分 调用函数quadl求定积分： format long; fx=inline(exp(-x)); [J,m]=quadl(fx,1,2.5,1e-10) J = 0.28579444254881 m = 18 数值积分 （3）．被积函数由一个表格定义 在MATLAB中，对由表格形式定义的函数关系的求定积分问题用trapz(X,Y)函数。其中向量X,Y定义函数关系Y=f(X)。 例8-4 用trapz函数计算定积分。I=∫12.5e-x命令如下： X=1:0.01:2.5; Y=exp(-X); %生成函数关系数据向量 trapz(X,Y) ans =