Mathematica4.0 ppt第三讲.pptVIP

用Mathematica做微积分 内容 1、极限的运算 2、函数求导运算 3、积分的运算 4、级数的运算 5、微分方程的求解 5.1极限的运算 命令表达式为：Limit[f[x],x-a] 例 求解 Limit[((2x+3)/(2x+1))^(x+1),x-Infinity] E （注：E即为常数e） 思考：如何验证结果的正确性？ 右极限 Limit[f[x],x-a,Direction--1] 左极限 Limit[f[x],x-a,Direction-+1] 5.2 函数求导运算 D[f[x],x] 表示函数对自变量x的一阶导数 D[f[x],{x,n}] 表示函数对自变量x的n阶导数 D[f[x,y],x] ,D[f[x,y],y]分别表示函数对自变量x,y的一阶偏导数,高阶导数类似。 例 求 的一阶和二阶导数。 5.2 函数求导运算 思考：对于隐函数和参数方程确定的函数如何求导？（提示：《高等数学》中的思路） 自己根据所学知识设计步骤，在上机课上完成相应题目。 5.3 积分的运算 Integrate[f[x],x] 不定积分（结果中不含任意常数C） Integrate[f[x],{x,xmin,xmax}] 定积分 例 求 5.3 积分的运算 Integrate[f[x,y,…],{x,a,b},{y,c,d},…] 多重积分 例 求 5.4 级数的运算 Series[f[x],{x,a,n}] 给出f[x]在a处的n次幂级数展开 例 将 展开成x和(x-1)的10次幂级数 Series[Exp[x],{x,0,10}] 5.4 微分方程的求解 DSolve[eqn,y[x],x] 求解微分方程y[x] DSolve[{eqn1,eqn2,…},{y1,y2,….},x] 求解微分方程组 DSolve[{eqn,y[a]＝＝b},y[x],x]求解微分方程y[x]在初始条件y[a]=b时的解。 例 求 的解。 y[0]=2时的解。 用Mathematica作向量、矩阵运算 在Mathematic中，有序数组被称为“表”（list）“表”既可以表示集合，又可以表示向量和矩阵。许多函数都可以作用在表上。 6.1 向量和矩阵的输入 6.2 获得表的元素 6.3 表的维数和加、减法 6.4 向量和矩阵的乘法 6.5 关于矩阵的几个常用函数 6.1 向量和矩阵的输入 从键盘输入一个表，用{ }将表的元素括起，元素之间用逗号分隔。 例1 输入数据列0，16，64，144，256。定义为变量data. data={0,16,64,144,256} 例2 输入矩阵M= M={{2,5,-1}，{0,-1,3}，{1,2,-2}} 矩阵的每一行用{ }括起。 对于某些有规律的表Mathematica提供了函数Table[ ],Nestlist[ ]。 例3 已知数列通项 ，给出前10项。 例4 给出30以内的奇数。 例5 给出 特殊矩阵的输入命令有 Table[f[i,j],{i,m}.{j,n}] 生成以f的计算值为元素的m行n列矩阵 Array[a,{m,n}] 生成以a[i,j]为元素的m行n列矩阵。 IdentityMatrix[n] 生成n阶单位阵。 DiagonaMatrix[List] 生成以表中元素为对角元的对角矩阵。 例6 生成三阶Hilbert矩阵 得到 例7 生成四阶单位阵 例8 生成一个以{1，2，3，4，5}为对角元的对角矩阵，并用矩阵形式表示 得到 6.2 获得表的元素 A是一个向量，则A[i]表示向量的第i个元素。 M是一个m行n列矩阵，则用M[[i]]表示矩阵的第i行； M[[i,j]] 表示第i行,第j列交叉点处的元素。 Transpose[m][[j]] 表示M的第j列. M[[{i1,i2},{j1,j2}]] 取M的第i1、i2行, j1、j2列构成子矩阵。 例9 构造一个3*3矩阵,再取出它的元素 . 取出第2行 取出第3行、第2列的元素 {a[1,3],a[2,3],a[3,3]}取出第3列 取出由1、3 行，2、3列构成子矩阵