quadl matlab,MATLAB求一元函数的数值积分(quad,quadl,quadv)，大师来详解

MATLAB提供了quad( ),quadl( )等函数命令计算一元函数的数值积分，其中quad( )采用遍历的自适应Simpson法计算函数的数值积分，quadl采用遍历的自适应Lobatto法计算函数的数值积分。其调用格式一般为[Q,Fcnt]=quad(function,a,b)，其中function为被积函数(形式为函数句柄/匿名函数)，a,b分别积分上限，[Q,Fcnt]分别返回数值积分的结果和函数计算的次数。

工具/材料

MATLAB

quad

quadl

操作方法

第一，计算下图两个一元函数的数值积分。

第二，启动MATLAB，新建脚本(Ctrl+N)，输入如下代码：

close all; clear all; clcformat compactfigure('Position',[50,50,600,500],'Name','sin(0.5x)')f1=@(x) sin(0.5*x);[Q1a,Fcnt1a]=quad(f1,0,2)[Q1b,Fcnt1b]=quadl(f1,0,2)fplot(f1,[0,2]);%-------------------%figure('Position',[100,100,600,500],'Name','e^x+1')f2=@(x) exp(x)+1;[Q2a,Fcnt2a]=quad(f2,-1,2)[Q2b,Fcnt2b]=quadl(f2,-1,2)fplot(f2,[-1,2]);

其中f1=@(x) sin(0.5*x)，f2=@(x) exp(x)+1是建立被积分的函数(形式为函数句柄/匿名函数)。

[Q,Fcnt]=quad(function,a,b)，quadl(function,a,b)是分别采用quad和quadl两种方法求数值积分。

其中被积函数function要为上述建立的函数句柄/匿名函数，a,b分别积分上限，[Q,Fcnt]分别返回数值积分的结果和函数计算的次数。

第三，保存和运行上述脚本，在命令行窗口(Command Window)得到如下结果：

Q1a = 0.9194Fcnt1a = 13Q1b = 0.9194Fcnt1b = 18Q2a = 10.0212Fcnt2a = 37Q2b = 10.0212Fcnt2b = 18

可以看出，采用quad( )和quadl两种方法计算函数(1)的数值积分都是0.9194，函数(2)的数值积分都是10.0212，只是计算次数不同。例如采用quad( )计算函数(1)的数值积分，计算次数为13，采用quadl( )计算次数为18。

第四，下图为被积分函数(1)的图像。

第五，下图为被积分函数(2)的图像，其中fplot( )绘图命令可以绘制函数句柄/匿名函数的图像。

好了，以上就是大致内容了，(END)

quad，quadl要求被积函数必须为函数句柄/匿名函数，即用f=@(x)创建的函数。

quad( )计算一元函数数值积分的速度快但准确性稍差，quadl( )计算速度慢但准确性强。

声明：该文观点仅代表作者本人，天晴资讯网系信息发布平台，仅提供信息存储空间服务。

quadl matlab,MATLAB求一元函数的数值积分(quad,quadl,quadv)，大师来详解相关推荐

java n%9==0_用C++实现求N!中末尾0的个数的方法详解
题目描述: 输入一个正整数n,求n!(即阶乘)末尾有多少个0? 比如: n = 10; n! = 3628800,所以答案为2 输入描述: 输入为1行,n(1≤n≤1000) 输出描述: 输出一个整数 ...
matlab输出李萨育图形,李萨如图形的详解与应用
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区. 您需要登录才可以下载或查看,没有帐号?注册 x 本帖最后由 Freeman_2016 于 2019-9-11 11:18 编辑李萨如图形 ...
matlab的NLP功能,pyhanlp 共性分析与短语提取内容详解
pyhanlp 共性分析与短语提取内容详解简介 HanLP中的词语提取是基于互信息与信息熵.想要计算互信息与信息熵有限要做的是文本分词进行共性分析.在作者的原文中,有几个问题,为了便于说明,这里首 ...
matlab梯形弦长公式,弦长公式证明及应用详解
弦长公式证明及应用详解公式为: |AB| 和:|AB|= 作用:应用弦长公式很方便,它所解决的问题是求直线与所有圆锥曲线所交弦的弦长,因为直线的斜率往往是已知的,这样再知道两个交点的横坐标或者纵坐标 ...
如何求matlab的in(2.0375),matlab2012年课后习题及答案详解
2012年12月7日星期五第2章 MATLAB 矩阵运算基础 2.1 在MA TLAB 中如何建立矩阵?? ? ???194375,并将其赋予变量a ? >> a=[5 7 3;4 9 ...
matlab函数表达式里分号_MatLab常见函数和运算符号使用详解
基本运算 convhull :凸壳函数 cumprod :累计积 cumsum :累计和 cumtrapz :累计梯形数值积分 delaunay :Delaunay三角化 dsearch :求最近点 ...
matlab 小波变换_连续小波变换实现方法的总结及其程序详解
在帖子"给大家分享我自己编的程序-连续小波变换" 中,pengzk版友给出了morlet小波变换的源代码,但其中的许多参数和语句意义不够明确,这就给一些希望了解连续小波变换实现方法 ...
【MATLAB】Parzen窗与K近邻算法原理与代码详解
文章目录 1.非参数估计原理 2.Parzen窗 2.1.算法原理 2.2.Matlab实现与参数探究 3.K近邻 3.1.算法原理 3.2.Matlab实现与参数探究 1.非参数估计原理 \qqua ...
matlab双目相机标定校正_双目相机的标定过程详解！-----MATLAB
基于双目视觉的测距.三维重建等过程中的第一步就是要进行标定.双目相机的标定过程在网上有很多资料,但是基本都没有matlab官方网址讲的好.所以请参考MATLAB官方文档:https://ww2.mat ...

quadl matlab,MATLAB求一元函数的数值积分(quad,quadl,quadv)，大师来详解

quadl matlab,MATLAB求一元函数的数值积分(quad,quadl,quadv)，大师来详解相关推荐

最新文章

热门文章