MATLAB求数值积分的方法

2008年12月第20卷第6期石家庄职业技术学院学报

JournalofShijiazhuangVocationalTechnologyInstituteDec.2008Vol.20 No.6

文章编号:100924873(2008)0620058203

用MATLAB求数值积分的方法

陈佩宁a, 刘 竞b

Ξ

(石家庄职业技术学院a.信息工程系;b.机电工程系,河北石家庄 050081)

摘 要:介绍了数值积分法的几种计算公式及相应的MATLAB命令,并给出了用MATLAB编程求数值积分的实例.

关键词:MATLAB;数值积分;矩形公式;梯形公式;辛普森公式中图分类号:O172   文献标识码:A

1 引言

(4)f(x),只有一些由

在一元微积分学中,若已知函数f(x)在闭区间[a,b]上连续且其原函数为F(x),求f(x)区间上的定积分可用牛顿∫

a

b

f(x)dx)|

a

(b)-F(a).而在

的方法求).2用数值积分的方法求一个函数在区间[a,b]上的定积分,可利用定积分的定义来求解:

I=

MATLABint求

a

b

f(x)dx,

该命令格式为:

int(f,x,a,b) %求函数f在区间[a,b]上的

a

b

n

n→∞k=1

f(x)dx=lim

n

ξf(k)∑

,n

设In=

k=1

ξf(k)∑

,则I=limIn.

n→∞n

定积分.

例1 求

sinxdx.

2

此时称In为数值积分.显然,数值积分In就是I的近似值,并且当n越大,In就越接近于精确值I.由

[2]

于ξk取值不同,数值积分In的结果会有所不同.数值积分的计算公式也有多种:

(1)矩形公式

解 输入命令:>>symsx>>I=int(sin(x),x,0,pi/2),结果显示为:I=1.

用牛顿-莱布尼兹公式计算定积分的方法在理论上和解决实际问题中起到了很大的作用,但它并不能解决定积分计算的所有问题.在工程技术领域常遇到十分复杂的情况而无法用牛顿-莱布尼兹公式求解.其可能出现的情况[1]有:

(1)某些被积函数f(x),其原函数无法用初等函数表示,如exdx,

将积分区间[a,b]n等分,每个小区间宽度均为h=(b-a)/n,h称为积分步长.

记a=x0

n-1

∫x

2

Ln=hRn=h

dx等.

k=0n

f(xk),h∑f(xk),h∑

==

nn

①②

(2)函数f(x)结构复杂,求其原函数非常困难.(3)函数f(x)的结构虽然简单且其原函数存

k=1

在,但其原函数的结构相对复杂.

Ξ收稿日期:2007212224

称公式①,②分别为左、右矩形公式,两个矩形面积分别小于和大于所求曲边梯形的面积

.

作者简介:陈佩宁(19712),女,河北望都人,石家庄职业技术学院讲师.

matlab利用数值 积分的方法求,MATLAB求数值积分的方法相关推荐

  1. matlab利用经纬度计算距离_【Matlab】根据经纬度计算两点间的球面距离

    做建模或者研究空间数据,可能会遇到"根据经纬度计算两点间的球面距离"的问题,网上的资料很多,都是各种公式推导,但是一旦按公式编程计算,很可能得不到正确的距离.根本原因是在" ...

  2. 用matlab计算稳态误差,利用Matlab求稳态误差的两种方法.

    利用Matlab求稳态误差的两种方法 摘要:稳态误差是系统控制精度或抗扰动能力的一种度量,它是稳态性能的一个重要指标.本文介绍利用Matlab的控制系统工具箱和Simulink工具箱求取系统误差稳态的 ...

  3. matlab绘制星形线并求面积,利用曲线积分,求星形线x=acos3t,y=asin3t所围图形面积...

    利用曲线积分计算曲线所围成图形的面积 :星形线x=acos³t,y=asin&a 由热心网友提供的答案1: 计算星形线:x=acos³t,y=asin³t (a>0 利用曲线积分,求星形 ...

  4. matlab 仿真光学实验报告,光学实验数值仿真的三种方法及MATLAB实现

    光学实验数值仿真的三种方法及 MATLAB实现 5 结 论 (1)数值模拟结果表明三种方法都能对光学 实验现象进行正确地仿 真,因此在课 堂教学 中适 当应用这种仿真模拟 ,将光学实验 中复杂的数学 ...

  5. 正交矩阵,(标准)正交基,正交投影,正交分解定理,最佳逼近定理,格拉姆-施密特方法求正交基(手算+MATLAB),QR分解(手算+MATLAB计算、分析)

    正交(orthogonality).正交集(orthogonal set).单位正交集 定义:如果向量uuu和向量vvv是相互正交的,则有:u⋅v=0u\cdot v=0u⋅v=0. 因为对于零向量, ...

  6. 微分方程求解matlab冲激信号,用Matlab求冲激响应的几种方法

    ·实验技术与方法· 用 Matlab 求冲激响应的几种方法 贺富堂,应柏青,张 锋,孙 敏 ( 西安交通大学 电工电子教学实验中心,陕西 西安 710049) 摘 要: 用 Matlab 软件求冲激响 ...

  7. matlab 微分命令 求导,Matlab微分和积分

    第六讲 Matlab 微分和积分 理论介绍:微分.有限差分.积分.离散求和 软件求解:函数及常见注意事项 一.一元函数导数与微分 Matlab 由命令函数diff 来完成求导运算,调用格式为:diff ...

  8. 【MATLAB】几种求矩阵伪逆的方法

    几种求矩阵伪逆的方法 A*x=y当A不可逆如何解x 求伪逆有五种方法,左右同时乘以A','chol()分解,qr() 分解,svd()分解和pinv()求伪逆,下面来比较一下哪个最优秀. 如果A是个3 ...

  9. 颜色量化 matlab,利用几种颜色量化方法提取图片颜色色调

    利用几种颜色量化方法提取图片颜色色调 这里主要用3种方法:包含两种k-means方法和一种matlab自带的最小方差量化方法. 为了方便对比,文中选用梵高的星空(The Starry Night)作为 ...

  10. matlab求递归问题,matlab利用递归求解差分方程

    导航:网站首页 > matlab利用递归求解差分方程 时间:2018-5-23 matlab利用递归求解差分方程 function y = recur(a,b,n,x,x0,y0); % % y ...

最新文章

  1. 64位的机器上VS2008不能使用SharePoint的workflow template
  2. Apache服务器的WSASocket failed to open the inherited socket错误
  3. linux alias命令参数及用法详解--linux定义命令别名alias
  4. svn 本地仓库使用
  5. 实现redis 手动_Redis精华所在,一口气说完Redis的主从复制和哨兵模式
  6. 用c语言求最大公约数的流程图,如何用c语言求最大公约数和最小公倍数
  7. Maven项目配置、检出、运行
  8. 分布式面试 - 集群部署时的分布式 session 如何实现?
  9. win10安装NET Framework 3.5提示0x800f0906原因及解决方法
  10. 阿里云云计算 50 云监控的使用
  11. appium环境搭建python_appium环境搭建python
  12. LOJ 6437 [PKUSC 2018] PKUSC
  13. TCP/IP协议都有哪几层协议
  14. HCIE学习笔记(2)之ISIS Overload
  15. 【实用网址】:平常用不到,想用的时候找不到
  16. [每日一氵] python创建简单 command-line interface (CLI) API
  17. PS|002自制夸张表情包
  18. 用硅光电池测量白炽灯频闪
  19. php微信网页登录,PHP 微信网页登录
  20. C#服务号推送微信公众号模板消息

热门文章

  1. MATLAB切换中文(无语言选项)
  2. 单模光纤和多模光纤的区别,及光纤收发机(光电收发器)的介绍
  3. VS安装VAssistX插件导致WPF-XAML文件输入中文出现乱码问题解决方案
  4. 发那科机器人圆弧指令怎么用_发那科机器人走弧线的指令是什么
  5. 通过ESP32 logging库来启发我们的思维吧~~
  6. 数学建模1---基本方法和步骤
  7. Wordnet 与 Hownet 比较
  8. PDMS Pipeline Tool 教程(二):通用功能
  9. 超轻量级Web安全漏洞扫描工具Netsparker使用教程介绍
  10. Hex Editor Neo(十六进制编辑器)官方正式版V6.54.1.6478 | 二进制编辑器hex下载 | 非hexeditor汉化版