Mathematica软件的基本用法_11数第一次练习.doc

第1章 Mathematica软件的基本用法1.1 启动与运行1.1.1 Mathematica的启动假设在Windows环境下已安装好Mathematica5.0，启动Windows后，在“开始”菜单的“程序”中单击，就启动了Mathematica5.0，在屏幕上显示如图的Notebook窗口，系统暂时为文件取名Untitled-1，直到用户保存时重新命名为止.1.1.2 Mathematica的运行输入1+1，然后按下小键盘上的Enter键(或Shif+Enter键. 注：直接按Enter键是表达式换行)，这时系统开始计算并输出计算结果，并给输入和输出附上次序标识In1和Out1，注意In1是计算后才出现的；再输入第二个表达式Expand(x+y)5，要求系统将一个二项式展开，按Shift+Enter输出计算结果后，系统分别将其标识为ln2和Out2. 如图.为了节省篇幅，本书后面不再把lni与Outi写出. 如将“In1:=1+1,Out1=2”写成“运行1+1得2”.说明：如果把本系统的输出结果复制到Word文档中，在没有安装Mathematica的计算机上无法正常显示，但只要安装了立即变得正常.1.2 基本用法及命令格式1.2.1 常量1数值类型在Mathematic中，基本的数值类型有四种：整数，有理数、实数和复数.如果你的计算机的内存足够大，Mathemateic可以表示任意长度的精确整数及实数，而不受所用的计算机字长的影响. 整数与整数进行四则运算的结果仍是精确的整数或是有理数.例如：2的100次方是一个31位的整数：运行2100得1267650600228229401496703205376.运行1/2+1/3得.当两个整数相除而又不能整除时，系统就用有理数来表示.实数既可以是精确值，也可以是用浮点数表示，Mathematica实数的有效位数可取任意位数，是一种具有任意精确度的近似实数，当然在计算的时候也可以控制实数的精度. 实数有两种表示方法：一种是小数形式表示，另外一种是用指数形式表示. 如：运行0.239998得0.239998.运行0.12*1011得1.21010.实数也可以与整数、有理数进行混合运算，结果还是一个实数.运行2+1/4+0.5得2.75.在Mathematica中，复数是由实部和虚部组成，实部和虚部可以用整数、有理数、实数表示，用I表示虚数单位如：运行3+0.7I得3+0.7 .2.不同类型数的转换在Mathematica的不同应用中，通常对数字的类型要求是不同的. 例如在公式推导中的数字常用整数或有理数表示，而在数值计算中的数字常用实数表示. 在一般情况下在输出行Outn中，系统根据输入行lnn的数字类型对计算结果做出相应的处理. 如果有一些特殊的要求，就要进行数据类型转换.在Mathematica中提供以下几个函数达到转换数据类型的目的：表1-1函数格式意义Nx将x转换成默认6位精度的实数Nx,n将x转换成近似实数，精度为nRationalizex给出实数x的有理数近似值Rationalizex,dx给出x的有理数近似值，误差小于dx举例如下：运行N5/3得1.66667，运行N5/3,20得1.66666666666666666667(怎么会多一个6呢？)1.6666666666666666667.注：在Mathematica4.0版本中，当时，系统自动以6位有效数字输出计算结果，在更高版本中不受此限制.Rationalizex给出x的有理数近似值，如运行Rationalize0.48得. Mathematica尽量保持计算的精确性，例如运行RationalizePi得，运行RationalizeSqrt2得.Rationalizex,dx给出x的误差小于dx的有理数近似值，如运行RationalizePi,1得3；运行RationalizePi,0.1得；运行RationalizePi,0.01得. 继续计算得的此后7个不同精度级的有理数依次为：.用同样的方法算得的前8个各种精度的近似值分别为.数学常数Mathematica中定义了一些常见的数学常数，这些数学常数都是精确数.表1-2符号意义值Pi圆周率=3.14159E自然对数的底e=2.71828Degree1度/180I虚数单位iInfinity无穷大Infinity负的无穷大-GoldenRatio黄金分割数0.61803的倒数黄金分割数应为：0.6180339887. 数学常数可用在公式推导和数值计算中，在数值计算中表示精确值. 如：运行Pi2得，运行NPi2得9.8696.1.2.2 变量1变量的命名变量的命名规则：以字母开头，后跟数字和字母的组合，不能跟具有特殊意义的符号如*,?,%,_,/等，长度不限，但Mathematica中内部函数和命令都是以大写字母开始的标识符. 为了不与它门混淆，读者在自定义变量时应该尽量以小写字母开始，另外在Mathematica中的变量也是区分大小写的. 例如：a12,ast,aST都是合法的变量名，而12a，z*a是非法的(提示与建议：给变量命名时，一方面不要用系统保留字，如Pi,E,C,Sin等；另一方面，尽可能与保留字相比至少有两个以上不一样的字符，不然在运行时系统会提示是否为拼写错误).符号%表示刚刚计算的结果，%2(或%)表示倒数第2个计算结果，%n表示倒数第n个计算结果，在进行交互式计算时可以使用，但在程序中尽量少用.2.给变量赋值在Mathmatica中用等号为变量赋值. 变量不必定义数值类型，可以表示一个数值，一个数组，一个表达式，甚至一个图形. 如：运行x=3得3；运行x2+2x得15；运行x=%+1得16.可同时对不同的变量赋不同的值，例如：运行u,v,w=1,2,3得1,2,3；运行2u+3v+w得11.对于已定义的变量，当你不再使用它时，为防止变量值的混淆，可以随时用“变量名.”清除它的值，如果变量本身也要清除则用函数Clearx，例如运行u=.后再运行2u+v得2+2u(前面的运行已有v=2，而u的值1已经清除).3.变量的替换在给定一个表达式时其中的变量可以取不同的值，这时可用变量替换来计算表达式的不同值. 格式为expr/.x-val，例如：运行Clearx; f=x/2+1得.运行f/.x-1得.运行f/.x-4得3.如果表达式中有多个变量也可以同时替换，例如有两个变量时格式为：expr/.x-xval,y-yval，例如：运行(x+y)*(x-y)2/.x-3,y-1-a得.1.2.3 内建函数在Mathematica的Notebook界面下，可以用交互方式完成各种运算，如绘制函数图象，求极限、解方程等，也可以用它编写象C语言那样的结构化程序. 在Mathematica系统中定义了许多功能强大的函数，我们称之为内建函数(built-in function)， 直接调用这些函数可以达到事半功倍的效果. 这些函数分为两类，第一类是数学意义上的函数，如：绝对值函数Absx，正弦函数Sinx，反正弦函数ArcSinx，余弦函数Cosx，指数函数Expx，以e为底的对数函数Logx，以a为底的对数函数Loga,x等；第二类是命令意义上的函数，如作函数图形的函数Plotfx,x,xmin,xmax，解方程函数Solveeqn,x，求导函数Dfx,x等.注意:(1)在Mathematica中，所有函数严格区分大小写，一般地，内建函数的首写字母必须大写，有时一个函数名是由几个单词构成，则每个单词的首写字母也必须大写，如：求局部极小值函数FindMinimumfx,x,x0等.(2)在Mathematica中，函数名和自变量之间的分隔符是用方括号“ ”，而不是一般数学书上用的圆括号“( )”，初学者很容易犯这类错误.如果输入了不合语法规则的表达式，系统会显示出错信息，并且不给出计算结果，例如：要画正弦函数在区间-10，10上的图形，输入plotSinx,x,-10,10，则系统提示：General:spell1: Possible spelling error: new symbol name plot is similar to existing symbol Plot. More同时原样输出原命令，表示“可能有拼写错误，新符号plot很象已经存在的符号Plot”，由于系统作图命令“Plot”第一个字母必须大写，错误出在”plot”中首字母没有大写. 再输入PlotSinx,x,-10,10，系统又提示“Syntax:bktmcp: Expression PlotSinx,x,-10,10 has no closing .”表示缺少右方括号，并且在命令中将不配对的括号用红色显示.一个表达式只有准确无误地输入后运行才能得出正确结果. 学会看系统出错信息能帮助我们较快找出错误，提高工作效率. 完成各种计算后，点击File-Exit退出，如果文件未存盘，系统提示用户存盘，文件名以“.nb”作为后缀，称为Notebook文件. 以后想使用本次保存的结果时可以通过File-Open菜单读入，也可以直接双击该文件图标，系统自动调用Mathematica将它打开.建议：最好在操作中途保存文件，第一次保存时，点击File-save出现“另存为”窗口，选定保存位置并输入文件名后单击“保存”按钮即可，以后保存直接击Ctrl+S即可.1.2.4 表达式的输入Mathematica提供了多种输入数学表达式的方法. 除了用键盘输入外，还可以使用工具栏或者快捷方式输入运算符、矩阵或数学表达式.1.一维格式输入形如x/(2+3x)+y*(x-w)的表达式称为一维格式，除特殊字符外，所有数学表达式均可用这种格式输入，其优点之一是只用键盘就能完成输入，优点之二是适合在编程中使用，优点之三是可以在Mathematica与Word之间相互复制命令表达式. 在利用该软件处理数学问题时，建议采用一维格式. 这种格式的缺点是有些数学表达式不直观.2.二维格式输入形如的表达式称为二维格式，如果为了特殊需要可以使用二维格式输入.方法一：先按一维格式输入：x/(2+3x)+y/(x-w)，在变量x,y,w未赋值的情况下运行的结果就是二维格式,将其复制到指定位置即可；方法二：可以使用快捷方式输入二维格式. 下面列出了用快捷方式输入二维格式的几种常见形式,如表1-3所示.表1-3数学运算数学表达式按键分式xCtrl+/2N次方x Ctrl+2下标x Ctrl+_2例1.1 如输入数学表达式可以按如下顺序输入按键：(,x,+,1,),Ctrl+,4,-,+,a,Ctrl+_,1,-,Ctrl+/,Ctrl+2,2,x,+,y,-,-.也可以用基本输入工具栏输入二维格式：从FILE菜单中激活Plaettes-Basic Input工具栏，对于常用的特殊字符(如下图)，只要单击这些字符按钮即可输入，使用工具栏可输入更复杂的数学表达式.3.特殊字符输入MathemMatica还提供了用以输入各种特殊符号的工具栏. 若要输入其它的特殊字符或运算符号，必须使用从File菜单中选取palettes，再选Complete Characters工具栏，如右图，单击对应的符号后即可输入.1.2.5 数值的输出形式在数值的输出中, 可以使用转换函数进行不同数据类型和精度的转换. 另外, 对一些特殊要求的格式还可以使用如下的表达式表示形式函数,如表1-4所示.表1-4函数格式意义NumberFormexpr,n以n位精度的实数形式输出实数exprScientificFormexpr以科学记数法输出实数exprEngineergFormexpr以工程记数法输出实数expr例1.2 数值的形式转换 运行A=NPi30,30得8.21289330402749581586503585434.运行NumberFormA,10得NumberForm=8.212893304.运行ScientificFormA的结果与此相同，但下面的命令输出幂值可被3整除的实数：运行EngineeringFormA得EngineeringForm=821.289330402749581586503585434.1.3 函数1.3.1 系统函数在Mathmatic中定义了大量的数学函数可以直接调用，这些函数其名称一般表达了一定的意义，可以帮助我们理解. 几个常用的函数如表1-5所示.表1-5函数格式意义Floorx不比x大的最大整数Ceilingx不比x小的最小整数Signx符号函数Roundx接近x的整数Absxx绝对值Maxx1,x2,x3, x1,x2,x3,中的最大值Minx1,x2,x3,x1,x2,x3,中的最小值Random01之间的随机函数RandomReal,xmax0xmax之间的随机函数RandomReal,xmin,xmaxxminxmax之间的随机函数Expx指数函数Logx自然对数函数lnxLogb,x以b为底的对数函数Sinx,Cosx,Tanx,Cotx,Secx,Cscx三角函数(变量以弧度为单位)ArcSinx, ArcCosx, ArcTanx, ArcCotx, ArcSecx, ArcCscx反三角函数Sinhx,Coshx,Tanhx,Cothx,Sechx,Cschx双曲函数ArcSinhx,ArcCoshx,ArcTanhx,ArcCothx,ArcSechx, ,ArcCschx反双曲函数Modm,nm被n整除的余数，余数与n的符号相同Quotientm,n整除函数(m/n的整数部分)GCDn1,n2,或GCDsn1,n2,的最大公约数，s为一数集合LCMn1,n2,或LCMsn1,n2,的最大公倍数，s为数据集合n!n的阶乘n!n的双阶乘Mathematica中的函数与数学上的函数有些不同的地方，Mathematica中函数是一个具有独立功能的程序模块，可以直接被调用. 同时每一函数也可以包括一个或多个参数，也可以没有参数. 参数的的数据类型也比较复杂. 更加详细的说明可以参看系统帮助，了解各个函数的功能和使用方法是学习Mathematica软件的基础.1.3.2 自定义函数的定义1函数的立即定义立即定义函数的语法是“fx_=expr”，其中f为函数名，x为自变量，expr为表达式. 在执行时会把expr中的x都换为f的自变量x. 函数的自变量具有局部性，只对所在的函数起作用. 函数执行结束后变量的值也就没有了，不会改变其它全局定义的同名变量的值.例1.3 定义函数f(x)=xsinx+x2，求函数值f(3)，绘制其图形.运行Clearx;fx_=Sinx+x2得运行f3得9+Sin3运行Plotft,t,0,2得如图1-4所示的图形.注意：如果运行前x已经赋了值，则不能得到上述函数.对于自定义函数我们可以使用命令Clearf清除掉，而Removef命令则从系统中删除该函数.2延迟定义函数从定义格式上延迟定义函数与即时定义的区别为“=”与“:=”，延迟定义的格式为“fx_:=expr”，其他操作与立即定义函数方式基本相同. 主要区别是：即时定义函数在输入函数后立即定义函数并存放在内存中并可直接调用，如果定义函数前变量已赋值将会把当前值直接代入函数表达式中计算出结果；延迟定义时，其中的变量只是形式变量，不受当前值的影响，只是在调用函数时，才把变量的值代入函数表达式进行计算.建议读者在编程时尽量用延迟定义函数方式定义函数.3多变量函数的定义也可