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…

element-ui+vue2实现粘贴上传

element-uivue2实现粘贴上传 <style scoped lang"scss">.img-upload{position: relative;display: inline-block;margin-right: 9px;}.image {width: 100px;height: 100px;margin-right: 9px;}.image:last-child{margin-right: 0;}.img-upload .el-upload--pic…

Java学习|初识热加载

文章目录 引言Java热加载概念热加载与传统部署的区别热加载的好处风险与挑战 技术背景类加载机制类加载器层次结构 实现方法1. 使用Java Agent示例代码 2. 利用JRebel (XRebel)配置使用场景 3. Spring Boot DevTools配置使用场景 4. 动态类加载示例代码 5. JIT编译器的热替换示…

【嵌入式英语教程--7】C语言中的数据类型与内存管理

C语言中的数据类型与内存管理 英文原文 In the C programming language, data types define the kind of data that can be stored in variables. Common data types include integers, floating-point numbers, characters, and strings. The choice of data type has impli…

Python面试整理-第三方库

Python社区提供了大量的第三方库,这些库扩展了Python的功能,覆盖了从数据科学到网络应用开发等多个领域。以下是一些非常流行和广泛使用的第三方库: 1. NumPy ● 用途:数值计算。 ● 特点:提供了一个强大的N维数组对象和大量用于数学运算的函数。 ● 应用场景:科学计算、…

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

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

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

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

shell语言中的 、;、|有什么不同?

在 shell 脚本语言中&#xff0c;&&、; 和 | 是三种不同的命令分隔符和控制结构符号&#xff0c;它们分别用于不同的目的&#xff1a; && (AND 运算符) 这个运算符用于连接两个命令&#xff0c;其中第二个命令只有在第一个命令成功执行后才会运行。 例如&am…

什么是线程安全?

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

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

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

Springboot 不同定时任务实现及场景

实现1、通过开启EnableScheduling 及注解Scheduled 实现定时执行任务 【完整示例】 package org.javatrip.springboottimer;import org.springframework.scheduling.annotation.EnableScheduling; import org.springframework.scheduling.annotation.Scheduled; import org.s…

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

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

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

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

平台数据脱敏方案

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

如何使用短链接生成接口

一、什么是短链接&#xff1f; 专业用于将长网址缩短&#xff0c;支持短域名定制&#xff0c;支持html5&#xff0c;l0S&#xff0c;Android&#xff0c;短信&#xff0c;微博分享&#xff0c;抖音分享&#xff0c;且短网址生成微信防屏蔽&#xff0c;360防屏蔽。智能域名分组…

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

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

M3U8流视频数据爬虫

M3U8流视频数据爬虫 HLS技术介绍 现在大部分视频客户端都采用HTTP Live Streaming&#xff08;HLS&#xff0c;Apple为了提高流播效率开发的技术&#xff09;&#xff0c;而不是直接播放MP4等视频文件。HLS技术的特点是将流媒体切分为若干【TS片段】&#xff08;比如几秒一段…

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

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

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)//控制行…