最新微积分的基本运算

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

微积分的基本运算

第4章微积分的基本运算

本章学习的主要目的:

1.复习高等数学中有关函数极限、导数、不定积分、定积分、二重积分、级数、方程近似求解、常微分方程求解的相关知识.

2.通过作图和计算加深对数学概念:极限、导数、积分的理解.

3.学会用MatLab软件进行有关函数极限、导数、不定积分、级数、常微分方程求解的符号运算;

4.了解数值积分理论,学会用MatLab软件进行数值积分;会用级数进行近似计算.

1 有关函数极限计算的MatLab命令

(1)limit(F,x,a) 执行后返回函数F在符号变量x趋于a的极限

(2)limit(F,a) 执行后返回函数F在符号变量findsym(F)趋于a的极限

(3)limit(F) 执行后返回函数F在符号变量findsym(F)趋于0的极限

仅供学习与交流,如有侵权请联系网站删除谢谢87

(4)limit(F,x,a,’left’) 执行后返回函数F在符号变量x趋于a的左极限

(5)limit(F,x,a,’right’) 执行后返回函数F在符号变量x趋于a的右极限

注:使用命令limit前,要用syms做相应符号变量说明.

例7 求下列极限

(1)«Skip Record If...»

在MatLab的命令窗口输入:

syms x

limit((cos(x)-exp(-x^2/2))/x^4,x,0)

运行结果为

ans =-1/12

理论上用洛必达法则或泰勒公式计算该极限:

方法1 «Skip Record If...»

«Skip Record If...»

方法2 «Skip Record If...»

«Skip Record If...»

(2) «Skip Record If...» %自变量趋于无穷大,带参数t

在MatLab的命令窗口输入:

仅供学习与交流,如有侵权请联系网站删除谢谢87

syms x t

limit((1+2*t/x)^(3*x),x,inf)

运行结果为

ans =exp(6*t)

理论上用重要极限计算:

«Skip Record If...»

(3) «Skip Record If...» %求右极限

在MatLab的命令窗口输入:

syms x

limit(1/x,x,0,’right’)

运行结果为

ans = inf

2 有关函数导数计算的MatLab命令

(1)diff(F,x) 表示表达式F对符号变量x求一阶导数,允许

表达式F含有其他符号变量,若x缺省,

则表示对由命令syms定义的变量求一阶

导数。

(2)diff(F,x,n) 表示表达式F对符号变量x求n阶导数。

例10求下列函数的导数

仅供学习与交流,如有侵权请联系网站删除谢谢87

(1)已知«Skip Record If...»,求«Skip Record If...»;

在MatLab的命令窗口输入如下命令序列:

syms x

y=x*asin(x/2)+sqrt(4-x^2)

diff(y,x) %执行结果ans = asin(1/2*x)与理论推导«Skip Record If...»完全吻合。

diff(y,x,3) %执行结果ans = 1/(4-x^2)^(3/2)*x与理论推导

«Skip Record If...»完全吻合。

(2)已知«Skip Record If...»,求«Skip Record If...»

在MatLab的命令窗口输入如下命令序列:

syms x y z

z=x^2*sin(2*y);

diff(z,x) %执行结果ans =2*x*sin(2*y)

diff(z,x,2) %执行结果ans =2*sin(2*y)

diff(diff(z,x),y) %执行结果ans =4*x*cos(2*y)

(3)已知«Skip Record If...»,求«Skip Record If...» (复合函数求导偏导数)

在MatLab的命令窗口输入如下命令序列:

syms x y z u

仅供学习与交流,如有侵权请联系网站删除谢谢87

z=x^2+y^2;

u=(x-y)^z;

diff(u,x) %执行结果«Skip Record If...»(x-y)^(x^2+y^2)*(2*x*log(x-y)+(x^2+y^2)/(x-y))

diff(u,y,2) %执行结果«Skip Record If...»(x-y)^(x^2+y^2)*(2*y*log(x-y)-(x^2+y^2)/(x-y))^2+

(x-y)^(x^2+y^2)*(2*log(x-y)-4*y/(x-y)-

(x^2+y^2)/(x-y)^2)

diff(diff(u,x),y) %执行结果«Skip Record If...»(x-y)^(x^2+y^2)*(2*y*log(x-y)-

(x^2+y^2)/(x-y))*(2*x*log(x-y)+(x^2+y^2)/(x-y))+(x-

y)^(x^2+y^2)*(-2*x/(x-y)

+2*y/(x-y)+(x^2+y^2)/(x-y)^2)

3 极值问题

MatLab软件提供了求一元和多元函数极值问题的命令:fmin(f,x1,x2) 求函数f(x)在x1

fmins('f',[x1,x2]),求二元函数在点(x1 x2)附近的极值

点。

仅供学习与交流,如有侵权请联系网站删除谢谢87

相关文档
最新文档