【matlab】数值积分公式的程序实现

这篇具有很好参考价值的文章主要介绍了【matlab】数值积分公式的程序实现。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

)专题实验(Newton-Cotes积分公式)

1、编写[a,b]上梯形积分公式、Simpson积分公式。

2、利用自己编写的程序计算定积分,计算一下数值解和精确解之间差的绝对值。

梯形积分:

function T=TX_int(f,a,b)

T=(b-a)/2*(f(a)+f(b));

>> TX_int(@(x)cos(x),0,pi/4)

ans =

0.6704

function T=TX_int(f,a,b)

T=(b-a)/2*(f(a)+f(b));

a=abs(sqrt(2)/2-T);

disp(a)

>> TX_int(@(x)cos(x),0,pi/4)

   0.0367

ans =

   0.6704

Simpson积分;

function T=TX_int(f,a,b)

T=(b-a)/6*(f(a)+4*f((a+b)/2)+f(b));

>> TX_int(@(x)cos(x),0,pi/4)

ans =

   0.7072

function T=TX_int(f,a,b)

T=(b-a)/6*(f(a)+4*f((a+b)/2)+f(b));

d=abs(sqrt(2)/2-T);

disp(d);

>> TX_int(@(x)cos(x),0,pi/4)

   9.5166e-05

ans =

0.7072文章来源地址https://www.toymoban.com/news/detail-445919.html

到了这里,关于【matlab】数值积分公式的程序实现的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 基于MATLAB的双重积分的数值求解

    目录 一. Gauss求积公式 1.1 数学理论 1.2 例题与MATLAB代码 例题1 二. 双重积分问题的数值解 2.1 数学理论 2.2 MATLAB代码与例题 例题2  例题3  例题4 例题5 例题6 三. 小结 我们知道数值的方法求解积分往往精度会差一点。这个时候 Gauss求积公式 可以让数值求解积分有较高的代数精度

    2023年04月26日
    浏览(27)
  • MATLAB 之 数值积分和离散傅里叶变换

    数值积分时研究定积分的数值求解方法,即借助于计算机,用数值逼近的方法近似计算定积分。 我们假设 I 1 = ∫ a b f ( x ) d x I_{1}=int_{a}^{b}f(x)mathrm{d}x I 1 ​ = ∫ a b ​ f ( x ) d x I 2 = ∫ a b p ( x ) d x I_{2}=int_{a}^{b}p(x)mathrm{d}x I 2 ​ = ∫ a b ​ p ( x ) d x 从高等数学中知道,当 ∣

    2024年02月08日
    浏览(31)
  • 10、MATLAB程序设计与应用刘卫国(第三版)课后实验十:方程数值求解

    目录  一、  二、  三、  四、  五、 分别用 3种不同的数值方法解线性方程组。   --------------------------------------- 示例代码 - -------------------------------------------- -------------------------------------- - 运行结果 --------------------------------------------- 求代数方程的数值解。 (1)3x +sin x-e

    2024年02月16日
    浏览(35)
  • Matlab绘制双纽线、莫比乌斯环,双纽线公式、重积分

    今天在做一道二重积分时遇到了这样一个问题,看题: 设 D = { ( x , y ) ∣ ( x 2 + y 2 ) 2 ≤ 4 ( x 2 − y 2 ) } D = {(x,y) | (x^2+y^2)^2 leq 4(x^2-y^2)} D = {( x , y ) ∣ ( x 2 + y 2 ) 2 ≤ 4 ( x 2 − y 2 )} ,则区域D的面积为? 当时我还不知道这是双纽线的公式,就直接这样做了: 算出来面积怎么等

    2024年02月06日
    浏览(24)
  • 实验6 Matlab数值计算

    实验目的: 掌握数据统计与分析的方法; 掌握数据插值和曲线拟合的方法及其应用; 掌握多项式的常用运算。 实验内容: 利用randn函数生成符合正态分布的10×5随机矩阵A,进行如下操作: 求A的最大元素和最小元素; 求A的每行元素的和以及全部元素的和; 分别对A的每列元

    2024年02月06日
    浏览(51)
  • 【数值分析实验】(八)常微分方程的数值解法(含matlab代码)

            科学技术中很多问题都可用常微分方程的定解问题来描述,主要有初值问题和边值问题两大类。常微分方程式描述连续变化的数学语言,微分方程的求解时确定满足给定方程的可微函数,要找出这类问题的解析解往往非常困难,甚至是不可能的。研究一阶常微分方

    2024年02月03日
    浏览(34)
  • 用C语言实现定积分求解的三种方法,梯形公式,辛普森公式,自适应辛普森公式

    1.梯形公式: 梯形公式(trapezoidal rule)是一种求定积分的方法。它假定函数在区间上是一条直线,因此可以通过计算梯形的面积来估计函数的定积分 可以用指针来初步优化这个代码: 2.辛普森公式: 辛普森公式(Simpson\\\'s rule)是一种求定积分的方法。它是由英国数学家 Tho

    2024年02月06日
    浏览(34)
  • 【数值分析实验】(一)插值法(含matlab代码)

            实际问题中许多变量的关系可以用数学函数概念进行刻画,但是在大多数情况下,这些函数的表达式是未知的,或者已知但十分复杂,需要我们将这个函数的未知解析式近似地构造出来,或者用一个简单的函数表达式来代替复杂的函数表达式。基于上述过程,我们

    2024年02月12日
    浏览(36)
  • 【数值分析实验】(七)特征值与特征向量(含matlab代码)

            利用已有的非线性方程的数值解法能够近似计算部分特征值,但要求出特征方程的所有根难度极大。幂法是一种计算矩阵主特征值及对应特征向量的迭代方法,特别适用于大型稀疏矩阵。反幂法是计算海森伯格阵或三对角阵的对应一个给定近似特征值的特征向量的

    2024年02月04日
    浏览(33)
  • 【数值分析实验】(五)线性方程组的迭代解法(含matlab代码)

            迭代法就是用某种极限过程去逐步逼近线性方程精确解的方法。迭代法具有需要计算机的存储单元较少、程序设计简单、原始系数矩阵在计算过程中始终不变等优点,但存在收敛性及收敛速度问题。 3.1.1 算法过程 3.1.2 代码 3.1.3 计算结果 3.2.1 算法过程 3.2.2 代码

    2024年02月03日
    浏览(33)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包