《圆周率的近似计算》课件.pptVIP

数值积分简介 在高等数学中有一类积不出的积分，如 （概率积分） 　　　（椭圆积分） 利用数值积分法 输入： x=0:0.1:1; y= exp(-x.^2); trapz(x,y) 输出： ans = 0.746211 MATLAB命令 梯形方法 ——trapz(x,y) 抛物线方法——quad(f,a,b) 如求积分的近似值 y=inline(exp(-x.^2)); quad(y,0,1) ans = 0.746826 梯形方法 抛物线方法 梯形方法 公式推导： 条件： 当区间划分为n等分时 ——trapz(x,y) 复化梯形方法 当区间等距划分为n个子区间时 现对每个子区间再二等分，得到2n个子区间 举例 利用复化梯形算法求Pi的近似值. 取 运用复化梯形算法 输入初值： 输出结果：STOP No Yes clear;a=0;b=1; f=inline(4/(1+x*x)); t1=(b-a)/2*(f(a)+f(b)); er=1;n=1; while er1.0e-6 h=(b-a)/n; s=0; for i=1:n s=s+f(a+i*h-h/2); end t2=(t1+h*s)/2