当前位置:首页 » 《随便一记》 » 正文

利用MATLAB进行求导

7 人参与  2023年04月09日 13:29  分类 : 《随便一记》  评论

点击全文阅读


导数是微积分中的一个重要的概念,又称为导函数值。在数学中,我们经常会利用导数来求切线、极值、拐点、渐近线等。在解决数学问题中,经常需要用到导数。MATLAB中提供了diff函数进行求导。本文对于MATLAB的求导操作进行简单介绍。

1、求导的定义

求导是数学计算中的一个计算方法,它的定义就是,当自变量的增量趋于零时,因变量的增量与自变量的增量之商的极限。在一个函数存在导数时,称这个函数可导或者可微分。可导的函数一定连续。不连续的函数一定不可导。本段定义来自百度百科(求导的定义)

2、diff函数

在MATLAB中提供了diff函数用于导数的求导,其中diff函数提供了四种调用方式。

(1)diff(y):这种调用方式没有指明需要求导的自变量以及导数的阶层数目,一般情况下,是对默认的变量进行一次的求导。

例如,在数学学习中,我们会知道一些比较常用的求导的公式有:

\small (x^u)'=ux^{x-1}

\small (e^x)'=e^x

\small (lnx)'=\frac{1}{x}

\small (sinx)'=cosx

\small (cosx)'=-sinx

(a^x)'=a^xlna

(log_{a}x)'=\frac{1}{xlna}

 下面利用MATLAB对其进行验证:

syms x u ay1=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/xy4 =    cos(x)y5 =    -sin(x)y6 =    a^x*log(a)y7 =    1/(x*log(a)) 

通过使用MATLAB代码对于上面的公式进行验证,可以看出MATLAB所计算的结果与公式相同。

在学习导数的时候,我们会学习到导数的乘法与除法的运算法则。

其中导数的乘法运算法则如下所示:

[f(x)g(x)]'=f'(x)g(x)+f(x)g'(x)

 导数的除法的运算法则如下所示:

[\frac{f(x)}{g(x)}]=\frac{f'(x)g(x)-f(x)g'(x)}{g^2(x)}

 下面利用MATLAB代码对其进行验证,例如下面这个式子对其进行求导操作:

y=xe^x

利用MATLAB代码对其进行验证得:

syms xy=x*exp(x);diff(y)

运行结果如下所示:

ans =    exp(x) + x*exp(x)

可以看出结果符合乘法的运算法则。

下面对于除法的运算法则进行验证:

syms xy=exp(x)/x;diff(y)

运行结果如下所示:

ans =    exp(x)/x - exp(x)/x^2

通过MATLAB所输出的结果看出符合除法的运算法则。

下面我们可以对于tan(x)的导数进行验证,tanx转换为sinx/cosx的形式,利用除法的运算法则可以对其进行计算:

syms xy=sin(x)/cos(x);diff(y)

运行结果如下所示:

ans =    sin(x)^2/cos(x)^2 + 1 

对于MATLAB所计算的结果进行验证:

\frac{sin^2x}{cos^2x}+1=\frac{sin^2x+cos^2x}{cos^2x}=\frac{1}{cos^2x}

 而对tanx进行求导的结果如下所示:

 通过对比可以看出MATLAB所求结果正确。

(2)当符号表达式中涉及到了多变量时,diff函数可以指定对于特定符号变量进行求导,调用格式为diff(y,'x')。其中diff(y,'x'):假设需要对其进行求导的变量是x,因为没有说明求导的阶数,这种情况下,MATLAB默认对其进行一阶求导。

例如对于下面的公式进行求导:

y=ax^3+bx^2+x+c

 对下面的式子进行求导:

syms x a b cy=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阶导数。

例如同样对下面的公式进行求导:

y=x^3+2x^2+4x+6

 MATLAB代码如下所示:

syms x a b cy=x^3+2*x^2+4*x+6;diff(y,2)

运行结果如下所示:

ans =    6*x + 4

(4)最后一种调用方式是对其进行求导的自变量、求导的阶数都进行说明,该方式的调用格式为diff(y,‘x’,n):该调用方式是以对于自变量x求其n阶导数。

例如,对于下面的式子的x变量求其二阶导数:

y=(x+1)lnx-x^2+x+1

 MATLAB代码如下所示:

syms xy=(x+1)*log(x)-x^2+x+1;diff(y,'x',2)

运行结果如下所示:

ans =    2/x - (x + 1)/x^2 - 2

3、利用MATLAB对于复杂的公式求导

通过使用MATLAB对于一些复杂的公式进行求导,方便我们操作。下面对其较为困难的公式进行求导,下面举几个例子:

y=\sqrt{x}sinx

y=\frac{1}{(3x-1)^2}

y=cos(sinx)

 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)^3y3 =    -sin(sin(x))*cos(x)

通过运行的结果可以看出,在MATLAB中对于较复杂的符号表达式都可以进行求导操作。


点击全文阅读


本文链接:http://zhangshiyu.com/post/58997.html

<< 上一篇 下一篇 >>

  • 评论(0)
  • 赞助本站

◎欢迎参与讨论,请在这里发表您的看法、交流您的观点。

最新文章

  • (此去经年无故人)南初陆南城:结局+番外精品选集起点章节+阅读即将发布预订
  • 沈凝夏叶晚怡附加完整在线阅读(归雁不栖故人枝)最近更新列表
  • 剧情人物是时初,白浩雄的玄幻言情小说《召诸神,踏万界,天命帝女逆乾坤》,由网络作家&ldquo;海鸥&rdquo;所著,情节扣人心弦,本站TXT全本,欢迎阅读!本书共计381345字,185章节,:结局+番外免费品鉴:结局+番外评价五颗星
  • 凤青禾,江明远,***枢小说(别人修仙我捡漏,卷王们破防了)最近更新(凤青禾,江明远,***枢)整本无套路阅读
  • 薛梨小说无删减+后续(曾经亲情似草芥)畅享阅读
  • 沈南栀小说(穿越时空,我要修补时空裂缝)章节目录+起点章节(沈南栀)全篇清爽版在线
  • 未婚妻被巨蟒缠身,我该吃就吃该喝就喝前言+后续_阿豪林月周然后续+番外_小说后续在线阅读_无删减免费完结_
  • 陆骁,陆本初小说(陆骁,陆本初)(癫!睁眼穿成老太太挥鞭***逆子)前传+阅读全新作品预订
  • 姐姐含冤而死后冥王另娶,我杀穿整个地府在线阅读_阎罗殿殷红别提一口气完结_小说后续在线阅读_无删减免费完结_
  • (书荒必看)毒后重生:疯王的神医小娇妻沈清歌,萧绝:+后续热血十足
  • 重生后我和太监联手灭了敌国喻辰,林雪续集(重生后我和太监联手灭了敌国)终极反转(喻辰,林雪)全篇一口气阅读
  • 我不做灵媒后,自称灵媒摆渡人的养妹害怕了内容精选_苏晓霍老阿姐无广告_小说后续在线阅读_无删减免费完结_

    关于我们 | 我要投稿 | 免责申明

    Copyright © 2020-2022 ZhangShiYu.com Rights Reserved.豫ICP备2022013469号-1