实验指导书2分析和总结.docxVIP

§1 Mathematica 系统简介 实验指导书 Mathematica 用于线性代数 1、Mathematica 系统是一个用计算机做数学的较小巧精干的系统，于 1987 年由Wolfram 创建。当时为Dos 版本 1.2，现已有Dos 与 Windows 版本，版本逐步升级至 4.0 以上。2、Mathematica 的安装（略） 3、Mathematica 的启动（略） 4、Mathematica 的功能，一般可以分为“高级对话式计算器” 与“程序设计”两类。在有了notebook 文件后，两类界限常不明显。 Mathematica 的对话计算器功能有加、减、乘、除、乘方等普通运算、符号运算、解方程、初等函数调用、微积分运算、解微分方程、线性代数计算等。 与其他高级语言（如Basic, Fortran, C, Pascal 等）相比，计算中明显的优点是： 可做整数与有理数的无误差计算； 可做有字符、变量、代数式的运算，如化简，因式分解等； 可做高精度计算，如N[Pi, 1000]可给出 的有 1000 个数的字的近似值； 可做求极限，求导数，求原函数，求定积分，求级数的展开，求解微分方程等高等数学中常见的符号运算与数值计算； 可定义矩阵，向量，数组等，完成线性代数中各种计算，如求行列式，求矩阵的乘积，求矩阵的逆，解线性方程组，求齐次线性方程组的基础解系，求矩阵的特征值与特征向 量等。 易作图。 易扩展。 §2.Mathematica 用于普通数学的最常见的命令 N（求近似值），Sqrt（求平方根），Exp[x]（指数函数 ex），Log[x]（对数函数 lnx），Sin[x] （正弦函数sinx），ArcTan[x（] 反正切函数arctanx），Limit（求极限），（D 求导数），Integrate （求定积分和不定积分），Series（求函数的 Taylor 展式），DSolve（求解微分方程）,Solve （解方程），FindRoot（求根），Plot（作（一元函数）图），Plot3D（作二元函数图），Print[expr] （屏显表达式expr）,Input[ ]（键盘输入）等等。 最后，向 Mathematica 系统询问，可用?A*，?Apart, ??Plot 等向系统查询，系统将（用英语）告诉你这些函数的名字或函数的格式、功能及使用方法。 也可点击Help 求得系统更详尽的帮助。 §3. Mathematica 用于线性代数的常用语句和函数。1、构造向量、矩阵 Table [f[i], {i ,n }]: 一维表或向量 Table [0, {m}, {n}]: 0 矩阵 Table [Random[ ],{m}, {n}]: 数值随机元素阵 Table [ If[I=j, 1, 0],{i,1,5},{j,1,5}]: 一个下三角阵例 1 做一个四阶Hilbert 矩阵. 解法 1 直接写出矩阵, 语句为： h={{1,1/2,1/3,1/4},{1/2,1/3,1/4,1/5},{1/3,1/4,1/5,1/6},{1/4,1/5,1/6,1/7}} 解法 2 利用函数Table, 语句为： h=Table [ 1/(i+j-1),{i,1,4},{j,1,4}] (2)Array[f,{m,n}] 例 2 作一个 4x4 Vandermonge 阵。解 如下几个语句可完成此任务。 x={x1,x2,x3,x4}; f[i_,j_]:=x[[j]]^(i-1); b=Array[f,{4,4}]; Print[b]; Print[MatrixForm[b]]; DiagonalMatrix[list]: 以中元素为对角元的矩阵例 3 DiagonalMatrix[{1,2,3}]//MatrixForm ? 1 0 0? 系统输出为? 0 2 0? ? ? ? ?? 0 0 ? ? IdentityMatrix[n]: n 阶单位阵 MatrixForm[m]这一函数在(2)中已出现，注意它只是形式上较好看，但不好计算。如要计算必须重新复制后方可使用。 2．求得矩阵的部分元素： m[[i,j]]: 矩阵m 的第(i,j)-元素m[[i]]: 矩阵的第i 行 m[[All,j]]: 矩阵的第j 列 Take[m,{i0,i1},{j0,j1}]: 由矩阵m 的第i0 至 i1 行，j0 至 j1 列所组成的子矩阵。m[[{i1,?,ir},{j1,?,js}]]由矩阵m 的第i1 行,?, 第 ir 行，第 j1 列,?，第 js 列组成的子矩阵。 Tr[m,List] 矩阵m 的对角元3．纯量（数量），向量及矩阵的判断与维数确定。 VectorQ[expr]: 表达式 expr 是向量吗