matlab能不能求不定积分,用MATLAB求定积分
一、符号积分
符号积分由函数int来实现。该函数的一般调用格式为:
int(s):没有指定积分变量和积分阶数时,系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分;
int(s,v):以v为自变量,对被积函数或符号表达式s求不定积分;
int(s,v,a,b):求定积分运算。a,b分别表示定积分的下限和上限。该函数求被积函数在区间[a,b]上的定积分。a和b可以是两个具体的数,也可以是一个符号表达式,还可以是无穷(inf)。当函数f关于变量x在闭区间[a,b]上可积时,函数返回一个定积分结果。当a,b中有一个是inf时,函数返回一个广义积分。当a,b中有一个符号表达式时,函数返回一个符号函数。
例:
求函数x^2+y^2+z^2的三重积分。内积分上下限都是函数,对z积分下限是sqrt(x*y),积分上限是x^2*y;对y积分下限是sqrt(x),积分上限是x^2;对x的积分下限1,上限是2,求解如下:
>>syms x y z %定义符号变量
>>F2=int(int(int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),x,1,2) %注意定积分的书写格式
F2 =
1610027357/6563700-6072064/348075*2^(1/2)+14912/4641*2^(1/4)+64/225*2^(3/4) %给出有理数解
>>VF2=vpa(F2) %给出默认精度的数值解
VF2 =
224.92153573331143159790710032805
二、数值积分
1.数值积分基本原理
求解定积分的数值方法多种多样,如简单的梯形法、辛普生(Simpson)法、牛顿-柯特斯(Newton-Cotes)法等都是经常采用的方法。它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1],i=1,2,…,n,其中x1=a,xn+1=b。这样求定积分问题就分解为求和问题。
2.数值积分的实现方法
基于变步长辛普生法,MATLAB给出了quad函数来求定积分。该函数的调用格式为:
[I,n]=quad('fname',a,b,tol,trace)
基于变步长、牛顿-柯特斯(Newton-Cotes)法,MATLAB给出了quadl函数来求定积分。该函数的调用格式为:
[I,n]=quadl('fname',a,b,tol,trace)
其中fname是被积函数名。a和b分别是定积分的下限和上限。tol用来控制积分精度,缺省时取tol=0.001。trace控制是否展现积分过程,若取非0则展现积分过程,取0则不展现,缺省时取trace=0。返回参数I即定积分值,n为被积函数的调用次数。
例:
求函数'exp(-x*x)的定积分,积分下限为0,积分上限为1。
>>fun=inline('exp(-x.*x)','x'); %用内联函数定义被积函数fname
>>Isim=quad(fun,0,1) %辛普生法
Isim =
0.746824180726425
IL=quadl(fun,0,1) %牛顿-柯特斯法
IL =
0.746824133988447
三、梯形法求向量积分
trapz(x,y)—梯形法沿列方向求函数Y关于自变量X的积分(向量形式,数值方法)。
>>d=0.001;
>>x=0:d:1;
>>S=d*trapz(exp(-x.^2))
S=
0.7468
或:
>>format long g
>>x=0:0.001:1; %x向量,也可以是不等间距
>>y=exp(-x.^2); %y向量,也可以不是由已知函数生成的向量
>>S=trapz(x,y); %求向量积分
S =
0.746824071499185
附:int与quad区别
int的积分可以是定积分,也可以是不定积分(即有没有积分上下限都可以积)可以得到解析的解,比如你对x^2积分,得到的结果是1/3*x^3,这是通过解析的方法来解的。如果int(x^2,x,1,2)得到的结果是7/3
quad是数值积分,它只能是定积分(就是有积分上下限的积分),它是通过simpson数值积分来求得的(并不是通过解析的方法得到解析解,再将上下限代入,而是用小梯形的面积求和得到的)。如果f=inline('x.^2');quad(f,1,2)得到的结果是2.333333,这个数并不是7/3
int是符号解,无任何误差,唯一问题是计算速度;quad是数值解,有计算精度限制,优点是总是能有一定的速度,即总能在一定时间内给出一个一定精度的解。
matlab能不能求不定积分,用MATLAB求定积分相关推荐
- matlab与数学实验王兵团,Matlab与数学实验(第3版)
第1章 Matlab基础知识 1.1 Matlab概述 1.1.1 Matlab简介 1.1.2 Matlab的启动与退出 1.1.3 Matlab操作的注意事项 1.2 Matlab的具体操作与操作 ...
- matlab实验求不定积分函数,实验四用matlab计算积分
实验四用matlab计算积分 实验四 用matlab计算积分 4.1积分的有关理论 定积分:积分是微分的无限和,函数在区间上的积分定义为 其中 从几何意义上说,对于上非负函数,记分值是曲线与直线及轴所 ...
- matlab 求不定积分与定积分
求不定积分与定积分 %%计算不定积分 clc; clear; syms x; % 定义一个符号 x. f = 'x'; % 定义一个原函数 F = int(f, x); % 计算不定积分 disp(F ...
- 运用数学软件matlab求无穷积分,matlab积分的计算及其简单应用论文.doc
积分的计算及其简单应用 摘要:本文简要的概述了MATLAB 在高等数学中积分的计算及应用:利用MATLAB 中符号积分和数值积分的命令,计算定积分和不定积分.同时,也可以通过这些命令来解决一些实际问题 ...
- MATLAB新手简明使用教程(七)——使用matlab建立多项式以及求导,商求导乘积求导等——新手来看,保证看懂。
前期回顾 上一期中,我们学了下面的知识: 定积分的基本概念和一些简单的几何意义. 使用 int 函数计算不定积分. 使用 int 函数计算定积分. 本期内容 本期我打算给大家介绍一下使用matlab对 ...
- matlab里面求线性方程组,利用matlab求线性方程组.doc
利用matlab求线性方程组.doc MATLAB语言课成论文利用MATLAB求线性方程组姓名郭亚兰学号12010245331专业通信工程班级2010级通信工程一班指导老师汤全武学院物电学院完成日期2 ...
- matlab怎么对sinx求导,用matlab程序求y=ln(sinx 1)的导数
哪位高手会画函数f(x,y)=(1+e^y)*cosx-y*e^y的图像啊?求真相,最好是用Matlab并附上程序! clearall;clc;%画函数f(x,y)=(1+e^y)*cosx-y*e^ ...
- matlab怎么求hadamard,hadamard matlab
3.图像Hadamard变换 数字图像处理 ? Matlab没有提供图像Hadamard变换功能,不过 提供了求Hadamard变换矩阵的功能.例如使用 命令hadamard(8)能够得到下面[8 8 ...
- matlab求微分方程精确解,matlab求微分方程精确解及近似解.ppt
matlab求微分方程精确解及近似解.ppt 还剩 24页未读, 继续阅读 下载文档到电脑,马上远离加班熬夜! 亲,喜欢就下载吧,价低环保! 内容要点: 求微分方程的解q 自牛顿发明微积分以来,微分方 ...
最新文章
- 循环斐波那契数列_剑指offer #10 斐波那契数列
- 人天生就会除法,宾大最新研究:儿童啥都不学也能算出来
- linux怎样查看内核参数,Linux 实例如何查看和修改 Linux 实例内核参数?
- java 技术 管理_java技术转型产品或者项目管理,请大神指导?
- Python操作数据库之 MySQL
- 使用Minify合并css和js减少http请求
- 0112作业 字符串练习
- linux配置usb主从_一种Linux下USB设备主从切换的实现
- SpringBoot最佳实践-dev-tools
- python模块(5)-Matplotlib 简易使用教程
- 前端技术基础(一):浏览器相关
- HTML 表格的样式属性
- 传智播客 python 私有化学习笔记
- Web笔记:jQuery的使用
- AutoCAD启动自动加载程序的研究
- win7 计算机右键没有属性,Win7系统右键点击计算机没有属性选项怎么办
- linux拷贝4g以上文件夹,32位Ubuntu中支持超过4G的大内存的方法
- html dashed 属性,css虚线样式dotted和dashed
- 第一次用示波器怎么使用(基础经验)
- Vue脚手架创建步骤
热门文章
- Element ui 组件中用键盘事件
- 两台计算机组装在一起,好看又实用 手把手教你把两台电脑装在一起
- 港科报道 | 汪扬教授:数字经济是不可逆转,不受地缘政治影响的全球化经济
- 反问疑问_反问疑问句的语法
- 自己总结的wireshark抓包技巧
- MP4 全介绍【转载】
- 企业为什么要选择人力外包?好处与坏处
- 字体图标iocnfont
- mysql远程连接数据库的二种方法_mysql 远程连接数据库的二种方法 | 速光网络博客...
- 查询出每一个雇员的姓名,工资,部门名称,工资在公司的等级及其领导的姓名,领导的工资,以及领导所相应的等级...