matlab实验求不定积分函数,实验四用matlab计算积分-程序员宅基地

技术标签: matlab实验求不定积分函数  

41528d3028836879cd698677c3999917.gif实验四用matlab计算积分

实验四 用matlab计算积分 4.1积分的有关理论 定积分:积分是微分的无限和,函数在区间上的积分定义为 其中 从几何意义上说,对于上非负函数,记分值是曲线与直线及轴所围的曲边梯形的面积。有界连续(或几何处处连续)函数的积分总是存在的。 微积分基本定理(Newton-Leibniz公式):在上连续,且,则有 这个公式表明导数与积分是一对互逆运算,它也提供了求积分的解析方法:为了求的定积分,需要找到一个函数,使的导数正好是,我们称是的原函数或不定积分。不定积分的求法有学多数学技巧,常用的有换元积分和分部积分法。从理论上讲,可积函数的原函数总是存在的,但很多被积函数的原函数不能用初等函数表示,也就是说这些积分不能用解析方法求解,需用数值积分法解决。 在应用问题中,常常是利用微分进行分析,而问题最终归结为微分的和(即积分)。一些更复杂的问题是含微分的方程,不能直接积分求解。 多元函数的积分称为多重积分。二重积分的定义为 当非负时,积分值表示曲顶柱体的体积。二重积分的计算主要是转换为两次单积分来解决,无论是解析方法还是数值方法,如何实现这种转换,是解决问题的关键。 4.2积分的数值方法 梯形法:将划分为若干小区间则 在每一小区间上近似为一直线,用弦代替,有 从而 称为梯形公式。通常将区间等分,, 可以证明,当时由上述公式给出的梯形法是收敛的。 重积分:重积分的数值计算可通过若干次单积分的组合实现,如对于二重积分 先化为二次计分 利用梯形法,先将区间等分,利用梯形积分公式可得 再将区间等分,利用梯形积分公式可得 4.3积分的MATLAB命令 MATLAB中主要用int进行符号积分,用trapz,dblquad,quad,quad8等进行数值积分。 R=int(s,v) %对符号表达式s中指定的符号变量v计算不定积分.表达式R只是表达式函数s的一个原函数,后面没有带任意常数C. R=int(s) %对符号表达式s中确定的符号变量计算计算不定积分. R=int(s,a,b) %符号表达式s的定积分,a,b分别为积分的上、下限 R=int(s,x,a,b) %符号表达式s关于变量x的定积分,a,b分别为积分的上、下限 trapz(x,y) 梯形积分法,x时表示积分区间的离散化向量,y是与x同维数的向量,表示被积函数,z返回积分值。 fblquad(‘fun’,a,b,c,d) 矩形区域二重数值积分,fun表示被积函数的M函数名,a,b分别为x的上、下限,c,d分别为y的上、下限. 可以用help int, help trapz, help quad等查阅有关这些命令的详细信息 例1 用符号积分命令int计算积分. MATLAB代码为: >>clear; syms x; >>int(x^2*sin(x)) 结果为 ans =-x^2*cos(x)+2*cos(x)+2*x*sin(x) 如果用微分命令diff验证积分正确性,MATLAB代码为: >>clear; syms x; >>diff(-x^2*cos(x)+2*cos(x)+2*x*sin(x)) 结果为 ans =x^2*sin(x) 例2 计算数值积分. 先用梯形积分法命令trapz计算积分,MATLAB代码为: >>clear; x=-2:0.1:2; y=x.^4; %积分步长为0.1 >>trapz(x,y) 结果为 ans = 12.8533 实际上,积分的精确值为。如果取积分步长为0.01, MATLAB代码为: >>clear; x=-2:0.01:2; y=x.^4; %积分步长为0.01 >>trapz(x,y) 结果为 ans =12.8005 可用不同的步长进行计算,考虑步长和精度之间的关系。一般说来,trapz是最基本的数值积分方法,精度低,适用于数值函数和光滑性不好的函数. 如果用符号积分法命令int计算积分, 输入MATLAB代码为: >>clear; syms x; >>int(x^4,x,-2,2) 结果为 ans =64/5 例3 计算数值积分,可将此二重积分转化为累次积分 输入MATLAB代码为: >>clear; syms x y; >>iy=int(1+x+y,y,-sqrt(1-x^2),sqrt(1-x^2)); >>int(iy,x,-1,1) 结果为 ans =pi 例4(广义积分) 计算广义积分。 输入MATLAB代码为: >>syms x; >>y=int(exp(sin(x)-x^2/50),-inf,inf); >>vpa(y,10) 结果为15.86778263。 习题16-4 1.(不定积分)用int计算下列不定积分,并用diff验证 ,,,, 2.(定积分)用trapz,int计算下列定积分 ,,, 3.(椭圆的周长) 用定积分的方法计算椭圆的周长 4.(二重积分)计算积分 5. (广义积分)计算广义积分 ,,

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_42524401/article/details/116412235

智能推荐

java后台利用POI将excle转换成html实现在线预览_javaexcel转html-程序员宅基地

文章浏览阅读6.6k次。 前一阵项目有office系列文件在线预览需求,所以查询了一些资料,参考其他一些博客,实现了通过POI将excle转化为html,后来需求被砍掉,没有继续深入研究,这里将查询到的一些资料做一个记录. 首先实现office系列文件在线预览主要查到有如下解决方案:1.flash 的flexpaper 将文档转换为swf格式,然后使用flash在网页中浏览 2.使用开源的软件openoffice+p..._javaexcel转html

html热点 位置不低,热点图 页面缩放不改变位置-程序员宅基地

文章浏览阅读281次。文章时间:2019年5月29日 23:42:02解决问题:热点图 随着页面的缩放比例而不改变其所在的位置推荐画图工具:Adobe Dreamweaver CC头部我们需要3个东西,请直接复制即可直接复制这段代码然后用设计模式画图就行了。任意地方加这段js代码。$("img[usemap]").each(function () {var img = $(this);var newImg = new ..._h5地图热点放大缩小如何不跑偏

maven打包报错_error locating assembly descriptor: src/main/assem-程序员宅基地

文章浏览阅读2k次。项目场景:maven打包报错:Failed to execute goal org.apache.maven.plugins:maven-assembly-plugin:2.4:single (make-assembly) on project hiveudf2: Error reading assemblies: Error locating assembly descriptor: assembly.xml问题描述:项目上执行目标org.apache.maven.plugins:maven-as_error locating assembly descriptor: src/main/assembly/assembly.xml

Chrome开发者工具使用教程-主要功能简介(一)_谷歌后台调测执行sql-程序员宅基地

文章浏览阅读1.4k次。Chrome开发者工具是谷歌浏览器自带的一款开发者工具,它可以给开发者带来很大的便利。主要作用有:快速定位问题查看站点页面加载信息在线页面样式调整前端代码调试前端日志分析浏览器存储查看在线性能分析等…1.如何调出开发者工具按F12键CTRL+SHIFT+I 或在页面上右键,点击“检查”按钮,如下图:2.面板功能简介元素(Elements):用于查看或修改H..._谷歌后台调测执行sql

C语言文件指针偏移的使用(点阵字库txt文件取字)_c语言指针偏移-程序员宅基地

文章浏览阅读3.5k次。从点阵字库txt取字,指向文件内容位置的文件指针是唯一的,无法多个副本分别定位(无法一个文件使用多个指针),这与传统上内存中的指针不同。_c语言指针偏移

Oracle Partition 分区详细总结_partition by range-程序员宅基地

文章浏览阅读651次。此文从以下几个方面来整理关于分区表的概念及操作: 1.表空间及分区表的概念 2.表分区的具体作用 3.表分区的优缺点 4.表分区的几种类型及操作方法 5.对表分区的维护性操作.(1.) 表空间及分区表的概念表空间: 是一个或多个数据文件的集合,所有的数据对象都存放在指定的表_partition by range

随便推点

transform,translation和animation-程序员宅基地

文章浏览阅读587次。css3在原来的基础上增加了变形和动画相关的属性,动画三兄弟:transform、transition和animation,通过使用这三个属性可以达到很炫酷的效果。需要注意的是这三个属性都是css3新增的属性,各大浏览器支持方面还不是特别好,使用时要特别注意浏览器的兼容性。Transform浏览器支持情况:Internet Explorer 10、Firefox、Opera_translation和animation

[redis] SpringBoot整合SpringDataRedis配置文件_spring boot 192.168.110.130-程序员宅基地

文章浏览阅读1.7k次。SpringBoot整合SpringDataRedis pom文件添加依赖 全局配置文件 配置类pom文件添加依赖<dependency> <groupId>org.springframework.boot</groupId> <artifactId>spring-boot-starter-data-redis<..._spring boot 192.168.110.130

JAVA安卓4.4.4,Could not find support-v4.jar (com.android.support:support-v4:24.0.0)-程序员宅基地

文章浏览阅读318次。Error:A problem occurred configuring project ':app'.Could not find support-v4.jar (com.android.support:support-v4:24.0.0). Searched in the following locations: https://jcenter.bintray.com/com/android/..._supportv4.jar

OpenLayers Map理解-程序员宅基地

文章浏览阅读165次。OpenLayers Map理解 1,视口坐标的原点在左上角,水平向右为x轴正向,垂直向下为y 轴正向;2,地图坐标原点为初始图层的中心点,水平向右为x轴正向,垂直向上为y轴正向;3,视口中心点永远与地图中心点重合,不一定与瓦片中心点重合;4,拖动图层的逻辑描述:地..._openlayer this.map.getsize()

vue运用 vue-qr 生成二维码-程序员宅基地

文章浏览阅读5.4k次,点赞4次,收藏18次。第一步,先安装 vue-qr 插件npm install vue-qr --save第二步,在想要生成vueQr 的Vue页面引入组件import vueQr from 'vue-qr'第三步,在components中引入VueQr组件components: { VueQr }最后就在在html上引用<VueQr :margin='8' :size='280' :whiteMargin="true" :logoMargin="3" :logoCornerRadius_vue-qr

各品牌电脑PE中找不到硬盘的解决方法/bios设置_联想ideapad710s进pe找不到硬盘-程序员宅基地

文章浏览阅读1.3w次。  大部分小伙伴遇到系统坏了之后,都会自己用U盘装系统,但是在u盘装系统过程中却容易遇到问题,特别是PE中找不到硬盘的情况,大家遇到这种情况该怎么解决呢,所以今天跟着快启动小编的脚步一起来详细了解一二吧。  首先,我们要了解硬盘的两种模式,一种是比较新的SATA模式,现在的新电脑大多采用此模式,而另一种则是比较旧的PATA模式,而这两种模式都需要配置驱动程序,但它们的驱动程序是不一样的。  然而,..._联想ideapad710s进pe找不到硬盘

推荐文章

热门文章

相关标签