精华内容
下载资源
问答
  • 分别用复化梯形公式、复化Simpson公式计算定积分分别验证结果(精确值I=4.006994)。复化梯形公式求定积分:function I=tquad(x,y)%复化梯形求积公式,其中,%x为向量,被积函数自变量的等距结点; %y为向量,被积函数...

    分别用复化梯形公式、复化Simpson公式计算定积分分别验证结果(精确值I=4.006994)。

    复化梯形公式求定积分:

    function I=tquad(x,y)

    %复化梯形求积公式,其中,

    %x为向量,被积函数自变量的等距结点; %y为向量,被积函数在结点处的函数值; n=length(x);

    m=length(y);

    %积分自变量的结点数应与它的函数值的个数相同 h=(x(n)-x(1))/(n-1);

    a=[1 2*ones(1,n-2) 1];

    I=h/2*sum(a.*y);

    复化Simpson公式求定积分:

    function I=squad(x,y)

    %复化Simpson求积公式,其中,

    %x为向量,被积函数自变量的等距结点; %y为向量,被积函数在结点处的函数值; n=length(x);

    m=length(y);

    %积分自变量的结点数应与它的函数值的个数相同 if rem(n-1,2)~=0

    I=tquad(x,y);

    return;

    end

    N=(n-1)/2;

    h=(x(n)-x(1))/N;

    a=zeros(1,n);

    for k=1:N

    a(2*k-1)=a(2*k-1)+1;

    a(2*k)=a(2*k)+4;

    a(2*k+1)=a(2*k+1)+1;

    end

    I=h/6*sum(a.*y);

    20 exdx,取n=2,4,8,16

    展开全文
  • 由误差余项公式可知区间过大,误差亦大;为避免可选取适当多的节点,即选取相对高阶的Newton-cotes公式,但由稳定性分析又知:当阶数过大时,...例如:复化Simpson公式的推导:复化Simpson公式误差分析:其中有个加...

    由误差余项公式可知区间过大,误差亦大;为避免可选取适当多的节点,即选取相对高阶的Newton-cotes公式,但由稳定性分析又知:当阶数过大时,会出现不稳定的Runge现象。

    复化求积法:将积分区间进行适当分段,在各分段子区间上采用低阶的Newton-Cotes求积公式,对各个小区间上的积分值进行一个近似,最后再累加起来。

    例如:复化Simpson公式的推导:

    ec6307b2f679dbedd971a8ce4ff0297a.png

    复化Simpson公式误差分析:

    6e31f35442c42f140533e1262cd73f84.png

    87ae764d7e9b0776f8c6987c86bd25db.png

    260b01aa4b765ebefa5c70fd031e2465.png

    其中有个加一项减一项

    区间逐次二分法

    由复化求积公式的截断误差可知,加密节点可以提高求积公式的精度,但困难在于:使用公式之前需给出合适的步长,h过大,满足不了精度;h过小,计算量过大,因而实用的方法是采用区间逐次二分法,反复利用求积公式计算,直至二分前后两次积分值的差满足精度为止。

    比如:对区间进行n等分,每个区间上先采用梯形公式,即复化梯形公式,若不能满足精度,则将每个小区间二等分,再分别采用梯形公式,不过端点处的值不用再算了,新算的只有新小区间上的二等分点处值,这样便可使计算量节约一半。

    龙贝格算法(自动调整等分数)

    ed0b29eddcb7bcb7fc9de9cef07e9568.png

    efe0ec50e0fbfaac5448264ab77a9926.png

    8842478db0792d6bfe16c8ed3cbc7cbf.png
    展开全文
  • 参考北京师范大学的《计算物理基础》第三章:数值微分积分计算物理基础_中国大学MOOC(慕课)​www.icourse163.org1.数值微分的计算数值微分计算公式:数值差分指令(diff):y=diff(x) y=diff(x,n) y=diff(x,n,dim)x...

    参考北京师范大学的《计算物理基础》

    第三章:数值微分积分

    计算物理基础_中国大学MOOC(慕课)www.icourse163.org
    bc0f562be0d632f8b766dab3c2dd16ca.png

    1.数值微分的计算

    数值微分计算公式:

    78a863e874f8bbe8b50ba850c8da4207.png

    数值差分指令(diff):

    y=diff(x)
    y=diff(x,n)
    y=diff(x,n,dim)
    x:代表矢量,矩阵,列阵
    n:计算n阶差分
    dim:矩阵或列阵的维度

    实例:矢量、矩阵(默认列向量)求微分

    d9be9e6a6ef970905ff82fdc23405491.png
    • 梯度(gradient):相邻两个差商值的平均
    内部元素:用一阶导数的中心差分公式计算
    端点元素:用前差公式或后差公式计算

    5bd62487e92c76db159093ae093e7715.png
    • 离散拉普拉斯算符(del2)

    拉普拉斯算子是n维欧几里得空间中的二阶微分算子,为梯度的散度。表示为

    4c9d57eb8d5f03f1676ac5e650edc5d1.png

    其中i表示x方向,j表示y方向。

    b14549acdbaf8218ac4b0a03e673fe65.png

    5de5d7bd5560d6faec8a2c1a540debca.png

    2.数值积分的算法

    定积分的梯形公式

    d053b7a9c26096ff4f582f2e627fa1f8.png

    d197634bae208d3d4de8f9739c7435a4.png
    • 梯形积分指令(trapz)
    trapz(y)            %y是函数值,可以是矢量、矩阵、列阵 
    trapz(x,y)          %x是变量值
    trapz(...,dim)      %dim是标明积分的维度。如果不写,则默认对列向量进行

    (1)对向量积分、对矩阵积分

    cf4c38fe03a12bd2fd282390fe7b5a56.png
    >> a=[1,5,7,2,3];
    >> trapz(a)
    
    ans =
    
        16

    计算过程:

    >> b=[0.95013 0.48598;0.23114 0.8913;0.60684 0.7621]
    
    b =
    
        0.9501    0.4860
        0.2311    0.8913
        0.6068    0.7621
    
    >> trapz(b)  %默认对列向量积分
    
    ans =
    
        1.0096    1.5153
    
    >> trapz(b,2) %指定对行向量积分
    
    ans =
    
        0.7181
        0.5612
        0.6845

    计算过程:

    (2)对于用列阵表示的数据也可作梯形积分,这种积分相当于多元函数的积分

    13e8fe416899bb3eae91051da2bd2839.png

    (3)累计梯形积分cumtrapz,相当于不定积分的数值计算

    结果为矢量,第n个元素是原矢量前n个元素的梯形积分
    矩阵是对列矢量进行,也可以指定行向量积分

    8cd1f884ff308c1df75fd5c6067b0a49.png
     %针对上面的b
    >> cumtrapz(b)
    
    ans =
    
             0         0
        0.5906    0.6886
        1.0096    1.5153 %其中最后一行就是trapz(b)的结果

    3.用函数做数值积分

    3.1 应用匿名函数

    对于matlab内部没有的函数,我们需要应用matlab的基本函数对我们的函数进行定义。

    主要内容:用函数句柄@来定义匿名函数。

    1.匿名函数的建立与使用

    用函数句柄符号@建立匿名符号
    函数句柄=@(变量)函数表达式
    要用矢量化运算表达式
    适用于函数表达式较简单的情况。

    6a4f698611c5f8e85216bdece73e7600.png

    指令是用函数句柄调用匿名函数

    2.多重匿名函数

    重复建立匿名函数就能构造多重匿名函数
    多重匿名函数必须按定义顺序分别对变量赋值

    a7b56951b538fb373db484a2c05aea39.png

    3.一元函数积分指令integral(用匿名函数定义的方式定义)

    算法:采用自适用积分法和默认的误差精度
    语法:
    integral(fun,xmin,xmax) %(积分函数,积分上限,积分下限)
    integral(fun,xmin,xmax,NameValue) %(...积分属性,积分属性值)

    cb3d54da4169094298e784c2341edcfd.png

    ee660ff5133f9c9f50ae18fd74706307.png

    RelTol:相对误差 AbsTol:绝对误差

    af2d1893769ccecc1b00c643bd7ca6a4.png

    4.二元函数积分

    integral2(fun,xmin,xmax,ymin,ymax,Name,Value) %(积分函数,积分上限,积分下限)

    a251750af582609e7f6838fa096aafcd.png

    ea8deb6e81da7abd14f8a16ea484df4d.png

    5.三元函数积分

    integral3(fun,xmin,xmax,ymin,ymax,zmin,zmax,Name,Value) %(积分函数,积分上限,积分下限)

    6f45971b63e41eb0c7261cfa6d160937.png

    a838522faf434b61404174cbb092e753.png

    3.2 应用函数文件

    1.函数文件的建立步骤

    在主页界面点击图标(新建->函数),打开程序编辑器,输入如下文件,完成后保存为函数文件test2.m

    037661ab3739c70a0fd230429f3a969f.png

    2.函数文件在指令窗口的调用

    01824fb4761fbc5028593a1d663dd138.png

    注:函数文件中所有的变量并不保存在内存空间中,只有函数值保存在内存空间中。表明:脚本文件和函数文件在变量存储上是有差别的。因此,在函数文件中一般都进行输出变量的设置,且在不同函数文件间,如果需要变量传递,一般需要设置全局变量。

    e4fb00a5d8fe55c8337aca3981cfe1f5.png

    3.函数与函数文件

    8ee7a5c0f65b2d6c5801e31049ec9ed9.png

    在matlab中有两类程序(M)文件,它们的扩展名都是m,一种是脚本文件,一种是函数文件。脚本文件和函数文件特点如下:

    a88ceae515773bceaf20e595aadd7cf7.png

    函数文件可完成脚本文件的功能,但是,脚本文件一般不具有函数文件的功能。

    脚本文件和函数文件比较如下:

    efffca4e7f208eb045e518ee0820c965.png

    4.用函数文件做积分

    函数文件建立的函数与匿名函数功能相同,如进行积分运算

    1ed4de8016b8811c0ff7a69b36accc1f.png

    1ad606b8b5d3ecabdb6400e1d1380067.png

    4.本地函数与嵌套函数

    1.主函数和子函数

    20784a640a7e190fed1a0de9733d2ef1.png
    从函数编写格式以及形态上来说,子函数和主函数并无区别,区别仅在于在函数文件中的位置以及调用关系。函数文件外部只能调用主函教文件,子函教文件只能在主函教文件内部调用。
    子函数是在同一函数文件中,主函数之后的由“function”引导的函数,一个函数文件可以有多个子函数。这些子函数在函数文件中的排列顺序可以随意,前提是都要位于主函数后。
    子函数文件分为本地函数和嵌套函数,二者的区别在于共享变量的方式不同。

    2.嵌套函数

    嵌套函数,即nested function,顾名思义,是嵌套在函数体内部的函数。以function声明,结束的时候加上end.

    d14e1d5db0225496697a0525fbfbfe0d.png

    d03a0b03302579ef07e1d4f9c60a11a2.png
    • 主函数文件与本地函数例子:

    7025aa3fbf5cd5983d5f0610f6b9d32d.png
    • 主函数文件与嵌套函数例子:可共享变量

    5fec4a5d934b6aaa897e19ff11f9bf8f.png

    如果主函数没有调用嵌套函数中某个变量,则这个变量仍然是嵌套函数的局部变量.例如,在下面的主函数中有两个嵌套函数,它们都有自己的变量x,两者并不相互影响。

    7023cfac7f53dd5332c353e3b6e997be.png
    • 嵌套函数小结
    优点:子函数与主函数可以共享变量。
    注意事项:
    1.在流程控制语句中如if,switch,for,while等不能定义嵌套函数。
    2.调用嵌套函数必须用函数名或者@建立的函数句柄,而不能用其它方法如feval,str2func.

    e422da8462e25a2added17f09a06af80.png

    9bbd7cbccf5516ec3a1cecd74c840c2a.png

    2020.12.4

    展开全文
  • 曲边梯形的面积(1)曲边梯形:由直线x=a、x=b(a≠b)、y=0和曲线 y=f(x)所围成的图形称为曲边梯形(如图①).(2)求曲边梯形面积的方法与步骤:①分割:把区间[a,b]分成许多小区间,进而把曲边梯形拆分为一些小曲边梯形...

    考纲原文

    (1)了解定积分的实际背景,了解定积分的基本思想,了解定积分的概念.

    (2)了解微积分基本定理的含义.

    知识点详解

    一、定积分

    1.曲边梯形的面积

    (1)曲边梯形:由直线x=ax=b(ab)、y=0和曲线 y=f(x)所围成的图形称为曲边梯形(如图①).

    (2)求曲边梯形面积的方法与步骤:

    ①分割:把区间[ab]分成许多小区间,进而把曲边梯形拆分为一些小曲边梯形(如图②);

    ②近似代替:对每个小曲边梯形“以值代曲”,即用矩形的面积近似代替小曲边梯形的面积,得到每个小曲边梯形面积的近似值(如图②);

    ③求和:把以近似代替得到的每个小曲边梯形面积的近似值求和;

    ④取极限:当小曲边梯形的个数趋向无穷时,各小曲边梯形的面积之和趋向一个定值,即为曲边梯形的面积.

    64d87b8741c8674190715938b1876017.png

    2.求变速直线运动的路程

    如果物体做变速直线运动,速度函数为v=v(t),那么也可以采用分割、近似代替、求和、取极限的方法,求出它在atb内所作的位移s.

    3.定积分的定义和相关概念

    aa6616d6a56a9ef72278994a5407f56f.png

    4.定积分的性质

    62229bd39b0d4f5fe32d5c9fd1368c36.png

    【注】定积分的性质(3)称为定积分对积分区间的可加性,其几何意义是曲边梯形ABCD的面积等于曲边梯形AEFD与曲边梯形EBCF的面积的和.

    6ff5cb4de13283162ccc292fedfc35b5.png

    5.定积分的几何意义

    83a2c4b6a061933b7d849c908b52a1a4.png
    bde5bdbdaac6fc4a657c5c340108b114.png

    6.定积分与曲边梯形的面积的关系(常用结论)

    定积分的概念是从曲边梯形面积引入的,但是定积分并不一定就是曲边梯形的面积.这要结合具体图形来确定:

    e33fb4332437084036e00bfa6aed71a3.png

    设阴影部分面积为S,则

    eda601994688f9c8a4c65115523777ea.png

    7.定积分的物理意义

    c820b35c62fe1f5d9012305ca3b8e774.png

    二、微积分基本定理

    1ec06591bf06d68205ecd163a2eb3937.png

    【注】常见的原函数与被积函数的关系

    e4492438afb598fdd18f9cba5a7c2e10.png

    考向分析

    考向一 定积分的计算

    1.求定积分的三种方法

    (1)利用定义求定积分(定义法),可操作性不强;

    (2)利用微积分基本定理求定积分;

    007f836c7528f6d1693908abd6a45e36.png

    2.用牛顿—莱布尼茨公式求定积分的步骤

    (1)把被积函数变形为幂函数、正弦函数、余弦函数、指数函数与常数的积的和或差;

    (2)把定积分用定积分性质变形为求被积函数为上述函数的定积分;

    (3)分别用求导公式找到一个相应的原函数;

    (4)利用牛顿—莱布尼茨公式求出各个定积分的值;

    (5)计算原始定积分的值.

    3.分段函数的定积分

    分段函数求定积分,可先把每一段函数的定积分求出后再相加.

    4.奇偶函数的定积分

    54ac22f1ebbdec375e767b5f7a1f5b7d.png

    考向二 利用定积分求平面图形的面积

    利用定积分求平面图形面积问题的常见类型及解题策略

    (1)利用定积分求平面图形面积的步骤

    ①根据题意画出图形;

    ②借助图形确定出被积函数,求出交点坐标,确定积分的上、下限;

    ③把曲边梯形的面积表示成若干个定积分的和;

    ④计算定积分,写出答案.

    (2)知图形的面积求参数

    求解此类题的突破口:画图,一般是先画出它的草图;然后确定积分的上、下限,确定被积函数,由定积分求出其面积,再由已知条件可找到关于参数的方程,从而可求出参数的值.

    (3)与概率相交汇问题

    解决此类问题应先利用定积分求出相应平面图形的面积,再用相应概率公式进行计算.

    考向三 定积分的物理意义

    利用定积分解决变速直线运动与变力做功问题

    利用定积分解决变速直线运动问题和变力做功问题时,关键是求出物体做变速直线运动的速度函数和变力与位移之间的函数关系,确定好积分区间,得到积分表达式,再利用微积分基本定理计算即得所求.

    【名师点睛】

    1、定积分的计算一般有三个方法:

    ①利用微积分基本定理求原函数;

    ②利用定积分的几何意义,即利用面积求定积分;

    ③利用奇偶性、对称性求定积分,如奇函数在对称区间的定积分值为0.

    2、由函数图象或曲线围成的曲边图形面积的计算及应用,一般转化为定积分的计算及应用, 但一定要找准积分上限、下限及被积函数,且当图形的边界不同时,要讨论解决.具体步骤如下:

    (1)画出图形,确定图形范围;

    (2)解方程组求出图形交点坐标,确定积分上、下限;

    (3)确定被积函数,注意分清函数图形的上、下位置;

    (4)计算定积分,求出平面图形的面积.

    3.由函数求其定积分,能用公式的利用公式计算,有些特殊函数可根据其几何意义,求出其围成的几何图形的面积,即其定积分.

    展开全文
  • 上章围绕曲边梯形的面积问题,引出了定积分的概念,并给出定积分的计算方法,刻画了可积条件,从理论上完善了定积分。而马克思再次告诉我们,理论要服务于实践,所以本章讨论定积分的应用!定积分的所有应用问题,都...
  • 一、为什么要数值积分传统的数值积分利用了Newton-Leibniz公式(牛顿-莱布尼茨公式),连续函数 在 上的定积分 其中 是 的原函数。对于大多数问题,牛顿-莱布尼茨公式不能使用,可能有这样几个原因。 找不到精确的原...
  • 哎呀,这位姑娘似乎认为解决微积分问题,只要套用背诵的公式就足够了。这就是那种在学校的考试中掌握了应试要领的典型人物。不过,对于如何看待微积分,还存在像上面这位博士一样的一类人,虽然会计算微积分更好,但...
  • 根据题意和复化梯形公式、复化辛普森公式的原理编辑程序求解代码如下: Matlab 代码 clc s=quad('sin(x)./x',0,1) p1=zeros(10,1); p2=......(xk1)] ; 根据题意和复化梯形公式、复化辛普森公式的原理编辑程序求解...
  • 复化梯形公式-辛普森公式的matlab程序》由会员分享,可在线阅读,更多相关《复化梯形公式-辛普森公式的matlab程序(2页珍藏版)》请在人人文库网上搜索。1、复化梯形公式与辛普森公式的matlab程序【程序代码】cclc;...
  • 《数值分析复化梯形公式,复化Simpson公式MATLAB程序》由会员分享,可在线阅读,更多相关《数值分析复化梯形公式,复化Simpson公式MATLAB程序(1页珍藏版)》请在人人文库网上搜索。1、分别用复化梯形公式、复化Simpson...
  • 包含代码和文档 采用复化梯形公式和复化辛普森公式求积分,并与精确值进行比较得下表。 采用复化梯形公式和复化辛普森公式求积分,并与精确值进行比较得下表。
  • matlab中利用复化梯形公式和复化simpson公式实现积分运算,对于数值计算类课程很有帮助。
  • 复化梯形公式与辛普森公式的matlab程序【程序代码】cclc;disp('1.复化梯形公式求解');disp('2.simpson公式求解');disp('请进行选择:');c=input(' ');if c==1clc;disp('复化梯形公式');disp('请输入积分下限 ');a=...
  • 2 分别用复化梯形公式复化 Simpson 公式计算定积分 1+ ex dx 取 n=2,4,8,16 0 分别验证结果精确值I=4.006994 复化梯形公式求定积分: function I=tquad(x,y) %复化梯形求积公式其中 %x 为向量被积函数自变量的等距...
  • PAGE PAGE 1 分别用复化梯形公式复化Simpson公式计算定积分取n=2,4,8,16分别验证结果精确值I=4.006994 复化梯形公式求定积分: function I=tquad(x,y) %复化梯形求积公式其中 %x为向量被积函数自变量的等距结点 %y为...
  • 分别用复化梯形公式、复化Simpson 公式计算定积分dx e x ?+201,取n=2,4,8,16分别验证结果(精确值I=)。复化梯形公式求定积分:function I=tquad(x,y)%复化梯形求积公式,其中,%x 为向量,被积函数自变量的等距结点;...
  • 一、复化梯形公式求积分 题目: 利用复化梯形公式计算 I(f)=∫15sinxxdxI(f) = \int_1^5{\frac{sinx}{x}dx}I(f)=∫15​xsinx​dx 的近似值,精确至7位有效数字。 Matlab程序: clc; clear all syms x %% 已知条件 ...
  • 复合梯形公式(compound trapezoidal formula)是1993年公布的数学名词。将积分区间分为若干份, 在每一个“小区间”上用低阶梯形求积公式可得复合梯形公式的收敛阶为2阶。matlab程序function I = ftrapz(fun,a,b,n)%...
  • 资源描述复化积分法(复化梯形求积,复化Simpson公式,变步长求积法)MATLAB编程实验报告一、 问题描述:编写函数实现复化积分法。二、 实验步骤(过程):(1) 复化求积法(1) 复化梯形求积:用复化梯形求积公式求解...
  • 复化梯形公式和复化辛普森公式 以及matlab源代码
  • 复化梯形公式求二重积分matlab源码

    千次阅读 2020-06-08 13:11:05
    复化梯形公式求二重积分matlab源码 %%%%%%%%%% 2020.6.5 %%%%%%%%% %%%%%%%%%%复化梯形公式求二重积分%%%%%%%%%%%%%%%%%%% Liu Deping %%%%%%%%%clear all;%%%被积函数及积分上下限导入;s=input('请输入函数表达式...
  • matlab代码
  • n T_n = C_t_r(f,n,a,b) %f为函数,n为精度,a,b为区间端点 h=(b-a)/n; %f=sqrt(x); for i=0:n x(i+1)=a+ih; end T_1=h/2(1+f(b)); T_2=0; for j=2:n F(j)=h*f(x(j)); T_2=T_2+F(j);......

空空如也

空空如也

1 2 3 4
收藏数 67
精华内容 26
关键字:

matlab复化梯形公式

matlab 订阅