Matlab编程资源库(15)数值积分

一、基本原理

       求解定积分的数值方法多种多样,如简单的梯形法、辛普生(Simpson)法、牛顿-柯特斯(Newton-Cotes)法等都是经常采用的方法。它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1] i=1,2,…,n其中 x 1=a xn+1=b。这样求定积分问题就分解为求和问题。

二、实现方法

1.变步长辛普生法

       基于变步长辛普生法,MATLAB给出了quad函数来求定积分。该函数的调用格式为:

[I,n]=quad('fname',a,b,tol,trace)

       其中fname是被积函数名。a 和 b分别是定积分的下限和上限。tol用来控制积分精度,缺省时取tol=0.001。trace控制是否展现积分过程,若取非0则展现积分过程,取0则不展现,缺省时取trace=0。返回参数I定积分值,n为被积函数的调用次数。

2.牛顿-柯特斯法
基于 牛顿-柯特斯法,MATLAB给出了 quad8函数来求定积分。该函数的调用格式为:
[I,n]=quad8('fname',a,b,tol,trace)
其中参数的含义和quad函数相似,只是tol的缺省值取10-6。该函数可以更精确地求出定积分的值,且一般情况下函数调用的步数明显小于quad函数,从而保证能以更高的效率求出所需的定积分值。
3 .被积函数由一个表格定义
       在 MATLAB 中,对由 表格形式 定义的函数关系的求定积分问 题用 trapz(X,Y) 函数。其中向量 X,Y 定义函数关系 Y=f(X)

8-1 求定积分。

(1) 建立被积函数文件fesin.m

function f=fesin(x)
f=exp(-0.5*x).*sin(x+pi/6);
(2) 调用数值积分函数 quad 求定积分。
[S,n]=quad('fesin',0,3*pi)
S =0.9008
n =77

运行结果 :

8-2 求定积分。
(1) 被积函数文件 fx.m
function f=fx(x)
f=x.*sin(x)./(1+cos(x).*cos(x));
(2) 调用函数 quad8 求定积分。
I=quad8('fx',0,pi)
I =2.4674

 运行结果 :

8-3 分别用 quad 函数和 quad8 函数求定积分的近 似值,并在相同的积分精度下,比较函数的调用
次数。 调用函数 quad 求定积分:
format long;
fx=inline('exp(-x)');
[I,n]=quad(fx,1,2.5,1e-10)
I =0.28579444254766
n =65

 运行结果 :

调用函数 quad8 求定积分:
format long;
fx=inline('exp(-x)');
[I,n]=quad8(fx,1,2.5,1e-10)
I =0.28579444254754
n =33
8-4 trapz 函数计算定积分。
命令如下:
X=1:0.01:2.5;
Y=exp(-X); %生成函数关系数据向量
trapz(X,Y)
ans =0.28579682416393

 运行结果 :

三、二重定积分的数值求解

       使用 MATLAB 提供的 dblquad 函数 就可以直接求出上述 二重定积分的数值解 该函数的调用格式为:
I=dblquad(f,a,b,c,d,tol,trace)
       该函数求 f(x,y) [a,b] × [c,d] 区域上 的二重定积分。参数 tol trace 用法与函数 quad 完全相同。
8-5 计算二重定积分
(1) 建立一个函数文件 fxy.m
function f=fxy(x,y)
global ki;
ki=ki+1; %ki用于统计被积函数的调用次数
f=exp(-x.^2/2).*sin(x.^2+y);
(2) 调用 dblquad 函数求解。
global ki;ki=0;
I=dblquad('fxy',-2,2,-1,1)
ki
I =1.57449318974494
ki =1038

运行结果 :

 结语 

别为不值得的人,少了微笑

别为鸡毛蒜皮的事,添了烦恼

努力做一个向阳而生的人

!!!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/877285.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

2024年PINN网络​还在火!发论文侧重点在哪儿?

2024年了,PINN网络依然火爆,各大顶会顶刊都能看见它的相关论文。 这是因为,AI交叉学科通常离不开求解偏微分方程PDE,而传统的求解方法受初始假设限制,一旦没设好就会导致很大的误差。 PINN作为一种新的思路&#xff…

气象水文耦合模WRF-Hydro建模技术

原文链接:气象水文耦合模WRF-Hydro建模技术https://mp.weixin.qq.com/s?__bizMzUzNTczMDMxMg&mid2247610398&idx4&sn34b4bbed4c74dcbbb0ac19ef8dcdaaff&chksmfa8271f9cdf5f8ef34ea6f721736a2fbbf8be896744ab7e46caa571c52a30628f056b4bd6964&t…

又一新AI搜索工具,OpenAI 推出新的搜索方式 SearchGPT

系列文章目录 每天推荐AI工具系列文章回顾: 选择 haiyi海艺图像生成、LoRA、模型的使用和训练网站 tusiart吐司艺术图像生成、LoRA 模型的使用和训练网站 解锁AI创造力的无限可能:探索Vivago.ai的革命性功能 文章目录 系列文章目录前言一、SearchGPT…

<数据集>手机识别数据集<目标检测>

数据集格式:VOCYOLO格式 图片数量:16172张 标注数量(xml文件个数):16172 标注数量(txt文件个数):16172 标注类别数:1 标注类别名称:[Phone] 使用标注工具:labelImg 标注规则:…

什么是线程安全?

什么是线程安全? 为什么需要线程安全?如何实现线程安全?1. 排队干活2. 自己带工具3. 用现成的安全工具 4、示例5、总结 💖The Begin💖点点关注,收藏不迷路💖 在编程里,特别是当程序能…

推荐一款专注批量推送消息的轻量工具,支持主流平台的消息推送,简单、高效、低成本(附源码)

前言 在数字化时代,企业和个人面临着日益增长的消息推送需求。然而,现有的推送处理方案往往存在一些挑战和不足,如cao作复杂、成本高昂、缺乏灵活性等。这些问题不仅影响了推送效率,也增加了用户的负担。此外,随着工作…

华为od 100问 持续分享10-华为OD的面试流程细说

我是一名软件开发培训机构老师,我的学生已经有上百人通过了华为OD机试,学生们每次考完试,会把题目拿出来一起交流分享。 重要:2024年5月份开始,考的都是OD统一考试(D卷),题库已经整…

Matlab编程资源库(16)数值微分

一、数值差分与差商 在Matlab中,数值差分与差商是数值分析中常用的概念,尤其在求解微分方程、插值、逼近等领域有广泛应用。下面简要介绍这两个概念及其在Matlab中的实现。 数值差分 数值差分是微分运算的离散化形式,用于近似求解导数。给定…

平台数据脱敏方案

在目前大环境下,这几年做事业政府单位的信息化项目,都特别强调安全,原因大伙都清楚。 安全包含两块,一是框架组件安全,二是业务信息安全。 框架组件安全一般就是漏洞修复,组件升级到对应没有漏洞的版本。 业…

数据结构【有头双向链表】

目录 实现双向链表 双向链表数据 创建双向链表 初始化双向链表创建(哨兵位) 尾插 打印双向链表 头插 布尔类型 尾删 头删 查询 指定位置后插入 指定位置删除数据 销毁 顺序表和链表的分析 代码 list.h list.c test.c 注意&#xff1a…

清华学姐熬夜肝了15天的软件测试面试题出炉(附答案)建议收藏!

一、Web自动化测试 1.Selenium中hidden或者是display = none的元素是否可以定位到? 不能,可以写JavaScript将标签中的hidden先改为0,再定位元素 2.Selenium中如何保证操作元素的成功率?也就是说如何保证我点击的元素一定是可以…

C:图案打印

引言 本篇文章讲了一些常见的图形编程题&#xff0c;并总结了一些规律。 1、打印空心正方形 1.1 代码展示&#xff1a; #include<stdio.h> int main() {int a 0;//边长初始化scanf("%d", &a);//输入边长的值{int i 0;for (i 0; i < a; i)//控制行…

知识图谱增强的RAG(KG-RAG)详细解析

转自&#xff1a;知识图谱科技 这是一个与任务无关的框架&#xff0c;它将知识图谱&#xff08;KG&#xff09;的显性知识与大型语言模型&#xff08;LLM&#xff09;的隐含知识结合起来。这是该工作的arXiv预印本 https://arxiv.org/abs/2311.17330 。 我们在这里利用一个名为…

自定义表格_可拖拽排序

在做后台管理系统的时候&#xff0c;经常需要表格里面的每行排序&#xff0c;自定义可拖拽表格&#xff0c;更改样式方便。 一、实现效果 进行拖拽演示&#xff1a; 可拖拽排序表格 无滚动条样式&#xff1a; 有滚动条样式&#xff1a; 二、代码 使用reactscssts,实现页面。 …

Linux(CentOS)ftp服务搭建

ftp服务器搭建 1. 下载ftp服务2. 查找ftp配置文件3. 查看配置文件信息4. Windows连接ftp服务1&#xff09;使用文件资源管理器连接2&#xff09;使用FlashFXP工具&#xff0c;比文件资源管理器好用&#xff0c;强烈推荐 5. Linux连接 1. 下载ftp服务 yum install -y vsftpd2. …

Docsify:快速用Markdown文档搭建网站的利器

Github官方地址&#xff1a;Docsify 什么是Docsify&#xff1f; 对于经常写博客的人来说&#xff0c;markdown大家都不陌生。今天介绍一个在最近需求中碰到的软件Docsify&#xff0c;通过它能够将Markdown直接转换为网页。话不多说&#xff0c;下面直接介绍它的快速用法。 D…

Apache、nginx

一、Web 1、概述 Web&#xff1a;为⽤户提供的⼀种在互联⽹上浏览信息的服务&#xff0c;Web 服务是动态的、可交互的、跨平台的和图形化的。 Web 服务为⽤户提供各种互联⽹服务&#xff0c;这些服务包括信息浏览服务&#xff0c;以及各种交互式服务&#xff0c;包括聊天、购物…

fastapi教程(五):中间件

一&#xff0c;什么是中间件 中间件是一种软件组件&#xff0c;它在请求到达应用程序处理程序之前和/或响应发送回客户端之前执行操作。 请求从客户端发出。 请求首先经过Middleware 1。 然后经过Middleware 2。 请求到达FastAPI路由处理器。 响应从路由处理器返回。 响应经过…

如何通过 CloudCanal 实现从 Kafka 到 AutoMQ 的数据迁移

01 引言 随着大数据技术的飞速发展&#xff0c;Apache Kafka 作为一种高吞吐量、低延迟的分布式消息系统&#xff0c;已经成为企业实时数据处理的核心组件。然而&#xff0c;随着业务的扩展和技术的发展&#xff0c;企业面临着不断增加的存储成本和运维复杂性问题。为了更好地…

《LeetCode热题100》---<双指针篇四道>

本篇博客讲解LeetCode热题100道双指针篇中的 第一道&#xff1a;移动零&#xff08;简单&#xff09; 第二道&#xff1a;盛最多水的容器&#xff08;中等&#xff09; 第一道&#xff1a;移动零&#xff08;简单&#xff09; class Solution {public void moveZeroes(int[] nu…