mathmatic用plot函数画多条曲线画不出曲线,坐标也没有出来

豆丁微信公众号
君,已阅读到文档的结尾了呢~~
mathematica_4.0_课件教程[整理版]
扫扫二维码,随身浏览文档
手机或平板扫扫即可继续访问
mathematica_4.0_课件教程[整理版]
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由:
将文档分享至:
分享完整地址
文档地址:
粘贴到BBS或博客
flash地址:
支持嵌入FLASH地址的网站使用
html代码:
&embed src='http://www.docin.com/DocinViewer--144.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布,请您等待!
3秒自动关闭窗口Mathematica教程ppt下载_爱问共享资料
Mathematica教程.ppt
Mathematica教程.ppt
Mathematica教程.ppt
简介:本文档为《Mathematica教程ppt》,可适用于高等教育领域,主题内容包含Mathematica教程Mathematica教程第章Mathematica概述第章Mathematica的基本量第章Mathematica的基本符等。
侵权或盗版
*若权利人发现爱问平台上用户上传内容侵犯了其作品的信息网络传播权等合法权益时,请按照平台要求书面通知爱问!
赌博犯罪类
67人已下载
在此可输入您对该资料的评论~
添加成功至
资料评价:mathematica-超级教程_福州大学数学与计算机科学学院:Mathematica课件(张碧霞)_pdf_大学课件预览_高等教育资讯网
福州大学数学与计算机科学学院:Mathematica课件(张碧霞):mathematica-超级教程
分类: 格式: 日期:日
Mathematica 是美国 Wolfram 研究公司生产的一种数学分析型的软件,以符号计算见长,也具有高精度的数值计算功能和强大的图形功能。假设在 Windows 环境下已安装好 Mathematica4.0,启 动 Windows 后,在,开始,菜单的,程序,中单击,就启动了 Mathematica4.0,在屏幕上显示如图的 Notebook 窗口,系统暂时取名 Untitled-1,直到用户保存时重新命名为止输入 1+1,然后按下 Shif+Enter 键,这时系统开始计算并输出计算结果,并给输入和输出附上次序标识 In[1]和 Out[1],注意 In[1]是计算后才出现的;再输入第二个表达式,要求系统将一个二项式展开,按 Shift+Enter 输出计算结果后,系统分别将其标识为 In[2]和 Out[2].如图在 Mathematica 的 Notebook 界面下,可以用这种交互方式完成各种运算,如函数作图,求极限、解方程等,也可以用它编写像 C 那样的结构化程序。在 Mathematica 系统中定义了许多功能强大的函数,我们称之为内建函数(built-in function ),直接调用这些函数可以取到事半功倍的效果。这些函数分为两类,一类是数学意义上的函数,如:绝对值函数 Abs[x],正弦函数 Sin[x],余弦函数 Cos[x],以 e 为底的对数函数 Log[x],以 a 为底的对数函数 Log[a,x]等;第二类是命令意义上的函数,如作函数图形的函数 Plot[f[x],{x,xmin,xmax}],解方程函数 Solve[eqn,x],求导函数 D[f[x],x]等。必须注意的是,Mathematica 严格区分大小写,一般地,内建函数的首写字母必须大写,有时一个函数名是由几个单词构成,则每个单词的首写字母也必须大写,如:求局部极小值函数 FindMinimum[f[x],{x,x0]等。第二点要注意的是,在 Mathematica 中,函数名和自变量之间的分隔符是用方括号,[ ]”,而不是一般数学书上用的圆括号,( ),,初学者很容易犯这类错误。如果输入了不合语法规则的表达式,系统会显示出错信息,并且不给出计算结果,例如:要画正弦函数在区间[-10,10] 上的图形,输入 plot[Sin[x],{x,-10,10}],则系统提示,可能有拼写错误,新符号 ‘plot’ 很像已经存在的符号 ‘Plot’”,实际上,系统作图命令,Plot”第一个字母必须大写,一般地,系统内建函数首写字母都要大写。再输入 Plot[Sin[x],{x,-10,10},系统又提示缺少右方括号,并且将不配对的括号用蓝色显示,如图一个表达式只有准确无误,方能得出正确结果。学会看系统出错信息能帮助我们较快找出错误,提高工作效率。 完成各种计算后,点击 File-&Exit 退出,如果文件未存盘,系统提示用户存盘,文件名以,.nb”作为后缀,称为 Notebook 文件。以后想使用本次保存的结果时可以通过 File-&Open 菜单读入,也可以直接双击它,系统自动调用 Mathematica 将它打开,1.1.2 表达式的输入Mathematica 提供了多种输入数学表达式的方法。除了用键盘输入外,还可以使用工具样或者快捷方式健入运算符、矩阵或数学表达式。1.数学表达式二维格式的输入Mathematic 担提供了两种格式的数学表达式。形如 x/(2+3x)+y*(x-w)的称为一维格式,形如的称为二维格式。你可以使用快捷方式输入二维格式,也可用基本输入工具栏 输入二维格式。下面列出了用快捷方式输入二维格式的方法数学运算 数学表达式 按键分式 x Ctrl+/ 2n 次方x Ctrl+^ 2开 n 次方 Ctrl +2 x下标 x Ctrl+_ 2如果要取消二维格式输入 按下 Ctrl+SPACE(空格) 例如输入数学表达式可以按如下顺序输入按键,(,x,+,1,),Ctrl+^,-&,+,a,Ctrl+_,1,-&,Ctrl+/,Ctrl+2,2,x,+,y,-&,-&另外也可从 FILE 菜单中激活 Plaettes-&Basic Input 工具栏,也可输入,并且使用工具栏可输入更复杂的数学表达式。如下图,2.特殊字符的输入MathemMatica 还提供了用以输入各种特殊符号的工具样。基本输入 工具样包含了常用的特殊字符(上图),只要单击这些字符按钮即可输入。若要输入其它的特殊字符或运算符号,必须使用从 FILE 菜单中选取 Complete Characters 工具栏,如下图,单击符号后即可输入。1.2.Mathematica 的联机帮助系统用 Mathematica 的过程中,常常需要了解一个命令的详细用法,或者想知系统中是否有完成某一计算的命令,联机帮助系统永远是最详细、最方便的资料库。1.获取函数和命令的帮助在 Notebook 界面下,用?或 可向系统查询运算符、函数和命令的定义和用法,获取简单而直接的帮助信息。 例如,向系统查询作图函数 Plot 命令的用法? Plot 系统将给出调用 Plot 的格式以及 Plot 命令的功能(如果用两个问号,”,则信息会更详细一些)。 Plot* 给出所有以 Plot 这四个字母开头的命令2.Help 菜单任何时候都可以通过按 F1 键或点击帮助菜单项 Help Browser,调出帮助菜单,如图所示,其中的各按钮用途如表所示Built-in Function 内建函数,按数值计算、代数计算、图形和编程分类存放Add-ons 有程序包(Standard Packages ) MathLink Library 等内容The Mathematica Book 一本完整的 Mathematica 使用手册Getting Started/Demos 初学者入门指南和多种演示Other Information 菜单命令的快捷键,二维输入格式等Master Index 按字母命令给出命令、函数和选项的索引表如果要查找 Mathematica 中具有某个功能的函数,可以通过帮助菜单中的 Mahematica 使用手册,通过其目录索引可以快速定位到自己要找的帮助信息。例如:需要查找 Mathematica 中有关解方程的命令,单击,The Mathematica Book”按钮,再单击,Contents”,在目录中找到有关解方程的节次,点击相应的超链接,有关内容的详细说明就马上调出来了。如果知道具体的函数名,但不知其详细使用说明,可以在命令按钮 Goto 右边的文本框中键入函数名,按回车键后就显示有关函数的定义、例题和相关联的章节。例如,要查找函数 Plot 的用法,只要在文本框中键入 Plot,按回车键后显示如图的窗口,再按回车键,则显示 Plot 函数的详细用法和例题。如果已经确知 Mathematica 中有具有某个功能的函数,但不知具体函数名,可以点击 Built-in Functions 按钮,再按功能分类从粗到细一步一步找到具体的函数,例如,要找画一元函数图形的函数,点击 Built-in Functions -&Graphics and Sound-&2DPlots-&Plot,找到 Plot 的帮助信息。如果知道具体的函数名,但不知其详细使用说明,可以在命令按钮 Goto 右边的文本框中键入函数名,按回车键后就显示有关函数的定义、例题和相关联的章节。例如,要查找函数 Plot 的用法,只要在文本框中键入 Plot,按回车键后显示如图 1-5 的窗口,再按回车键,则显示 Plot 函数的详细用法和例题。如果已经确知 Mathematica 中有具有某个功能的函数,但不知具体函数名,可以点击 Built-inFunctions 按钮,再按功能分类从粗到细一步一步找到具体的函数,例如,要找画一元函数图形的函数,点击 Built-in Functions -&Graphics and Sound-&2D Plots-&Plot,找到 Plot 的帮助信息。2.1 数据类型和常数1 数值类型在 Mathematic 中,基本的数值类型有四种:整数,有理数、实数和复数如果你的计算机的内存足够大,Mathemateic 可以表示任意长度的精确实数,而不受所用的计算机字长的影响。整数与整数的计算结果仍是精确的整数或是 有理数。例如:2 的 100 次方是一个 31位的整数,ln[1]:=2^100Out[1]=在 Mathematica 中允许使用分数,也就是用有理数表示化简过的分数。当两个整数相除而又不能整除时,系统就用有理数来表示,即有理数是由两个整数的比来组成如,Ln[2]:=Out[2]=实数是用浮点数表示的,Mathematica 实数的有效位可取任意位数,是一种具有任意精确度的近似实数,当然在计算的时候也可以控制实数的精度。实数有两种表示方法,一种是小数点另外一种是用指数方法表示的。 如,ln[3]:=0.239998Out[3]=0.23998ln[4]:=0.12*10^11Out[4]=0.12*10^11实数也可以与整数,有理数进行混合运算,结果还是一个实数。Ln[5]:=2+1/4+0.5Out[5]=2.75复数是由实部和虚部组成。实部和虚部可以用整数,实数,有理数表示。在 Mathematica 中,用i 表示虚数单位如,Ln[6]:=3+0.7iOut[6]:=3+0.7i2.不同类型数的转换在 Mathematica 的不同应用中,通常对数字的类型要求是不同的。例如在公式推导中的数字常用整数或有理数表示,而在数值计算中的数字常用实数表示。在一般情况下在输出行 Out[n]中,系统根据输入行 ln[n]的数字类型对计算结果做出相应的处理。如果有一些特殊的要求,就要进行数据类型转换。在 Mathematica 中的提供以下几个函数达到转换的目的,N[x] 将 x 转换成实数N[x,n] 将 x 转换成近似实数,精度为 nRationalize[x] 给出 x 的有理数近似值Rationalize[x,dx] 给出 x 的有理数近似值,误差小于 dx[举例 ]ln[1]=N[5.3,20]Out[1]=1.ln[2]:=N[%,10]Out[2]=1.二行输出是把上面计算的结果变为 10 位精度的数字。%表示上一输出结果。Ln[3]=Rationalize[%]Out[3]=5/33.数学常数Mathematica 中定义了一些常见的数学常数,这些数学常数都是精确数,例如表示圆周率。Pi 表示 =3.14159……E 自然对数的底,e=2.71828……,Degree /180i 虚数单位Infinity 无穷大- infinity 负的无穷大 -GondenRatio 黄金分割数 0.61803数学常数可用在公式推导和数值计算中。在数值计算中表示精确值:如,n[1]:=Pi^2Out[1]=ln[2]:=Pi^2//NOut[2]=9.8696l4.数的输出形式在数的输出中可以使用转换函数进行不同数据类型和精度的转换。另外对一些特殊要求的格式还可以使用如下的格式函数,NumberForm[expr,n] 以 n 位精度的实数形式输出实数 exprScientificFormat[expr] 以科学记数法输出实数 exprEngineergForm[expr] 以工程记数法输出实数 expr例如,ln[1]:=N[Pi^30,30]ln[2]:=NumberForm[%,10]Out[2]//NumberForm=8.下面的函数输出幂值可被 3 整除的实数Ln[3]=EngineeringForm[%%]Out[6]//EngineeringForm=1.变量的命名Mathematica 中内部函数和命令都是以大写字母开始的标示符。为了不会与它门混淆,我们自定义的变量应该是以小写字母开始,后跟数字和字母的组合,长度不限。 例如,a12,ast,aST 都是合法的,而 12a,z*a 是非法的。另外在 Mathematica 中的变量是区分大小写的 在 Mathematica 中,变量不仅可以存放一个数值,还可以存放表达式或复杂的算式。2.给变量赋值在 Mathmatica 中用等号=为变量赋值。同一个变量可以表示一个数值,一个数组,一个表达式,甚至一个图形。如,Ln[1]:=x=3Out[1]=3Ln[2]:=x^2+2xOut[2]=15Ln[3]:=x=%+1Out[3]=16对不同的变量可同时赋不同的值例如,Ln[4]:={u,v,w}={1,2,3}Out[4]={1,2,3}Ln[5]:=2u+3v+wOut[5]=11对于已定义的变量,当你不再使用它是,为防止变量值的混淆,可以随时用=,清除他的值,如果变量本身也要清除用函数 Clear[x]例如ln[6]:=u=,ln[7]:=2u+vOut[7]=2+2u3.变量的替换在给定一个表达式时其中的变量可能取不同的值,这是可用变量替换来计算表达式的不同值。方法为用 expr/.例如,Ln[1]:=f=x/2+1Out[1]=Ln[2]:=f/.x-&1Out[2]=Ln[3]:=f/.-&2Out[3]=3如果表达式中有多个变量也可以同时替换方法为例如有两个,expr/.{x-&xval,y-&val}Ln[4]:=(x+y)(x-y)^2/.{x-&3,y-&1-a}Out[4]=2.3 函数1.系统函数在 Mathmatic 中定义了大量的数学函数可以直接调用,这些函数其名称一般表达了一定的意义,可以帮助我们理解。下面是几个常用的函数,Floor[x] 不比 x 大的最大 整数Ceiling[x] 不比 x 小的最小整数Sign[x] 符号函数Round[x] 接近 x 的整数Abs[x] x 绝对值Max[x1,x2,x3……..] x1,x2,x3…….中的最大值Min[x1,x2,x3……..] x1,x2,x3…….中的最小值Random[] 0~1 之 间的随机函数Random[Real,xmax] 0~xmax 之间的随机函数Random[Real,{xmin,xmax}] xmin~xmax 之间的随机函数Exp[x] 指数函数Log[x] 自然对数函数 lnxLog[b,x] 以 b 为底的对数函数Sin[x],Cos[x],Tan[x],Csc[x],Sec[x],Cot[x] 三角函数(变量是以弧度为单位的)Sinh[x],Cosh[x],Tanhx[x],Csch[x],Sech[x],Coth[x] 双曲函数ArcSech[x],ArcCoth[x] 双曲函数Mod[m,n] m 被 n 整除的余数,余数与 n 的符相同Quotient[m,n] m/n 的整数部分GCD[n1,n2,n3……]或 GCD[s] n1,n2,…的最大公约数,s 为一数集合LCM[n1,n2……]或 LCM[s] n1,n2…….的最大公倍数,s 为数据集合N! n的阶程N!! n 的双阶程Mathematica 中的函数与数学上的函数有些不同的地方,Mathematica 中函数是一个具有独立功能的程序模块,可以直接被调用。同时每一函数也可以包括一个,或多个参数,也可以没有参数。参数的的数据类型也比较复杂。更加详细的可以参看系统的帮助,了解各个函数的功能和使用方法是学习 Mathematica 软件的基础2.函数的定义( 1)函数的立即定义立即定义函数的语法如下 f[x_]=expr 函数名为 f,自变量为 x,expr 是表达式。在执行时会把 expr 中的 x 都换为 f 的自变量 x(不是 x_)。函数的自变量具有局部性,只对所在的函数起作用。函数执行结束后也就没有了,不会改变其它全局定义的同名变量的值。请看下面的例子定义函数f(x)=x*Sinx+x2 对定义的函数我们可以求函数值,也可绘制它的图形。对于定义的函数我们可以使用命令 Clear[f]清除掉而 Remove[f]则从系统中删除该函数。( 2).多变量函数的定义也可以定义多个变量的函数,格式为 f[x_,y_,z_,…]=expr自变量为 x,y,z….,相应的 expr中的自变量会被替换。例如定义函数f(x,y)=xy+ycosx( 3),延迟定义函数延迟定义函数从定义方法上与即时定义的区别为,=”与,,=”延迟定义的格式为 f[x_],=expr其他操作基本相同。那么延迟定义和即时定义的主要区别是什么?即时定义函数在输入函数后立即定义函数并存放在内存中并可直接调用。延时定义只是在调用函数时才真正定义函数。( 4),使用条件运算符定义和 If命令定义函数如果要定义如,这样的分段函数应该如何定义,显然要根据 x 的不同值给出不同的表达式。一种办法是使用条件运算符,基本格式为 f[x_]:=expr/;condition,当 condition条件满足时才把 expr赋给 f.下面定义方法,通过图形可以验证所定义函数的正确性当然使用 If命令也可以定义上面的函数,If语句的格式为 If[条件,值 1,值 2]如果条件成立取,值 1”,否则取,值 2”,下面用 If语句的定义结果可以看出用 If定义的函数 g(x)和前面函数 f(x)相同,这里使用了两个 If嵌套。逻辑性比较强关于其他的条件命令的进一步讨论请看后面的章 节2.4 表将一些相互关联的元素放在一起,使它们成为一个整体。既可以对整体操作,也可以对整体中的一个元素单独进行操作。在 Mathematica 中这样的数据结构就称作表( List)。表主要有三个用法:表{ a,b,c}可以表示一个向量;表 {{a,b},{c,d}}可表示一个矩阵。1.建 表在表中元素较少时,可以采取直接列表的方式列出表中的元素,如{1,2,3 },请看下面的操作Ln[1]:={1,2,3}Out[1]={1,2,3}下面是符号表达式的列表Ln[2]:=1+%x+x^%Out[2]={1+2x,1+2x+x2,1+3x+x2}下面是对列表中的表达式对 x求 导Ln[3]:=D[%,x]Out[3]={2,2+2x,3+2x}Ln[4]:=%/.x-&1Out[4]={2,4,5}如果表中的元素较多时,可以用建表函数进行建表。Table[f,{I,min,max,step}] 以 step为步长给出 f的数值表,i由 min变到 max,Table[f,{min,max}] 给出 f的数值表,I由 min变到 max 步长为 1Table[f,max] 给出 max个 f的表Table[f,{I,imin,imax},{j,jmin,jmax},….] 生成一个多维表TableForm[list] 以表格格式显示一个表Range[n] 生成一个{ 1,2,……..}的列表Range[n1,n2,d] 生成{ n1,n1+d,n1+d,….,n2}的列 表下面给出 x乘 i的值的表,i的变化范围为 [2,6]:Ln[1]:=Table[x*i,{i,2,6}]Out[1]={2x,3x,4x,5x,6x}Ln[2]:=Table[x^2,{4}]Out[2]={x2,x2,x2,x2}用 Range函数生成一个序列数Ln[3]:=Range[10]Out[3]={1,2,3,4,5,6,7,8,9,10}下面这个序列是以步长为 2,范围从 8 到 20Ln[4]:=Range[8,20,2]Out[4]={8,10,12,14,16,18,20}上面的参数变化都是只有一个,也可制成包括多个参数的表,下面生成一个多维表:Ln[5]:=Table[2i+j,{i,1,3},{j,3,5}]Out[5]={{5,6,7},{7,8,9},{9,10,11}}使用函数 TableForm可以以表格的方式输出Ln[6]:=%//TableFormOut[6]//TableForm=5 6 77 8 99 10 112.表的元素的操作当 t表示一个表时,t[[i]]表示 t中的第 i个子表。如果 t={1,2,a,b}那么 t[[3]]表示,a”。如:ln[1]:=t=Table[I+2,j{I,1,3},{j,3,5}]Out[1]={{7,9,11},{8,10,12},{9,11,13}}ln[2]:=t[[2]]Out[2]={8,10,12}对于表的操作 Mathematica提供了丰富的函数,详细的可以查阅后面的附录或者系统帮助。3.1 多项式的表示形式可认为多项式是表达式的一种特殊的形式,所以多项式的运算与表达式的运算基本一样,表达式中的各种输出形式也可用于多项式的输出。 Mathematica 提供一组按不同形式表示代数式的函数。Expand[ploy] 按幂次展开多项式 ployExpand[ploy] 全部展开多项式 ployExpandAll[ploy] 全部展开多项式 ployFactor[ploy] 对多项式 poly 进行因式分解FactorTerms[ploy,{x,y,…}] 按变量 x,y,… 进行分解Simplify[poly] 把多项式化为最简形式FullSimplify[ploy] 把多项式展开并化简Collect[ploy,x] 把多项式 poly 按 x 幂展开Collect[poly,{x,y…}] 把多项式 poly 按 x,y….的幂次展开1.下面是一些例子(1).对 x8-1 进行分解(2).展开多项式( 1+x) ^5(3).展开多项式( 1+x+3y) 4(4).化简( 2+x) ^4(1+x)^4(3+x)^32.多项式的代数运算多项式的运算有加、减、乘、除运算,+,-,*,/ 下面通过例子说明。(1).多项式的加运算 a2+3a+2 与 a+1 相加(后面例子中也使用这两个多项式运算(2).多项式相减(3).多项式相乘(4).多项式相除(5).另外使用 Cancel函数可以约去公因式两个多项式相除,总能写成一个多项式和一个有理式相加 Mathematic中提供两个函数 PolynomialQuotient和 PolynomialRemainder分别返商式和余式。例如,3,2 方程及其根的表示因为Mathematica 把方程看作逻辑语句。在数学方程式表示为形如,x2-2x+1=0”的形式。在 Mathematica中,=”用作赋值语句,这样在 Mathematica中用,==”表示逻辑等号,则方程应表示为,x^2-2x+1==0” 。方程的解同原方程一样被看作是逻辑语句。例如用 Roots求方程 x^2-3x+2 的根显示为这种表示形式说明 x取 1 或 2 均可。而用 Solve[]可得解集形式。1 求解一元代数方程下面是常用的一些方程求解函数Solve[lhs==rhs,vars] 给出方程的解集NSolve[lhs==rhs,vars] 直接给出方程的数值解集Roots[lhs==rhs,vars]求表达式的根FindRoot[lhs==rhs,{x,x0}求 x=x0 时,方程的解值先看 Solve函数例子Solve函数可处理的主要方程是多项式方程。 Mathematica总能对不高于四次的方程进行精确求解,对于三次或四次方程,解的形式可能很复杂。例如求 x3+5x+3=0这时可用 N函数近似数值解,当方程中有一些复杂的函数时,Mathematica可能无法直接给出解来。在这种情况下我们可用 FindRoot[]来求解,但要给出起始条件。例如:求 3Cosx=logx的解但只能求出 x=1 附近的解,如果方程有几个不同的解,当给定不同的条件时,将给出不同的解。 如上例若求 x=10 附近的解命令为,因此确定解的起始位置是比较关键,一种常用的方法是,先绘制图形观察后再解如上例通过图形可断定在 x=5 附近有另一根2.求方程组的根使用 Solve 和 NSolve,FindRoot 也可求方程组的解,只是使用时格式略有不同下面给出一个 Solve 函数的例子,求解,3 求方程的全解如果我们求 ax2+bx+c=0的根我们用 Solve函数解的结果是,这显然是不合理的,因为对不同的 a,b,c方程的解有不同的情况,而上面只是给出部分解如果要解决这个问题可用 Reduce命令,它可根据,a,b,c的取值给出全部值。因此 Solve,Roots只给出方程的一般解,而 Reduce函数数可以给出方程的全部可能解。4.解条件方程在作方程计算时,可以把一个方程看作你要处理的主要方程,而把其他方程作为必须满足的辅助条件,你将会发现这样处理很方便。譬如在求解像 这样的方程时,通常我们采用 的代换方法使求解方程得到简化。在 Mahematica中,我们通常是首先命名辅助条件组,然后用名字把辅助条件包含在你要用函数 Solve[] 求解的方程组中。用 Sc定义方程:,在这种条件下,求解方程。3.3 求和与求积在 Mathematica 中,数学上的各式符号 用 Sum 表示,连乘 用 Product 表示。下面列出求各与求积函数的形式和意义,Sum[f,{i,imin,imax}] 求和Sum[f,{i,imin,imax,di}] 以步长 di 增加 i 求和Sum[f,{i,imin,imax},{j,jmin,jmax}] 嵌套求和Product[f,{i,imain,imax}] 求积Product [f,{i,imin,imax,di}] 以步长 di 增加 i 求和Product[f,{I,imin,imax},{j,jmin,jmax}] 嵌套求积Nsum[f,{i,imin,Infinity}] 求 近似值NProduct[f,{i,imin,Infinity}] 求 近似值一些例子.1 基本的二维图形Mathematica 在直角坐标系中作一元函数图形用下列基本命令。Plot[f,{x,xmin,xmax},option-&value] 在指定区间上按选项定义值画出函数在直角坐标系中的图形,Plot[{f1,f2,f3,…},{x,xmin,xmax},option-&&value] 在指定区间上按选项定义值同时画出多个函数在直角坐标系中的图形Mathematica 绘图时允许用户设置选项值对绘制图形的细节提出各种要求。例如:要设置图形的高宽比,给图形加标题等。每个选项都有一个确定的名字,以,选项名 -&选项值,的形式放在 Plot 中的最右边位置,一次可设置多个选项,选项依次排列,用逗号隔开,也可以不设置选项,采用系统的默认值。选项 说明 默认值AspectRatio 图形的高、宽比 1/0.618AxesLabel 给坐标轴加上名字 不加PlotLabel 给图形加上标题 不加PlotRange 指定函数因变量的区间 计算的结果PlotStyle 用什么样方式作图(颜色,粗细等) 值是一个表PlotPoint 画图时计算的点数 251.举例(1).例如绘制 的图形。(2).如果要取消刻度可以使用 Ticks 选项(3).如果要标注坐标名称 x 轴为,Time”,y 轴为,Height”(4).将坐标交点(3,0),并标注图形名称。(5).修改 x 方向的刻度,y 轴方向的刻度则用默认值。(6).定义 y 轴的绘图范围(7).另外我们也可以将图形结果定义给变量,但不显示图形,后用 Show 命令显示。2.数据集合的图形Mathematica 用于绘数字集合的图形的命令与前而介绍的绘函数图形的命令是相似的。如下,ListPlot[{y1,y2,…..}] 绘出在 x 的值为 1,2…时 y1,y2,…的图形ListPlot[{{x1,y1},{x2,y2},…..}] 绘出离散点(xi,yi )ListPlot[List,PlotJoined-&True] 把离散点连成曲线(1).下面举例说明下面是一个离散数据的集合的图形3.二维参数作图前面我们使用 Plot 命令可以绘出直角坐标系下的函数图形,使用 ParametrecPlot 可以绘制参数曲线下面给出 ParametricPlot 的常用形式ParametricPlot[{fx,fy},{t,tmin,tmax}] 绘出参数图ParametricPlot[{fx,fy},{gx,gy},….{t,tmin,tmax}] 绘出一组参数图ParametricPlot[{fx,fy},{t,tmin,tmax},AspectRatio-&Automatic] 设法保持曲线的形(1).绘制参数方程 的图形(2).下面将一个园与上面参数绘在同一个坐标下,并保证图形的形状正确。4,1 二维图形元素用图形元素绘图适合于绘制结构复杂的图形。Mathematica 中还提供了各种如绘制点、线段、圆弧等函数。同样我们可先用 Grahpics 作出平面图形的表达式,再用 Show 显示守成的图形。下面给出在 Mathematica 中常用的二维图形元素。____________________________________________________________________Point[[x,y]] 点Line[{{x1,y1},{x2,y2},…}] 线段Rectangle[{xmin,ymin},{xmax,ymax}] 填充矩阵Polygon[{{x1,y1},{x2,y2},…….}] 填充多边形Circle[{x,y},r] 圆Circle[{x,y},{rx,ry}] 半轴分别为 rx,ry 的椭圆Circle[{x,y},r,{theta1,thata2}] 圆弧Circle[{x,y},{rx,ry},{theta1,theta2}] 椭圆弧Disk[{x,y},r] 填充圆Raster[{{a11,a12,…..},{a21,……},….}] 灰度在 0 到 1 之间的灰层组Text[Expr,{x,y}] 文本大小下图绘出一个有颜色和大小的点,且在图形四周插入文本下面绘制一些有线条组成的图形当然也可以添加坐标轴下面的例子,说明了这一点。下面的例子,是说明了 Retangle 的图形绘制,例子中用一些小矩形逼近正弦曲线与 x 轴所成面积。程序中生成一个图形集合并显示出来。4,3 图形的样式我们称图形的颜色、曲线的形状和宽度等特性为图形样式。在本节中,我们就图形的各种样式,尤其是曲线的样式进行学习。下面给出选项用于设置图形样式。Graykvel[] 灰度介于 0(黑 )到 l(白 )之间RGBColor[r,g,b]由红、绿,蓝组成的颜色,每种色彩取 0 到 1 之间的数Hue[A] 取 0 到 1 之间的色彩Hue[h,s,b] 指定色调,位置和亮度的颜色,每项介于 0 到 1 之间PointSize[d] 给出半径为 d 的点,单位是 Plot 的一个分数AbsolutePointSize[d] 给出半径为 d 的点 (以绝对单位量取 )Thickness[w] 给所有线的宽度 w,单位是 Plot 的分数AbsoluteThickness[w] 给所有线的宽度 w,(以绝对单位量取 )Dashing[wl,w2,….] 给所有线为一系列虚线,虚线段的长度为 wl,w2,…Absolutedashing[{w1,w2,…..}] 以绝对单位给出虚线长度PlotStyle-&style 设立 Plot 中所有曲线的风格PlotStyle-&{{Style1},{Style2}…….} 设立 Plot 中一些列曲线的风格MeshStyle-&Style 设立宽度和表面网格的风格1.图形颜色的设置在 Mathematicaa 提供各种图形指令中,对图形元素颜色的设置是一个很重要的设置,。下面给出三条不同颜色的正弦曲线,此处以灰度表示,即颜色深浅不同。下面用不同的色调对三个菱形进行着色。2.图形大小下面是一些点,注意点大小的控制。下面的点的控制是用绝对单位3.线段的控制下面的例子是控制线段的宽度,使用的是绝对控制。Mathematica 提供的虚线指令可生成多种不同的复杂虚线。4.4 图形的重绘和组合每次绘制图形后,Mathematica 保存了图形的所有信息,所以用户可以重绘 这些图形。我们在重绘图形的时候,还可以改变一些使用。下面是常用重绘图形的函数。Show[plot] 重绘图形Show[plot,option-&value] 改变方案重绘图形Show[plot1,plot2,plot3…] 多个图形的绘制Show[GraphcisArray[{{plot1,plot2,…}…}]] 绘制图形矩阵InputForm[plot] 给出所有的图形信息1.使用 Show 显示图形下面绘制函数 Sin[x^2] 的图形。重绘图形时,可以改变命令的设置,下面改变 y 的比例同时给 图边框2.使用 Show 命令进行组合也可使用 Show 进行图形组合。图形组合与图形是否有相同的比例无关,这是 Mathematica 会自动选择新的比例来绘制图形。下面绘制函数 -xsin(2x+Pi)的图形和 xcos(2x) 然后绘制在一张图时。3,将多个图形组合为一个图形我们也可把图形组合为一个图形,我们还可以用 GraphicsArray 把多个图形绘制在一个图形矩阵中如下图。4,5 基本三维图形绘制函数 f(x,y)在平面区域上的三维立体图形的基本命令是 Plot3D,Plot3D 和 Plot 的工作方式和选项基本相同。ListPlot3D 可以用来绘制三维数字集合的三维图形,其用法也类似于 listPlot,下面给出这两个函数的常用形式。Plot3D[f,(x,xmin,xmax),(y,ymin,ymax)] 绘制以 x 和 y 为变量的三维函数/的图形ListPlot3D[{Z11,Z12,…},{Z21,Z22,…},…..]] 绘出高度为 Zvx 数组的三维图形Plot3D 同平面图形一样,也有许多输出选项,你可通过多次试验找出你所需的最佳图形样式。选项 取值 意义Axes True 是否包括坐标轴AxesLabel None在轴上加上标志,zlabel 规定 z 轴的标志,{xlabel,ylabel,zlabel}规定所有轴的标志Boxed True 是否在曲面周围加上立方体ColorFunction Automatic使用什么颜色的明暗度; Hue 表示使用一系列颜色TextStyle STextStyle 用于图形文本的缺省类型ormatType StandardForm 用于图形文本的缺省格式类型DisplayFunction SdlisplayFunction 如何绘制图形,Indentity 表示不显示FaceGrids None如何在立体界面上绘上网格; All 表示在每个界面上绘上网格HiddenSurface True 是否以立体的形式绘出曲面Lighdng True 是否用明暗分布米给表面加色Mesh True 是否在表面上绘出 xy 网格PlotRange Automatic图中坐标的范围;可以规定为 All,{zmin,zmax}或 {xminn,xmax},{ymin,ymax},{zmin,zmax}Shading True 表面是用阴影还是留空白ViewPoint {1,3,-2.4,2} 表面的空间观察点1.三维绘图举例(1).函数 sin(x+y)cos(x+y)的立体图(2),对于三维图形中 Axes,Axeslabel,Boxed 等操作同二维图形的一些操作很相似。用 PlotRange 设定曲线的表面的变化范围。(3).图形轴上加上标记,且在每个平面上画上网格。(4).视图的改变学习过画法几何或工程制图的都知道,制图时通常用三视图来表示一个物体的具体形状特性。我们在生活中也知道从不同观察点观察物体,其效果是很不一样的。 Mathematica 在绘制立体图形时,在系统默认的情况下,观察点在 (1,3,-2,4,2)处。这个参考点选择是具有一般性的,因此偶尔把图形的不同部分重在一起也不会发生视觉混乱。下面例子改变观察视点。从上面我们可以看出,观察点位于曲面的上方有利于看清对于图形全貌。对于较复杂的图形,我们在所绘的图形上包括尽可能多的曲线对于我们观察很有帮助。同时,在曲面的周围直接绘出立方体盒子也有利于我们认清曲面的方位。(6).下面是没有网格和立体盒子的曲面图,它看起来就不如前面的图形清晰明了。(7).下图给出没有阴影的曲面带有阴影和网格的图形对于理解曲面的形状是很有好处的。在有些矢量图形的输出装置中,你可能得不到阴影,但是当有阴影时,输出装置可能要花很长时间来输出它。(8).给空间立体曲面着色通常情况下,Mathematica 为了使图形更加逼真而用明暗分布的形式给空间立体曲面着色。在这种情况下,Mathematica 假定在图形的右上方有三种光源照在物体上。但有时这种方法会造成混乱,此时你可用 Lighting-&False 来采取根据高度在表面上涂以不同灰度的阴影的方法。2 用数据来进行绘图同二维绘图一样,三维图形也可用数据来进行绘图。下面给出数据矩阵,因其较大未表示其结果。3.三维空间的参数方程绘图三维空间中的参数绘图函数 ParametricPlot3D[{fx,fv,fz},{t,tmin,tmax}]和二维空间中的 ParametricPlot很相仿。在这种情况下,Mathematica 实际上都是根据参数 t 来产生系列胡点,然后再连接起来。三维参数作图的基本形式为,ParametricPlot3D[{fx,fv,fz},{t,tmin,tmax}] 给出空间曲线的参数图ParametricPlot3D[{fx,fv,fz},{t,tmin,tmax},{u,umin,umax}] 给出空间曲面的参数图ParametricPlot3D[{fx,fv,fz,s}…..]按照函数关系 s 绘出参数图的阴影部分ParametricPlot3D[{fx,fv,fz},{gx,gy,gz}…..]] 把一些图形绘制在一起下面是一些空间曲线的例子结果为命令 ParametricPlot3D[{fx,fv,fz},{t,tmin,tmax},{u,umin,umax}] 产生一个见曲面而不是一条曲线曲面是由四边形组成。下面这个图形也很漂亮5.1 极限Mathematica 计算极限的命令是 Limit 它的使用方法主要有Limit[expr,x-&x0] 当 x 趋向于 x0 时求 expr 的极限Limit[expr,x-&x0,Direction-&1] 当 x 趋向于 x0 时求 expr的左极限Limit[expr,x-&x0,Direction-&-1] 当 x 趋向于 x0 时求 expr的右极限趋向的点可以是常数,也可以是 +∞,-∞ 例如1.求2.求3.求5.2 微分1.函数的微分在 Mathematica 中,计算函数的微分或是非常方便的,命令为 D[f,x],表示对 x 求函数 f 的导数或偏导数。该函数的常用格式有以下几种D[f,x]计算微分D[f,x1,x2,…]计算多重偏微分D[f,{x,n}]计算 n 阶微分D[f,x,NonConstants-&{v1,v2,----}]计算微分 其中 v1,v2…依赖于 x例如1.求函数 sinx的导数2.求函数 exsinx的 2 阶导数3.假设 a是常数可以对 sinax求导4.如果对二元函数 f(x,y)=x^2*y+y^2 求对 x,y 求一阶和二阶偏导Mathematica可以求函数式未知的函数微分,通常结果使用数学上的表示法例如,对链导法则同样可用如果要得到函数在某一点的导数值可以把这点代入导数如:2.全微分在 Mathematica 中,D[f,x]给出 f 的偏导数,其中假定 f 中的其他变量与 x 无关。当 f 为单变量时,D[f,x]计算 f 对 x 的导数。函数 Dt[f,x]给出 f 的全微分形式,并假定 f 中所有变量依赖于 x.下面是 Dt 命令的常用形及意义Dt[f] 求全微分 dfDt[f,x]求全微分Dt[f,x1,x2,…]求多重全微分Dt[f,x,Constants-&{c1,c2,….}] 求全微分其中 c1,c2..是常数下面我们求 x^2+y^2 的偏微分和全微分可以看出第一种情况 y 与 x 没有关系,第二种情况 y 是 x 的函数。再看下列求多项式 x^2+xy^3+yz 的全微分并假定 z 保持不变是常数。如果 y 是 x 的函数哪么,y 被看成是常数5.3 计算积分1.不定积分在 Mathematica 中计算不 定积分命令为 Integerate[f,x],当然也可使用工具栏直接输入不定积分式。来求函数的不定积分。当然并不是所有的不定积分都能求出来。例如若求 Mathematica 就无能为力。但对于一些手工计算相当复杂的不定积分,MatheMatica 还是能轻易求得,例如求积分变量的形式也可以是一函数,例如输入命令也可求得正确结果。对于在函数中出现的除积分变量外的函数,统统当作常数处理,请看下面例子。2.定积分定积分的求解主要命令也是用 Integrate 只是要在命令中加入积分限 Integrate[f,{x,min,max}]或者使用式具栏输入也可以。例如求显然这条命令也可以求广义积分例如:求求无穷积也可以例如如果广义积发散也能给出结果,例如如果无法判定敛散性,就用给出一个提示,例如如果广义积分敛散性与某个符号的取值有关,它也能给出在不同情况下的积分结果例如结果的意义是当|p|&1 时,积分值为 1/1-p,否则不收敛。在 Integrate 中可加两个参数 Assumptions 和GenerateConditions 例如上例中,只要用 Assumptions-&{Re[p]&1}就可以得到收敛情况的解3.数值积分数值积分是解决求定积分的另一种有效的方法,它可以给出一个近似解。特别是对于用 Integrate 命令无法求出的定积分,数值积分更是可以发挥巨大作用。它的命令格式为Nintegrate[f,{x,a,b}] 在[a,b] 上求 f 数值积分Nintegrate[f,{x,a,x1,x2,…,b}] 以 x1,x2….为分割求[a,b] 上的数值积分Nintegrate[f,{x,a,b},MaxRecursion-&n] 求数值积分时指定迭代次数 n,下面我们求 Sinsinx 在[0,Pi] 上的积分值,由于这个函数的不定积分求不出,因此使用 Integrate 命令无法得到具体结果,但可以用数值积分求如果积分函数存在不连续点,或存在奇点我们可对积分进行分段求解。例如函数 在[-1,1]上,显然 x=0 点是一个无穷间断点。因此若要求其数值积分,必须在其中插入点 0对无穷积分,也可求数值积分,例如。多变量函数的微分下面是计算多变量函数的偏导数及全微分的命令与单变量基本相同,通过分析下面的例子我们可以我们可以轻松掌握。( I ) 计算偏导数下面是实际的例子,( II )中依赖于 x,下面是实际的例子,( III ) 计算全微分 df,下面是实际的例子,( IV ) 求隐函数的导数下面是实际的例子,( V ) 计算全微分 df,其中是常数,下面是实际的例子,( VI ) 计算多重全微分,下面是实际的例子,多变量函数的积分 ( 重积分 )多变量函数的积分类似于一元函数的积分,可以利用 Integrate 函数来完成,命令如下,计算重积分,数值积分或重积分的数值解,下面是具体的例子,以下是数值积分的例子,6,1 微分方程解在 Mathematica 中使用 Dsolove[]可以求解线性和非线性微分方程,以及联立的微分分方程组。在没有给定方程的初值条件下,我们所得到的解包括 C[1],C[2]是待定系数。求解微分方程就是寻找未知的函数的表达式,在 Mathematica 中,未稳中有降函数用 y[x]表示,其微分用 y'[x],y''[x]等表示。下面给 出微分方程(组)的求解函数。Dsolve[eqn,y[x],x] 求解微分方程 y[x]Dsolve[eqn,y,x] 求解微分方程函数 yDsolve[{eqn1,eqn2,…},{y1,y2,….},x] 求解微分方程组1.用 Dsolve 求解微分方程 y[x]解 y[x]仅适合其本身,并不适合于 y[x]的其它形式,如 y’[x],y[0]等,也就是说 y[x]不是函数,例如我们如果有如下操作,y’[x],y[0]并没有发生变化,2.解的纯函数形式使用 Dsolve 命令可以给出解的纯函数形式,即 y,请分析下面的例子这里 y 适合 y 的所有情况下面的例子可以说明这一点在标准数学表达式中,直接引入亚变量表示函数自变量,用此方法可以生成微分方程的解。如果需要的只是解的符号形式,引入这样来变量很方便。然而,如果想在其他的的计算中使用该结果,那么最好使用不带亚变量的纯函数形式的结果。3.求微分方程组请分析下面的例子当然微分方程组也有纯函数形式。4.带初始条件的微分方程的解当给定一个微分方程的初始条件可以确定一个待定系数。请看下面的例子第二个例子由于给出一个初始条件所以只能确定 C[1],5.进一步讨论对于简单的微分方程的解比较简单,对一些微分方程它的解就复杂的多。特别是对一些微分方程组或高阶微分方程,不一定能得具体的解,其解中可能含有一些特殊函数。并且很多特殊函数的提出就是为了解这些方程的如,上面三个方程中分别使用了三种类型的函数,可以查看系统帮助了解他们的性质和含义。对于非线性微分方程,仅有一些特殊的情况可用标准数学函数得到解。 Dsolve 能够处理所有在标准数学手册有解的非线性微分方程。例如,可以看出第二个方程的解已经非常复杂。6.2 微分方程的数值解在 Mathematica 中用函数 DSolve[]得到微分方程的准确解,用函数 NDSolve 得到微分方程的数值解,当然在此处要给出求解区间(x,xmin,xmax )。NDSolve 也是既能计算单个的微分方程,也能计算联立微分方程组。它能对大多数的常微分方程和部分偏微分方程求解。在常微分可能有一些未知函数 yi,但这些未知函数都依赖于一个单变量 x。NDSolve[{eqn1,eqn2,…},y,{x,xmin,xmax}]求函数 y 的数值解,x 属于[xmin,xmax]NDSolve[{eqnl,eqn2,…},{y1,y2,…}{x,xmin,xmax}]求多个函数yi的数值解NDSolve以InterpolatingFunction目标生成函数yi的解,InterpolatingFunction目标提供在独立变量x的xmin到xmax范围内求少的近似值。NDSolve用迭代法求解,它以某一个x值开始,尽可能覆盖从xmin到xmax的全区间。为使迭代开始,NDSolve指定yi及其导数为初始条件。初始条件给定某定点x处的yi[x]及尽可能的导数y'i[x],一般情况下,初始条件可在任意x处,NDSolve将以此为起点自动覆盖xmin到xmax的全区域。下面对初始条件y[0]=0和y[1]=0分别求出x从0到1的范围内y’[x]=y[x]的解。再看下面的微分方程的数值解使用Mathematica页可以很容易的得到解的图形。这儿给出如何观察微商的逆函数的近似值图形。我们使用命令Evaluate代替InterpolatingFunction能够节省时间。例如:7.1 模块和块中的变量前面我们学习了有关 Mathematica 的各种基本运算及操作,为了使 Mathematica 更有效的工作,我们可对 Mathematica 进行模块化运算。在模块内部通过编写一系列表达式语句,使其实现一定的功能。在Mathematica 内部也提供了很多程序包,我们将学习如何调用它们。一般情况下,Mathematica 假设所有变量都为全局变量。也就是说无论何时你使用一个你定义的变量,Mathematica 都假设你指的是同一个目标。然而在编制程序时,你则不会想把所有的变量当作全局变量,因为如果这样程序可能就不具有通用性,你也可能在调用程序时陷入混乱状态。给出定义模块或块和局部变量的常用形式,Module[{x,y,...},body] 具有局部变量 x,y…的模块Module[{x=x0,y=y0,…},body] 具有初始值的局部变量的模块lhs:=Module[vars,rhs/:cond] rhs 和 cond 共享局部变量Block[{x,y,..,},body]运用局部值 x,y,…计算 bodyBlock[{x=x0,y=y0,…},bddy] 给 x,y,..赋初始值Mathematica 中的模块工作很简单,每当使用模块时,就产生一个新的符号来表示它的每一个局部变量。产生的新符号具有唯一的名字,互不冲突,有效的保护了模块内外的每个变量的作用范围。首先我们来看 Module 函数,这个函数的第一部分参数,里说明的 变量只在 Module 内起作用,body 执行体,包含合法的 Mathematica 语句,多个语句之间可用,;”分割下面定义有初值的变量 t,Mathematica 默认它为全局变量,ln[1]:=t=lOOut[1]=10模块中的 t 为局部变量,因此它独立于全局变量 tln[2]:=Module[{t},t=8;Print[t]]全局变量 t 的值仍为 10,ln[3]=t=lOOut[3]=10下面定义函数中的中间变量 t 为局部变量并调用 f,全局变量 t 的值仍为 10,ln[6]:=t=10Out[6]=10我们可以对模块中的任意局部变量进行初始化,这些初始值总是在模块执行前就被计算出来。下面给局部变量 t 赋初值 u:调用函数Mathematica 中的模块允许你把某变量名看作局部变量名。然而又存在有时你又希望它们为全局变量时,但变量值为局部的矛盾,这时我们可以用 Block[]函数。下面是一个含有全局变量 x 表达式,使用 x 的局部值计算上面的表达式,在 Mathematica 中编制程序时,必须使程序中的各个部分尽可能的独立,这样程序才便于读懂、维护和修改。确保程序各部分不相干的主要方法是设置具有一定作用域的变量。在 Mathematica 中有两种限制变量作用域的基本方法,模块 (Module)和块 (Block)。我们在书写实际程序中,模块比块更具普遍性。然而在交互式计算中需要定义作用域时,块更实用。Module[vars,body]所要做的是把执行模块时表达式 body 的形式看成 Mathematica 程序的,代码” 。然而当,代码” 中直接出现变量 vats 时,这些 vars 都将被看作局部的。 Block[vats,body]并不查看表达式body 的形式,而在整个计算 Body 的过程中,实用 vars 的局部值。下例中我们根据 i 定义 m,ln[12]:=m=i^2 ’Out[12]:=i2在计算 i+m的整个过程中使用块中 i的局部值,h[13]:=BlockI(i; a),!+m]Out[13]=a+a2而对于下面的例子,只有直接出现在 i+m中的 i,才被看作局部变量,In[14]:=Module[{i=a},i+m]Out[14]=a+i‘8.2.1 条件结构我们在用计算机语言进行编程时,常用到条件语句。在 Mathematica 中也提供了多种设置条件的方法,并规定只有在该条件满足时才计算表达式。下面条件结构的常用形式。lhs:=rhsl/:test 当 test 为真时使用定义If[test,then,else] 如 test 为真计算 then,反之计算 elsewhich[test1,value1,test2,...] 依次计算 testl,给出对应的第一个为真的值Switch[expr,forml,value1,form2,...]expr 与每一个 formi 相比较,给出第一个相匹配的值Switch[expr,form1,value1,form2,…,_,def] 用 def 为系统默认值1.If 命令下面的 test 为真,故返回第一表达式的值,ln[1]:=If [1&0,1+2,2+3]Out[1]=3用 Mathematica 编程时,不可避免的要在单个或多个定义之间进行选择。单个定义的右边包含多个由 If函数控制的分支,多个定义是用 /;condition 来表示的。运用多个定义进行编程你常能得到结构很好的程序。下面定义了一个阶跃函数,即当 x&O 时值为 1,反之值为 -1,ln[2]:=If[x&0,1,-1]下面运用 /;condition 形式分别定义阶跃函数的正数和负数部分,ln[3]:g=1/:x&0ln[4]:g=-1/:x&0用,?”显示用 /:condition 定义的函数 g 的完整信息,ln[5]:=?gGlobal`gg[x_]:=1/:x&Og[x_]:=-1/,x&O我们用函数 If 时,还可以用 if(test,expr)结构,即当 test 真时,计算表达式 expr,表达式 expr 的值就是整个If 结构的值,反之返回空值。2,Which 命令对于一般情况函数 If 提供一个两者择一的方法。然而,有时条件多于两个,在这种情况下可用 If 函数的嵌套方式来处理,但在这种情况下使用 Whitch 或 Switch 函数将更合适。下面用 Which 定义具有三个条件的函数,调用这个函数,用 Switch 定义一个与模的余数有关的函数,Mod[17,3]=2,因此运用了 Switch 中的第三种情况,3 符号条件在 Mathemahca 中,有一种可能的情况就是你给出的条件结果既不是真也不为假。下面测试的结果既不是真也不是假,因此 If 的两个分支保持不变,ln[1]:=If[x==y,a,b]Out[1]:If[x==y,a,b]你可以给 If 加上第三个条件结果,这允许你测试的结果既不是真也不是假的情况下使用它,ln[2]:=If[x==y,a,b,c]Out[2]=c下面给出处理符号条件的函数。因 x=x,故 Mathematica 给出结果为真,但 Mathematica 在下面情况下以符号等式输出,ln[4]:=x==yOut[4]:=x==y除非表达式能得出真,否则都被假设为假,ln[5]:=TrueQ[x==x]Out[5]=Trueln[6]:=TrueQ[x==y]Out[6]=false我们用,===”可直接测试两个表达式的等同性,In[7]:x===yOut[7]:=False一般情况下,“===”返回值为真 (Truee)或假 (False),而,==”为符号形式输出,表示一个符号等式。在特殊情况下可用,===”测试一个表达式的结构,而用,==”测试数学上的等同性。下例用,===”来测试表达式的结构,下面给出一个无用结果,在建立条件时,你常需要运用组合条件,如 testl&&test2&&…。对于这些组合条件,如果其中有一个为假,则最后结果为假。,Mathematica 依次对 test 进行计算,直到其中有一个为假为止。4.是逻辑表达式的运算形式。逻辑表达式expr1&&expr2&&expr3 计算 expri,直到其中有一个为假为止exprl||expr2||expr3 计算 expri,直到其中有一个为真为止下面的函数包括两个组合条件,In[l0]:=t[x_]:=(x!=0&&1/x&3)对这两个测试条件进行计算,下面的第一次测试得出为假,因此不进行第二个条件的测试,第二测试结果可能为 1 或 0,因此输出结果为假,ln[12]:=t[0]Out[12]=FalseMathematica 处理逻辑表达式的方法允许你组合一系列的测试条件,且只有当前面条件满足时才处理后面的条件。8.3 循环结构Mathematica 程序的执行包括对一系列 Mathematica 表达式的计算。对简单程序,表达式的计算可用分号,;” 来隔开,然后一个接一个地进行计算。然而,有时你需要对同一表达式进行多次计算,即循环计算。1 Do 循环结构简单地 Do 循环结构形式,Do[expr,{i,imax}] 循环计算 expr,以步长 1,i 从 1 增加到 imaxDo[expr,{i,imin,imax,di}] 循环计算 expr,以步长 di,i 从 imin 增加到 imaxDo[expr,{n}] 循环计算 expr n 次计算 Print[i+i^2],i 从 1 增加到 3,Do 中的定义的循环方式与函数 Table 和 Sum 中的定义一样。在函数 Do 中,你同样能建立重循环。下面给出的 i 从 1 到 4 进行循环,而对于每个 i,j 又从 1 到 i-1 进行循环,我们还可把一个过程放入 Do 函数中,2 While 与 For 结构在 Mathematica 程序中,Do 是以结构方式进行循环的,然而有时你需要生成非结构循环。此时,运用函数 While 和 For 是合适的。下面是 While 和 For 函数的循环结构形式,While[test,body] 只要 test 为真,就重复计算 bodyFor[start,test,incr,body]以为 start 起始值,重复计算 body 和 incr,直到 test 为假为止当条件满足时,While 循环一直进行,因此为了防止死循环,在 While 中应包括命令能改变 test 的值。下面给出 For 循环的例子,i++表示 i 的值加 1(在本节的最后我们给出在编程时常会 用到的赋值方法 ),下面再给出一个较复杂的 For 循环的例子,一旦 I^2&lO 不成立,就中止循环,Mathematica 中的函数 While 和 For 循环总是在执行循环体前对循环条件进行测试。 一旦测试结果为假。就中止 While 和 For 循环。因此,循环体的计算总是在测试结果为真 的情况下进行的3.一些特殊的赋值方式一些赋值方式在循环结构中有时能带来一些方便。i++ 变量 i 加 1i-- 变量 i 减 1++i,变量 i 先加 1--i 变量 i 先减 1i+=di i 加 dii-=di i 减 dix*=C x 乘以 Cx/ =c x 除以 c{x,y}={y,x} 交换 x 和 y 值4 重复运用函数我们除了可用 Do,While,For 等进行循环计算外,我们还可以运用函数进行编程。运用函数编程结构你能得出非常有效的程序。例如 Nest[f,x,n]允许你对某一表达式重复运用函数 f给出重复运用函数的方式。Nest[expr,n] 对表达式 expr 重复调用函数 fn 次FixedPoint[y,expr] 对表达式 expr 重复调用函数 fn 次,直到结果不变为止NestWhile[f,expr,test]对表达式 expr 重复调用函数 f,直到产生的结果为假时为止下面对函数 f 迭代 5 次,对纯函数进行迭代,你能得出与运用 Do 函数得出的结果一样,Nest 函数允许你重复运用某函数。然而,有时你想在结果不再发生变化的情况下就中止对函数的使用,此时立刻使用函数 FixPoint[f,x]。 FixPoint 函数重复运用某一函数直到结果不再发生变化,8,3.5 流程控制函数程序结构的流程控制一般来说比较简单,但是在应用 While 或 For 等循环时就比较复杂了,这是因为他们的流程控制依赖于表达式的值。而且在这样的循环中,流程的控 制并不依赖于循环体中表达式的值。有时你在编制 Mathematica 程序时,在该程序中,流程控制受某一过程或循环体执行结果的影响。这时,我们可用 Ma~ematka 提供的流程控制函数来控制流程。这些函数的 工作过程与 C语言中的很相似。常用的流程控制函数。Break[] 退出本层的循环Continue[] 转入当前循环的下一步Return[expr] 退出函数中的所有过程及循环,并返回 expr 值Goto[name] 转入当前过程中的元素 Label[name]Throw[value] 返回 expr 值当 t&20 时,Break[] 就引起循环体的中断,当 k&3 时,Continue[]继续执行循环,下面给出 Return 的一个例子,Return[]允许你退出一函数,并返回一个值。 Mathematica 可以进行局部返回,这可 允许你退出一列迭代函数。非局部返回在错误处理时是很有用的。 下面给出的例子中如函数变量小于 0 则输出 error~ln[6]不产生 error,且出示 Catch 的结果无效,当变量小于 0 时输出 error运算符及特殊符号 系统常量1.运算符及特殊符号Linel 执行 Line,不显示结果Linel,line2 顺次执行 Line1,Line2,并显示结果name 关于系统变量 name 的信息name 关于系统变量 name 的全部信息!command 执行 Dos 命令N! N 的阶乘!!filename 显示文件内容&&filename 读入文件并执行Expr,&&filename 打开文件写Expr&&&filename 打开文件从文件末写() 结合率[] 函数{} 一个表&*MathFun*& 在 c 语言中使用 math 的函数(*Note*) 程序的注释#n 第 n 今参数## 所有参数Rule& 把 role 作用于后面的式子% 前一次的输出%% 倒数第二次的输出Var::mote 变量 var 的注释“Astring” 字符串Context 上下文A+b 加a-b 减A*b 或 ab 乘A/b 除2.系统常量Pi 3.1415 的无限精度数值E 2.17828 的无限精度数值Catalan 0,915966Catalan 常数EulerGamma 0,5772Euler 常数Khinchin 2,68545KhinchinGlaisher 0,915966GlaisherGoldenRatio 1,61803 黄金分割数Degree Pi/l80 角度弧度换算I 复数单位Infinity 无穷大-Infinity 负无穷大Complexlnfinity 复无穷大Indeterminate 不定式代数计算Expand[expr] 展开表达式Factor[expr] 展开表达式Simplify[expr] 化简表达式FullSimplify[expr] 将特殊函数也进行化简PowerExpand[expr] 展开所有的幂次形式ComplexExpand[expr,{x1,x2…)] 按复数实部虚部展开FunctionExpand[expr] 化简表达式中的特殊函数Collect[expr,x] 合并同次项Collect[expr,{x1,x2,…)] 合并 x1,x2,..,的同次项Together[expr] 通分Apart[expr] 部分分式展开Apart[expr,var) 对 var 的部分分式展开Cancel[expr] 约分xpandAll[expr] 展开表达式ExpandAll[expr,patt] 展开表达式FactorTermsrpoly] 提出共有的数字因子FactorTerms[poly,x] 提出与 x 无关的数字因子FactorTerms[poly,(x1,x2…)] 提出与 xi 无关的数字因子Coefficient[expr,form] 多项式 expr 中 form 的系数Coefficient[expr,form,n) 多项式 expr 中 form^n 的系数Exponent[expr,form] 表达式 expr 中 form 的最高指数Numerator[expr] 表达式 expr 的分子Denominator[expr] 表达式 expr 的分母ExpandNumerator[expr] 展开 expr 的分子部分Solve[eqns,vats] 从方程组 eqns 中解出 VatsSolve[eqns,vats,elims] 从方程组 eqns 中削去变量 elims,解出 vatsDSolve[eqn,y,x] 解微分方程,其中,y 是 x 的函数DSolve[{eqnl,eqn2,…},{y1,y2…},] 解微分方程组,其中 yi 是 x 的函数DSolve[eqn,y,{x1,x2…}] 解偏微分方程Eliminate[eqns,Vats] 把方程组 eqns 中变量 vars 约去SolveAlways[eqns,vars] 给出等式成立的所有参数满足的条件Reduce[eqns,Vats] 化简并给出所有可能解的条件LogicalExpand[expr] 用 &&和,,将逻辑表达式展开InverseFunctionI 刀 求函数 f 的反函数Root[f,k1 求多项式函数的第 k 个根Roots[1hs==rhs,var] 得到多项式方程的所有根微积分D[f,x] 求 f[x]的微分D[f,{x,n}] 求 f[x]的 n 阶微分D[f,x1,x2…] 求 f[x]x1,x2…偏微分Dt[f,x] 求 f[x]的全微分 df/dxDt(f) 求 f[x]的全微分 dfDt[f,{x,n}] n 阶全微分 df^n/dx^nDt[f,x1,x2..] 对 x1,x2..的偏微分Integrate[f,x] f[x]对 x 在的不定积分Integrate[f,{x,xmin,xmax}] f[x]对 x 在区间 (xmin,xmax)的定积分Integrate[f,{x,xmin,xmax},{y,ymin,ymax}] f[x,y]的二重积分Limit[expr,x-&x0] x 趋近于 x0 时 expr 的极限Residue[expr,{x,x0}] expr 在 x0 处的留数Series[f,{x,x0,n)] 给出 f[x]在 x0 处的幂级数展开Series[f,{x,xO,nx},{y,y0,ny}] 先对 y 幂级数展开,再对 x 幂级数展开Normal[expr] 化简并给出最常见的表达式SeriesCoefficient[series,nJ 给出级数中第 n 次项的系数SeriesCoefficient[series,{n1,n2…}] 一阶导数InverseSeries[s,x] 给出逆函数的级数ComposeSeries[seriel,serie2…] 给出两个基数的组合SeriesData[x,x0,{a0,a1,..},nmin,nmax,den] 表示一个 x0 处 x 的幂级数O[x]^n n 阶小量 x^n多项式函数Variableslpoly] 给出多项式 poly 中独立变量的列表CoefficientList[poly,var] 给出多项式 poly 中变量 var 的系数CoefficientList[poly,{varl,var2…)) 给出多项式 poly 中变量 var(i)的系数列PolynomialMod[poly,m] poly 中各系数 mod m 同余后得到的多项式,m 可为整式PolynomialQuotient[p,q,x] 以 x 为自变量的两个多项式之商式 p/PolynomialRemainder[p,q,x] 以 x 为自变量的两个多项式之余式P01ynomialGCD[polyl,poly2,...] poly(i)的最大公因式PolynomialLCM[polyl,poly2...,] poly(i)的最小公倍式PolynomialReduce[pjoly,{polyl,Poly2,..},{xl,x2…})得到一个表 I(a1,a2,·),b) 其中 Sumld*polyi]+b=polyResultant[polyl,poly2,var] 约去 polyl,poly2 中的 varFactor[poly] 因式分解 (在整式范围内 )FactorTerms[poly] 提出 poly 中的数字公因子FactorTermslpoly,{x1,x2…)) 提出 poly 中与 xi 无关项的数字公因子FactorList(poly],FactorSquareFreeList[p01y],FactorTermsList[poly,{x1,x2…}]给出各个因式列表Cyclotomic[n,x] n 阶柱函数Decomposet[poly,x] 迭代分解,给出 {p1,p2,...),其中 P1(p2(…))=polyInterpolafinSPolynomial[data,Var] 在数据 data 上的插值多项式RootSum[f,form] 得到 f[x]=0 的所有根,并求得 Sum[form[xi]]随机函数RandomCtype,range] 产生 type 类型且在 range 范围内的均匀分布随机数Random[] 0-1 上的随机实数SeedRandom[n1 以 n 为 seed 产生伪随机数Randomldistribution] 可以产生各种分布数值函数N[expr] 表达式的机器精度近似值N[expr,n) 表达式的 n 位近似值,n 为任意正整数NSolve[1hs=---rhs,val] 求方程数值解NSolveleqn,Var,n1 求方程数值解,结果精度到 n 位NDSolve[eqns,y,{x,xmirl,xmax}] 微分方程数值解NDSolve[eqns,{y1,y2,·,1,{x,xmin,xmax}} 微分方程组数值解FindRoot[1hs==rhs,{x,x0)1 以 x0 为初值,寻找方程数值解FindRoot[1hs=--rhs,{x,xstart,xmin,xmax}] 以 xstart 为初值,在[xmin,xmax] 范围内寻找方程数值解NSum[f,{imin,imax,di}] 数值求和,出为步长NSum[f,{imin,imax,di},{j,..},..] 多维函数求和NProduct[f,{i,imin,imax,di}] 函数求积NIntegrate[f,{x,xmin,xmax}] 函数数值积分FindMinimum[f,{x,xO}] 以 x0 为初值,寻找函数最小值FindMinimum[f,{x,xstart,xmin,xmax}] 以 xstart 为初值,在[xmin,xmax] 范围内寻找方程解ConstrainedMin[f,{inequ},{x,y,..}]inequ 为线性不等式组,f 为 x,y,.之线性函数,得到最小值及此时的 x,y..取值ConstrainedMax[f,{inequ),{x,y,..}] 同上LinearProgramming[C,m,b]解线性组合 c.x 在 m,x&=b&&x&=0 约束下的最小值,x,为向量,m 为矩阵LatticeReduce[{v1,v2...}] 向量组 Vi 的极小无关组Fit[data,funs,vats] 用指定函数组对数据进行最小二乘拟合Interpolation[data] 对数据进行插值Lisfinterpolation[array] 对离散数据插值,array 可为 n 维ListInterpolafion[array,{{xmin,xmax},{min,ymax},..}] 在特定网格上进行插值FunctionInterpolation[expr,{x,xmin,xmax},{y,ymin,ymax},..] 以对应 expr[xi,yi]的数值为数据进行插值Fourier[list] 对复数数据进行傅氏变换InverseFourieroiSt) 对复数数据进行傅氏逆变换制表函数 元素操作 表的操作1.制表函数{e1,e2,...} 一个表,元素可以为任意表达式,无穷嵌套Table[expr,{imax}] 生成一个表,共 imax 个元素Table[expr,{i,imax}} 生成一个表,共 imax 个元素 expr 间Table(expr,{i,imin,imax},{j,jmin,jmax},..] 多维表Range[imax] 简单数表 f1,2+,imax)Range[imin,imax,di] 以 di 为步长的数表Array[f,n] 一维表,元素为 fI¨(i 从 1 到 n)Array[f,{n1,n2..}] 多维表,元素为玎 i 小.1 ( 各自从 1 到 ni)IdentityMatrix[n] n 阶单位阵DiagonalMatrix[list] 对角阵2、元素操作Part[expr,i]或 expr[[i]] 第 i 个元素expr[[-i]] 倒数第 i 个元素expr[{i,j,..}] 多维表的元素expr[{i1,i2,..}] 返回由第 i(n)的元素组成的子表FirstCexpr] 第一个元素Last[expr] 最后一个元素Head[expr] 函数头,等于 expr[[0]]Extract[expr,list] 取出由表 list 指定位置上 expr 的元素值Take[list,n] 取出表 list 前 n 个元素组成的表Take[list,{m,n}] 取出表 list 从 m 到 n 的元素组成的表Drop[list,n] 去掉表 list 前 n 个元素组下的表Rest[expr] 去掉表 list 第一个元素剩下的表Select[USt,crit] 把 crit 作用到每一个 list 的元素上,为 True 的所有元素组成的表Length[expr] expr 第一层元素的个数Dimensions[expr] 表的维数返回 (n1,n2..},expr 为一个 nl*n2…的阵TensorRank[expr) 秩Depth[expr] expr 最大深度Level[expr,n] 给出 expr 中第 n 层子表达式的列表Count[USt,paUem] 满足模式的 list 中元素的个数MembefQ[1ist,form] list 中是否有匹配 form 的元素FreeQ[expr,form] MemberQ 的反函数FreeQ[expr,form] 表中匹配模式 pattern 的元素的位置列表Cases[{e1,e2…},pattem} 匹配模式 pattem 的所有元素 ei 的表3.表的操作Append[exp[elem] 返回在表 expr 的最后追加 elem 元素后的表Prepend[expr,elem) 返回在表 expr 的最前添加 elem 元素后的表Insert[1ist,elem,n] 在第 n 元素前插入 elemlnsert[expr,elem,{i,j,...}] 在元素 expr[[{i,j,..}]]前插入 elemDelete[expr,{i,j,..}] 删除元素 expr[[{i,j,..}]]后剩下的表DeleteCases[expr,pattem] 删除匹配 pattern 的所有元素后剩下的表ReplacePart[expr,new,n] 将 expr 的第 n 元素替换为 newSort[list] 返回 list 按顺序排列的表Reverse[expr] 把表 expr 倒过来RotateLeft[expr,n] 把表 expr 循环左移 n 次RotateRight[expr,n] 把表 expr 循环右移 n 次Partition[list,n] 把 list 按每 n 个元素为一个子表分割后再组成的大表Flatten[1istl 抹平所有子表后得到的一维大表Flatten[1ist,n] 抹平到第 n 层Split[1ist] 把相同的元素组成一个子表,再合成的大表二维绘图 二维绘图设置三维绘图 三维绘图设置 密度图 图形显示 图元函数 着色及其他1.二维绘图Plot[f,{x,xmin,xmax}] 一维函数 f[x]在区间[xmin,xmax] 上的函数曲线Plot[{fl,f2..},{x,xmin,xmax}] 在同一图形上画几条曲线ListPlot[{y1,y2,..}] 绘出由离散点对 (n,yn)组成的图ListPlot[{{x1,y1},{x2,y2},}} 绘出由离散点对 (xrl,yrl)组成的图ParametricPlot[{fx,fy},{t,tmin,tmax}] 由参数方程在参数变化范围内产生的曲线ParametricPlot[{fx,fy},{gx,gy},...],{t,tmin,truax}]2.二维设置PlotRange-&{0,1} 作图显示的值域范围AspectRatio-&1/GoldenRatio 生成图形的纵横比PlotLabel-&label 标题文字Axes-&{false,True} 分别制定是否画 x,y 轴AxesLabel-&{xlabel,ylabel} x,y 轴上的说明文字Ticks-&None,Automatic,fun 用什么方式画轴的刻度AxesOrigin-&{x,y} 坐标轴原点位置AxesStyle-&{{xstyle},{ystyle}} 设置轴线的线性颜色等属性Frame-&True,False 是否画边框FrameLabel-&{xmlabel,ymlabel,xplabel,yplabel} 边框四边上的文字FrameTicks 同 Ticks 边框上是否画刻度GridLines 同 Ticks 图上是否画栅格线Framestyle-&{{xmstyle},{ymstyle}} 设置边框线的线性颜色等属性ListPlot[data,PlotJoined-&True] 把离散点按顺序连线Plotsytle-&{{style1},{style2},..} 曲线的线性颜色等属性PlotPoints-&15 曲线取样点,越大越细致3.三维绘图Plot3D[f,{x,xmin,xmax},{y,ymin,ymax}] 二维函数 flx,y]的空间曲面Plot3D[{f,s},{x,xmin,xmax},{y,ymin,ymax}] 同上,曲面的染色由 slx,y1 值决定ListPlot3D[array] 二维数据阵 array 的立体高度图ListPlot3D[array,shades] 同上,曲面的染色由 shades[数据 ]值决定ParametricPlot3D[{fx,fy,fz},{t,tmin,tmax}] 三维参数图形ContourPlot[f,{x,xmin,xmax},{y,ymin,ymax}] 二维函数 flx,刀在指定区间上的等高线图ListContourPlot[array] 二维函数 fix,yJ 在指定区间上的等高线图4.三维设置Contours-&n 画 n 条等高线Contours-&{z1,z2,...} 在 zi 处画等高线ContourShading-&False 是否用深浅染色ContourLines-&True 是否画等高线ContourStyle-&{{stylel},{style2},..} 等高线线性颜色等属性5.密度图DensityPlot[f,{x,xmin,xmax},{y,ymin,ymax)] 二维函数 f[x,y]在指定区间上的密度图ListDensityPlot[array] 二维函数 f[x,y]在指定区间上的密度图5.图形显示Show[graphics,options] 显示一组图形对象,options 为选项设置Show[g1,g2…] 在一个图上叠加显示一组图形对象GraphicsArray[{g1,g2,...}] 在一个图上分块显示一组图形对象SelectionAnimate[notebook,t] 把选中的 notebook 中的图画循环放映6.图元函数Graphics[prim,options] prim 为下面各种函数组成的表,表示一个二维图形对象Graphics3D[prim,options] prim 为下面各种函数组成的表,表示一个三维图形对象SurfaceGraphics[array,shades] 表示一个由 array 和 shade 决定的曲面对象ContourGraphics[array] 表示一个由 array 决定的等高线图对象DensityGraphics[array] 表示一个由 array 决定的密度图对象Point[p] p={x,y}或 {x,y,2},在指定位置画点Line[{p1,p2,,.}] 经由 Pi 点连线Rectangle[{xmin,ymin),{xmax,ymax}] 画矩形Cuboid[{xmin,ymin,zmin},{xmaxffmax,zmax}] 由对角线指定的长方体Polygon[{p1,p2,,..}] 封闭多边形Circle[{x,y},r] 画圆Circle[{x,y},{rx,ry}] 画椭圆,rx,ry 为半长短轴Circle[{x,y},r,{a1,a2}] 从角度 al-a2 的圆弧Disk[{x,y},r] 填充的园、椭圆、圆弧等参数同上Raster[array,ColorFunction-&f] 颜色栅格Text[expr,coords] 在坐标 coords 上输出表达式PostScrip[&string&] 直接用 Postscript 图元语言写Scaled[{x,y...}] 返回点的坐标,且均大于 0 小于 17.着色及其他GrayLevel[level] 灰度 level 为 0~1 间的实数RGBColor[red,green,blue] RGB 颜色,均为 0~I 间的实数Hue[h,s,b] 亮度,饱和度等,均为 0~1 间的实数CMYKColor[cyan,magenta,yellow,block] CMYK 颜色Thicknessr[r] 设置线宽为 rPointSize[d] 设置绘点的大小Dashing[{r1,r2,...}] 画一个单元的间隔长度的虚线ImageSize-&{x,y} 显示图形大小 (单位为像素 )流程控制If[condition,t,f] 如果 condition 为 True,执行 t,否则执行 f 段if[condition,t,f,u]如果 condition 为 Tme,执行 t,为 False 执行 f,既非 True又非 False,则执行 u 段Which[test1,blockl,test2,block2..] 执行第一为 True 的 tesfi 对应的 blockiSwitch[expr,forml,blockl,form2,block2...] 重复执行 expr imax 次Do[expr,{imax}] 重复执行 expr imax 次Do[expr,{i,imin,imax},{j,jmin,jmax}] 多重循环While[test,body] 循环执行 body 直到 test 为 FalseFor[start,test,incr,body] 循环执行 body 直到 test 为 FalseThrow[value] 停止计算,把 value 返回给最近一个 Catch 处理Throw[value,tag] 停止计算,把 value 返回给最近一个 Catch 处理Catch[exprl 计算 expr,遇到 Throw 返回的值则停止Catch[expr,form] 当 Throw[value,tag]中 Tag 匹配 form 时停止Return[expr) 从函数返回,返回值为 exprReturn[] 返回值 NullBreakl[] 结束最近的一重循环Continuel[] 停止本次循环,进行下一次循环Goto[tag] 无条件转向 Label[Tag]处Label[tag] 设置一个断点Check[expr,fmlexpr] 计算 expr,如果有出错信息产生,则返回 failexpr 的值Check[expr,failexpr,s1::t1,s2::t2,...] 当特定信息产生时则返回 failexprCheckAbort[expr,failexpr] 当产生 abort 信息时返回 failexprInterrupt[] 中断运行Abort[] 中断运行TimeConstrained[expr,t] 计算 expr,当耗时超过 t 秒时终止MemoryConstrained[expr,b] 计算 expr,当耗用内存超过 b 字节时终止运算Print[exprl,expr2,...] 顺次输出 expri 的值Input[] 产生一个输入对话框,返回所输入的任意表达式Input[&prompt&] 同上,prompt 为对话框的提示Pause[n] 运行暂停 n 秒
课件名称:课件分类:数学课件类型:电子教案文件大小:1.88MB下载次数:19评论次数:5用户评分:6.4
1. 2. 3. 4. 5.

我要回帖

更多关于 qwtplot3d画三维曲线 的文章

 

随机推荐