数学软件选讲29278.pptVIP

数 学 软 件 选 讲 Mathematica Matlab SAS 第 一 篇Mathematica 基础知识 作为一门新的编程语言 图形处理（二维、三维及其参数方程的形式） 极限、微分与积分 求解方程（组）、微分方程（组） 在线形代数方面的应用 数值处理 文件及其它高级操作 1· 条件语句 ◆ 逻辑判断符 == = = != === =!= ◆ 逻辑运算符 ! || ◆ /；运算符 x = a /;test 仅当test为True时才执行赋值语句 ◆ If 语句 语法：If [test, then, else] 若test为 True，则执行then，若test为False，则执行else. 第 二 篇Matlab 截取矩阵块 M[[i]] 取矩阵M的第 i 行 Map[#[[i]], M] 取矩阵M的第 i 列 M[[i, j ]] 取矩阵M的i, j 位置的元素 M[[{i1,…,ir}, {j1,…,js}]] 矩阵M的r×s子 矩阵，元素行标为ik，列标为jk M[[Range{i0,i1}, Range{j0,j1}]]矩阵M的从 i0到i1行， j0到j1列元素组成的子矩阵 3. 矩阵及向量的运算 M.N 对M、N做矩阵乘法（向量内积） M*N 将M、N的对应位置元素相乘 Outer[Times,M,N] 求M、N的外积 Dimensions[ M ] 给出矩阵M的维数 Transpose[ M ] 转置 Inverse[ M ] 求逆 Det[ M ] 方阵M的行列式值 MatrixPower[M,n] n阶矩阵幂 MatrixExp[M] 矩阵指数 Eigenvalues[ M ] M的特征值 Eigenvectors[M] M的特征向量 第四章 求解方程（组）、微分方程（组） 1. 求解多项式方程（组） Solve[ eqns ,vars] 求解多项式方程 Solve[{eqn1,…eqnn}, {var1,…varn}]求解多项式方程组 注：Solve只能给出多项式方程（组）的解，因此它们只适用于幂次不高、规模不大的多项式方程（组）。 NSolve[ eqns ,vars]求多项式方程的数值解 NSolve[{eqn1,…eqnn}, {var1,…varn}]求多项式方程组的数值解 对于数值解，可以直接用NSolve求解 例：求解以下方程（组） x2+ax=2 x3+34x+1=0 x5-1331x+11= 0 2. 求解微分方程（组） DSolve[ eqns ,y[x], x] 求解y[x]的微分方程 DSolve[ eqns ,y, x]以纯函数的形式给出y的解 DSolve[{eqn1, eqn2,…}, {y1, y2, …}, x]求解微分方程组 例：求解以下微分方程（组） y’ = yy’’－ k y =1 第五章 数值处理 1. 数值积分 NIntegrate[expr , {x,xmin,xmax}] 注意，NIntegrate直接计算数值积分，不先给出符号结果，而Integrate[…]//N会尽可能的先求精确解的形式。 数值根求解 FindRoot[lhs=rhs , {x, x0}]以x0为初始点求方程的数值解 FindRoot[lhs=rhs , {x, {x0 ,x1}}] 给出两个 初值求数值根（方程的符号导数无法求出 时，必须使用此形式） FindRoot[{eqn1, eqn2,…}, {x, x0},{y, y0 }, …] 对联立方程 eqni 求数值解 例：求解下列方程（组） cosx =x x600+5x+3=0 3. 微分方程数值解 NDSolve[{eqn1, eqn2,…}, y,{x,xmin,xmax}] 求函数y的数值解，x的范围为{xmin,xmax} NDSolve[{eqn1, eqn2,…},{y1, y2, …}, {x,xmin,xmax}] 求函数yi的数值解 注：以上两种形式用于求解常微分方程（组） NDSolve以InterpolatingFunction目标生成函数yi的解。 InterpolatingFunction目标提供独立变量x在xmin到xmax范围内yi的近似值。 例：求解以下微分方程（组）并画出函数y的图形 NDSolve[{eqn1, eqn2,…}, y,{x,xmin,xmax}, {t,tmin,tmax}] 求由函数y构成的偏微分方 程的数值解 NDSolve[{eqn1, eqn2,…