实验4_数值积分与数值微分.pdfVIP

数值分析实验报告四 数值积分与数值微分实验（ 2 学时） 班级专业 统计学 1 班 姓名 何斌 学号 06 日期 2011-5-10 一 实验目的 1．掌握复化的梯形公式、辛扑生公式等牛顿 - 柯特斯公式计算积分。 2. 会用高斯公式计算积分。 3. 掌握数值微分的计算方法。 二 实验内容 1．分别用复化的梯形公式和辛扑生公式计算积分。 9 x dx M=4 0 2．用高斯公式（ n=3）计算积分 3 1 dx 1 x 3．给定下列表格值 X 50 55 60 65 Y 1.6990 1.7404 1.7782 1.8129 利用四点式（ n=3）求 f (50)和f (50) 的值。 三 实验步骤（算法）与结果 9 1. 用复化的梯形公式计算 x dx M=4 0 function I=tx(fun,a,b,M) %复合梯形公式计算定积分 %b,a 分别为上下限 %M为面积段数 %n为总结点数 n=M+1; 1 h=(b-a)/(n-1); %步长 x=a:h:b; f=feval(fun,x); I=h*(0.5*f(1)+sum(f(2:n-1))+0.5*f(n)); f=inline(x.^(1/2)); jf_1=tx(f,0,9,4) jf_1 = 17.3686 9 用复化的辛扑生公式计算积分 x dx M=4 0 function I=simpson(fun,a,b,M) %复合辛扑生公式计算定积分 %b,a 分别为上下限 %M为面积段数 %n为总结点数 n=2*M+1; h=(b-a)/(n-1); %步长 x=a:h:b; f=feval(fun,x); I=(h/3)*(f(1)+4*sum(f(2:2:n-1))+2*sum(f(3:2:n-2))+f(n)); f=inline(x.^(1/2)); jf_2=simpson(f,0,9,4) jf_2 = 17.9031 3 1 2．用高斯公式（ n=3）计算积分 dx 1 x function I=gsquad(fun,a,b)