一、数值积分问题的数学表达

$\int_{1}^{2}f(x)dx\approx \sum_{k=0}^{n}A_{k}f(x_{k}))$

$x_{k}$ 称为积分节点， $A_{k}$ 称为求积系数。右端公式称为左端定积分的某个数值积分。MATLAB支持三重以下的积分运算。分别为：

$\int_{x_{min}}^{x_{max}}f(x)dx$

$\int_{x_{min}}^{x_{max}}\int_{y_{min}}^{y_{max}}f(x,y)dxdy$

$\int_{x_{min}}^{x_{max}}\int_{y_{min}}^{y_{max}}\int_{z_{min}}^{z_{max}}f(x,y,z)dxdydz$

在计算积分时，要求积分区间是确定的。

二、一元函数的数值积分

1、quad函数：采用遍历的自适应辛普森法计算函数的数值积分。

适用于精度要求低，被积函数平滑性差的数值积分。

Q = quad(FUN,A,B)
Q = quad(FUN,A,B,TOL)
Q = quad(FUN,A,B,TOL,TRACE)
[Q,FCNT] = quad(...)

FUN为被积函数的句柄，FUN应该接收向量输入，并输出相同长度的向量。

A，B分别是积分的起始值和结束值。

TOL用于控制自适应辛普森法的误差，增大TOL可以加快计算速度，但是计算精度下降。

FCNT为函数计算的次数。

2、quadl函数：采用遍历的自适应Lobatto法计算函数的数值积分。

适用于精度要求高，被积函数曲线比较平滑的数值积分。

用法与quad相同。

Q = quadl(FUN,A,B)
Q = quadl(FUN,A,B,TOL)
Q = quadl(FUN,A,B,TOL,TRACE)
[Q,FCNT] = quadl(...)

quad有较快的计算速度，但是准确性较差；而quadl函数需要更多的计算时间，但具有较高的准确性。

三、quadv函数：

一次可以计算多个一元函数的数值积分。k取不同的值，积分结果不尽相同。

$\int_{0}^{1}x^{_k} dx,k=1,2,\cdot \cdot \cdot ,n$

Q = quadv(FUN,A,B)
Q = quadv(FUN,A,B,TOL)
Q = quadv(FUN,A,B,TOL,TRACE)
[Q,FCNT] = quadv(...)

quad ，quadl，quadv三个函数都要求被积函数FUN必须是函数句柄，同时积分限必须是有限的。

例题：计算 $\int_{0}^{2}\frac{1}{x^{3}-2x-5}dx$

clear all;clc;
F = @(x) 1./(x.^3-2*x-5);
Q = quadl(F,0,2)Q =-0.4605

二、二重数值积分

二重数值积分计算函数dblquad

dbquad函数可以用来计算被积函数在积分区域 $x\in \begin{bmatrix}x_{min}, &x_{max} \end{bmatrix},y\in \begin{bmatrix}y_{min}, &y_{max} \end{bmatrix}$ 内的积分数值。该函数先计算内积分值，然后利用内积分的中间结果来计算二重积分。

dxdy：所以x为内积分变量，y为外积分变量。

dblquad函数的调用格式：

Q=dblquad(FUN,XMIN,XMAX,YMIN,YMAX)
Q=dblquad(FUN,XMIN,XMAX,YMIN,YMAX,TOL)
Q=dblquad(FUN,XMIN,XMAX,YMIN,YMAX,TOL,@QUADL)
Q=dblquad(FUN,XMIN,XMAX,YMIN,YMAX,TOL,MYQUADF)

FUN是被积函数的句柄。

XMIN,XMAX,YMIN,YMAX分别为x，y两个方向上的积分限。

TOL是绝对计算精度。

@QUADL和MYQUADF指定计算以为积分时使用的函数。

例题：计算积分 $\int_{0}^{\pi }\int_{\pi}^{2\pi}(ysinx+3cosy-1)dxdy$

clear all;clc;
F = @(x,y) y*sin(x)+3*cos(y)-1;
xmin=pi;
xmax=2*pi;
ymin=0;
ymax=pi;
q=dblquad(F,xmin,xmax,ymin,ymax)q =-19.7392

三、三重数值积分

triplequad函数可以用来计算被积函数在空间区域 $x\in \begin{bmatrix}x_{min}, &x_{max} \end{bmatrix},y\in \begin{bmatrix}y_{min}, &y_{max} \end{bmatrix},z\in \begin{bmatrix}z_{min}, &z_{max} \end{bmatrix}$ 内的数值积分。

其调用格式如下：

Q=triplequad(FUN,XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX)
Q=triplequad(FUN,XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX,TOL)
Q=triplequad(FUN,XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX,TOL,MYQUADF)
Q=triplequad(FUN,XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX,TOL,@QUADL)

上式中相关参数的含义同dblquad函数。

【MATLAB】数据分析之函数数值积分相关推荐

MATLAB数据分析2
MATLAB数据分析2 %特殊矩阵%通用的特殊矩阵 %zeros函数:产生全0矩阵,即零矩阵 %格式: %zeros(m):产生mxm零矩阵 %zeros(m,n):产生mxn零矩阵 %zeros(s ...
MATLAB数据分析
MATLAB数据分析 %求解x^2 - 3x + 1 = 0 的根%方法一: % 求根 % p = [1,-3,1]; % x = roots(p);%绘制图像 % x = -5 : 0.1 : 5; ...
MATLAB 数据分析方法（第2版）1.2　MATLAB基础概述
1.2 MATLAB基础概述 1.2.1 MATLAB的影响 MATLAB源于Matrix Laboratory,即矩阵实验室,是由美国Mathworks公司发布的主要面对科学计算.数据可视化.系统仿 ...
MATLAB数据分析（插值运算和曲线拟合）
MATLAB数据分析(插值运算和曲线拟合) MATLAB数据分析在数学研究和工程应用中,我们常常需要对某些环节进行反复的测验或测试,以获取更多的数据支持最优决策,但是由于环境条件的限制,我们往往只能 ...
matlab多项式计算题目,MATLAB数据分析与多项式计算-习题答案
<MATLAB数据分析与多项式计算-习题答案>由会员分享,可在线阅读,更多相关<MATLAB数据分析与多项式计算-习题答案(4页珍藏版)>请在人人文库网上搜索. 1.第6章 M ...
matlab exp(),matlab的exp函数
matlab中exp用法的问题请问下面式子中的.有什么用,这里为什么要 x=(exp(-1*tx)).*rectpuls((texp(a),如a是一标量(即一常数),返回自然常数的a次方:如a是以矩 ...
MATLAB 数据分析
MATLAB 数据分析多项式的估值,可用polyval函数: x=-1:0.01:1; %生成自变量 g=[1 3 5 7 9]; %函数对应的多项式 h=polyval(g,x); %进行估值运算 ...
【MATLAB数据分析】02数据的中位数和分位数
引言在[MATLAB数据分析]01数据的均值.变异度.偏度和峰度一文中我们介绍了均值.方差.标准差等数字特征,它们都是总体相应特征值的一种矩估计,更适合来自正态分布的数据的分析.但若总体的分布未知, ...
matlab离散型数据求和,MATLAB程序设计教程(6)---MATLAB数据分析与多项式计算
第6章 MATLAB数据分析与多项式计算 6.1 数据统计处理 6.2 数据插值 6.3 曲线拟合 6.4 离散傅立叶变换 6.5 多项式计算 6.1 数据统计处理 6.1.1 最大值和最小值 MAT ...

【MATLAB】数据分析之函数数值积分