06新 数值积分+常微分方程求解.pdfVIP

目录CONTENTS 数值积分 常微分方程求解 数值积分 b y a f (x )dx 的几何意义 y=f(x) 三种方法： x  矩形法 a b 用小矩形面积代替小曲边梯形的面积，然后求和  复合梯形公式 用小梯形面积代替小曲边梯形的面积，然后求和  复合辛普生公式 用抛物线代替小曲边梯形的曲边计算小面积，再求和 数值积分 矩形法 1 4 例7-28 计算定积分0 1  x2 dx , 与精确值 比较。 解：MATLAB命令为: h=0.01;x=0:h:1; y=4./(1+x.^2); format long z1 = 3.151575986923129 z1=sum(y(1:length(x)-1))*h %左矩形公式 z2 = 3.131575986923129 z2=sum(y(2:length(x)))*h %右矩形公式 数值积分 复合梯形公式 1 4 例7-28 计算定积分0 1  x2 dx , 与精确值 比较。 解：MATLAB命令为: h=0.01;x=0:h:1; y=4./(1+x.^2); format long z1 = 3.151575986923129 z1=sum(y(1:length(x)-1))*h %左矩形公式 z2 = 3.131575986923129 z2=sum(y(2:length(x)))*h %右矩形公式 z3 = 3.141575986923129 z3=trapz(x,y) %复合梯形公式 数值积分 复合辛普生公式 命令：quad(‘fun’,a,b,tol)  fun是被积函数表达式字符串或者是M 函数文件；  a,b是积分的下限与上限；  tol代表精度，可以缺省（tol=0.001 ）；  trace=1时用展示积分过程，省略时不展示。 数值积分 复合辛普生公式 1 4 例7-28 计算定积分0 1  x2 dx , 与精确值 比较。 解： 方法1: z4=quad(‘4./(1+x.^2)’,0,1) %复合辛普生公式 z1 = 3.151575986923129 方法2： z2 = 3.131575986923129 function y=fun(x) z3 = 3.141575986923129 函数文件 y=4./(1+x.^2); z4 = 3.141592682924567 z5 = 3.