利用MATLAB进行求导
导数是微积分中的一个重要的概念,又称为导函数值。在数学中,我们经常会利用导数来求切线、极值、拐点、渐近线等。在解决数学问题中,经常需要用到导数。MATLAB中提供了diff函数进行求导。本文对于MATLAB的求导操作进行简单介绍。
1、求导的定义
求导是数学计算中的一个计算方法,它的定义就是,当自变量的增量趋于零时,因变量的增量与自变量的增量之商的极限。在一个函数存在导数时,称这个函数可导或者可微分。可导的函数一定连续。不连续的函数一定不可导。本段定义来自百度百科(求导的定义)
2、diff函数
在MATLAB中提供了diff函数用于导数的求导,其中diff函数提供了四种调用方式。
(1)diff(y):这种调用方式没有指明需要求导的自变量以及导数的阶层数目,一般情况下,是对默认的变量进行一次的求导。
例如,在数学学习中,我们会知道一些比较常用的求导的公式有:
下面利用MATLAB对其进行验证:
syms x u a
y1=diff(x^u,x)
y2=diff(exp(x))
y3=diff(log(x))
y4=diff(sin(x))
y5=diff(cos(x))
y6=diff(a^x,x)
y7=diff(log(x)/log(a),x)
(注:因为在MALTAB中仅有以自然数e、2和10为底的对数公式,这里转换为log(x)/log(a)的形式)。
运行结果如下所示:
y1 = u*x^(u - 1)
y2 =exp(x)
y3 =1/x
y4 =cos(x)
y5 =-sin(x)
y6 =a^x*log(a)
y7 =1/(x*log(a))
通过使用MATLAB代码对于上面的公式进行验证,可以看出MATLAB所计算的结果与公式相同。
在学习导数的时候,我们会学习到导数的乘法与除法的运算法则。
其中导数的乘法运算法则如下所示:
导数的除法的运算法则如下所示:
下面利用MATLAB代码对其进行验证,例如下面这个式子对其进行求导操作:
利用MATLAB代码对其进行验证得:
syms x
y=x*exp(x);
diff(y)
运行结果如下所示:
ans =exp(x) + x*exp(x)
可以看出结果符合乘法的运算法则。
下面对于除法的运算法则进行验证:
syms x
y=exp(x)/x;
diff(y)
运行结果如下所示:
ans =exp(x)/x - exp(x)/x^2
通过MATLAB所输出的结果看出符合除法的运算法则。
下面我们可以对于tan(x)的导数进行验证,tanx转换为sinx/cosx的形式,利用除法的运算法则可以对其进行计算:
syms x
y=sin(x)/cos(x);
diff(y)
运行结果如下所示:
ans =sin(x)^2/cos(x)^2 + 1
对于MATLAB所计算的结果进行验证:
而对tanx进行求导的结果如下所示:
通过对比可以看出MATLAB所求结果正确。
(2)当符号表达式中涉及到了多变量时,diff函数可以指定对于特定符号变量进行求导,调用格式为diff(y,'x')。其中diff(y,'x'):假设需要对其进行求导的变量是x,因为没有说明求导的阶数,这种情况下,MATLAB默认对其进行一阶求导。
例如对于下面的公式进行求导:
对下面的式子进行求导:
syms x a b c
y=a*x^3+b*x^2+x+c;
diff(y,'x')
运行结果如下所示:
ans =3*a*x^2 + 2*b*x + 1
通过上面的例子可以看到,当存在多个符号变量的时候,diff函数可以指定特定的自变量对其进行求导。
(3)另外一种调用方式是指定特定的阶数,未说明对其求导的自变量,调用方式是diff(y,n),其中对于MATLAB指定的自变量求其n阶导数。
例如同样对下面的公式进行求导:
MATLAB代码如下所示:
syms x a b c
y=x^3+2*x^2+4*x+6;
diff(y,2)
运行结果如下所示:
ans =6*x + 4
(4)最后一种调用方式是对其进行求导的自变量、求导的阶数都进行说明,该方式的调用格式为diff(y,‘x’,n):该调用方式是以对于自变量x求其n阶导数。
例如,对于下面的式子的x变量求其二阶导数:
MATLAB代码如下所示:
syms x
y=(x+1)*log(x)-x^2+x+1;
diff(y,'x',2)
运行结果如下所示:
ans =2/x - (x + 1)/x^2 - 2
3、利用MATLAB对于复杂的公式求导
通过使用MATLAB对于一些复杂的公式进行求导,方便我们操作。下面对其较为困难的公式进行求导,下面举几个例子:
MATLAB代码如下所示:
syms x
y1=diff(sqrt(x)*sin(x),'x',1)
y2=diff(1/((3*x-1)^2),'x',1)
y3=diff(cos(sin(x)),'x',1)
运行结果如下所示:
y1 = x^(1/2)*cos(x) + sin(x)/(2*x^(1/2))
y2 =-6/(3*x - 1)^3
y3 =-sin(sin(x))*cos(x)
通过运行的结果可以看出,在MATLAB中对于较复杂的符号表达式都可以进行求导操作。
利用MATLAB进行求导相关推荐
- 如何利用MATLAB来求二元函数的一阶偏微分
如何利用MATLAB来求二元函数的一阶偏微分 在上一篇文章中我写到了如何用差分法来求像y=f(x)这样的一元函数的一阶二阶导数,那如果是像f=f(x,y)这样的二元函数,怎么利用差分法求它们的偏微分呢 ...
- matlab 曲线拟合求导,基于matlab曲线拟合的数据预测分析
龙源期刊网 http://www.doczj.com/doc/2f22f716cbaedd3383c4bb4cf7ec4afe05a1b148.html 基于matlab曲线拟合的数据预测分析 作者: ...
- Matlab之求导篇
matlab求导命令diff调用格式: diff(函数) , 求的一阶导数: diff(函数, n) , 求的n阶导数(n是具体整数): diff(函数,变量名), 求对的偏导数: diff(函数, ...
- Matlab中求导操作
文章目录 前言 一.求导命令 二.使用方法 总结 前言 正文 一.求导命令 diff % 差分和近似导数%此 MATLAB 函数 计算沿大小不等于 1 的第一个数组维度的 X 相邻元素之间的差分: 二 ...
- [MATLAB]多项式求导/加减/乘除
函数公式: 多项式的加减运算 多项式的加减运算非常简单,即相应向量相加减 多项式乘法 conv(p1,p2):多项式相乘函数.在这里,P1.P2是两个多项式系数向量 多项式除法 [Q,r]=decon ...
- Matlab符号计算求导与化简
依托于Maple的符号计算引擎,Matlab可以完成高等数学中几乎一切符号计算.这里简单记录一下蔡自兴<机器人学>第三版中第四章的双连杆的动能推导过程. 在书本上我们已知的信息是: x2= ...
- 利用matlab软件求数列的中值,matlab求两个数列的方差
matlab怎么求矩阵所有元素的期望和方差? 标准差s=std(X(1:end),flag)flag=0,采用1/(N-1)的系数,flag=1,采用1/(N)的系数 matlab用var函数算出的方 ...
- 【matlab】求导
一:matlab求一阶导 diff(y,x) 其中y是x的函数 二:matlab求二阶导 diff(y,x,2) 三:举例 syms x y y=x^3; dy=diff(y,x) ddy=diff( ...
- matlab 离散点求导_Matlab的离散点曲线导数曲率数值模拟方法
模拟方法: 复制内容到剪贴板 代码 : clc; clear all; close all; x0 = linspace(0, 1); y0 = sin(x0).*cos(x0); h = abs(d ...
最新文章
- ssm 异常捕获 统一处理_SSM 统一异常处理
- java object finalize_Java Object finalize() 方法
- 文件下载类型__response
- MySQL设置数据库及表的字符集
- boost::fusion::map用法的测试程序
- 窗口分析函数_1_生成不重复排名序号
- cannot import caffe
- mysql+io+参数_MySQL IO线程及相关参数调优
- VMware内虚拟机自适应及最大化窗口调整方式
- Oracle 隐含参数的查询
- 机器学习入门(二):工具与框架的选择
- linux之VMware安装Centos7
- 摸鱼摸累了,看几个有趣的网站
- php做后端的优势,【后端开发】php和java的优势
- php支付宝聚合支付源码,ThinkPHP5.1聚合支付源码
- linux 可视化分区,可视化linux块设备的工具(分区,LVM PV,LV,mdadm设备……)
- IE系列浏览器把application/json响应视为文件,并尝试下载
- TeamViewer打开后没有密码的解决办法
- 编程实现输入两个正整数m和n,求其最大公约数和最小公倍数
- 项目管理实践【五】自动编译和发布网站【Using Visual Studio with Source Control System to build and publish website autom
热门文章
- 解决死锁问题的基本方法
- MMX,SSE,SSE2扫盲
- 基于哈工大生物信息学专业培养方案提炼的学程
- 蒙特卡洛 股票 matlab,风险管理matlab蒙特卡洛模拟股票价格
- 怎么把备忘录里的照片存储到相册
- 算法与数据结构(三) 时间复杂度分析 [例题]
- 【渝粤教育】电大中专办公设备使用与维护作业 题库
- stm32工程模板的创建
- Collection coll = new ArrayList() 为什么可以这么写!!!
- JavaScript -- 对象与原型链 [[Prototype]]