南通大学计算机类课程实验报告,南通大学实验报告

南通大学实验报告

(15页)

本资源提供全文预览，点击全文预览即可全文预览,如果喜欢文档就下载吧，查找使用更方便哦！

11.9 积分

南通大学实验报告定积分与定积分的近似计算 学院： 理学院 班级： 数师153班 学号： 1502012072 姓名： 顾阳 第一部分实验报告书解读一、实验目的实验主要是分析用矩阵公式，梯形公式，辛普森公式求定积分的近似值，并比较它们与定积分的近似情况。可以先学习定积分的数值计算方法，理解定积分的定义，掌握牛顿-莱布尼茨公式。 二、实验材料1.1定积分的数值计算 计算定积分的近似值，可将积分区间等分而得矩形公式程序为或也可用梯形公式近似计算如果要准确些，可用辛普森公式 对于，矩形公式、梯形公式、辛普森公式的Mathematica程序为 a=0;b=1;k=10; f[x_]:=Sin[x]; d=N[Integrate[f[x],{x,a,b}],k];(计算精确值) s1[m_]:=N[Sum[f[a+i*(b-a)/m]*(b-a)/m,{i,0,m-1}],k];(取小区间左端点的矩形公式) s2[m_]:=N[Sum[f[a+(i+1/2)*(b-a)/m]*(b-a)/m,{i,0,m-1}],k](取小区间中点的矩形公式) s3[m_]:=N[Sum[f[a+i*(b-a)/m]*(b-a)/m,{i,1,m}],k]; (取小区间右端点的矩形公式)s4[m_]:=N[Sum[(f[a+i*(b-a)/m]+f[a+(i+1)*(b-a)/m])/2*(b- a)/m,{i,0,m-1}],k]; (梯形公式) s5[m_]:=N[(b-a)/m/6*((f[a]+f[b])+2*Sum[f[a+i*(b-a)/m],{i,1,m-1}] +4*Sum[f[a+(i-1/2)*(b-a)/m],{i,1,m}]),k];(辛普森公式) t=Table[{s1[m],r1[m],s2[m],r2[m],s3[m],r3[m],s4[m],r4[m],s5[m],r5[m]}, {m,100,1000,100}]1.2可积的条件设f(x)=sinx,取a=0,b=1对于，矩形公式、梯形公式、辛普森公式的Mathematica程序为 a=0;b=1;k=10; f[x_]:=Sin[x]; d=N[Integrate[f[x],{x,a,b}],k];(计算精确值) s1[m_]:=N[Sum[f[a+i*(b-a)/m]*(b-a)/m,{i,0,m-1}],k];(取小区间左端点的矩形公式) s2[m_]:=N[Sum[f[a+(i+1/2)*(b-a)/m]*(b-a)/m,{i,0,m-1}],k](取小区间中点的矩形公式) s3[m_]:=N[Sum[f[a+i*(b-a)/m]*(b-a)/m,{i,1,m}],k]; (取小区间右端点的矩形公式)s4[m_]:=N[Sum[(f[a+i*(b-a)/m]+f[a+(i+1)*(b-a)/m])/2*(b- a)/m,{i,0,m-1}],k]; (梯形公式) s5[m_]:=N[(b-a)/m/6*((f[a]+f[b])+2*Sum[f[a+i*(b-a)/m],{i,1,m-1}] +4*Sum[f[a+(i-1/2)*(b-a)/m],{i,1,m}]),k];(辛普森公式) r1[m_]:=d-s1[m];r2[m_]:=d-s2[m];r3[m_]:=d-s3[m];r4[m_]:= t=Table[{s1[m],r1[m],s2[m],r2[m],s3[m],r3[m],s4[m],r4[m],s5[m],r5[m]}, {m,100,1000,100}]1.3牛顿-莱布尼茨公式设函数在上连续，而且是的一个原函数，则有牛顿-莱布尼兹公式。函数在不连续、不存在原函数，但在上可积；函数在不连续，但在上可积。此外函数处处不连续、不存在原函数，在任意区间(长度大于0)上不可积。求原函数并验证牛顿-莱布尼兹公式的Mathematica程序 f[x_]:=Sin[x]; Integrate[f(x),x](求不定积分) F[x_]:=%(定义原函数) d=NIntegrate[f(x),{x,a,b}](求定积分) df=F[b]-F[a] (计算原函数的增量)三、实验所用软件及版本 Mathematica 5.0第二部分 实验计划(一)定积分的数值计算 1.程序修改a=0;b=1;k=10; f[x_]:=Sin[x]; d=N[Integrate[f[x],{x,a,b}],k]; s1[m_]:=N[Sum[f[a+i*(b-a)/m]*(b-a)/m,{i,0,m-1}],k]; s2[m_]:=N[Sum[f[a+(i+1/2)*(b-a)/m]*(b-a)/m,{i,0,m-1}],k] s3[m_]:=N[Sum[f[a+i*(b-a)/m]*(b-a)/m,{i,1,m}],k]; s4[m_]:=N[Sum[(f[a+i*(b-a)/m]+f[a+(i+1)*(b-a)/m])/2*(b- a)/m,{i,0,m-1}],k]; s5[m_]:=N[(b-a)/m/6*((f[a]+f[b])+2*Sum[f[a+i*(b-a)/m],{i,1,m-1}] +4*Sum[f[a+(i-1/2)*(b-a)/m],{i,1,m}]),k]; r1[m_]:=d-s1[m];r2[m_]:=d-s2[m];r3[m_]:=d-s3[m];r4[m_]:=d-s4[m];r5[m_]:=d-s5[m]; t=Table[{s1[m],r1[m],s2[m],r2[m],s3[m],r3[m],s4[m],r4[m],s5[m],r5[m]}, {m,100,1000,100}]利用以上程序计算，，，，，并对几个公式比较。2.实验思路对以上程序，分别将sinx的x替换成1，x，x2，ex，In(1+x)(二)可积的条件1.实验思路：(1)如果函数f(x)在区间[a,b]上连续，则f(x)在区间[a,b]上可积，反之亦然。(2)设一连续函数，判断其是否可积。2.程序修改 a=0;b=1;k=10; f[x_]:=Sin[x]; d=N[Integrate[f[x],{x,a,b}],k 关 键 词： 南通 大学 实验 报告

 天天文库所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。