实验3 插值和数值积分.pptVIP

将f(x)代入计算得 用n个节点，Gn的代数精度可达2n-1, 但是需解 复杂的非线性方程组，实用价值不大。 常 用 的 高 斯 公 式 将(a,b)分小，把小区间变换为(-1，1), 再用G2 代数精度为3 节点加密时，原计算信息无法利用 思路：将积分区间分小，在小区间上用n不太 大的 。而在节点加密一倍时能够利用原节点的函数值，可以把区间的端点作为固定节点。 改进的高斯公式 Gauss-Lobatto求积公式 其中a, b为小区间的端点， 为2n-2个参数， 代数精度可达到2n-3 注意：实际计算中一般采用自适应方法确定步长 用MATLAB 作数值积分 矩形 公式 Sum(x) 输入数组x(即fk),输出x的和(数) cumsum(x) 输入数组x,输出x的依次累加和(数组) 梯形 公式 trapz(x) 输入数组x,输出按梯形公式x的积分(单位步长) trapz(x,y) 输入同长度数组 x,y,输出按梯形公式 y对x的积分(步长不一定相等) 用MATLAB 作数值积分 辛普森公式 quad(@fun,a,b,tol,trace) [I,fn]=quad(…) 用自适应辛普森公式计算 tol为绝对误差，缺省时为10-6 Gauss-Lobatto公式 quadl(@fun,a,b,tol,trace) [I,fn]=quadl(…) 用自适应Gauss-Lobatt