如何用matlab计算定积分

定积分的方法总结

定积分的方法总结 定积分是新课标的新增内容，其中定积分的计算是重点考查的考点之一，下面例析定积分计算的几种常用方法． 一、定义法 例1、求 s i n b a x d x ? ， （b a <） 解:因为函数s i n x 在],[b a 上连续，所以函数sin x 在],[b a 上可积，采用特殊的 方法作积分和．取h = n a b -，将],[b a 等分成n 个小区间, 分点坐标依次为 ?=+<<+<+

用递推公式计算定积分(matlab版)

用递推公式计算定积分 实验目的： 1．充分理解不稳定的计算方法会造成误差的积累，在计算过程中会导致误差的迅速增加，从而使结果产生较大的误差。 2．在选择数值计算公式来进行近似计算时，应学会选用那些在计算过程中不会导致误差迅速增长的计算公式。 3．理解不稳定的计算公式造成误差积累的来源及具体过程； 4．掌握简单的matlab语言进行数值计算的方法。 实验题目： 对n=0,1,2,…,20，计算定积分： 实验原理： 由于y(n)= = – 在计算时有两种迭代方法，如下： 方法一： y(n)=– 5*y(n-1),n=1,2,3, (20) 取y(0)= = ln6-ln5 ≈ 0.182322 方法二： 利用递推公式：y(n-1)=-*y(n),n=20,19, (1) 而且，由 = * ≤≤* =

可取：y(20)≈*()≈0.008730. 实验容： 对算法一,程序代码如下： function [y,n]=funa() syms k n t; t=0.182322; n=0; y=zeros(1,20); y(1)=t; for k=2:20 y(k)=1/k-5*y(k-1); n=n+1; end y(1:6) y(7:11) 对算法二，程序代码如下： %计算定积分； %n--表示迭代次数； %y用来存储结果； function [y,n]=f(); syms k y_20;

y=zeros(21,1); n=1; y_20=(1/105+1/126)/2; y(21)=y_20; for k=21:-1:2 y(k-1)=1/(5*(k-1))-y(k)/5; n=n+1; end 实验结果： 由于计算过程中，前11个数字太小，后9个数字比较大，造成前面几个数字只显示0.0000的现象，所以先输出前6个，再输出7—11个，这样就能全部显示出来了。 算法一结果： [y,n]=funa %先显示一y(1)—y(6) ans = 0.1823 -0.4116 2.3914 -11.7069 58.7346

七大积分总结

七大积分总结 一． 定积分 1. 定积分的定义：设函数f(x)在[a,b]上有界，在区间[a,b]中任意插入n －1个分点： a=x 0

? ??==b a b a b a du u f dt t f dx x f )()()(。 （2） 定义中区间的分法与ξi 的取法是任意的。 （3） 定义中涉及的极限过程中要求λ→0，表示对区间[a,b]无限细分的过程，随λ →0必有n →∞，反之n →∞并不能保证λ→0，定积分的实质是求某种特殊合式的极限： 例：∑?=∞→=n i n n i f dx x f 1 1 0n 1 )()(lim （此特殊合式在计算中可以作为公式使用） 2. 定积分的存在定理 定理一 若函数f(x)在区间[a,b]上连续，则f(x)在[a,b]上可积。 定理二 若函数f(x)在区间[a,b]上有界，且只有有限个间断点，则f(x)在区间上可积。 3. 定积分的几何意义 对于定义在区间[a,b]上连续函数f(x)，当f(x)≥0时，定积分 ? b a dx x f )(在几何上表示由曲线y=f(x),x=a,x=b 及x 轴所围成的曲边梯形的面积；当f(x) 小于0时，围成的曲边梯形位于x 轴下方，定积分?b a dx x f )(在几何意义上表示曲边梯形面积的负值。若f(x)在区间上既取得正值又取得负值时，定积分的几何意义是：它是介于x 轴，曲线y=f(x),x=a,x=b 之间的各部分曲边梯形的代数和。 4．定积分的性质 线性性质（性质一、性质二）

定积分总结

定积分讲义总结 内容一 定积分概念 一般地，设函数()f x 在区间[,]a b 上连续，用分点0121i i n a x x x x x x b -=<<<<<<<=L L 将区间[,]a b 等分成n 个小区间，每个小区间长度为x ?（b a x n -?= ），在每个小区间[]1,i i x x -上取一点()1,2,,i i n ξ=L ，作和式：1 1 ()()n n n i i i i b a S f x f n ξξ==-=?=∑∑ 如果x ?无限接近于0（亦即n →+∞）时，上述和式n S 无限趋近于常数S ，那么称该常数S 为函数()f x 在区间[,]a b 上的定积分。记为：()b a S f x dx = ? 其中()f x 成为被积函数，x 叫做积分变量，[,]a b 为积分区间，b 积分上限，a 积分下限。 说明：（1）定积分 ()b a f x dx ? 是一个常数，即n S 无限趋近的常数S （n →+∞时）称为()b a f x dx ?，而不是n S ． （2）用定义求定积分的一般方法是：①分割：n 等分区间[],a b ；②近似代替：取点[]1,i i i x x ξ-∈；③求和： 1()n i i b a f n ξ=-∑；④取极限：()1()lim n b i a n i b a f x dx f n ξ→∞=-=∑? 例1．弹簧在拉伸的过程中，力与伸长量成正比，即力()F x kx =（k 为常数，x 是伸长量），求弹簧从平衡位置拉长b 所作的功． 分析：利用“以不变代变”的思想，采用分割、近似代替、求和、取极限的方法求解． 解： 将物体用常力F 沿力的方向移动距离x ，则所作的功为W F x =?． 1．分割 在区间[]0,b 上等间隔地插入1n -个点，将区间[]0,1等分成n 个小区间： 0,b n ??????，2,b b n n ?? ????，…，()1,n b b n -?????? 记第i 个区间为()1,(1,2,,)i b i b i n n n -???=? ? ??L ，其长度为()1i b i b b x n n n -??=-= 把在分段0, b n ? ???? ?，2,b b n n ?? ????，…，()1,n b b n -?????? 上所作的功分别记作：1W ?，2W ?，…，n W ? （2）近似代替 有条件知：()()11i i b i b b W F x k n n n --???=??=?? ? ?? (1,2,,)i n =L （3）求和 ()1 1 1n n n i i i i b b W W k n n ==-=?=??∑∑ =()()22222 110121122n n kb kb kb n n n n -?? ++++-==-?? ?? ??? L

数值积分算法与MATLAB实现陈悦5133201讲解

东北大学秦皇岛分校 数值计算课程设计报告 数值积分算法及MATLAB实现 学院数学与统计学院 专业信息与计算科学 学号5133201 姓名陈悦 指导教师姜玉山张建波 成绩 教师评语： 指导教师签字： 2015年07月14日

1 绪论 数值分析是计算数学的一个主要部分，计算数学是数学科学的一个分支，它研究用计算机求解各种数学问题的数值检索方其理论与软件的实现.而数值分析主要研究数值计算. 现科学技术的发展与进步提出了越来越多的复杂的数值计算问题，这些问题的圆满解决已远人工手算所能胜任，必须依靠电子计算机快速准确的数据处理能力.这种用计算机处理数值问题的方法，成为科学计算.今天，科学计算的应用范围非常广泛，天气预报、工程设计、流体计算、经济规划和预测以及国防尖端的一些科研项目，如核武器的研制、导弹和火箭的发射等，始终是科学计算最为活跃的领域. 1.1 数值积分介绍 数值积分是数值分析的重要环节，实际问题当中常常需要计算积分，有些数值方法，如微分方程和积分方程的求解，也都和积分计算相联系. 求某函数的定积分时，在多数情况下，被积函数的原函数很难用初等函数表达出来，因此能够借助微积分学的牛顿-莱布尼兹公式计算定积分的机会是不多的.另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数，对这类函数的定积分，也不能用不定积分方法求解.由于以上原因，数值积分的理论与方法一直是计算数学研究的基本课题.对微积分学做出杰出贡献的数学大师，如I.牛顿、L.欧拉、C.F.高斯、拉格朗日等人都在数值积分这个领域作出了各自的贡献，并奠定了这个分支的理论基础. 构造数值积分公式最通常的方法是用积分区间上的n 次插值多项式代替被积函数，由此导出的求积公式称为插值型求积公式.特别在节点分布等距的情形称为牛顿-科特斯公式，例如梯形公式(Trapezoidal Approximations)与抛物线公式(Approximations Using Parabolas)就是最基本的近似公式.但它们的精度较差.龙贝格算法是在区间逐次分半过程中，对梯形公式的近似值进行加权平均获得准确程度较高的积分近似值的一种方法，它具有公式简练、计算结果准确、使用方便、稳定性好等优点，因此在等距情形宜采用龙贝格求积公式(Rhomberg Integration).当用不等距节点进行计算时，常用高斯型求积公式计算，它在节点数目相同情况下，准确程度较高，稳定性好，而且还可以计算无穷积分.数值积分还是微分方程数值解法的重要依据.许多重要公式都可以用数值积分方程导出.现探讨数值积分算法以及运用MATLAB软件的具体实现

定积分的性质与计算方法

定积分的性质与计算方法 摘要: 定积分是微积分学中的一个重要组成部分,其计算方法和技巧非常 丰富。本文主要给出定积分的定义及讨论定积分的性质和计算方法,并通过一些很有代表性的例题说明了其计算方法在简化定积分计算中的强大功能。 关键词:定积分 性质 计算方法 定积分的定义 设函数f(x) 在区间[a,b]上连续，将区间[a,b]分成n 个子区间[x 0,x 1], (x 1,x 2], (x 2,x 3], …, (x n-1,x n ]，其中x 0=a ，x n =b 。可知各区间的长度依次是：△x 1=x 1-x 0, △x 2=x 2-x 1, …, △x n =x n -x n-1。在每个子区间(x i-1,x i ]中任取一点i ξ（1,2,...,n ），作和式1()n i i f x ι=ξ?∑。设λ=max{△x 1, △x 2, …, △x n }（即λ是 最大的区间长度），则当λ→0时，该和式无限接近于某个常数，这个常数叫做函数f(x) 在区间[a,b]的定积分，记为: ()b a f x dx ?。 其中：a 叫做积分下限，b 叫做积分上限，区间[a, b]叫做积分区间，函数f(x)叫做被积函数，x 叫做积分变量，f(x)dx 叫做被积表达式，∫ 叫做积分号。 对于定积分，有这样一个重要问题：函数()f x 在[a,b]上满足怎样的条件， ()f x 在[a,b]上一定可积？下面给出两个充分条件： 定理1： 设()f x 在区间[a,b]上连续，则()f x 在[a,b]上可积。 定理2： 设()f x 在区间[a,b]上有界，且只有有限个间断点，则 ()f x 在[a,b]上可积。 例：利用定义计算定积分1 20x dx ?. 解：因为被积函数2()f x x =在积分区间[0,1]上连续，而连续函数是可积的，所以积分与区间[0,1]的分法及点i ξ的取法无关。因此，为了 便于计算，不妨把区间[0,1]分成n 等份，分点为i i x n = ，1,2,,1i n =?-;这样，

matlab求定积分之实例说明

一、符号积分 符号积分由函数int来实现。该函数的一般调用格式为： int(s)：没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分； int(s,v)：以v为自变量，对被积函数或符号表达式s求不定积分； int(s,v,a,b)：求定积分运算。a,b分别表示定积分的下限和上限。该函数求被积函数在区间[a,b]上的定积分。a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷(inf)。当函数f关于变量x在闭区间[a,b]上可积时，函数返回一个定积分结果。当a,b中有一个是inf时，函数返回一个广义积分。当a,b中有一个符号表达式时，函数返回一个符号函数。 例： 求函数x^2+y^2+z^2的三重积分。内积分上下限都是函数，对z积分下限是sqrt(x*y)，积分上限是x^2*y；对y积分下限是sqrt(x)，积分上限是x^2；对x的积分下限1，上限是2，求解如下： >>syms x y z %定义符号变量 >>F2=int(int(int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),x,1,2) %注意定积分的书写格式 F2 = 1610027357/6563700-6072064/348075*2^(1/2)+14912/4641*2^(1/4)+64/225*2 ^(3/4) %给出有理数解 >>VF2=vpa(F2) %给出默认精度的数值解 VF2 = 224.92153573331143159790710032805 二、数值积分 1.数值积分基本原理 求解定积分的数值方法多种多样，如简单的梯形法、辛普生(Simpson)?法、牛顿－柯特斯(Newton-Cotes)法等都是经常采用的方法。它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b。这样求定积分问题就分解为求和问题。 2.数值积分的实现方法 基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为： [I,n]=quad('fname',a,b,tol,trace) 基于变步长、牛顿－柯特斯(Newton-Cotes)法，MATLAB给出了quadl函数来求定积分。该函数的调用格式为： [I,n]=quadl('fname',a,b,tol,trace) 其中fname是被积函数名。a和b分别是定积分的下限和上限。tol用来控制积分精度，缺省时取tol=0.001。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。返回参数I即定积分值，n为被积函数的调用次数。 例： 求函数'exp(-x*x)的定积分，积分下限为0，积分上限为1。 >>fun=inline('exp(-x.*x)','x'); %用内联函数定义被积函数fname

定积分计算的总结论文

定积分计算的总结论文公司内部档案编码：[OPPTR-OPPT28-OPPTL98-OPPNN08]

定积分计算的总结 闫佳丽 摘 要：本文主要考虑定积分的计算,对一些常用的方法和技巧进行了归纳和总结.在定积分的计算中,常用的计算方法有四种：（1）定义法、（2）牛顿—莱布尼茨公式、（3）定积分的分部积分法、（4）定积分的换元积分法. 关键词：定义、牛顿—莱布尼茨公式、分部积分、换元. 1前言 17世纪后期，出现了一个崭新的数学分支—数学分析.它在数学领域中占据着主导地位.这种新数学思想的特点是非常成功地运用了无限过程的运算即极限运算.而其中的微分和积分这两个过程，则构成系统微积分的核心.并奠定了全部分析学的基础.而定积分是微积分学中的一个重要组成部分. 2正文 那么，究竟什么是定积分呢我们给定积分下一个定义：设函数()f x 在[],a b 有定义，任给[],a b 一个分法T 和一组{}k ξξ=，有积分和 1 (,)()n k k k T f x σξξ==?∑，若当()0l T →时，积分和(,)T σξ存在有限极限， 设()0()0 1 lim (,)lim ()n k k l T l T k T f x I σξξ→→==?=∑，且数I 与分法T 无关，也与k ξ在[]1,k k x x -的取法无关，即{}0,0,:(),k T l T εδδξξ?>?>?

用MATLAB算多元函数积分

用MATLAB 计算多元函数的积分 三重积分的计算最终是化成累次积分来完成的，因此只要能正确的得出各累次积分的积分限，便可在MA TLAB 中通过多次使用int 命令来求得计算结果。但三重积分的积分域Ω是一个三维空间区域，当其形状较复杂时，要确定各累次积分的积分限会遇到一定困难，此时，可以借助MATLAB 的三维绘图命令，先在屏幕上绘出Ω的三维立体图，然后执行命令 rotate3d on ↙ 便可拖动鼠标使Ω的图形在屏幕上作任意的三维旋转，并且可用下述命令将Ω的图形向三个坐标平面进行投影： view(0,0),向XOZ 平面投影； view(90,0),向YOZ 平面投影； view(0,90),向XOY 平面投影. 综合运用上述方法，一般应能正确得出各累次积分的积分限。 例11.6.1计算zdv Ω ???，其中Ω是由圆锥曲面222z x y =+与平面z=1围成的闭区域 解 首先用MA TLAB 来绘制Ω的三维图形，画圆锥曲面的命令可以是： syms x y z ↙ z=sqrt(x^2+y^2); ↙ ezsurf(z,[-1.5,1.5]) ↙ 画第二个曲面之前，为保持先画的图形不会被清除，需要执行命令 hold on ↙ 然后用下述命令就可以将平面z=1与圆锥面的图形画在一个图形窗口内： [x1,y1]=meshgrid(-1.5:1/4:1.5); ↙ z1=ones(size(x1)); ↙ surf(x1,y1,z1) ↙ 于是得到Ω的三维图形如图:

由该图很容易将原三重积分化成累次积分： 111zdv dy -Ω=???? 于是可用下述命令求解此三重积分： clear all ↙ syms x y z ↙ f=z; ↙ f1=int(f,z.,sqrt(x^2+ y^2),1); ↙ f2=int(f1,x,-sqrt(1- y^2), sqrt(1- y^2)); ↙ int(f2,y,-1,1) ↙ ans= 1/4*pi 计算结果为4 π 对于第一类曲线积分和第一类曲面积分，其计算都归结为求解特定形式的定积分和二重积分，因此可完全类似的使用int 命令进行计算，并可用diff 命令求解中间所需的各偏导数。 例11.6.2用MATLAB 求解教材例11.3.1 解 求解过程如下 syms a b t ↙ x=a*cos(t); ↙ y=a*sin(t); ↙ z=b*t; ↙ f=x^2 +y^2+z^2; ↙ xt=diff(x,t); ↙ yt=diff(y,t); ↙ zt=diff(z,t); ↙ int(f*sqrt(xt^2 +yt^2+zt^2),t,0,2*pi) ↙ ans= 2/3*( a^2 +b^2)^1/2*a^2*pi+8/3*( a^2 +b^2)^1/2*b^2*pi^3 对此结果可用factor 命令进行合并化简： factor （ans ） ans= 2/3*( a^2 +b^2)^1/2*pi*(3* a^2 +4*b^2*pi^2) 例11.6.3用MATLAB 求解教材例11.4.1 解 求解过程如下 syms x y z1 z2↙ f= x^2 +y^2; ↙ z1=sqrt(x^2 +y^2); ↙ z2=1; ↙ z1x=diff(z1,x); ↙ z1y=diff(z1,y); ↙ z2x=diff(z2,x); ↙ z2y=diff(z2,y); ↙

定积分应用方法总结(经典题型归纳).docx

精品文档 定积分复习重点 定积分的考查频率不是很高，本讲复习主要掌握定积分的概念和几何意义，使 用微积分基本定理计算定积分，使用定积分求曲边图形的面积和解决一些简单的物 理问题等． 1. 定积分的运算性质 (1) b b kf (x)dx k f (x)dx(k 为常数 ). a a (2) b b f 1 ( x)dx b 2 ( x)dx. [ f 1 ( x) f 2 ( x)]dx f a a a b c b 其中 a

详解Matlab求积分的各种方法

详解Matlab求积分的各种方法 一、符号积分由函数int来实现。 该函数的一般调用格式为： int(s)： 没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分；int(s,v)： 以v为自变量，对被积函数或符号表达式s求不定积分；int(s,v,a,b)： 求定积分运算。 a,b分别表示定积分的下限和上限。 该函数求被积函数在区间[a,b]上的定积分。 a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷(inf)。 当函数f关于变量x在闭区间[a,b]上可积时，函数返回一个定积分结果。 当a,b中有一个是inf时，函数返回一个广义积分。 当a,b中有一个符号表达式时，函数返回一个符号函数。 例： 求函数x^2+y^2+z^2的三重积分。 内积分上下限都是函数，对z积分下限是sqrt(x*y)，积分上限是x^2*y；对y积分下限是sqrt(x)，积分上限是x^2；对x的积分下限1，上限是2，求解如下： >>syms x y z %定义符号变 量>>F2=int(int(int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),x,1,2) %注意定积分的书写格式F2 =57/-

/348075*2^(1/2)+14912/4641*2^(1/4)+64/225*2^(3/4) %给出有理数 解>>VF2=vpa(F2) %给出默认精度的数值解VF2 = 224.9 232805二、数值积分 1.数值积分基本原理求解定积分的数值方法多种多样，如简单的梯形法、辛普生(Simpson)?法、牛顿－柯特斯(Newton-Cotes)法等都是经常采用的方法。 它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]， i=1,2,…,n，其中x1=a，xn+1=b。 这样求定积分问题就分解为求和问题。 2.数值积分的实现方法基于变步长辛普生法，MATLAB给出了quad函数来求定积分。 该函数的调用格式为： [I,n]=quad('fname',a,b,tol,trace)基于变步长、牛顿－柯特斯(Newton-Cotes)法，MATLAB给出了quadl函数来求定积分。 该函数的调用格式为： [I,n]=quadl('fname',a,b,tol,trace)其中fname是被积函数名。 a和b分别是定积分的下限和上限。 tol用来控制积分精度，缺省时取tol= 0.0 01。 trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace= 0。

[全]高等数学之不定积分的计算方法总结[下载全]

高等数学之不定积分的计算方法总结不定积分中有关有理函数、三角函数有理式、简单无理函数的求法,是考研中重点考察的内容，也是考研中的难点。不定积分是计算定积分和求解一阶线性微分方程的基础，所以拿握不定积分的计算方法很重要。不定积分考查的函数特点是三角函数、简单无理函数、有理函数综合考查,考查方法是换元积分法、分部积分法的综合应用。不定积分的求法的理解和应用要多做习题,尤其是综合性的习题,才能真正掌握知识点,并应用于考研。 不定积分的计算方法主要有以下三种： (1)第一换元积分法,即不定积分的凑微分求积分法; (2)第二换元积分法 (3)分部积分法常见的几种典型类型的换元法:

樂，Q? o 金J犷- / .乍治阳必厶二如皿盒.「宀丄" 名% =a仏 找.』x二a沁沁r 年”十I '九久二严詈严妬5inx八ic5兄厶 整 I—炉 叶严 山二启虫? 常见的几种典型类型的换元法 题型一:利用第一换元积分法求不定积分

分析: 1-3 ? - IK ）-忑.旦r x 二）祝成）；网＞＜可久切 二2氐化如（長）寸 a 花不直押、朱 J 、 解： 2少弋協“尤十C__

-辿迪牆H JS m 弟 R Eff 洱 ->1和弟r 直 - —7朮呻' g 丄 U P A J 齐—系卩￡.§计 一 H a8~t ' J 乂 u D y " ?朮?

p o r t v 卩 J (r 4 5*〉J" 卩?对渎 t-k )+c p T + T d ? g T + c m -辿」

当积分j/O心(X)不好计算容易计算时［使用分部私jf(A-)Jg(.v)二f(x)g(x)- J g(x)df(x).常见能使用分部积分法的类型： ⑴卩"“dx J x n srn xdx J尢"cos皿等，方法是把。',sin-t, cosx 稽是降低X的次数 是化夫In 尢9 arcsine arctanx. 例11: J (1 + 6-r )arctanAz/.r ：解：arctan f xdx等，方法是把疋； Jx" arcsm11xdx

数值积分的matlab实现

实验10 数值积分 实验目的： 1．了解数值积分的基本原理； 2．熟练掌握数值积分的MATLAB 实现； 3．会用数值积分方法解决一些实际问题。 实验内容： 积分是数学中的一个基本概念，在实际问题中也有很广泛的应用。同微分一样，在《微积分》中，它也是通过极限定义的，由于实际问题中遇到的函数一般都以列表形式给出，所以常常不能用来直接进行积分。此外有些函数虽然有解析式，但其原函数不是初等函数，所以仍然得不到积分的精确值，如不定积分?1 0 d sin x x x 。这时我们一般考虑用数值方法计算其 近似值，称为数值积分。 10.1 数值微分简介 设函数()y f x =在* x 可导，则其导数为 h x f h x f x f h ) ()(lim )(**0* -+='→ （10.1） 如果函数()y f x =以列表形式给出（见表10-1），则其精确值无法求得，但可由下式求得其近似值 h x f h x f x f ) ()()(*** -+≈' （10.2） 表 10-1 一般的，步长h 越小，所得结果越精确。（10.2）式右端项的分子称为函数()y f x =在 *x 的差分，分母称为自变量在*x 的差分，所以右端项又称为差商。数值微分即用差商近似 代替微商。常用的差商公式为： 000()() ()2f x h f x h f x h +--'≈ （10.3） h y y y x f 243)(2 100-+-≈ ' （10.4）

h y y y x f n n n n 234)(12+-≈ '-- （10.5） 其误差均为2 ()O h ，称为统称三点公式。 10.2 数值微分的MATLAB 实现 MATLAB 提供了一个指令求解一阶向前差分，其使用格式为： dx=diff(x) 其中x 是n 维数组，dx 为1n -维数组[]21321,, ,n x x x x x x ---,这样基于两点的数值导 数可通过指令diff(x)/h 实现。对于三点公式，读者可参考例1的M 函数文件diff3.m 。 例1 用三点公式计算()y f x =在=x 1.0,1.2,1.4处的导数值，()f x 的值由下表给 解：建立三点公式的M 函数文件diff3.m 如下： function f=diff3(x,y) n=length(x);h=x(2)-x(1); f(1)=(-3*y(1)+4*y(2)-y(3))/(2*h); for j=2:n-1 f(j)=(y(j+1)-y(j-1))/(2*h); end f(n)=(y(n-2)-4*y(n-1)+3*y(n))/(2*h); 在MATLAB 指令窗中输入指令： x=[1.0,1.1,1.2,1.3,1.4];y=[0.2500,0.2268,0.2066,0.1890,0.1736];diff3(x,y) 运行得各点的导数值为：-0.2470，-0.2170，-0.1890，-0.1650，-0.0014。所以()y f x =在=x 1.0,1.2,1.4处的导数值分别为-0.2470，-0.1890和-0.0014。 对于高阶导数，MATLAB 提供了几个指令借助于样条函数进行求导，详细使用步骤如下： step1：对给定数据点（x,y ），利用指令pp=spline(x,y)，获得三次样条函数数据pp ，供后面ppval 等指令使用。其中，pp 是一个分段多项式所对应的行向量，它包含此多项式的阶数、段数、节点的横坐标值和各段多项式的系数。 step2：对于上面所求的数据向量pp ，利用指令[breaks,coefs,m,n]=unmkpp(pp)进行处理，生成几个有序的分段多项式pp 。 step3：对各个分段多项式pp 的系数，利用函数ppval 生成其相应导数分段多项式的系数，再利用指令mkpp 生成相应的导数分段多项式 step4：将待求点xx 代入此导数多项式，即得样条导数值。 上述过程可建立M 函数文件ppd.m 实现如下： function dy=ppd(pp) [breaks,coefs,m]=unmkpp(pp);

定积分计算公式和性质

第二节 定积分计算公式和性质 一、变上限函数 设函数在区间上连续，并且设x 为上的任一点， 于是， 在区间 上的定积分为 这里x 既是积分上限，又是积分变量，由于定积分与积分变量无关，故可将此改为 如果上限x 在区 间上任意变动，则对 于每一个取定的x 值，定积分有一个确定值与之对应，所以定积分在 上定义了一个以x 为自变量的函数，我们把 称为函数 在区间 上 变上限函数 记为 从几何上看，也很显然。因为X 是上一个动点， 从而以线段 为底的曲边梯形的面积，必然随着底数 端点的变化而变化，所以阴影部分的面积是端点x 的函数（见图5-10） 图 5-10

定积分计算公式 利用定义计算定积分的值是十分麻烦的，有时甚至无法计算。因此，必须寻求计算定积分的简便方法。 我们知道：如果物体以速度作直线运动，那么在时间区间上所经过的路程s 为 另一方面，如果物体经过的路程s 是时间t 的函数，那么物体 从t=a 到t=b 所经过的路程应该是（见图5-11） 即 由导数的物理意义可知：即 是 一个原函数，因此，为了求出定积分，应先求出被积函数 的原函数 ， 再求 在区间 上的增量 即可。 如果抛开上面物理意义，便可得出计算定积分的一般 方法： 设函数在闭区间上连续， 是 的一个原函数， 即 ，则 图 5-11

这个公式叫做牛顿-莱布尼兹公式。 为了使用方便，将公式写成 牛顿-莱布尼兹公式通常也叫做微积分基本公式。它表示一个函数定积分等于这个函数的原函数在积分上、下限处函数值之差。它揭示了定积分和不定积分的内在联系，提供了计算定积分有效而简便的方法，从而使定积分得到了广泛的应用。 例1 计算 因为是的一个原函数所以 例 2 求曲线 和直线x=0、x= 及y=0所围成图形面积A(5-12) 解 这个图形的面积为 二、定积分的性质 设 、 在相应区间上连续，利用前面学过的知识，可以 得到定积分以下几个简单性质： 图 5-12

大学微积分1方法总结

第一章 函数、极限、连续 注 “★”表示方法常用重要. 一、求函数极限的方法 ★1.极限的四则运算；★2.等价量替换；★3.变量代换；★4.洛比达法则；★5.重要极限；★6.初等函数的连续性；7.导数的定义；8. 利用带有佩亚诺余项的麦克劳林公式；9.夹逼定理；10利用带有拉格朗日余项的泰勒公式；11.拉格朗日定理；★12. 无穷小量乘以有界量仍是无穷小量等. ★二、已知函数极限且函数表达式中含有字母常数，确定字母常数数值的方法 运用无穷小量阶的比较、洛必达法则或带有佩亚诺余项的麦克劳林公式去分析问题，解决问题。 三、无穷小量阶的比较的方法 利用等价无穷小量替换或利用洛必达法则，无穷小量的等价代换或利用带有皮亚诺余项的佩亚诺余项公式展开 四、函数的连续与间断点的讨论的方法 如果是)(x f 初等函数，若)(x f 在0x x =处没有定义，但在0x 一侧或两侧有定义，则0x x =是间断点，再根据在0x x =处左右极限来确定是第几类间断点。如果)(x f 是分段函数，分界点是间断点的怀疑点和所给范围表达式没有定义的点是间断点。

五、求数列极限的方法 ★1.极限的四则运算；★2. 夹逼定理；★3. 单调有界定理； 4. )()(lim )()(lim ∞=?∞=∞ →+∞→A n f A x f n x ；5. 数列的重要极限；6.用定积分的定义求数列极限；7. 利用若∑∞ =1n n a 收敛，则0lim =∞→n n a ；8. 无穷小量乘以有界量 仍是无穷小量；9.等价量替换等. 【评注】1. 数列的项有多项相加或相乘式或∞→n 时，有无穷项相加或相乘，且不能化简，不能利用极限的四则运算， 2.如果数列的项用递推关系式给出的数列的收敛性或证明数列极限存在，并求极限.用单调有界定理 3.对数列极限的未定式不能用洛比达法则。因为数列作为函数不连续，更不可导，故对数列极限不能用洛比达法则. 4.由数列{}n a 中的通项是n 的表达式，即).(n f a n =而)(lim )(lim x f n f x n ∞ →∞→与是特殊与一般的关系，由归结原则知 ★5. 有lim 1011()()n n i i f f x dx n n →∞ ==?∑或1lim 1001()()n n i i f f x dx n n -→∞==?∑ 第二章 一元函数微分学 ★一、求一点导数或给处在一点可导推导某个结论的方法： 利用导数定义，经常用第三种形式 二、研究导函数的连续性的方法：

MATLAB计算积分

函数的积分和椭圆的周长 1．正弦函数的积分 [问题]求正弦函数从0到π的积分 y = sin x 当x = 0时，积分为0，画出积分的函数曲线。 [数学模型] 定积分的结果为 ππ00 sin d cos 2S x x x ==-=? 不定积分的结果为 sin d cos I x x x C ==-+? 其中C 是积分常量，由初始条件决定。当x = 0时，积分为I = 0，必有C = 1。结果为 I = -cos x + 1 [算法]根据积分的基本概念，将积分区域分为多份，用矩形法求曲线下的近似面积表示积分的近似值 1()n i i S f x x ==?∑ 矩形法的函数是sum(f)。 用梯形法求曲线下的近似面积表示积分的近似值 1 101[()()]2 n i i i S f x f x x -+==+?∑ 梯形法的函数是trapz(f)。 用数值积分的函数是quad 和quadl ，常用使用格式是 S = quad(f,a,b) 其中，f 表示被积函数，a 表示积分的下限，b 表示积分的下限。 用符号的函数是int ，常用使用格式是 S = int(f,a,b) [程序]zqy4_1.m 如下。 %正弦函数的积分 clear %清除变量 x=linspace(0,pi); %自变量向量 dx=x(2); %间隔 y=sin(x); %被积函数 s1=sum(y)*dx %矩形法积分 s2=trapz(y)*dx %梯形法积分 f=inline('sin(x)'); %被积的内线函数 s3=quad(f,0,pi) %数值定积分

s4=int('sin(x)',0,pi) %符号积分 sc1=cumsum(y)*dx; %矩形法累积积分(精度稍差) sc2=cumtrapz(y)*dx; %梯形法累积积分 figure %创建图形窗口 plot(x,-cos(x)+1,x,sc1,'.',x,sc2,'o') %画解析式和矩阵法以及梯形法积分曲线 s=int('sin(x)') %符号积分 sc3=subs(s,'x',x); %替换数值求符号积分的值 C=-sc3(1) %求积分常数 hold on %保持图像 plot(x,sc3+C,'c*') %画符号法积分曲线 grid on %加网格 fs=16; %字体大小 xlabel('\itx','FontSize',fs) %横坐标 ylabel('\intsin\itx\rmd\itx','FontSize',fs)%纵坐标 title('正弦函数的积分','FontSize',fs) %标题 legend('解析解','矩形法','梯形法','符号法')%图例 zqy4.1图 zqy4.2图 2.三角函数和指数的积分 [问题]求如下函数的积分 y = e ax sin bx 其中a = 0.5，b = 2。积分下限为0。画出积分的函数曲线。 [数学模型] 设 11e sin d sin de {e sin e cos d }ax ax ax ax I bx x bx bx b bx x a a == =-??? 11{e sin cos de }{e sin [e cos e sin d ]}ax ax ax ax ax b b bx bx bx bx b bx x a a a a =-=-+?? 因此不定积分为 221e (sin cos )ax I a bx b bx C a b =-++ 当x = 0时，I 应该为零，所以 22b C a b =+