matlab实验求不定积分函数,实验四用matlab计算积分
实验四用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. (广义积分)计算广义积分 ,,
matlab实验求不定积分函数,实验四用matlab计算积分相关推荐
- matlab中求分段函数的分段点,matlab求解分段函数问题是如何用下面的算法求解下面的分段函数 爱问知识人...
问题是如何用下面的算法求解下面的分段函数,目前出现的问题是:??? Error using ==> at 11Function 'gt' is not implemented for MuPAD ...
- 计算机函数的实验报告,excel函数实验报告.doc
excel函数实验报告 excel函数实验报告 篇一:Excel常用函数与公式应用实验报告 学生实验报告 实验课名称:计算机应用基础 实验项目名称:Excel常用函数与公式应用 专业名称: 班 学级: ...
- 如何利用MATLAB来求二元函数的一阶偏微分
如何利用MATLAB来求二元函数的一阶偏微分 在上一篇文章中我写到了如何用差分法来求像y=f(x)这样的一元函数的一阶二阶导数,那如果是像f=f(x,y)这样的二元函数,怎么利用差分法求它们的偏微分呢 ...
- 【优化求解】基于matlab禁忌搜索算法求解函数极值问题【含Matlab源码 1204期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源: [优化求解]基于matlab禁忌搜索算法求解函数极值问题[含Matlab源码 1204期] 点击上面蓝色字体,直接付费下载,即可. 获取 ...
- matlab中求包络线函数,matlab包络线函数
的包络 解调方法 一直 是机 械故 障诊 线 ,只 需要经 过大 脑 的取舍 ,...(c )可 以看 出,经过 MATLAB 中的 Hilbert 函数 后也 出现 了皱 纹 ,...... yma ...
- matlab怎么求ln函数,用matlab求函数e的1+ln[erf(1/t)]+t次方对t的积分求详细编程过程...
答:syms t; f = exp(1+log(1/t)+t); int(f,t) 希望对你有帮助! 答:没有精确解,试试数值解吧 答:不可以,这个没有原函数的. 答:1.双击matlab软件图标,打 ...
- matlab求二元函数极值算法_最优化计算与matlab实现(3)——进退法
参考资料 <精通MATLAB最优化计算(第二版)> 数值实现 Matlab 2019a 目录 石中居士:最优化计算与Matlab实现--目录zhuanlan.zhihu.com 进退法 ...
- matlab 导数表达式,matlab语言怎样表示x的导数,matlab 如何求自定义函数的导数
matlab的求导命令与求导法 建立符号变量命令sym和syms调用格式: x=sym('x'), 建立符号变量x: syms x y z , 建立多个符号变量x,y,z: matlab求导命令dif ...
- matlab如何输出一个函数的图像,在matlab中如何输出一个函数的图像
点击查看在matlab中如何输出一个函数的图像具体信息 答:A. imread imread函数用于读入各种图像文件,其一般的用法为 [X,MAP]=imread('filename','fmt')其 ...
最新文章
- 决策树的C++实现(CART)
- 怎么检测JDK环境变量是否配置正确
- 关于oracle分组后组外排序的问题
- Spring MVC请求处理流程分析
- php微信浏览器自带后退功能禁用,jquery如何实现禁用浏览器后退
- Linux chapter test 10
- 浅析:通过自定义DSL实现一个序列号生成器
- 软件测试工程师调研报告,软件测试工程师竞聘报告范文.docx
- 软件测试主要种类大全
- 模电基础讲解02:晶体二极管的结构类型
- react-custom-scrollbars滚动组件
- elasticsearch两年学习资料整理分享
- c语言元音字母辅音字母,英语字母中的5个元音字母和21个辅音字母
- 爬虫实战——爬取电影天堂的电影详情页信息
- 【Java实现】南京地铁导航系统的简单实现(一)—— 存储站点信息
- SpringMVC视图及如何在域对象中共享数据
- 【C生万物】 函数递归篇
- Mac下浏览器安装证书
- QuickBI 使用记录
- 【0109】Linux系统监测工具sysstat介绍
热门文章
- 制表符(tab 键)的概念
- 为什么能在子线程通过setText进行更新UI
- 石家庄地区十大有影响力的调查研究咨询公司
- (转)如何有效使用Project-编制进度计划、保存基准
- Latex 之table
- html.action 访问分部视图,MVC+EF 随笔小计——分部视图(Partial View)及Html.Partial和Html.Action差异...
- android 应用之间相互跳转,拉起 app跳转app
- 「设计模式」六大原则之三:里氏替换原则小结
- chrome 64 bit 用于最新版本下载
- htmltd文本自动换行,完美解决table中td里面的内容自动换行