数学分析实验-Mathematica_软件的应用(高等数学实验课讲义)
Mathematica软件在高等数学教学中的应用
Mathematica软件在高等数学教学中的应用【摘要】Mathematica软件是一套专门进行数学计算的软件，具有形象性、直观性、互动性和时效性。
在高等数中引入Mathematica软件进行辅助教学，有利于激发学生的学习兴趣，培养学生基本计算及用数学方法解决实际应用问题的能力。
【关键词】高等数学；Mathematica软件；应用0.引言数学软件Mathematica是处理数学问题的一种应用软件，它的功能非常强大，不仅可以用于符号运算和数值计算，还可以用来方便地绘制一元和二元函数的图形。
运用数学软件作为平台，进行数学实验，把数学软件作为认知工具、计算工具和应用工具，不仅丰富了教学手段，也给传统的数学教学注入了活力。
教师利用数学软件将传统教学中的粉笔加黑板、教师一言堂这种枯燥的教学过程制作成课件，通过计算机呈现给学生，使学生多种感观并用，增强学生的好奇心，吸引学生的注意力，提高对信息的吸收率。
同时，在高等数学教学中适当使用数学软件，能对抽象的数学概念和烦琐的数学运算加以几何解释和简化，通过图形的动态变化使“死”的知识活起来，真正将形与数有机地结合起来，把运动和变化呈现在学生面前，加深对知识的理解，充分调动学生学习的积极性和主动性。
1.Mathematica软件在数学教学中的运用高等数学内容十分丰富，包括微积分、空间解析几何和微分方程等。
高等数学涉及大量的数学计算，Mathematica强大的计算功能，能很好地解决高等数学中的计算问题。
利用Mathematica软件的计算功能，可以提高学生的计算能力。
1.1用Mathematica求极限计算极限命令格式：Limt[函数表达式，自变量→定值]如计算，只要利用Mathematica“基本输入工具栏”，在其“工作窗口”中输入：Limit[Log[1+x]/2x，x->0]，运行软件，即可得到计算结果。
1.2用Mathematica求导数计算导数的命令格式：D[f[x]，x] （求f（x）一阶导数）D[f[x]，{x，n}] （求f（x）n阶导数）D[f[x，y]，x，y] （求f（x，y）的二阶混合偏导数f’’xy（x，y））如已知y=xex，求y’。
mathematica实验报告
mathematica实验报告《使用Mathematica进行实验报告：探索数学的奥秘》Mathematica是一款强大的数学软件，它不仅可以进行数学计算和图形绘制，还可以进行数据分析和模拟实验。
在本实验报告中，我们将使用Mathematica来探索数学的奥秘，展示其强大的功能和应用。
首先，我们将使用Mathematica进行数学计算。
通过输入数学表达式和方程式，我们可以快速地进行数值计算和符号运算。
Mathematica还提供了丰富的数学函数和算法，可以帮助我们解决复杂的数学问题，如微积分、线性代数和离散数学等。
其次，我们将利用Mathematica进行图形绘制。
通过输入函数表达式和参数设置，我们可以绘制出各种数学图形，如函数图像、曲线图和三维图形等。
Mathematica还提供了丰富的绘图工具和选项，可以帮助我们定制和美化图形，使其更加直观和具有艺术感。
接下来，我们将利用Mathematica进行数据分析。
通过输入数据集和统计方法，我们可以进行数据的可视化和分析，帮助我们发现数据的规律和趋势。
Mathematica还提供了丰富的数据处理和建模工具，可以帮助我们进行数据挖掘和预测分析，为决策和规划提供有力的支持。
最后，我们将利用Mathematica进行模拟实验。
通过输入模型和参数设置，我们可以进行各种科学和工程问题的模拟实验，帮助我们理解和预测实际现象。
Mathematica还提供了丰富的模拟工具和仿真方法，可以帮助我们进行虚拟实验和验证假设，为科学研究和工程设计提供有力的工具支持。
总之，Mathematica是一款强大的数学软件，它可以帮助我们探索数学的奥秘，解决数学问题，展示数学图形，分析数学数据，进行数学模拟实验，为科学研究和工程应用提供有力的支持。
希望本实验报告可以激发更多人对数学和科学的兴趣，让我们一起来探索数学的奥秘吧！。
数学软件Mathematica的应用数学实验一例
] 用 >?@ABC?@DE?验证结合律
已 知 在 集 合 ^_ ‘"=$=!=ab上 定 义 了 c 运 算 G’Md5O"K=ef =g h是 否 构 成 半 群i若 按 ’Md5O$规定其c运算=ef =gh是否构成半群i
’Md5O] g "$!a ""$!a $$$aa !!$!a aa$aa
J 用 KLMNOPLMQRL找出所有 S阶么半群
对于么半群)用 (表示单位元T即运算表的第 (行U第 (列均必为 (U5U&V)只需考虑结合 律W执行下述程序 -可得 ((个结果)其中第 (U第 5个表对应的半群同构)第 &U第 (软件 UD@BGVD@6WD的应用 XX 数学实验一例
889
889
8888889
数学软件 >?@ABC?@DE?的应用 FF 数学实验一例
张小红
G宁波大学理学院=浙江 宁波 !"H$""I 西北工业大学计算机科学与工程系=陕西 西安 J"##J$K
摘 要 L 本文介绍作者编写的几个用于研究有限代数系统的 %M;NOPM;Q:M程序=它是数学实验的极好素材< 关 键 词 L 数学实验I%M;NOPM;Q:MI代数系统
程序 J* +,- ..&!/!)!"0!./!&!"!"0!.)!)!&!&0!."!"!/!&001 2345?- &!?7 &’’’’’!?8 8 !
37O
数学的实践与认识
77卷
!"##$%$&&’()*+,-#.*/0!01%23%45&6!"##$%7&&’()*+,-#.*/0!01%23%45&6 !"##$%4&&’()*+,-#.*/0!01%23%45&6!"##7%$&&’()*+,-#.*/0!01%23%45&6 !"##7%7&&’()*+,-#.*/0!01%23%45&6!"##4%$&&’()*+,-#.*/0!01%23%45&6 !"##4%7&&’()*+,-#.*/0!01%23%45&6 !"##4%4&&’()*+,-#.*/0!01%23%45&6 "’ 869,1#:3’ 3%:3; <%:3= = %
数学软件Mathematica的应用
数学软件Mathematica的应用一、数学软件Mathematica简介★Mathematica是由美国Wolfram公司研究开发的一款著名的数学软件；★Mathematica能够完成符号运算、数学图形的绘制等，功能非常强大；★Mathematica能够做精确计算；★Mathematica的界面操作非常友好；★Mathematica是数学建模常用的数学软件之一。
二、利用模板进行微积分运算File（文件）→Palettes（模板）→BasicInput（基本输入）File（文件）→Palettes（模板）→BasicCalculations（基本计算）三、Mathematica中一些常用的函数（1（2（3（（5（6（8）数值分析函数在Mathematica 中，一个逻辑表达式的值有三个：真（True ）、假（False ）和“非真非假”。
条件控制函数If（1） If 语句的结构与一般的程序设计语言中的If 的结构类似。
它有三种情况：If[逻辑表达式，表达式1]当逻辑表达式的值为真时则计算表达式1，表达式1的值就是整个If 结构的值；If[逻辑表达式，表达式1，表达式2]当逻辑表达式的值为真时则计算表达式1，为假时则计算表达式2； If[逻辑表达式，表达式1，表达式2，表达式3]当逻辑表达式的值为真时则计算表达式1，为假时则计算表达式2，其它情况则计算表达式3。
循环控制语句Mathematica 中有3种描述循环的语句，它们是Do,While 和For 语句。
下面是其一般形式：For[初值，条件，修正，循环体] While[条件，循环体] Do[循环体，{循环围}]四、结合图形进行分析1．作出函数xx f y 1sin )(==在区间]1,1[-上的图像，观察当0→x 时函数的变化情况；作出函数xx x f y 1sin)(==在区间]1,1[-上的图像，观察当0→x 时函数的变化情况；2．作出双曲抛物面xy z =的图形； 3．作weierstracs 函数)13cos(21)(1x x f n n nπ∑∞==（处处连续但处处不可导）的图像；4．x ∈(-5,5), y ∈(-5,5)的所有根；五、验证与探索1．x sin 的泰勒级数2．x sin 的无穷乘积猜想六、算法与程序1．分形图（迭代）2．将矩阵化为行最简形（步骤）七、实际问题的Mathematica 求解1．椭圆弧长的计算问题计算椭圆βα≤≤⎩⎨⎧==t t b y ta x ,sin cos 的弧长及近似值。
Mathematica软件介绍及其应用
自然对数的 －Infinity 负无穷大 底e 虚数单位 ComplexIn 复平面上无 i 1 finity 穷远点
11
常数可以直接参与计算
（1）2*E （2）N[Pi,30] 给出30位有效数字
在函数中使用
（1）Sin[Pi/6] （2）Sqrt[3.3+1.5I] （3）Log[10,2.0]
8
1:算术运算
Mathematica软件的算术运算是指加减乘除以及 乘方、开方运算。 例如： （1）3*(5-2)+4^(6-3)/2 （2）3^(1/3) （3）1/3+2/5
精确计算与近似计算： 在Mathematica软件中，若输入的数据是精确的 ，计算结果保留精确数字。若计算近似值，可以 采用下面的方法： （4）N[3^(1/3)] 函数N[x]表示x的近似值
26
2.集合的运算
命令 First[list] Last[list] Part[list,i](Part[list,-i]) Take[list,n](Take[list,n]) Take[list,{m,n}] 取list的第m到第n个元素 Drop[list,{m,n}] 去掉list的第m到第n个元素 （1）First[t2] (*取t2的第一个元素*) （2）Part[t4,4]
(2): ?u (3): ?v 可以看出，上例中函数u[x]被定义为Sin[x]，而v[x]被定 义为常函数0.909297（=Sin2.0） 变量与函数在定义后可以用Clear命令清除其内容 (4): Clear[f,g,x];?f 在编写程序时,一般在每个程序段的开始将程序中所用的 21 符号的内容都用Clear命令加以清除,以防出错.
19
如果需要计算函数f[x]在某一点x=x0的值: (1) f[Pi/2] Mathematica软件中的函数可以用递归的方法进行定义 h[0]=1;h[n_]:=n*h[n-1];h[5] 在这里定义了h函数的初始值,以及一般的函数定义,需要 计算h[5]的具体数值
Mathematica使用入门数学软件Mathematica课件
16
内置函数
Mathematica 具有超过 3000 个内置函数，具体见 “帮助” “参考资料中心” 的左下角处的 “函数索引”
Mathematica 函数命名规则 ① 第一个字母大写，后面跟小写字母，如 Sin[x], Log[x] ② 大多数函数名与数学中的名称相同 ③ 当函数名分为几段时，每一段的头一个字母大写，后面 的用小写字母，如：ArcSin[x]
In[1]:= Clear[x,y]; In[2]:= f=2*x+y; In[3]:= f./{x->2,y->3} (* f(2,3) 的值 *) In[3]:= f./{2->5} (*把 2 替换成 5*)
14
数的基本运算
Mathematica 中的实数分精确数和双精度数
N[x,n] N[x] IntegerPart[x]
可以运行单个命令或语句 也可以运行多个语句：输入全部语句后再按执行键 命令（语句）分隔符：回车或分号 如果不需要显示运行结果：在语句后面加分号
输入和输出标识符：In[n]: 和 Out[n]
In[1]:= Sin[Pi/4]+Cos[Pi/4]
Out[1]= 2
8
数学公式的输入
数学公式的输入
lim(1 1
n
2
黄金分割数：1 1 5 2
虚部单位
1 ln n) 0.577215 n
无穷大
圆周率
11
基本运算符
基本运算：
算术运算 +-*/ ^ !
加减乘除 幂 阶乘（运算级别比加减乘除和幂运算高）
比较运算与逻辑运算 == > < >= <= != && || ! Xor
高数mathematica上机实验报告
15-16-2《高等数学》数学实验报告学号：姓名：得分： . ..实验：已知函数f(x)=（5≤x≤4）,作出并比较当c分别取1,0,1,2,3时的图形，并从图上观察极值点、驻点、单调区间、凹凸区间以及渐进线。
当c=时In[1]= f[x]:=1/(-1+2*x+x^2)Plot[f[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0]]In[3]= f'[x]:=Dt[f[x],x]Plot[f'[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f'[x]"]In[5]= f''[x]:=Dt[f[x],{x,2}]Plot[f''[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f''[x]"]In[7]=Solve[f'[x]==0,x]Solve[f''[x]==0,x]Out[7]=Out[8]=极大值点为x=，驻点为x=，单调递增区间为,单调递减区间为，下凸区间为，上凸区间为，渐进线为x=,x=当c=0时In[9]= f[x]:=1/(2*x+x^2)Plot[f[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0]]In[11]= f'[x]:=Dt[f[x],x]Plot[f'[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f'[x]"]In[13]= f''[x]:=Dt[f[x],{x,2}]Plot[f''[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f''[x]"]In[15]=Solve[f'[x]==0,x]Solve[f''[x]==0,x]Out[15]=Out[16]=极大值点为x=，驻点为，单调递增区间为,单调递减区间为，下凸区间为，上凸区间为，渐进线为x=,x=0当c=1时In[17]= f[x]:=1/(1+2*x+x^2)Plot[f[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0]]In[19]= f'[x]:=Dt[f[x],x]Plot[f'[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f'[x]"]In[21]= f''[x]:=Dt[f[x],{x,2}]Plot[f''[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBor[1,0,0],PlotLabel→"A Graph of f''[x]"]In[23]=Solve[f'[x]==0,x]Solve[f''[x]==0,x]Out[23]=Out[24]=无极值点，无驻点，单调递增区间为,单调递减区间为，下凸区间为，无上凸区间，渐进线为x=当c=2时In[25]= f[x]:=1/(2+2*x+x^2)Plot[f[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0]]In[27]= f'[x]:=Dt[f[x],x]Plot[f'[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f'[x]"]In[29]= f''[x]:=Dt[f[x],{x,2}]Plot[f''[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f''[x]"]In[31]=Solve[f'[x]==0,x]Solve[f''[x]==0,x]Out[31]=Out[32]=极大值点为x=，驻点为，单调递增区间为单调递减区间，下凸区间为，上凸区间为，无渐进线当c=3时In[33]= f[x]:=1/(3+2*x+x^2)Plot[f[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0]]In[35]= f'[x]:=Dt[f[x],x]Plot[f'[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f'[x]"]In[37]= f''[x]:=Dt[f[x],{x,2}]Plot[f''[x]//Evaluate,{x,-5,4},GridLines→Automatic,Frame→True,PlotStyle→RGBColor[1,0,0],PlotLabel→"A Graph of f''[x]"]In[39]=Solve[f'[x]==0,x]Solve[f''[x]==0,x]Out[39]=Out[40]=极大值点为x=，驻点为，单调递增区间为单调递减区间，下凸区间为，上凸区间为，无渐进线。
数学软件Mathematica的应用
数学软件Mathematica的应用一、数学软件Mathematica简介★ Mathematica是由美国Wolfram公司研究开发的一款著名的数学软件；★ Mathematica能够完成符号运算、数学图形的绘制等，功能非常强大；★ Mathematica能够做精确计算；★ Mathematica的界面操作非常友好；★ Mathematica是数学建模常用的数学软件之一。
二、利用模板进行微积分运算File（文件）→Palettes（模板）→BasicInput（基本输入）File（文件）→Palettes（模板）→BasicCalculations（基本计算）三、Mathematica中一些常用的函数（1）数学常数（2）数学函数（3）数值函数（4）表操作函数（5）代数函数（6）微积分运算函数（7）作图函数（8）数值分析函数下面是其常用的几种形式：（9）编程相关在Mathematica中，一个逻辑表达式的值有三个：真（True）、假（False）和“非真非假”。
条件控制函数If（1）If语句的结构与一般的程序设计语言中的If的结构类似。
它有三种情况： If[逻辑表达式，表达式1]当逻辑表达式的值为真时则计算表达式1，表达式1的值就是整个If结构的值；If[逻辑表达式，表达式1，表达式2]当逻辑表达式的值为真时则计算表达式1，为假时则计算表达式2； If[逻辑表达式，表达式1，表达式2，表达式3]当逻辑表达式的值为真时则计算表达式1，为假时则计算表达式2，其它情况则计算表达式3。
循环控制语句Mathematica 中有3种描述循环的语句，它们是Do,While 和For 语句。
下面是其一般形式：For[初值，条件，修正，循环体] While[条件，循环体] Do[循环体，{循环范围}]四、结合图形进行分析1．作出函数xx f y 1sin )(==在区间]1,1[-上的图像，观察当0→x 时函数的变化情况；作出函数xx x f y 1sin )(==在区间]1,1[-上的图像，观察当0→x 时函数的变化情况；2．作出双曲抛物面xy z =的图形； 3．作weierstracs 函数)13cos(21)(1x x f n n nπ∑∞==（处处连续但处处不可导）的图像；4x ∈(-5,5), y ∈(-5,5)内的所有根；五、验证与探索1．x sin 的泰勒级数 2．x sin 的无穷乘积猜想六、算法与程序1．分形图（迭代）2．将矩阵化为行最简形（步骤）七、实际问题的Mathematica 求解1．椭圆弧长的计算问题计算椭圆βα≤≤⎩⎨⎧==t t b y ta x ,sin cos 的弧长及近似值。
Mathematica软件在高等数学教学中的应用
容 易 出错 的地 方 。Ｍａ ｅ ａｉ ｔ ｍ ｔ ａ对 计 算 函数 的 导 数 也 是 相 当容 易 ｈ ｅ
的。
Ｍ ｔｅ ａｉ 计 算 函数 导 数 命 令 ： 【 数表 达式 ， ａｈｍ ｔ ａ ｅ Ｄ函 自变量 】 如 已 知 Ｙ ｌ ２ 求 。 要 利 用 Ｍａ ｅ ａ ｅ“ ＝ ＿ｘ ， 只 ｔ ｍ ｔ ａ 基本 输 入工 ｈ ｉ
ｂ ｌ ｆ ｅｔ ｎ ｔ ｓｃ ｍ ｐ ｅ ｅ ｓ ｅ ｐ ｏ ｌｍｓ ｉ ｔ ｏ ｔ ｉ ｇ ｍａｈ ｏ ｒ ｈ ｎ ｉ ｒ ｂ ｅ ． ｉ ｙ ｓｌ ｖ Ｋｅ ｙ ｗｏｒ ｓ Ｍ ａｈ ｍａ ｃ ｏ ｔ ａｅａ ｖ ｎ ｅ ｔ ｅ ｔ ｓ ｏ ｕ ａ ｏ ；ｌ ｔ ｏ ｒｈ ｎ ｉｅ ｐ ｏ ｌｍｓ ｄ ： ｔ ｅ ｔ ａｓ ｆ ｒ ；ｄ ａ ｃ ｄ ｍａ ｍａｉ ； ｍｐ ｔｔ ｎ ｐ ｏ ； ｍｐ ｅ ｅ ｓ ｒ ｂ ｅ ｉ ｗ ｈ ｃｃ ｉ ｃ ｖ
４ ｘ
结 果 一 （ ２ｘ ） ３ — ２￣ ｔｅ ｔ ａ软 件简 介 ｃ
Ｍａｈｍ ｔ ａ是 由 位 于 美 国 伊 利 诺 州 的 伊 利 诺 大 学 香 槟 ｔｅ ａｉ ｃ （ ｈｍ ａ ｎ Ｃ ａ ｐ ｉ ）分 校 附近 的 ＷｏｒｍＲ ｓａｅ ｇ ｌａ ｅｅｒｈ公 司 开 发 的 一 套 专 门 ｆ 进行数学计算 的软件 ， 目前 最 新 的发 布 版 是 Ｍａ ｅ ａｅ ．。该 ｔ ｍ ｔ ａ５ ｈ ｉ ０
Ａｂ ｔ ｃ ： ａｈ ｍａｉａ ｓ ｆ ｒ ｓａｓ ｉ ｏ ｅ ｉ ｔ ｔｓｃ ｍｐ ｔｔ ｎ ｓ ｆ ｒ ． ａｈ ｍａｉａ ｓ ｆ ａ ｅ ｉ ｉｔｏ ｕ ｅ Ｏ ａｓ ｔｎ ｅ ｃ ｉ ｇ ｉ ｓｒ ｔ ａ Ｍ ｔ ｅ ｔ ｏｔ ｅ ｉ ｕ ｔ ｆｓ ｃａ ｙ ｍａｈ ｏ ｕ ａｉ ｏｔ ｅ Ｍ ｔ ｅ ｔ ｏ ｔ ｒ ｓ ｎ ｒ ｄ ｃ ｄ ｔ ｓｉ ａ ｔ ａ ｈ ｎ ｃ ｗａ ｐ ｌ ｏ ｗａ ｃ ｗ ｓ ｔ ｎ
《数学实验》实验报告——用Mathematica软件解微分方程
Out[1]= In[2]：= s1=LaplaceTransform[1,t,s]
1
Out[2]= s In[3]：= x''[0]=x'[0]=x[0]=0;Solve[f1 s1,LaplaceTransform[x[t],t,s]]
３
Out[3]= In[4]：= Out[4]= 上例中的 LaplaceTransform[x[t]，t，s]就是教材中的 X（s） ，In[3]解出 X（s） ， 其余过程与教科书完全相同。
例 5 求微分方程ቤተ መጻሕፍቲ ባይዱ：
x 2 x y 2 y 0 t x y 2 x 2e 满足条件 x（0）=3，x′（0）=2，y（0）=0 的特解。
In[1]：= f1=LaplaceTransform[ {x''[t] - 2x'[t] - y'[t] + 2y[t],x'[t] + y'[t] - 2x[t]},t,s]
例2 求常微分方程 y′= x2 + y2，满足初始条件 y（0）= 0 的数值解 例3 求函数 t 5 和 et sint 的拉氏变换 例 4 用拉氏变换解微分方程：
x + 3x″+ 3x′+ x = 1 满足条件 x″(0) = x′(0) = x(0) = 0 的解。
例 5 求微分方程组：
《数学实验》实验报告
班级 试验 内容 **** 学号 **** 姓名 试验 类别 **** 成绩 试验 时间
用 Mathematica 软件解微分方程
自选试验
试验问题：
例1 求解下列微分方程： 1) 2) 3)
y 2 (1 y) (2 y) 2
Mathematica软件使用入门讲解
Mathematica软件使用入门目录第一章基本知识与基本操作 (3)1.1 Mathematica的基本语法特征 (3)1.2 Mathematica的启动、基本操作 (4)1.3 操作小技巧 (7)1.4 数值计算 (8)1.5 赋值与替换 (9)1.6 自定义函数 (10)1.7 方程与方程组解 (11)1.8 解不等式与不等式组 (12)1.9 由递推式求数列的通项公式 (13)1.10 作函数图像 (14)第二章运用Mathematica实现高等数学中的基本运算 (16)2.1 求极限运算 (16)2.2 求导数与微分 (18)2.3 求不定积分 (25)2.4 求定积分 (25)第三章实验练习题 (28)Mathematica是当今世界上最为流行的计算机代数系统之一．Mathematica系统是美国物理学家Stephen.Wolfram领导的一个小组开发的,后来他们成立了Wolfram研究公司．1987年推出了系统的1.0版；现在的最新版本是8.0版．Mathematica可以做：●符号计算和数值计算问题,如：能做多项式的计算、因式分解和展开等；●做各种有理式计算，求多项式、有理式方程和超越方程的精确解和近似解；●做向量、矩阵的各种计算；●求极限、导数、积分，做幂级数展开，求解某些微分方程等；●做任意位数的整数或分子分母为任意大整数的有理数的精确计算，做具有任意位精度的数值（实、复数值）的计算．●可以很方便地画出用各种方式表示的一元和二元函数的图形,通过图形,可以立即形象地掌握函数的某些特性，而这些特性一般是很难从函数的符号表达式中看清楚．第一章 基本知识与基本操作1.1 Mathematica 的基本语法特征使用Mathematica ，一定要牢牢记住：● Mathematica 中大写小写是有区别的，如Name 、name 、NAME 等是不同的变量名或函数名；● 系统所提供的功能大部分以系统函数的形式给出， 内部函数一般写全称， 而且一定是以大写英文字母开头， 如Sin[x], Cos[z]等；● 乘法即可以用*，又可以用空格表示，如2 3＝2*3＝6 , 2 Sin[x]＝2* Sin[x] ● 乘幂可以用“^”表示，如x^0.5 表示： Tan[x]^y 表示： ● 自定义的变量可以取几乎任意的名称，长度不限，但不可以数字开头． ● 当你赋予变量任何一个值，除非你：明显地改变该值或 使用Clear[变量名] 或 使用“变量名=.”取消该值，否则它将始终保持原值不变．● 一定要注意四种括号的用法：0.5x yTan[x]( )：表示项的结合顺序，如： (x+(y^x+1/(2x)))；[ ]：表示函数，如：Log[x], Sin[x]；{ }：表示一个“表”(即是一组数字、或任意表达式、或函数等的一个有序集合)，如：{2x,Sin[12 Pi]，A，1}, {1+A,y*x，1，2}；[[ ]]：双方括号表示“表”或“表达式”的下标，如：a； {3,5,7}[[2]]=5．a[[2,3]]表示：23●Mathematica的语句书写十分方便，一个语句可以分为多行写，同一行可以写多个语句（但要以分号间隔）．●当语句以分号结束时，语句计算后不做输出（输出语句除外），否则将输出计算的结果．●Mathematica命令中的标点符号必须是英文的．1.2 Mathematica的启动、基本操作1.2.1 启动“Mathematica”：在windows操作系统中安装了Mathematica后，与其他的常用软件一样，可从“开始”→“程序”→“Mathematica5” Mathematica的主窗口并出现第一个notebook窗口（Untitled-1）:1.2.2 简单使用：例1.1 计算＋33的值①在“Ｕntitled－１”窗口中输入：329/412+3^3②按下“Shift＋Enter”（或数字键盘上的Enter键）,就得到计算结果：其中“In[1]:=”是Mathematica自动加上的，表示第一个输入；“Out[1]:=”表示第一个输出．一般地：In[n]:= 表示第n个输入Out[n]:=表示第n个输出．注意：“In[n]:=”自动加上的，不能人工输入!1.2.3 保存结果：保存方法同一般的Windows软件：“文件”→“保存”⇒“另存为”窗口→在“查找范围”内找到目标文件夹→输入文件名（比如输入“1”）→“”．Mathematica 4或Mathematica 5的文件的后缀是“nb”，当输入“1”时，即产生文件“1.nb”．1.2.4打开文件1.nb启动Mathematica →“文件”→“打开”⇒打开”窗口：→在“查找范围”内找到文件“1.nb”→“”即可．1.2.5 退出Mathematica ：与一般应用软件一样，单击右上方的“ ”按钮（或用菜单：“文件”→“退出”）．1.3 操作小技巧1.3.1 Ctrl+K 的用途如果只知道命令的首写字母， 可在输入该首写字母（要大写），再按下“Ctrl+K ”组合键， 则所有以该字母为首的命令都列出来，只要用鼠标双击命令名就输入了该命令． 1.3.2 使用前面已有的结果 举例如下：例1.2 做如下操作：① 输入：Integrate[x^2*(11-Sin[x]),{x,-1,1}]按：“Shift ＋Enter ”； ② 输入：%＋１，按：“Shift ＋Enter ”； ③ 输入：%＋１，按：“Shift ＋Enter ”； ④ 输入：%1＋1，按：“Shift ＋Enter ”； ⑤ 输入：%3＋1，按：“Shift ＋Enter ”， 计算结果如下：Integrate[f,x]是求：()f x dx ⎰ Integrate[f,{x,xmin,xmax}]是求：maxmin()x x f x dx ⎰可见，“％”表示前一个计算结果；“％n”表示第n个计算结果.1.3.3删除行：见下图示1.4 数值计算请看下例：只要选定且删除此即可系统默认的计算结果,是精确的．N[]，取近似值函数，默认输出6位有效数字．N[]，取近似值函数，指定输出3位有效数字．1.5 赋值与替换X=. 或Clear[x] 清除赋给x的值expr/.{x->xval,y->yval} 用xval、yval分别替换expr中的x、y．例1.3输入：x=3；y=4；w=x+y 计算输入：Clear[x,y]；计算输入：z=(x+y)^2 计算输入：z/.x->5 计算输入：Clear[x,y]；计算输入：u=x+y 计算将(x+y)^2赋给z清除变量的定义和值变量替换：用5代替表达式z中的变量变量替换：输入：u/.{x->5,y->6} 计算 计算结果如下：1.6 自定义函数用户可以自行定义函数，一个函数一旦被定义好之后就可以象系的内部函数一样使用．例1.4 如要定义函数f(x)=x 2＋3x-2只要键入：f[x_]:=x^2+3x-2即可．又如要定义分段函数2+1 < 0()= 2sin 0x x g x x x ⎧⎨≥⎩“:=”是定义符．左边f 是函数名,方括号内x 是自变量，其后的下划线“_”不能少．右边是函数的表达式．可键入：g[x_]:= Which[x<0,x^2+1,x>=0,2Sin[x]]或g[x_]:=If[x<0,x^2+1,2Sin[x]]请见以下计算结果：1.7 方程与方程组解例1.5 ① 解方程:0652=+-x x输入：Solve[x^2-5x+6==0,x]即可．② 解方程组 输入：Solve[{x+y==1,3x^2-y^2==0},{x,y}] 即可（结果见下图）．Solve 是解方程或方程组的函数．其格式为：Solve[eqns,vars]其中方程用exp==0的形式（其中exp 为未知元的表达式，“= =”必须是2个等号）；方程列表 2213x y x y +=⎧⎨-=⎩未知数列表1.8 解不等式与不等式组例1.6 ① 解不等式组⎪⎩⎪⎨⎧>-<--0101222x x x输入: <<Algebra`InequalitySolve`InequalitySolve[{x^2-5x-6<0,x^2-1>0}, x]即可．② 解不等式3)3(12>--x x输入： <<Algebra`InequalitySolve`即可（结果见下图）不等式列表 变量列表加载解不等式的程序包，这是必须的，可谓是固定的格式， “< ”为键盘上的小于号, “`”为数字键1的左侧的Algebra —— 代数类InequalitySolve —— 解不等式程序包绝对值函数注： Mathematica 系统有内部函数.还有一些系统扩展的功能但不是作为内部函数的、以文件的形式存储在磁盘上的文件，要使用它们，必须用一定的方式来调用这些文件，这些文件我们称之为程序包. 调用方式之一如上所述：<<Algebra`InequalitySolve`或用：Needs["Algebra`InequalitySolve`"]1.9 由递推式求数列的通项公式例1.7 设 求数列的通项公式只要输入：<<DiscreteMath`RSolve` RSolve[{a[n]==n ＊a[n-1], a[1]==1}, a[n], n]即可（结果见下图）11,1,n n a na a -==函数名 递推关系 初始条件调用程序包 类名，此处是函数类 函数类中的这个函数离散类 离散类中的这个函数1.10 作函数图像例1.8在同一坐标系中作出2-1y x 和y=sinx在[－2,2]内的图像．输入： Plot[{x^2-1,Sin[x]},{x,-2,2}] 结果见下图例1.9作出sinxcosy的三维图形输入：Plot3D[Sin[x]*Cos[y],{x,－2Pi,2Pi},{y,－即可（结果见下图）增加取样点提高光滑度第二章 运用Mathematica 实现高等数学中的基本运算极限、导数和积分是高等数学中的主要概念和基本运算，如果你在科研中遇到较复杂的求极限、求导数或求积分问题，Mathematica 可以帮你快速解决这些问题。
mathematica数学实验报告
高等数学实验报告
实验一
一、实验题目
1：作出各种标准二次曲面的图形
2：作出曲面所围的图形
二、实验目的和意义
方法的理论意义和实用价值。
如利用数形结合的方法观察数列的极限，可以从点图上看出数列的收敛性，以及近似地观察出数列的收敛值；通过编程可以输出数列的任意多项值，以此来得到数列的收敛性。
通过此实验对数列极限概念的理解形象化、具体化。
三、计算公式
请写出在程序中所需要的计算公式。
比如定积分的数值计算中，如用梯形法计算的，请描述梯形法的公式。
四、程序设计
五、程序运行结果
六、结果的讨论和分析
如初值对结果的影响；不同方法的比较；该方法的特点和改进；整个实验过程中（包括程序编写，上机调试等）出现的问题及其处理等广泛的问题，以此扩大知识面和对实验环节的认识。
Mathematica软件在高等数学教学中应用
Mathematica软件在高等数学教学中应用摘要：本文通过一些具体的例子，介绍了Mathematica 软件在高等数学教学中的应用。
说明在高等数学教学中融入软件的学习，不仅使得抽象概念变得形象生动，而且能避免冗长繁杂的计算，从而激发学生学习高等数学的兴趣。
关键字：Mathematica软件高等数学教学应用一、引言极限、导数、定积分等概念，可以说是高等数学中最重要、最具有代表性的概念，它们体现了应用微积分的思想和方法，其应用几乎涵盖了所有的自然学科。
但上述概念对于学生来说也是最难理解的，因为从本质上来说它们有三种表示形态：逻辑形态、算法形态和直观形态。
大学老师呈现最多的是前两种形态，因此造成大部分学生觉得高等数学的学习抽象枯燥，运算繁琐冗长。
为了帮助学生解决认知中的困难，首先通过数学软件的直观演示，加深学生对一些重要概念的理解，然后再详细地介绍它们的逻辑形态和算法形态，这样使得抽象概念的学习更加形象生动。
下面就Mathematica软件在教学中的具体应用谈谈心得体会。
二、Mathematica软件在高等数学教学中的应用1.运用软件演绎极限的概念在同济版的高等数学教材中，数列极限的引入借用的是刘徽的割圆术，即利用圆内接正多边形来推算圆的面积，具体过程如下：设有半径为r的圆，首先作内接正六边形，把它的面积记为A■；再作内接正十二边形，其面积为A■；循此下去，每次边数加倍，一般的把内接正6×2■边形的面积记为A■。
当n越大，内接正n边形与圆的差别就越小，从而用其内接正n边形的面积A■逼近圆面积S，由图1经过计算可知A■=nr■sin■cos■ （n=3，4，5，…），当n无限增大时，A■无限逼近S。
上述的文字叙述过程在课本中非常繁琐，如果我们只用语言表达，学生理解起来会比较吃力，因为他们看不到n无限增大时，A■与S逼近的程度。
如果用Mathematica 软件，在图1中用动画的方式将上述过程演示出来，学生就会更加直观地看到上述逼近的过程，从而对极限概念有一个更直接的感官认识。
数学软件Mathematica的应用
数学软件Mathematica的应用一、数学软件Mathematica简介★ Mathematica是由美国Wolfram公司研究开发的一款着名的数学软件；★ Mathematica能够完成符号运算、数学图形的绘制等，功能非常强大；★ Mathematica能够做精确计算；★ Mathematica的界面操作非常友好；★ Mathematica是数学建模常用的数学软件之一。
二、利用模板进行微积分运算File（文件）→Palettes（模板）→BasicInput（基本输入）File（文件）→Palettes（模板）→BasicCalculations（基本计算）三、Mathematica中一些常用的函数（1（2（3（（5（6（8）数值分析函数在Mathematica中，一个逻辑表达式的值有三个：真（True）、假（False）和“非真非假”。
条件控制函数If（1）If语句的结构与一般的程序设计语言中的If的结构类似。
它有三种情况： If[逻辑表达式，表达式1]当逻辑表达式的值为真时则计算表达式1，表达式1的值就是整个If结构的值；If[逻辑表达式，表达式1，表达式2]当逻辑表达式的值为真时则计算表达式1，为假时则计算表达式2；If[逻辑表达式，表达式1，表达式2，表达式3]当逻辑表达式的值为真时则计算表达式1，为假时则计算表达式2，其它情况则计算表达式3。
循环控制语句Mathematica中有3种描述循环的语句，它们是Do,While和For语句。
下面是其一般形式：For[初值，条件，修正，循环体]While[条件，循环体]Do[循环体，{循环范围}]四、结合图形进行分析1．作出函数xx f y 1sin )(==在区间]1,1[-上的图像，观察当0→x 时函数的变化情况；作出函数xx x f y 1sin )(==在区间]1,1[-上的图像，观察当0→x 时函数的变化情况；2．作出双曲抛物面xy z =的图形； 3．作weierstracs 函数)13cos(21)(1x x f nn nπ∑∞==（处处连续但处处不可导）的图像；4x ∈(-5,5), y ∈(-5,5)内的所有根；五、验证与探索1．x sin 的泰勒级数2．x sin 的无穷乘积猜想六、算法与程序1．分形图（迭代）2．将矩阵化为行最简形（步骤）七、实际问题的Mathematica 求解1．椭圆弧长的计算问题计算椭圆βα≤≤⎩⎨⎧==t t b y ta x ,sin cos 的弧长及近似值。
Mathematica4.0使用方法(数学实验课讲义)
Mathematica4.0使用方法数学实验课教材首钢工学院Mathematica数学实验Mathematica 是一个交互式的计算系统．这里说的交互是指：在使用Mathematica 系统的时候，计算是在使用者（用户）和Mathematica 互相交换、转递信息数据的过程中完成的．用户通过输入设备（一般指计算机键盘）给系统发出计算的指令（命令），Mathematica 完成给定的计算工作后把计算结果告诉用户（一般通过计算机显示器）．Mathematica 是一个集成化的计算机软件系统．它的主要功能包括三个方面：符号演算、数值计算和图形绘制．例如，它可以完成多项式的各种计算（四则运算、展开、因式分解）；可以求多项式方程、有理式方程和超越方程的精确解和近似解；做数值的或一般表达式的向量和矩阵的各种计算；求一般函数表达式的极限、导函数、积分、幂级数展开，求解微分方程等等．根据教学大纲的要求及学校的课时安排（共12课时，内含2课时考试），我们将Mathematica数学软件的学习缩编成下面的四个实验，以期在短时间内使同学们掌握该软件的基本使用方法，学会用它解决高等数学中的一些常见问题．目录第一篇微积分 (1)实验一……………………………………………………实验二……………………………………………………实验三……………………………………………………实验四……………………………………………………第二篇线性代数……………………………………………………实验一……………………………………………………实验二……………………………………………………第三篇概率统计……………………………………………………第四篇复数与积分变换……………………………………………附录Mathematiac一部分函数及意义……………………第一篇微积分实验一一、实验目的1.学习在Windows下Mathematica 4.0软件的启动与退出，并熟悉其界面；2.建立文件与保存文件；3.学习用基本运算符号和模板进行加、减、乘、除、乘方、开方等常用的算术运算；4.学习表示计算结果的近似结果；5.会用符号或模板进行常见函数的输入及多项式的变换；6.会给变量赋值．二、内容与步骤1.Mathematica 4.0的启动与退出启动计算机，屏幕上显示Windows界面，单击“开始”进入主菜单，将鼠标移向“程序”，找到包含Mathematica 4.0的程序组，单击可执行程序Mathematica 4.0就进入了该系统，此时系统已进入交互状态，在等待用户输入命令．当软件使用完毕后，需要退出Mathematica系统时，只须单击工作窗口右上方的“File”菜单中选用命令“Exit”，或者按“Alt+F4”键均可退出系统，回到操作系统状态．例如：输入2+3后，按Enter+Shift组合键或右边小键盘上的Enter键运行，屏幕上就显示出In[1]:=2+3Out[1]=5其中In[1]：= 表示第一个输入，Out[1] = 表示第一个输出，它们是在运行后由系统自动显示的，用户不必输入．注意：若直接按左边的Enter键，只是在输入的组合命令中起换行的作用．2.建立文件与保存文件在工作窗口做好的某些内容，如果想要保留以供今后多次使用，通常是建立一个文件，将做好的内容保存在文件中．单击File/ Save as，在文件名N一栏内键入一个文件名,然后左击保存S．3.算术运算与模板的使用a)：输入基本运算符号加+减-乘*（或用一个空格表示相乘）除/幂乘yx^优先运算：用圆括号，并可重复多次使用．b)：模板的调出与运用方法一：在Mathematica 3.0以上版本的输入中，可以使用工具按钮输入各种运算，其步骤如下：①单击菜单栏中的文件File选项；②在下拉菜单中选择调色板Palettes选项；③在下一级菜单中单击基本计算BasicaCalculations选项，将会另外出现一个工具窗口；④在其窗口中单击计算与数值Arithmetic and Numbers选项前的符号“”，使其符号变成“▽”，将出现加、减、乘、除、乘方、开方等工具按钮；⑤单击需要的按钮，在原Notebook窗口中将会出现相应输入格式，将光标移到标有“□”的位置上，输入数值或表达式，就可以完成输入格式．方法二：在第③步，在下一级菜单中单击基本计算BasicInput 选项，出现一个常用的含有多种运算的模板（加、减可以直接从键盘输入+、-号） 4. 近似与精确 a ） 命令输入：N[表达式，n] 精确到n 位有效数字；N[表达式] 近似值按计算机默认的数位（6位）处理； [表达式]// N 同上；% 表示最近一次计算机运行后的输出结果；注意：1）当输出结果是610以下的数字，近似值按计算机默认的6位有效数字处理；610及610以上的近似值计算机按科学计数法处理．2）N[表达式，n] 表示精确到n 位有效数字（注：当n=1~16时，结果都按计算机默认的6位处理）． b) 模板调出：与上述算术运算模板调出的方法一相同． 例1 1）输入： N310,结果显示：0.0141592653589792）输入：N ,结果显示： 3.1（按计算机默认的6位处理） 3）输入：N %, 表示对当前结果取18位有效数字近似 4）输入：4566000.66777777777777结果显示：4.5665.Mathematica中的常数、数学函数与常见的多项式变换a)直接从键盘输入（在英文状态下）Mathematica的常数：Pi 表示πE 表示eDegree （π/180）表示度I 表示虚数iInfinity 表示无穷大∞Mathematica中常用的数学函数：幂函数Sqrt[x] (求平方根) ；指数函数Exp[x] (以e为底的指数函数)；对数函数Log[x] (以e为底的对数函数)；Log[a，x] (以a为底的对数函数)；三角函数Sin[x]，Cos[x]，Tan[x]，Cot[x]，Sec[x]，Csc[x]；反三角函数ArcSin[x]，ArcCos [x]，……；双曲函数Sinh[x]，Cosh[x]，Tanh[x]，Coth[x]，……；反双曲函数ArcSinh[x]，……．Mathematica中常见的多项式变换：Factor[表达式] 将表达式分解因式Expand[表达式] 将表达式展开成多项式和的形式Simplify[表达式] 将表达式化简成最简形式Apart[表达式] 将表达式分解为部分分式之和函数表达式的运算规则有：1）.它们都以大写字母开头，后面用小写字母．当函数名可以分成几段时，每一个段的头一个字母用大写，后面的字母用小写．例如，ArcSin[x]．2）.函数的名字是一个字符串，其中不能有空格．3）.函数的自变量表用方括号括起来，不能用圆括号．4）.多元函数的自变量之间用逗号分隔．b）模板介绍在Mathematica3.0以上版本的输入中，可以使用工具按钮输入各种函数，其步骤如下：①击菜单栏中的文件File选项；②在下拉菜单中选择调色板Palettes选项；③在下一级菜单中单击基本计算BasicaCalculations选项，将会另外出现一个工具窗口；在其窗口中单击三角与指数函数Trigonometric and Exponential Finctions选项前的符号“”，使其符号变成“▽”并列出子选项的清单；在此清单中单击三角Trigonometric选项前的符号“”，使其符号变成“▽”，将会出现一些三角函数和反三角函数工具按钮；单击需要的按钮，在原Notebook窗口中将会出现相应三角函数或反三角函数输入格式，将光标移到标有“□”的位置上，输入数值或表达式，就可以完成输入格式；在此清单中单击指数与对数Exponential and Logarithmic选项前的符号“”，使其符号变成“▽”，将会出现一些指数与对数函数工具按钮；单击需要的按钮，在原Notebook窗口中将会出现相应指数或对数函数输入格式，将光标移到标有“□”的位置上，输入数值或表达式，就可以完成输入格式；在此清单中单击双曲函数Hyperbolic选项前的符号“”，使其符号变成“▽”，将会出现一些双曲函数和反双曲函数工具按钮；单击需要的按钮，在原Notebook窗口中将会出现相应双曲函数和反双曲函数输入格式，将光标移到标有“□”的位置上，输入数值或表达式，就可以完成输入格式． 在其窗口中单击计算与数值Algebra 选项前的符号“”，使其符号变成“▽”，将出现Polynomial Manipulation ，Simplifyication 等工具按钮进行相关选择即可完成多项式的变换； 例2In[13]：＝Log[2，3.256] Out[13]：＝1.7031 例3：已知1 ,1232221-=-+=x p x x p ，计算2121 ,p p p p ⨯+，21p p ÷并将2121 ,p p p p ⨯+的结果分解因式、展开多项式，将21p p ÷的结果分解为部分分式 输入：p1 3x^22 结果显示： 12xp2 x^1P122x p11x212xp1Factor p11x1x21Expand p1p212x 4x 22x3Apart p 136. 变量赋值：命令格式：x= a 将值a 赋给变量xu=v=a 将值a 赋给变量u 、v （给多个变量赋值）f[x]/. x->a 变量x 赋值为a （求函数f[x]在x=a 时的值） u := 延迟赋值，按Shift+Enter 键没有结果输出，待给变量赋值运行后才有结果u= 直接赋值，按Shift+Enter 键后有结果输出 u=. 清除变量u 的值Clear[x] 清除变量x 的值，多用作清除函数注意：应随时将以后不再使用的变量的值清除掉，以免影响后面某些计算结果的正确性．习题一１． 计算1)62456log 3e -+并保留15位有效数字.2) sin(30)+tan(6π)并精确到小数点后7位.3)7lg 21arctan 1arcsin ++2. 给变量赋值并计算1) 若x=6,y=e,z=x+3y ，计算3z-5y 2+6(x-7)52)x=3,y=5π,计算（lgx ）⨯arcos(2y)- 9并保留18位有效数字.3．设p1=2x-1, p2=3x-7, 求 p1×p2, 并展开它，再分解因式，最后将 1／（p1×p2）分解为部分分式. 练习过程及答案N 34Log 2,566,316．8．1.0z x 3y . x 6,y3z 5y 26 x 75. x 6,y ,z 665 23 69.000000000000000000.33490675722196522x 1 3x 12x 73Expandy9．实 验 二一、实验目的1、学习使用自定义函数，会求函数值；2、学习用绘图语句作函数图形；3、学习用解方程的语句解方程、方程组；4、会建立表，进行表的基本运算. 二、内容与步骤 1、自定义函数:一般函数： f[x_]= 表达式 定义的规则x 可以被替代 f[x_]：= 表达式 延迟赋值 f[x_]=. 清除f[x_]的定义Clear[f] 清除所有以f 为函数名的函数定义 分段函数：Which[条件1，表达式1，条件2，表达式2，…条件n ，表达式n]Which 语句是表示分段函数的常用语句． 例1：定义函数：x x x x f cos )(2++=，并求f (2)的值输入命令：显示输出： 4.9输入显示结果注意：f[2.]表示求自变量为2时函数的近似值；f[2]表示为精确值．．10．例2:定义函数....0()0.. 0....0x x g x x x x >⎧⎪==⎨⎪-<⎩并求：)0(),3(),1(g g g -的值 输入命令g x_: Which x 0,x,x 0,0,x0,（将分段函数自定义成一个函数）显示结果 1 输入显示结果 3 输入显示结果 0注意：中括号内的等号要输成双等号 2．作图：1）基本作图命令格式（a ）只规定自变量范围的作图命令：Plot[f(x),{x,x1,x2}](b) 不仅规定自变量范围，还规定因变量范围的作图命令Plot[f(x),{x,x1,x2}，PlotRange->{y1,y2}](c) 不仅规定自变量范围，还可以加标注（函数名称，坐标轴） Plot[f(x),{x,x1,x2}，PlotLabel->“表达式 ”，AxesLabel ->{“x ”,“y ”}11．2）观察函数图形的叠加情况设)...(),(21x f y x f y==，若在一个坐标系里观察这几个函数图像命令格式为：Plot[{ )(),(21x f x f },{x,x1,x2}]注意：不要将“ )(),(21x f x f ”写成“ )(),(21x f y x f y ==”例3：做出y=sinx 在[-4π4π]之间的图像Plot S in x , x ,4Pi,4例4：做出y=tanx 在[0,4π]，y ∈[-5,5]之间的图像PlotT an x , x ,0,4 ,PlotRange 5,．12．例5：做出y=sinx,sin2x,sin3x 在[0,2π]内的标出坐标轴的且用三种不同颜色标示的图像．3) 分段函数的作图先利用条件语句Which 自定义分段函数，然后用Plot 语句画出分段函数的图形格式步骤：首先输入 f [x _]：= Which[条件1，表达式1，条件2，表达式2，…条件n ，表达式n]再输入 Plot[f(x),{x,x1,x2}] 例6 作出....0()0.. 0....0x x g x x x x >⎧⎪==⎨⎪-<⎩的图像g x _ : Which x 0,x,x 0,0,x0, Plot g x , x ,2,13 ．4）参数方程作图使用 ParametricPlot 函数可以画参数形式的图形，格式如下： ParametricPlot[{x(t)，y(t)}，{t ，a ，b}，可选项]ParametricPlot[{{x1(t)，y1(t)，{x2(t)，y2(t)}，...}，{t ，a ，b}，可选项]例7 画出圆的参数方程的⎩⎨⎧==ty tx sin cos ，0＜t ＜2π曲线图形解 In[5]：＝ParametricPlot[{Sin[t]，Cos[t]}，{t ，0，2Pi}，AspectRatio －>Automatic] Out[5]：＝AspectRatio ：指定作图的纵横比例．系统默认值约0.618：1．可以为 AspectRatio 指定任何一个其他数值．如果希望系统按实际情况作图即纵横比例为1：1，则需要将这个可选项设置为Automatic ． 5）二元函数的图像命令格式：首先定义二元函数: z[x_,y_]：=表达式 然后作图Plot3D[z[x,y],{x,x1,x2},{y,y1,y2}]．14．例8 做出222y x z +=的图像输入： 输出：上述命令大多可以通过模板调出 ① 左击菜单栏中的文件File 选项； ②在下拉菜单中选择调色板Palettes 选项；③在下一级菜单中单击基本计算BasicaCalculations 选项，将会另外出现一个工具窗口；④ 在其窗口中单击图形Graphics 选项前的符号“”，使其符号变成“▽”并列出子选项的清单进行选择3．解方程： 1）解方程命令格式：Solve[f(x)= =0，x] 2) 解方程组命令格式：Solve[{f (x)= =0，g (y)= =0,…},{x,y,…}]15．上述命令可以通过模板调出 ①左击菜单栏中的文件File 选项； ②在下拉菜单中选择调色板Palettes 选项；③在下一级菜单中单击基本计算BasicaCalculations 选项，将会另外出现一个工具窗口；④在其窗口中单击图形Algebra 选项前的符号“”，使其符号变成“▽”并列出Solving Equations 选项的清单进行选择 例9 求方程063523=++-x x x 的根． 解： 输入Solve x 35x 23x 60 输出：例10 求方程组⎩⎨⎧=+=-ny x m y x 2的根 解： 输入Solvex 2y m,x y n , x ,输出：例11求解方程b x x =++-11 解： 输入输出：4．表的操作 1）表的生成．16．一维表：{a,b,c…}二维表(表中表)：{{一维表1},{一维表2},{一维表n}…} 如：一维表{1，2，3}，二维表{{1，2}，{5，2}，{6}}2）表中元素的提取一维表b 的第i 个元素: b[[i]] 或Part[[b,i]] 二维表b 的第i 个分表：b[[i]] 或Part[[b,i]] 二维表b 的第i 个分表中的第j 个元素: b[[i,j]] 如：b={{1，2}，{5，2}，{6}} b[[2]]-----显示 {5，2} b[[2,1]]----- -显示53）表的运算设b1，b2表示结构完全相同的两个表，表b1，b2的和、差、积、商等于对应元素的相应运算（分母不为零）b1={{1，2}，{5，2}，{6}}，b2={{3，1}，{0，2}，{2}} b1+ b2={{4，3}，{5，4}，{8}}习题21． 解方程组⎪⎩⎪⎨⎧=+=342y x x y 2．f(x)=2x 2+5x-8, 求f （1) f （3）f( 2)作出图像3．作出⎪⎪⎩⎪⎪⎨⎧>-≤≤-<=2 (22)0........20................sin )(32x x x x x x x x f 的图像，并求f(0.3)17．4．作出y=cosx,cos2x,cos3x 在[0,2π],标出坐标轴并带有三种不同颜色的图像 答案：Solvey 2 4x,x y 3 , x ,yx 1,y 2 , x 9,y 6f x _2x^25x 85x 2fPlot f x , x ,5,5Graphicsg 0.0.79895Plot g x , x ,5,5GraphicsPlot C os x ,Cos 2x ,Cos 3x , x ,0,2Pi ,AxesLabel "x","y" PlotStyle R GBColor 1,0,1 ,RGBColor 0,1,0 ,RGBColor 0,0,1GraphicsSurfaceGraphics实 验 三一、 实验目的1．学习用软件计算极限，判断函数的连续性；2．学习用软件计算一元函数的导数、多元函数的偏导数；3．学习用软件计算隐函数、参数式函数的导数及函数的微分、全微分； 4．学习用软件计算微分方程的解； 5．导数的简单应用． 二、 内容与步骤 1． 极限、连续：1）求一元函数的极限的命令格式是：Limit[f[x]，x －>x 0] 表示求函数x →x 0 的极限；Limit[f[x]，x －>x 0，Direction －>1] 表示求函数x →x 0-的极限（左极限）； Limit[f[x],x －>x 0，Direction －>-1] 表示求函数x →x 0+的极限（右极限）．2）若x 趋于无穷，即 x → ∞,则格式为Limit[f[x],x → ∞] x 趋于负无穷或正无穷格式为：Limit[f[x],x → - ∞] , Limit[f[x],x → + ∞]3）注：->∞ 也可由File → Palettes → BasicInput 中的符号输入 例1 求下列函数的极限：（1）443lim 24---→x x x x输入： Limit[4 ,4432→---x x x x ]输出：5 （2）xxx 3arctan lim+∞→输入：Limit[ArcTan[x]3x,x→+∞]输出：0 （3）x x x 2)4751(lim -+∞→ 输入：Limit[x x 2)4751(-+,x→∞] 输出：例2 求 x x e --→133lim 及x x e +-→133lim输入：Limit[，x→3,Direction→1]Limit[，x→3,Direction→-1] （e 为BasicInput 符号栏中的 ）输出：0 输出：∞还有一些函数没有极限,此时系统会进行相应的处理,返回一些特殊的结果.例3 求当x →0时，y ＝sinx1的极限． 解：输入：Limit[Sin[1/x]，x→0]输出：Interval[{－1，1}]上面这个例子表示当x →0时，函数sin x1在－1与1之间无穷震荡，所以没有确定的极限．例4 判定函数⎪⎩⎪⎨⎧≤+>=02302sin )(x x x xxx f 在 x=0点是否连续．解：输入：Limit[Sin[2x]x，x →0，Direction→-1] 右极限 输出：2输入： Limit[3x +2, x →0, Direction→1]] 左极限 输出：2输入：3x+2/.x→0 计算函数值 输出：2∴ 函数在x =0这一点连续． 2． 导数、偏导数1）一阶导数)(x f '的命令格式为： D[f ，x] （f 为函数表达式，x 为自变量） 2）n 阶导数)()(x f n 的命令格式为： D[f,{x,n}] （n 为导数的阶数） 3）用BasicInput 工具栏输入： （函数表达式变量∂ 此时的函数表达式可以是一元或多元函数，变量可有一个或多个，使用灵活．如输入： x x 3（求一元函数x 3对x 的一阶导数） 输出：8x输入： x,x x 3（求一元函数x 3对x 的二阶导数） 输出：输入： x x 3y 4x （求二元函数x 3y 4对x 的一阶偏导数）输出：3x 2y 8输入： y x 3y 4x （求二元函数x 3y 4对y 的一阶偏导数）输出：x 38输入： x,x x 3y 4x （求二元函数x 3y 4对x 的二阶偏导数）输出：6x y输入：x,y x 3y 4x （求二元函数x 3y 4先对x 后对y 的二阶偏导数）输出：3x 21输入： y,y x 3y 4x（求二元函数x 3y 4对y 的二阶偏导数）输出：例1 求下列显函数的导数：（1）3532x x y += （2）x e x y 2= （3）12ln +=x x y 解:（1）输入： D[2 x 5+3 x 3，x]输出： 9x 2+10x 4（2）输入：x x输出：2 xx（3）输入：x Log x2x输出：例2 求函数22ln ),(y x y x f +=的偏导数x f ∂∂，y f∂∂，y x f ∂∂∂2解： 输入：输出：x输入： 输出：x输入：输出：例3 求函数5-202Q Q R =，当Q=15和Q=20时的()20)15(R R ''及 解：求函数在一点x 0处的导数值,只需在输入表达式后面再继续输入“/.x→x 0”即可.方法一：输入：D[Q ,5Q -Q 202]/．Q →15 输出：14 （即 (15)14R '=）输入：D[Q ,5Q -Q 202]/．Q →20 输出：12． （即 (20)12R '=）方法二：（函数表达式）变量∂/.x->a输入：输出：输入：输出：12例4 求函数f （x ）＝sin ax cos bx 的一阶导数dx df ，并求ba x dxdf+=1．解: 输入：x S in a x Cos b x.x a输出：例5 求下列函数的高阶导数：（1）5x y = 求：y ''' （2）x xe y 3= 求：y '' （3）xx xy cos sin sin += 求：y ''解：（1）输入：D[x ^5，{x ，3}]输出：60x 2（2）输入：D[x Exp[3 x]，{x ，2}]输出：6 3x9输入：Simplify[D[Sin[x]/(Sin[x]+Cos[x]),{x,2}]] 输出：Cos[x]-Sin[x]Cos[x]+Sin[x]-22()()3． ㈠求隐函数的导数由方程F (x , y )=0 确定的函数)(x f y =，称为隐函数．方法：1）自定义一个导函数G[x_]对F (x ,y ）求导，但必须将变量y 输入成y[x],即y 是x 的函数．2）用Solve 函数将y [x]'解出即可．即先求导再解方程．例6 求由方程12222=+by a x 所确定的隐函数的导数．解：方法一输入：D[2222x y[x]+a b-1，x ]（先自定义一个导函数G[x]，这里表达式中的y 应写成y[x]）输出：22b [x]2y[x]y'a 2x + 输入：Solve[G[x]＝＝0，y'[x]]（用解方程Solve 命令，从导函数的方程G[x ]＝0 中解出y'[x]，这里方程必须使用双等号“＝＝” ）输出：{{y'[x] → －y[x]a xb 22}}方法二：利用工具栏与解方程语句：输入：输出：例7 已知方程0=-y xe xy 确定一个y 是x 的函数)(x f ，求 )(x f '． 解: 输入：Solve x xx y xy x0,y'输出：例8.设函数满足方程sin x x y ye +=0，求 ()y x '. 解：输入：Solve x x Sin y xy xx 0,y'输出：㈡ 求函数的微分、全微分求函数的微分dy ，其形式为Dt[f(x)].输出的表达式中所含的Dt[x],这里可以视为dx .求函数f (x, y )的全微分dz , 其形式为 Dt[f[x ，y]] 例9 求y ＝sin2x 的微分dy . 解: 输入：Dt[Sin[2x]]输出：2 Cos[2 x] Dt[x]例10 求函数x e x x y 23ln +=的微分dy . 解: 输入：Dt[x ∧ 3 Log[x]+Exp[2 x]]输出：2 e 2 x Dt[x]+x 2Dt[x]+3x 2Dt[x]Log[x] 再化简一下输入：Simplify[%]输出：Dt[x](2 e 2 x +x 2+3x 2 Log[x]) 即 dx x x x e dy x )ln 32(222++= 例11 求函数u xy z =23的全微分． 解: 输入：Dt[x y^2 z^3]输出：y 2 z 3 Dt[x] ＋ 2 x y z 3 Dt[y] ＋ 3 x y 2 z 2 Dt[z] ㈢ 参数式函数的求导形如 ⎩⎨⎧==)()(t x t y ψϕ 的函数为参数式函数，其导数 t t x x y y ''='． 其输入方式为：例12．设 ⎪⎩⎪⎨⎧==ta y ta x 33sin cos ，求 dx dy解： 输入： 输出：Ta例13.求椭圆⎩⎨⎧==tb y t a x sin cos 在 4π=t 处的导数解： 输入：输出：4．用 Mathematica 解微分方程其格式为： DSolve[微分方程，y[x]，x] 注意要将y 输入成y[x] 例14 解微分方程 ()()y x y x '+=1解: 输入：DSolve[y'[x]＋y[x]＝＝1，y[x]，x]输出：{{y[x]－>1＋xE C[1]}} 例15 求微分方程（x 2＋y 2）dx －xydy ＝0的通解．解: 输入：DSolve[(x^2＋y[x]^2)Dt[x]－x y[x] Dt[y[x]]＝＝0，y[x]，x]输出：{{y[x]－>－Sqrt[x 2 (C[1]＋2 Log[x])]}，y[x]－> Sqrt[x 2 (C[1]＋2 Log[x])]}}例16 求微分方程 ()x y xy '''+=212满足初始条件10==x y ，3'0==x y 的特解． 解: 输入：DSolve[{(x^2＋1)y''[x]＝＝2x y'[x], y[0]＝＝1, y'[0]＝＝3}, y[x]，x]输出：{{y[x]－> 1＋3 x ＋x 3}}5．导数的简单应用 （1）求函数的单调区间例17 求函数123+-=x x y 的单调区间解：函数的单调区间需要用到一阶导函数的图像、一阶导函数为零的驻点．输入:f x _ : x 32（建立函数） Plotf x ,f' x, x ,3,3 ,PlotStyle G rayLevel 0.01 ,Dashing0.01（画函数与导函数图像，其中虚线为导函数图像）输出：输入：Solve f ' x（求函数的驻点） 输出：观察图像，两个驻点将定义域分成三个区间，可看出函数在 ),32,(--∞),32(+∞内为增函数，在)32,32(-内为减函数．（2）求函数的极值例18 求函数21xxy +=的极值 解： 输入：g x _ : 1 Plotg x ,g' x, x ,3,3 ,PlotStyle G rayLevel 0.01 ,Dashing0.01输出：输入：Solve g ' x输出： x 1 , x（从图中可看出两驻点分别是极小值点和极大值点）输入： g输出：2（3）求极值的近似值 例19 求函数)2(cos 25)2(sin 222xx x y +=位于),0(π内的极值的近似值． 解：输入：Plot f x , x ,0,输出：观察图形，函数约在x=0.8、x=2.3处有极大值，在x=1.6处有极小值，可用命令FindMinimum 直接求极值的近似值，其格式为：FindMinimum[f[x]，{x ，x 0}]，求以x 0为初始点的局部极小值．FindMinimum 只可求极小值的近似值，欲求极大值的近似值，须将函数换成相反函数．输入： FindMinimum f x , x ,1 输出：1.94461, x 1.623即同时得到极小值1.94461和极小值点1.62391 输入：FindMinimum f x , x ,0输出： 3.73233,x 0.8641输入： FindMinimum f x , x ,2输出：2.95708,x 2.244即函数-y 的两个极小值和两个极小值点，从而得到函数y 的两个极大值和极大值点．（4）最大、最小值的应用例20 要制造一个容积为2,上端为半球形，下端为圆柱的粮仓，问：当圆柱的高和底半径为何值时，粮仓的表面积最小？ 解: 设粮仓的表面积为S ，圆柱的高为h>0, 底半径为r>0.由题意，粮仓的容积2=323421 r h r ππ⋅+，则 )31 1(2 322223r r r r h -=-=πππ ∴粮仓的表面积S=⋅r 2π)31 1(22r r -π+324 42122r r r ππ+=⋅. 输入： FindMinimum[4/r+2πr 2/3,{r,10}] 输出：{6.09295,{r →0.984745}}.（5）微分方程的应用例21 一质量为m 千克的物体从高处下落，所受空气阻力与速度成正比，设物体开始下落时（t=0）的速度为零，求物体下落速度与时间的函数关系v (t). 解：设物体所受空气阻力为f ，由题得 kv f =（k 为比例系数），下落时所受重力为mg ，根据牛顿第二定律有 v m ma kv mg f mg '==-=- 输入：DSolvem v' t m g k v t ,v 0 0 ,v t输出： 输入：Simplif输出：习题3 （每小题中括号内为该题答案）1. 求导数：（1）tan )2xy =[（2）1124=y （3）sin cos cos x y y y -+=220 求 .y '（4）cos()sin ,y xy x =223求 .y ' 326s i n [3x ]c o s [3x ]+y s i n [x y ][]2ycos[xy]-xy sin[xy]（5）,6x e y x ⋅= 求 )1()5(y [4051e] （6）x y z cos = ， 求 y x z z '' , [y Si ，Co]（7）xy e z =，求 y x z z '' , [，]（8）求 y e z x cos sin = 的二阶偏导数 [SinxCos x2Cos ySin xCos y Si，SinxCos x Si，Sin xCo] （9） 求函数 ⎪⎩⎪⎨⎧==-tt ey tex 的导数[（10）求函数 ⎩⎨⎧-=-=)cos 1()sin (t a y t t a x 的导数[2．求微分及全微分：（1）674335+-+=x x x y [7Dt x 12x 2Dt x 15x 4D] （2）32cot(ln )=x y ex[（3）xxx y ++=1sin ln [（4）y e z x sin = [ xCos y Dt yxDt x Si] （5）)cos(y x x z += [Cos x y Dt x x Dt x Dt ySinx] 3．解微分方程 （1）求微分方程yxdx dy -=的通解． [（2）求微分方程0)1(22=++dy x dx xy 的通解．[（3）求微分方程x yx y dx dy tan +=的通解． [{y xx ArcSin x}] （4）求微分方程x x x y dxdysin 2cot =-的通解．[y x x 2Sin x C 1 Sin] （5）求微分方程42x y y x =+'满足初始条件61)1(=y 的特解． [ y x]4．求下列极限 （1）1lim1-+→x xx [∞]（2）11lim31++-→x x x [31 ] （3）121lim +-∞→⎪⎭⎫⎝⎛+x x x x [ 2e ]（4）判断函数 ⎪⎩⎪⎨⎧>≤-=0,0,1)(2x x x x x f 在 0=x 处是否连续？ [ 不连续 ]实 验 四一、 实验目的1、 学习用软件计算不定积分；2、 学习用软件计算定积分、二重积分和广义积分；3、 定积分的简单应用，求平面面积和旋转体体积． 二、 内容与步骤 1．不定积分输入格式： BasicInput 符号栏中的符号注意：输出结果均不带积分常数． 例1 求下列不定积分 ⎰dx x5解：输入：x输出：6x 62． 定积分输入格式： BasicInput 符号栏中的符号例2 求下列定积分 ⎰-212 1dx x x解：输入：输出：3 例3 计算广义积分⎰+∞∞-+dx x 211解：输入：输出：例4 计算由抛物线2x y =和直线x y =所围成的平面图形的面积及该图形绕x 轴旋转一周所得的旋转体体积（表示出必要的步骤）解：（1）求交点输入：Solvey x,y x^2 , x ,输出： y 0,x 0 , y 1,x（2）作图 输入：Plotx ,x^2 , x ,2,输出：GraphClea Clea（3）定积分求面积输入：1 x x^2输出：6（4）定积分求体积输入：1x 2x输出：13．二重积分用Mathematica 计算二重积分的命令格式是：输入方法：先输入一元定积分符号，在中间积分变量的位置再输入一次定积分符号，作为累次积分的第一次积分．括号内为第一次积分，括号外为第二次积分． 例5 计算⎰⎰+1212x xxydy dx解: 输入：012xx 21x y输出：121 例6 计算⎰⎰+=Ddxdy y x I )(22， 其中D 由2 ,21,===y x y x y 围成解：①画平面区域图输入：输出：② Y - 型区域输入： 02y2yx 2y 2输出：3习题4（每小题中括号内为该题答案）求下列积分：(1)⎰-dx x x x)11(2[x+x](2)⎰+dx xsin 11[x 2Sin[]2x x Cos[]+Sin[]22] （3）⎰+dx x x 3)cos (sin [1(-9Cos[x]-Cos[3x]+9Sin[x]-Sin[3x])6]（4）2ln(sin )sin x dx x⎰[ -x-Cot[x]-Cot[x] Log[Sin[x]]] （5）⎰xdx x arctan 2 [（6）21sin 1cos x xdx x++⎰ [(7)⎰--1145dx xx [6](8)⎰∞--02dx xex []（9）[2](10) 计算由曲线282yx =-和x 轴所围成的平面图形的面积及该图形绕x 轴旋转一周所得的旋转体体积（表示出必要的步骤）[3，（11）计算由曲线21yx =-和22y x =+所围成的平面图形的面积及该图形绕x 轴旋转一周所得的旋转体体积（表示出必要的步骤） 过程： Solvey x 21,y 2x 2 , x ,y 0,x 1 , y 8,xPlotx 21,2x 2 , x ,2,Graph132x 2 x 213132x 2 ^2 x 21 ^2（12）计算二重积分3y Ded σ-⎰⎰，其中D 由20,1,x y y x ===围成过程：①画平面区域图② Y - 型区域1y 2 y3第二篇 线性代数实验 一一、实验目的6．掌握Mathmatica 中矩阵的输入方法； 7．学习用Mathmatica 软件计算行列式；8．学习用Mathmatica 软件进行矩阵的基本运算； 9．学习用Mathmatica 求逆矩阵及矩阵的秩.二、内容与步骤1．Mathmatica 中矩阵的输入方法 (1)按表的格式输入： （一般方法）}}{},{},{{1212222111211mn ,m ,m n ,,n ,,a a a a a a a a a A ，生成m 行n 列的矩阵(2)菜单输入：(适用于大矩阵) a)打开主菜单Input 项；b)单击Create Table/Matrix 项，输入行数及列数，填数即可。
数学分析实验-Mathematica 软件的应用
七、自定义函数
格式 f[x_]:=… f[x_,y_]:=… f[x_,y_,…]:= … Clear[f] 意义 定义一元函数 f(x) 定义二元函数 f(x,y) 定义多元函数 f(x,y,…) 取消对 f 的定义
(1) 出现在 f[x_]中的x是一类实体，用来表示函数定 义中的变量。x_可以用任何形式的参数来替代，右端定 义式中的x将会随之变化。看下面的例子：
1
三、函数与变量的命名规则
在Mathematica系统中，变量名和函数名遵从如下命 名规则： （1） 以字母开头的任意长度的字符或数字串； （2） 区分大小写； （3） 为与系统函数相区别，通常以小写字母开头； （4） 函数的形式为f[x].
四、变量赋值与变量替换
举例如下： x=5; （ 赋值） x=. （取消赋值） /. x ->3 （变量替换，计算时用3暂时替换x的值） Clear[x] （清除x的定义及其赋值）
H(n)-ln(n) 0.5777155816 0.5774656441 0.5773823223 0.5773406597 0.577315661 0.5772989959
7000
8000 9000 10000
9.430952520
9.564474984 9.682251076 9.787606036
5 4 3 2 1
20
40
60
80ቤተ መጻሕፍቲ ባይዱ
100
猜测1 调和数列的前n项和H(n)是发散数列，它的数值与 ln(n)+C 很接近。
猜测2 数列H(n)- ln(n)可能是收敛的。
Step5 用计算数据作印证 对充分大的n，计算H(n)-ln(n)的值： t2=Table[N[{n,H[n],Log[n],H[n]-Log[n]},10], {n,1000,10000,1000}]
Mathematica软件介绍及其应用
(1) 2<3<4<3
(2) 3!=2!=3!=4 (3) x>=y
2:逻辑运算
!p p||q 否 或
23
p&&q
且
(5.8>4.1)&&!(3.2==2.0)
集合运算
Mathematica软件中的集合是形如{a,b,c,….}的结构。 其中的元素有多种形式，可以是常数、变量、函数、方 程 图形等。在Mathematica软件中向量、矩阵都可以看作 为集合。
19
如果需要计算函数f[x]在某一点x=x0的值: (1) f[Pi/2] Mathematica软件中的函数可以用递归的方法进行定义 h[0]=1;h[n_]:=n*h[n-1];h[5] 在这里定义了h函数的初始值,以及一般的函数定义,需要 计算h[5]的具体数值
20
(1): x=2.0;u[x_]:=Sin[x];v[x_]=Siatica软件的算术运算是指加减乘除以及 乘方、开方运算。 例如： （1）3*(5-2)+4^(6-3)/2 （2）3^(1/3) （3）1/3+2/5
精确计算与近似计算： 在Mathematica软件中，若输入的数据是精确的 ，计算结果保留精确数字。若计算近似值，可以 采用下面的方法： （4）N[3^(1/3)] 函数N[x]表示x的近似值
3
本章内容
1. 2. 3. 4. 5. 7. 8. Mathematica软件概述 Mathematica软件的使用 数列极限 函数的最值与导数 定积分计算 级数 线性方程组和矩阵的应用
4
0 引言
Mathematica是什么 ?
Mathematica是什么? Mathematica是一位数学家，她能做一位物理系大 学毕业生能做的几乎所有的数学题，从简单的代数运算 直到解偏微分方程。 Mathematica不能做什么? Mathematica没有脑袋，所以她不能思考问题，也 不能主动解决问题。她只能按照人类的指令一步一步地 去解决某一特定问题。Mathematica是研究人员和工程 技术人员的好朋友和最佳助手之一。
数学实验-Mathematic应用实验
目录实验01 基本语法 (1)实验02 一元函数极限与导数运算 (8)实验03 一元函数微分学及其应用 (20)实验04 一元函数积分学及其应用 (33)实验05 绘制空间图形 (43)实验06 多元函数微分学 (61)实验07 多元函数积分学 (72)实验08 无穷级数及其应用 (82)实验09 常微分方程及其应用 (94)实验10 编程 (106)实验01 基本语法实验内容：Mathematica软件在数值计算、符号计算、编程方面的基本语法数据类型在Mathematic中，基本的数据类型有四种：整数、有理数、实数和复数。
整数与整数的计算结果是精确的整数或有理数。
例如2的100次方是一个31位的整数：ln[1]:=2^100Out[1]=1267650600228228229401496703205376有理数是由两个整数的比来组成如：In[2]:=12345/5555Out[2]=2469 1111实数有两种表示形式：(1)用数学表达式精确表示,例如：2(2)用浮点数近似表示，包括小数形式和指数形式。
例如：In[3]:=0.239998In[4]:=1.23*^12复数是由实部和虚部组成，实部和虚部可以用整数、实数、有理数表示。
用I表示虚数单位。
如：In[6]:=3+0.7I数值类型转换在Mathematica中的提供以下几个函数达到转换的目的：函数功能N[x] 将x转换成实数（有效位一般为6位）N[x,n] 将x转换成近似实数，精度为nRationalize[x] 给出x的有理数近似值Rationalize[x,dx] 给出x的有理数近似值，误差小于dx 举例:In[1]:=N[5/3,20]Out[1]=1.6666666666666666667In[2]:=Rationalize[%]Out[2]=5 3数学常数Mathematica定义了一些常见的数学常数，这些数学常数都是精确数。
常数意义Pi 表示π＝3.14159……E 自然对数的底e＝2.71828……Degree 1度，π/180弧度I 虚数单位iInfinity 无穷大∞数学常数表示精确值。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3、研究数列H(n)-ln(n)的收敛性
Step1 令C(n)=H(n)-ln(n)，通过图象观察其特性： Cup[n_] :=H[n]-Log[n] tup=Table[ {n, N[Cup[n],6]},{n,1,100}] ph4=ListPlot[tup,PlotStyle->RGBColor[0,0,1]]
我们将 lnx 的图象向上平移C个单位后再进行观察。 c1=H[100]-Log[100] ph3=Plot[Log[x]+c1,{x,1,100}] Show[ph1,ph3]
5 4 3 2 1
20
40
60
80
100
猜测1 调和数列的前n项和H(n)是发散数列，它的数值与 ln(n)+C 很接近。
1、运行和退出Mathematica系统。 2、界面介绍。
3、输入和计算表达式。 4、保存和打开文件。 5、使用帮助系统。
变量与函数
一、Mathematica中的数据类型
Mathematica系统中，数值分成四种类型：整数、有 理数、实数和复数。
整数型数据可以表示任意长度的精确整数，不受计 算机字长的限制。
Integer Rational
Real Complex
整数型 有理数(分数)型
(近似)实数型 复数型
精确运算 精确运算 近似计算
二、系统中的数据常数
Pi E Degree GoldenRatio Infinity I
Step4 与对数函数 y=lnx 作比较 ph2=Plot[Log[x],{x,1,100}] Show[ph1,ph2]
5 4 3 2 1
20
40
60
80
100
根据图象比较的结果可以看出，当n很大时，H(n)的 图象与ln(n)的图象非常相似，但它们大致相差一个常数。 这个常数约为
C=H(100)-ln100≈0.5822.
Step2 令c(n)=H(n)-ln(n+1)，通过图象观察其特性： Clow[n_] :=H[n]-Log[n+1] tlow=Table[ {n, N[Clow[n],6]},{n,1,100}] ph5=ListPlot[tlow,PlotStyle->RGBColor[1,0,0]]
Step3 比较C(n)和c(n)，在同一坐标系中作出它们的图象。 Show[ph4,ph5]
8.006367568 8.2940496401 8.517193191 8.699514748 8.853665428 8.9871968207 9.104979856 9.210340372
H(n)-ln(n) 0.5777155816
0.5774656441
0.5773823223 0.5773406597 0.577315661 0.5772989959 0.5772870918 0.5772781636 0.5772712194 0.5772656641
数学分析实验
－ Mathematica 软件的应用
数学实验实例一 调和数列研究
1、调和数列
自然数的倒数组成的数列 1, 1 , 1 ,, 1 , 23 n
称为调和数列。它的前n项和数列 n 1 记作H(n)。 k 1 k
2、提出问题：H(n)是否收敛？
我们借助于数学软件Mathematica 对H(n)的收敛性进 行观察。
如：264； 最大素数 21257787-1等。
有理型数据可以精确表示任意的既约分数，当两个 整数相除而又不能整除时，就用有理数型表示。
如：3/39.
实数型数据可以表示任意精度的近似实数。 如：Pi.
复数型数据可以表示复数，其实部和虚部可以是整 数型、有理型或实数型。
如：I2；Arg[1+I].
Mathematica中的数据类型
n 2 3
n
把这个极限值记为C，C ≈0.5772，称为欧拉（Euler） 常数。
数学实验工具
常见的数学工具软件： 1、Mathematica； 2、Matlab； 3、Maple； 4、MathCad.
我们主要使用Mathematica这一数学工具软件。
Mathematica介绍
Mathematica是一个功能强大的数学工具软件，具有 数值计算、符号演算、图象制作、公式编辑和编程等各 项功能。
0.62
0.58 0.56 0.54 0.52
20
40
60
80
100
通过观察可知如下事实： 1、C(n)是单调递减数列； 2、c(n)是单调递增数列； 3、c(n) ≤ C(n)； 4、c(n)，C(n)都是收敛数列，而且它们有相同的极限。
4、结论与证明
结论：
极限 lim (1 1 1 1 ln n) 存在。
猜测2 数列H(n)- ln(n)可能是收敛的。
Step5 用计算数据作印证 对充分大的n，计算H(n)-ln(n)的值： t2=Table[N[{n,H[n],Log[n],H[n]-Log[n]},10], {n,1000,10000,1000}]
可以得到如下的数据表：
n 1000
2000
3000 4000 5000 6000 7000 8000 9000 10000
…… 5.17738 5.18738 7.48547 9.78761 14.39273
Step3 根据数据表画出H(n)的图形 ph1=ListPlot[t]
5 4 3 2
20
40
60
80
100
Байду номын сангаас
通过对所得图象的观察和分析，我们发现它很接近对 数函数的图象。我们把它与对数函数 y=lnx 的图象一起比 较一下。
H(n) 7.485470861
8.178368104
8.583749890 8.871390300 9.094508853 9.276813744 9.430952520 9.564474984 9.682251076 9.787606036
ln(n) 6.907755279
7.6009024595
Step1 定义前n项和H(n) H[n_]:=Sum[1/k,{k,1,n}]
Step2 列出H(n)随n变化的数据表 t=Table[{n,N[H[n],6]},{n,1,100}]
n 1 2 3 4 …… 99 100 1000 10000 1000000
H[n] 1.00000 1.50000 1.83333 2.08333