Mathlab编程-微积分在Matlab中的解法

这一章节将介绍一系列典型的微积分问题(求极限、级数、定积分、导数、重积分等)在Matlab中的求解。

首先关于极限：

(1) 数列极限：

给出下面三段例程。

求解数列极限的limit函数参数说明：可以看到该函数可以有三个参数也可以有两个参数，对于三个变量(比如说第二个例程)，第一个参数是数列的通项，第二个参数是确认离散变量，因为在通项中有两个字母a、n，第三个参数表示这个离散变量趋于某个范围，这个函数将返回在离散变量趋于的那个范围(第三个参数)时的极限。

当然，借助Matlab自身强大的图像功能，我们对于求解数列极限，我们还有一个可行的方法就是得到这个数列的散点图，通过观察这个数列趋近于一个怎样的数，但是散点的数量要始终，而且得到的结果很可能不准确。

那第一个例子来说：

>> n=1:50;

>> plot(n,n.^(1./n),'k.') % k.这里表示制自变量整数位置的散点图

得到如下的大致图像：

导数和微分在Matlab中的求解：

对于最简单的一阶导数，基于之前我们介绍的极限在Matlab中的写法，可根据定义直接求出导函数。

定义法求导函数：

diff函数直接求导函数：

基于定义法和左右导数的概念所编写的DerivativeDefinition函数：

function df=DerivativeDefinition(fun,x,x0,type)
%DERIVATIVEDEFINITION   根据导数的定义求函数的导函数或在某点处导数值
% DF=DERIVATIVEDEFINITION(FUN,X)或
% DF=DERIVATIVEDEFINITION(FUN,X,[])  求函数FUN关于X的导函数
% DF=DERIVATIVEDEFINITION(FUN,X,X0)  求函数FUN在点X0处的导函数
% DF=DERIVATIVEDEFINITION(FUN,X,X0,TYPE)  根据TYPE指定导数类型求函数在点X0处的导数，
%                                                 TYPE有以下取值：
%                                                 1.'double'或0：双侧导数值，此为缺省值
%                                                 2.'left'或-1：左导数
%                                                 3.'right'或1：右导数
% DF=DERIVATIVEDEFINITION(FUN,X,[],TYPE)  根据TYPE指定导数类型求函数的导函数
%
% 输入参数：
%     ---FUN：符号函数表达式
%     ---X：符号自变量
%     ---X0：求导点
%     ---TYPE：导数类型
% 输出参数：
%     ---DF：返回的导函数或导数值
%
% See also limit, diffif nargin<4type=0;
end
if nargin==2 || isempty(x0)x0=x;
end
syms h
delta_y=subs(fun,x,x0+h)-subs(fun,x,x0);
switch typecase {0,'double'}df=limit(delta_y/h,h,0);  % 求导数case {-1,'left'}df=limit(delta_y/h,h,0,'left');  % 求左导数case {1,'right'}df=limit(delta_y/h,h,0,'right');  % 求右导数otherwiseerror('The Style of Derivative is Illegal.')
end

结合一个实际问题对这段代码进行应用。

根据结果可知该点导数并不存在。

这里值得提的一点是，为了实现Matlab自定义函数功能，需要在Matlab编辑器中定义函数及其功能，然后在命令行窗口进行编程计算时直接调用即可。

转载于:https://www.cnblogs.com/rhythmic/p/5814391.html

Mathlab编程-微积分在Matlab中的解法相关推荐

matlab用双精度编程,如何在Matlab中得到一个单、双精度数的整数部分和小数部分？...
对于一个单.双精度数,获得其单独的整数部分或小数部分在一些程序中有特殊的作用.如对一个数进行去模运算.对一幅图像进行剪切使其高.宽是某个数的整数倍等等. 对于取整,在C语言中有专门的取整符号[],在M ...
matlab编程求平均,matlab中的分组平均函数grpstats的用法
群组函数grpstats 前面讨论到之平均值求法,通常应用于整个数组之值,若要应用到比较复杂的分组平均问题,则必须使用不同的函数才能达成.此项指令之格式如下: means = grpstats(X, ...
matlab 符号运算数值,MATLAB中的微积分运算（数值符号）
显然这个函数是单词differential(微分)的简写,用于计算微分.实际上准确来说计算的是差商. 如果输入一个长度为n的一维向量,则该函数将会返回长度为n-1的向量,向量的值是原向量相邻元素的差, ...
matlab 死区,matlab中怎么产生死区
主要是针对非线性系统死区环节进行SIMULINK仿真,通过仿真分析来实现控制系统的优化.本文介绍了MATLAB中SIULINK工具的应用,对控制系统非线性环节的死区环节的分析,并用...... (ex ...
matlab累加数组的前n行,有一组m行n列的数组，在matlab中如何编程，每三个数求一次平均数（行与列分别求）？...
有一组m行n列的数组,在matlab中如何编程,每三个数求一次平均数(行与列分别求)? 关注:216 答案:3 手机版解决时间 2021-03-05 21:45 提问者繁华落尽 2021-03- ...
LabVIEW与MATLAB混合编程——调用Matlab中.m的函数
LabVIEW与MATLAB混合编程--调用Matlab中.m的函数 1.调用labview中的matlab script 通过Labview中的Matlab script可以调用Matlab的.m程 ...
labview中java联合编程_LabVIEW与MATLAB混合编程的实现
本栏目责任编辑:谢媛媛系统软件与软件工程 ComputerKnowledgeand Technology 电脑知识与技术第 6 卷第 15 期 (2010 年 5 月) LabVIEW 与 MATL ...
在matlab中编程n,Matlab与C++混合编程 1--在C++中调用自己写的matlab函数
在Visual Studio中使用C++中调用MATLAB程序在matlab中可以通过mbuild工具将.m文件编译成dll文件供外部的C++程序调用,这样就可以实现matlab和C++混合编程的目 ...
matlab中p是什意思,如下的用matlab进行编程的勒让德函数，求解释，看不懂什么意思。如:p和t分别代表什么?p(1,:)是什么意思...
答:找本数学物理方法的书,弄好两件事: 1.勒让德函数与缔合勒让德函数,MATLAB的legendre函数是后者. 2.顺便找一下勒让德函数的递推公式. 我刚才粗看了一下,可能P(2,:),t(2,: ...

Mathlab编程-微积分在Matlab中的解法

Mathlab编程-微积分在Matlab中的解法相关推荐

最新文章

热门文章