用MATLAB计算某些区域上的二重积分

摘要:本文研究某些区域上二重积分数值积分公式的构造,数学软件MATLAB 实

现所构造的数值积分的计算,通过MATLAB 用所得公式计算某些典型的二重积分。

主要工作包括:将定积分数值计算的几个公式推广到二重积分,编制MATLAB 程

序,最后通过具体的数值算例进行精度比较,从中选出精度高的二重积分的计算

公式,并用公式计算一些典型的二重积分。

关键词:二重积分;数值计算;插值多项式;求积公式

1.引言

二重积分的计算在科学计算中起着重要的作用,关于矩形区域上的二重积分

的计算一般都是化重积分为累次积分,然后借助定积分已有的数值积分计算公式

推导出,MATLAB 已经有这些计算公式的相应的函数,但是往往我们建模得到的

二重积分的积分区域都不是矩形区域,对于一般的非矩形区域的二重积分,直接

用MATLAB 是无法计算的。又当被积函数比较复杂,无法用初等函数表示或求其

原函数很困难时,就只能求积分的数值解。

若 在 上连续,二重积分 f (x , y )dxdy 存在且为一确定的常数,这

f (x , y ) D ∫∫

D

个数值与 的结构、 的几何形状有关,二重积分计算的基本途径是在一

f (x , y ) D

定条件下化为二次积分,本文研究的某些区域的二重积分,要求二重积分在该区

域上能化为二次积分。

[1]

二重积分的存在性 : 在闭区域 上连续,则 f (x ,y )dx dy 必存在。

f (x , y ) D ∫∫

D

[1]

定理 :若f (x , y ) 在闭区域D {a ≤x ≤b ;c (x ) ≤ y ≤d (x )}上连续,且c(x) 、

d (x) 在[a,b]上连续,则:

b d ( x)

( , ) ( , )

I ∫∫f x y dxdy ∫ ∫( ) f x y dxdy

a c x

D

上式右端是一个先对 后对 的二次积分:先把 看作 的函数,在区

x y f (x , y ) x

间 上对 计算定积分(这时 看作常数),把得到的结果(是 的函数)

[c(x),d (x)] x y y

再在 [a,b] 上对 计算定积分即为二重积分。具体处理办法是设:

y

d ( x ) b

ϕ(x ) ∫ f (x , y )dy (1),I ∫ϕ(x )dx (2)

c ( x ) a

计算二重积分转化为计算单次积分 (1)、(2)。

2.定积分的数值积分公式[2]

1

关于定积分,有Newton-Leibniz 公式:

b

∫a f (x )dx F (b ) =−F (a )

在实际问题中,往往会遇到被积函数的原函数无法用初等函数来表示,或有

的虽然能用初等函数表示,但过分复杂,这些情形都需要去建立定积分的近似计

算公式,也就是求积分的数值解。

定义数值积分如下:I n ( f ) 是离散点上的函数值的线性组合

matlab二重积分有奇异点,用MATLAB计算某些区域上二重积分.pdf相关推荐

  1. matlab trapz二重积分函数_如何使用 MATLAB 求解定积分、不定积分和多重积分问题...

    介绍几种 MATLAB 中求解积分的方法,首先是采用符号积分的方法,求解积分的符号函数为 \(int\),使用方法如下: int(f,x,a,b) 例如计算如下积分函数: $$\int \frac{1 ...

  2. matlab trapz二重积分函数_科学网—MATLAB中的数值积分方法 - 王福昌的博文

    实际应用中在MATLAB里面都有开发好的命令可以使用,如  quad(), quadl(),quad2d(),triplequad() .需要掌握这些命令的用法. 1. 定积分 trapz(),qua ...

  3. MATLAB实战系列(三十七)-MATLAB基于PQ解耦风电场并网潮流计算

    前言 文中涉及代码请参见 电力系统仿真-MATLAB基于PQ解耦风电场并网潮流计算 IEEE30节点.14节点.4节点标准算例,潮流计算的功能:风力发电机组并网潮流计算,并网对大电网的影响. 以下是我 ...

  4. 2021-01-07 matlab数值分析  矩阵特征值与特征向量的计算 改进乘幂法 反幂法

    matlab数值分析  矩阵特征值与特征向量的计算 1改进乘幂法 function [t,y]=eigIPower(A,v0,ep) [tv,ti]=max(abs(v0)); lam0=v0(ti) ...

  5. matlab鬼成像,一种基于光计算的可视化计算鬼成像系统及成像方法与流程

    本发明涉及成像探测领域,具体为一种基于光计算的可视化计算鬼成像系统及成像方法. 背景技术: 光学鬼成像技术是通过双路光信号符合探测恢复待测物体空间信息实现的.其特点是包含物体信息的信号光(signal ...

  6. matlab广泛用于什么,目前,Matlab已成为世界上最受欢迎的科学和工程计算软件工具....

    目前,Matlab已成为世界上最受欢迎的用于科学和工程计算的软件工具. 它是一种新型的计算机高级编程语言,具有广阔的应用前景. 有人将其称为"第四代"计算机语言. 就影响而言,仍然 ...

  7. matlab 电力系统分析 毕设,基于MATLAB的P-Q分解法电力系统潮流计算毕业设计

    基于MATLAB的P-Q分解法电力系统潮流计算毕业设计 摘 要 电力系统潮流计算是研究电力系统稳定运行情况的一种重要的计算,在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用它来定量地分析比 ...

  8. matlab线圈磁场,利用Matlab计算螺线管内磁场分布研究.pdf

    利用Matlab计算螺线管内磁场分布研究.pdf ( ) 2010 12 688 Journal of Beijing Normal University ( Natural Science) 46( ...

  9. matlab p q的确定,基于MATLAB软件的P-Q分解法潮流计算 (1)

    基于MATLAB软件的P-Q分解法潮流计算 摘要 电力系统潮流计算是研究电力系统稳态运行情况的一种重要的分析计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态:各母线的电压,各元 ...

  10. matlab如果不想立即在,科学计算与MATLAB 1.5

    单元测验已完成成绩:100.0分 1 [单选题] MATLAB一词来自()的缩写. 窗体顶端 A. MathematicaLaboratory B. MatrixLaboratory C. MathW ...

最新文章

  1. python子进程通信_python执行子进程实现进程间通信的方法
  2. js判断ipad还是安卓_还考虑iPad?荣耀平板V6麒麟985支持5G才3000多元
  3. adb push和pull使用
  4. 华为3com交换路由命令详解
  5. 下载Nacos源码并运行
  6. JPA /休眠刷新策略初学者指南
  7. linux脚本获取usb设备,Linux基于USB端口执行脚本
  8. 复制release文件到另一台电脑.exe文件无法运行_电脑技巧:电脑版微信双开(或微信多开)?用start指令可以解决...
  9. Java常用接口与类——String类、StringBuffer类、StringBuilder类
  10. javaweb学习总结五(内省、beanUtils工具包)
  11. 安卓listview点击空白事件_要权限才给用?安卓毒瘤APP滚蛋吧!
  12. Oracle 安装报错 [INS-06101] IP address of localhost...
  13. java是什么类型语言_java属于什么类型语言
  14. t检验自由度的意义_统计学常用概念:T检验、F检验、卡方检验、P值、自由度...
  15. Thinkpad T410加内存
  16. Ubuntu解压安装包及make命令相关操作(针对make: *** [btest] Error 127以及make: *** [btest] Error 1)
  17. nginx openresty DNS resolver配置实例,通过配置resolver解决proxy_pass中使用变量参数,高性能负载均衡 NGINX Plus 中 RESTful API
  18. 乒乓球比赛赛程_2018国际乒联乒乓球重大赛事详细赛程表,看点爆点早知道
  19. 帝都某外企裁员,补偿673252.62元 !准备回家躺平了 !
  20. 毕业一周年总结-不忘初心,砥砺前行

热门文章

  1. python工程师简历项目经验怎么写_班长项目经验简历范文
  2. 物联网(IoT)课程
  3. 转自博客园好文---关于程序员今后的道路规划,借鉴并勉励一下
  4. Tomcat 端口被javaw.exe占用 有效解决方法
  5. python数据挖掘电影评分分析_豆瓣电影数据与票房数据分析
  6. DTcms-【知识点】-知识点缴获
  7. 求极限的若干方法经验和教训
  8. 龙果学院mysql分布式集群代码_MySQL分布式集群搭建
  9. 2011年上半年国内优秀初创企业产品汇总
  10. div嵌套的div水平垂直居中