mathematica 数学实验报告材料 实验一

天水师范学院数学与统计学院
实验报告
实验项目名称一元函数的图形
所属课程名称数学实验
实验类型微积分实验
实验日期2011.9.21
班级
学号
姓名
成绩
附录1：源程序
附录2：实验报告填写说明
1．实验项目名称：要求与实验教学大纲一致.
2．实验目的：目的要明确，要抓住重点，符合实验教学大纲要求.
3．实验原理：简要说明本实验项目所涉及的理论知识.
4．实验环境：实验用的软、硬件环境.
5．实验方案（思路、步骤和方法等）：这是实验报告极其重要的内容.概括整个实验过程.
对于验证性实验，要写明依据何种原理、操作方法进行实验，要写明需要经过哪几个步骤来实现其操作.对于设计性和综合性实验，在上述内容基础上还应该画出流程图、设计思路和设计方法，再配以相应的文字说明.对于创新性实验，应注明其创新点、特色. 6．实验过程（实验中涉及的记录、数据、分析）：写明具体实验方案的具体实施步骤，包括实验过程中的记录、数据和相应的分析.
7．实验结论（结果）：根据实验过程中得到的结果，做出结论.
8．实验小结：本次实验心得体会、思考和建议.
9．指导教师评语及成绩：指导教师依据学生的实际报告内容，给出本次实验报告的评价.。

mathematica数学实验报告本次实验使用Mathematica进行数学建模实验，主要包括以下内容：三角函数、极限和导数、积分和微分方程。

一、三角函数1. 三角函数的绘制使用Mathematica的Plot函数绘制正弦函数和余弦函数的图像。

代码：Plot[{Sin[x], Cos[x]}, {x, -2 Pi, 2 Pi},PlotStyle -> {Blue, Red}, PlotTheme -> "Web"]结果：![trigonometric_functions.png](2. 求三角函数的值使用Mathematica的N函数计算正弦函数和余弦函数在不同角度下的取值。

代码：N[Sin[Pi/6]]N[Cos[Pi/6]]N[Sin[Pi]]N[Cos[Pi]]结果：0.50.8660251.22465*10^-16-1.二、极限和导数1. 求函数的极限使用Mathematica的Limit函数计算函数x^2/(4-x)在x趋近于4时的极限。

代码：Limit[x^2/(4 - x), x -> 4]结果：82. 求函数的导数使用Mathematica的D函数计算函数x^3 - 3x的导数。

代码：D[x^3 - 3x, x]结果：3 x^2 - 3三、积分和微分方程1. 求定积分使用Mathematica的Integrate函数计算函数e^x * cos(x)在0到π/2之间的定积分。

代码：Integrate[E^x * Cos[x], {x, 0, Pi/2}]结果：1/2 (1 + E^(π/2))2. 解微分方程使用Mathematica的DSolve函数求解微分方程y''(x) + 4y(x) = 0。

代码：DSolve[y''[x] + 4 y[x] == 0, y[x], x]结果：y[x] -> C[1] Cos[2 x] + C[2] Sin[2 x]本次实验使用Mathematica进行数学建模实验，主要包括三角函数的绘制、求三角函数的值，函数的极限、导数，积分和微分方程等内容。

实验一 熟悉Mathematica 的基本使用1、 写出圆周率π的前50位小数，看看它的前100位，1000位小数，能不能发现什么规律? In[ ]:= N[Pi,50]Out[ ]= 3.1415926535897932384626433832795028841971693993751In[ ]:= N[π,100]Out[ ]= 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068In[ ]:= N[Pi,1000]（结果略）说明：100位，1000位的都可以不记录，看看有没有规律就可以。

2、 第1234个素数是什么？15485863是素数吗？、In[ ]:= Prime[1234]Out[ ]= 10061In[ ]:= PrimeQ[15485863]Out[ ]= True (*不要看成或写作Ture*)3、26(1)π+位于哪两个整数之间？In[ ]:= 26Floor[(1)]π+ (*取整数*)Out[ ]= 1649234In[ ]:= 26Round[(1)]π+ (*四舍五入取整数,不小于x 的最小整数*)Out[ ]= 16492354、构造一个表，由不超过100的所有素数组成；In[ ]:= Table[Prime[i],{i,1,25}]Out[ ]= {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97} 但是，这个25是怎么来的呢？可以先观察一下：In[ ]:= Table[If[PrimeQ[i],i],{i,1,100}]Out[ ]= {Null,2,3,Null,5,Null,7,Null,Null,Null,11,Null,13,Null,Null,Null,17,Null ,19,Null,Null,Null,23,Null,Null,Null,Null,Null,29,Null,31,Null,Null,Null,Null,Nu ll,37,Null,Null,Null,41,Null,43,Null,Null,Null,47,Null,Null,Null,Null,Null,53,Nu ll,Null,Null,Null,Null,59,Null,61,Null,Null,Null,Null,Null,67,Null,Null,Null,71,Null,73,Null,Null,Null,Null,Null,79,Null,Null,Null,83,Null,Null,Null,Null,Null,89,Null,Null,Null,Null,Null,Null,Null,97,Null,Null,Null}或者：In[ ]:= Table[If[Prime[i]<100,Prime[i]],{i,1,100}]Out[ ]= {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97, Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null, Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null, Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null, Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null, Null,Null,Null,Null,Null,Null,Null,Null,Null,Null,Null}经过观察，即可得到25个。

mathematica实验报告Mathematica 实验报告一、实验目的本实验旨在深入了解和掌握 Mathematica 软件的基本功能和操作方法，通过实际的案例和问题解决，提升运用 Mathematica 进行数学计算、数据分析、图形绘制以及编程的能力。

二、实验环境操作系统：Windows 10Mathematica 版本：121三、实验内容与步骤（一）数学计算1、基本运算在 Mathematica 中，直接输入数学表达式进行计算，例如：计算 2＋ 3 4 的结果，输入｀2 ＋ 3 4` ，得到结果 14。

2、函数计算使用内置函数进行复杂的数学运算，如计算正弦函数｀SinPi ／ 6`的值，结果为 05。

（二）数据分析1、数据导入通过｀Import` 函数导入外部数据文件，如 CSV 格式的数据文件。

假设我们有一个名为｀datacsv` 的文件，包含两列数据｀x` 和｀y` ，使用｀data ＝ Import"datacsv"｀即可将数据导入。

2、数据处理对导入的数据进行处理，如计算平均值、方差等统计量。

可以使用｀Meandata` 计算平均值，｀Variancedata` 计算方差。

（三）图形绘制1、二维图形绘制简单的函数图形，如｀PlotSinx, ｛x, 0, 2 Pi}｀绘制正弦函数在｀0` 到｀2 Pi` 区间的图形。

2、三维图形绘制三维图形，如｀Plot3Dx^2 ＋ y^2, ｛x, －2, 2}，｛y, －2, 2}｀绘制一个抛物面。

（四）编程实践1、定义函数使用｀Function` 关键字定义自己的函数，例如定义一个计算阶乘的函数｀factorialn_ ：＝ Ifn ＝＝ 0, 1, n factorialn 1` 。

2、循环结构使用｀For` 循环和｀While` 循环实现重复操作，例如使用｀For`循环计算 1 到 10 的和，｀sum ＝ 0; Fori ＝ 1, i ＜＝ 10, i+＋， sum ＋＝ i; sum` 。

mathematica-数学实验报告-实验一————————————————————————————————作者：————————————————————————————————日期：数学实验报告实验一数学与统计学院信息与计算科学(1)班郝玉霞201171020107数学实验一一、 实验名：微积分基础二、实验目的：学习使用Mathematica 的一些基本功能来验证或观察得出微积分 学的几个基本理论。

三、实验环境：学校机房，工具：计算机,软件：Mathematica 。

四、实验的基本理论和方法：利用Mathematica 作图来验证高中数学知识与大学数学内容。

五、实验的内容和步骤及结果内容一、验证定积分dtt s x⎰=11与自然对数x b ln =是相等的。

步骤1、作积分dtt s x⎰=11的图象； 语句：S[x_]:=NIntegrate[1/t,{t,1,x}] Plot[S[x],{x,0.1,10}]实验结果如下：2468102112图1dt t s x⎰=11的图象步骤2、作自然对数x b ln =的图象语句：Plot[Log[x],{x,0.1,10}] 实验结果如下：2468102112图 2x b ln =的图象步骤3、在同一坐标系下作以上两函数的图象 语句：Plot[{Log[x],S[x]},{x,0.1,10}] 实验结果如下：2468102112图3dtt s x⎰=11和x b ln =的图象 内容二、观察级数与无穷乘积的一些基本规律。

（1）在同一坐标系里作出函数x y sin =和它的Taylor 展开式的前几项构成的多项式函数3!3xx y -=，!5!353x x x y +-=，⋅⋅⋅的图象，观察这些多项式函数的图象向x y sin =的图像逼近的情况。

语句1：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]Plot[{Sin[x],s[x,2]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}] 实验结果如下：6422464224图4x y sin =和它的二阶Taylor 展开式的图象语句2：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]Plot[{Sin[x],s[x,3]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,1,1]}] 实验结果如下：6422463211234图5x y sin =和它的三阶Taylor 展开式的图象语句3：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]Plot[{Sin[x],s[x,4]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,1,0]}] 实验结果如下：642246321123图6x y sin =和它的四阶Taylor 展开式的图象语句4：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]Plot[{Sin[x],s[x,5]},{x,-2Pi,2Pi},PlotStyle->{RGB[1,0,0]}] 实验结果如下：642246321123图7x y sin =和它的五阶Taylor 展开式的图象语句5：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}] Plot[{Sin[x],s[x,2],s[x,3],s[x,4],s[x,5] },{x,-2Pi,2Pi}] 实验结果如下： 642246224图8xy sin=和它的二、三、四、五阶Taylor展开式的图象（2）分别取n=10，20，100，画出函数xkkynk)12sin(1211--=∑=在区间[-3π，3π]上的图像，当n→∞时，这个函数趋向于什么函数？语句1：f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}]Plot[f[x,10],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}]实验结果如下：6422460.50.5图9 n=10时，xkkynk)12sin(1211--=∑=的图像语句2：f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}]Plot[f[x,20],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}]实验结果如下：6422460.50.5图10 n=20时，xk k y nk )12sin(1211--=∑=的图像语句3：f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}]Plot[f[x,100],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}] 实验结果如下：6422460.50.5图11 n=100时，xk k y nk )12sin(1211--=∑=的图像（3）分别取5,15,100,，在同一坐标系里作出函数x x f sin )(=与∏=-⋅=nk k x x x p 1222)1()(π在区间[-2π，2π]上的图像。

一、 实验目的：学习使用Mathematical 的一些基本功能来验证或观察得出微积分的几个基本结论.二、 实验的内容1、调和级数自然数倒数组成的数列称为调和级数.我们把它的前n 项和11nk k=∑记作：H （n ）.2、泰勒级数在同一坐标系内作出区间[]ππ,-∈x 上正弦函数x y sin =及多项式函数!7!5!3,1206,6753533x x x x y x x x y x x y -+-=+-=-=的图像。

观察这些多项式函数的图像逼近正弦曲线的情况。

3、求自然对数e在同一坐标系中画出下面三个函数的图像：xx y 101011⎪⎭⎫ ⎝⎛+=，1101011+⎪⎭⎫ ⎝⎛+=x x y ，e y =，观察当x 增大时图像的走向。

4、正弦函数的叠加在自变量区间]2,2[ππ-∈x 上使用Mathematica 8.0绘制出函数nnxx x y sin 22sin sin +⋅⋅⋅++=的图像。

观察当n 增大时图像形状的变化趋势。

三、 实验过程1.步骤1将坐标为（n ，H （n ））（n=1,2,…,100）的点依次连接成光滑曲线.观察曲线的形状.它与什么函数的图像形状类似. H[n_]:=NSum[1/k,{k,1,n}]; t =Table[{n,H[n],{n,1,100}}];pic1=ListPlot[t]好像与对数函数图像类似步骤2为了验证这些店连成的曲线是否接近对数函数的图像，可以将对数函数的图像与上述p1画在同一坐标系中进行比较.pic2=Plot[Log[x],{x,1,100},PlotStyle->{RGBColor[0,0,1]}];Show[pic1,pic2, AxesLabel->{“x”, “y”}]步骤3c=H[100]-Log[100]pic3=Plot[Log[x]+c,{x,1,100}, PlotStyle->{RGBColor[0,0,1]}]Show[pic1,pic2,pic3, AxesLabel->{“x ”, “y ”}]有图像可以看出，pic1与pic2基本重合。

实验报告1 函数与极限院系 班号姓名学号成绩一、实验内容函数图形的显示，极限的运算，最值的计算．二、预期目标1.熟悉Mathematica 软件的基本操作．2.掌握函数与极限的有关操作命令．3.学会利用Mathematica 软件对函数进行分析研究．三、常用命令1． 作图命令： 2． 参数作图命令： 3． 图形显示命令： 4． 求极限命令： 5． 求极值名命令：四、练习内容1.画出下列函数的图形： （1） y=cos3x作图命令：（2） f （x ）=x 5+3e x+log 3（3－x ） x ∈[－2,2]作图命令：（3）⎪⎩⎪⎨⎧=+=ty t t x 2sin作图命令：（4）⎪⎩⎪⎨⎧==ty t x 33sin cos t ∈[0,2π]作图命令：2.求下列极限：（1）110002lim+∞>-n nn （2）113)2(3)2(lim ++∞>-+-+-n n n n n （3）35)3)(2)(1(limnn n n n +++∞>- （4）3522lim -+>-x x x （5）131lim +->-x x x（6）x e xx arctan lim -+∞>-（7）156182221lim +-->-x x x x （8）)sin 11sin (lim x x x x x -∞>-计算结果：（1）（2）（3）（4）（5）（6）（7）（8）3.讨论函数f（x）=2x3－6x2－18x+7在点2.8附近的极值．命令：结果：五、思考与提高1.怎样对隐函数的图形进行显示?2.怎样利用软件对函数极限存在性进行判断?3.如何利用软件对函数的连续性进行判断?4.如何求函数的最大（小）值?实验报告2 微分及其应用院系 班号姓名学号成绩一、实验内容导数的运算法则,复合函数求导法及参数方程求导法等．二、预期目标1.进一步理解导数及其几何应用．2.学习Mathematica 的求导命令与求导法．三、常用命令1．求导命令： 2．求微分命令： 3．隐函数求导命令： 4． 参数方程所确定的函数求导命令：四、练习内容1.求下列函数的导数： （1）x y 2ln 1+=求导命令： 求导结果： （2）21121xx y +++=求导命令：求导结果：（3） y=cos 2（cos2x ）求导命令： 求导结果：（4）y=2x/lnx求导命令： 求导结果： （5）y=ln[ln(lnt)]求导命令：求导结果： （6）xxy arccos arcsin =求导命令： 求导结果：（7）y=e arcsinx +arctane x求导命令： 求导结果：（8）xey 1sin 2-=求导命令：求导结果：2.求下列函数的二阶导数：（1） y=tanx 计算结果：y ” =（2）y=（1+x 2）arctanx 计算结果：y ” =（3）y=xtanx －cscx 计算结果：y ” =（4）y=21ln （x －1）－21ln （x+1） 计算结果：y ” = （5）⎪⎩⎪⎨⎧-==21arcsin ty t x 计算结果：y ” =（6）⎪⎩⎪⎨⎧==tb y t a x sin cos计算结果：y ” =3.求下列方程所确定的隐函数y=y （x ）的导数xyd d ： （1） sin （xy ）+cosy=0 计算结果：xyd d =（2）arctan x y =ln 22y x + 计算结果：xyd d =（3）x y =y x计算结果：xy d d =4.验证参数方程⎪⎩⎪⎨⎧==te y t e x tt cos sin 所确定的函数y 满足关系：)d d (2)(d d 222y x y x y x xy -=+ 程序：五、思考与提高1. 如何利用函数的导数判定函数的单调性、凹凸性?2.如何求由方程F(x,y,z)=0确定的隐函数z=z(x,y)的偏导数?实验报告3 积分及其应用院系 班号姓名学号成绩一、实验内容一元函数的不定积分与定积分二、预期目标1.加深理解积分理论中分割、近似、求和、取极限的思想方法．2.学习求积分的命令Integrate 与NIntegrate ．3.熟悉Mathematica 软件在积分运算的重要作用．三、常用命令1．求和命令： 2．求不定积分命令： 3．求定积分命令：四、练习内容1.求下列函数的一个原函数：（1）41x （2）212x +积分命令： 积分命令： 积分结果： 积分结果：（3）)1()1(22x x x ++ （4）4211xx -+ 积分命令： 积分命令： 积分结果：积分结果：（5）x x 22sin cos 2cos （6）xxe e +1积分命令： 积分命令： 积分结果： 积分结果：（7）)tan 1(cos 12x x + （8）x e x 32 积分命令： 积分命令： 积分结果：积分结果：（9）x cos 1+ （10）)34cos()23sin(+⋅+x x 积分命令： 积分命令： 积分结果： 积分结果：2.计算下列定积分：（1）⎰2/6/2d cos ππx x （2）⎰+-4/02sin 12sin 1πxxdx计算结果： 计算结果：（3）⎰-2/0d cos 351πx x（4）⎰30d cot arc x x x计算结果： 计算结果：（5）⎰---222d 11x x （6）t t e td cos 2/02⎰π计算结果：计算结果：（7）⎰+12/3d 1x xx （8）⎰π222d sin x x 计算结果： 计算结果：3.计算下列积分,并求其结果关于变量x 的导数：（1）⎰+02d 1x t t （2）⎰-xt t te 0d 2积分结果： 积分结果： 关于x 的导数：关于x 的导数：（3）⎰0sin 2d )cos(x t t （4）⎰+203d 11x t t 积分结果： 积分结果： 关于x 的导数： 关于x 的导数： 4.判定广义积分⎰∞++12)1(1x x dx 及⎰--2022)2(x exdx 的敛散性，收敛时计算出积分值． ⎰∞++12)1(1x x dx ⎰--2022)2(x exdx 程序： 程序： 结果： 结果： 5．求积分⎰-102)43(x x dx 具有6位、10位有效数的近似值． 命令： 五、思考与提高1. Mathematica 系统对分段函数的积分能否进行自动处理?2.《高等数学》中所学的积分换元法在软件系统里如何应用?3.怎样用Mathematica 中动画来演示定积分的定义?实验报告4 三角级数院系 班号 姓名 学号 成绩一、实验内容级数敛散性的判定．二、预期目标1.掌握级数的展开与求和命令．2. 学习使用Mathematica 进行级数敛散性的判定．三、常用命令1．求taylor 展式命令：四、练习内容1.求下列泰勒展开式，并在同一坐标系下画出函数图形及展开式图形． （1） ln （1+x ） 在x0=0点的8阶Taylor 展开． 程序：（2） P （x ）=x 4－5x 3+x 2－3x+4 在x0=4点的4阶Taylor 展开． 程序：（3） f （x ）=x1在x0=－1点的n 阶Taylor 展开． 程序：2.求下列级数的和函数：（1）∑∞=--112121n n x n （2）∑∞=+1)1(1n n x n n （－1≤x ≤1） 命令：命令：结果： 结果：（3）∑∞=-+112)1(n n x n n 命令： 结果：3.判定下列级数的敛散性：（1）∑∞=12n n n（2）∑∞=++13211n nn 结论：结论：（3）∑∞=1!2n nn （4）∑∞=1)(sin n n n n n π结论： 结论：（5）∑∞=+112tann n n π（6）∑∞=12)!(n nn n 结论： 结论：4.判定下列级数是否收敛，收敛时请指出是绝对收敛，还是条件收敛? （1）∑∞=---11121)1(n n n （2）∑∞=+-122)1(n n n 结论：结论：（3）∑∞=--1ln )1(n nn n （4）∑∞=12sin n n na （a 为常数） 结论： 结论：五、思考与提高用判别法可以判别级数的敛散性，但在实际应用时，往往要使用其和，原则上可用Sum 语句求和，但许多数项级数仅仅使Sum 语句求不出其和，而另－Mathematica 命令NSum 可与判别结果一起用来求出其近似值，问：是否对任一级数均可用NSum 来求其近似值?试以∑∞=-1)1(n n 为例观察．实验报告5 空间解析几何院系 班号 姓名 学号 成绩一、实验内容空间图形的显示，简单动画的制作．二、预期目标1.能正确显示空间图形．2.能用Mathematica 制作简单的动画．三、常用命令1．三维作图命令： 2．参数方程作图命令（三维曲线）： （曲面）： 3．动画命令：四、练习内容1.显示下列函数图形：（1） 椭球面⎪⎩⎪⎨⎧===v z v u y v u x cos 3sin sin 5sin cos 2，),0(),2,0(ππ∈∈v u作图命令：（2） 椭球抛物面⎪⎩⎪⎨⎧===23sin 3cos 3u z v u y v u x ，其中)2,0(),2,0(π∈∈v u作图命令：（3） 双曲抛物面⎪⎩⎪⎨⎧-===3/)(22v u z v y u x ，其中)4,4(,-∈v u作图命令：（4） 圆柱螺线⎪⎩⎪⎨⎧===t z t y t x 4sin 34cos 3，其中)5,0(∈t作图命令：3. 制作平面振动动画（利用函数y x y x f 3sin 2cos ),(=，其中x,y 均属于（-1,1））．程序：五、思考与提高用参量函数与直接函数显示图形有什么区别？比较谁更容易作出图形？实验报告6 多元微分学院系 班号 姓名 学号 成绩一、 实验题目隐函数的导数，函数的偏导数，函数的极值．二、 预期目标1.求隐函数的导数．2.求函数的偏导数和全微分．3.用微分知识求函数的极值．三、常用命令1．求偏导命令： 2．求全微分命令： 3．解方程（组）命令：四、练习内容1. 设xx xy tan =，求dxdy ． 命令： 结果：2. 设),(y x f z =由方程02=+--z xye z e所确定，求xz ∂∂． 程序： 结果：3. 设0ln 2=--xyz xy xz 确定函数),(y x f z =，求z 的全微分． 程序： 结果：4.求下列函数的偏导数：（1）yz x z y y x y y x z ∂∂∂∂-=,sin cos sin cos 2323，求结果：（2）yzx z v u y v u x y x z ∂∂∂∂+=-==,2,22，求，其中结果：4. 求函数22y x z +=在平面x+y=1上的最小值．程序： 结果：五、思考与提高1． 隐函数的二阶(偏)导数应如何求？2．函数的方向导数怎样求？实验报告7 多元积分学院系 班号 姓名 学号 成绩一、 实验题目空间立体体积和表面积．二、 预期目标1.用Mathematica 软件计算重积分．2.能解决空间立体体积和表面积的计算．三、常用命令1．求二重积分命令：四、练习内容1.计算下列重积分：（1）⎰⎰1D dxdy x y，其中D 1是由y=2x ，y=x ，x=4，x=2所围成的区域 ． 积分命令：计算结果：（2）⎰⎰+2)(22D dxdy y x，其中D 2是由y=x ，y=x+2，y=2，y=6所围成的区域．积分命令：计算结果：（3）⎰⎰++3)1ln(22D dxdy y x ，其中D 3：0,0,122≥≥≤+y x y x ． 积分命令：计算结果：（4）⎰⎰⎰Ω++3)(z y x dxdydz，其中Ω：21≤≤x ，21≤≤y ，21≤≤z ． 积分命令：计算结果：（5）⎰⎰⎰Ω++222zy x dxdydz ，其中Ω是由222z y x =+及1=z 所围成的区域． 积分命令：计算结果：2.求抛物面x y x y 2,==及平面z=0，z+x=6所围成的物体（密度为1）的质量．程序： 结果：五、思考与练习1.在实验步骤1中{x,0,1}与{y,2*x,x^2+1}能不能交换次序？为什么？2.在重积分中，如果可以用换元法，也可以用Integrate直接积分时，用哪一种方法好，为什么？3.曲线积分和曲面积分如何计算？实验报告8 常微分方程院系 班号 姓名 学号 成绩一、 实验题目常微分方程(组)的精确解．二、 预期目标1.求一阶常微分方程的精确解．2.求解简单的微分方程组和高阶方程．三、常用命令1．求解微分方程命令： 2．求解微分方程组命令： 3．求微分方程数值解命令：四、练习内容1. 求x y x y tan cos '2=+的通解．命令：结果：2. 求13232=-+y xx dx dy ，且满足y(1)=0的特解． 命令：结果：3. 求⎩⎨⎧=--=++03'5'y x y e y x x t ，满足⎩⎨⎧==0)0(1)0(y x 的特解．命令：结果：五、思考与提高如果遇见无法直接用DSolve 求解的常微分方程，如22112'x y y +=+，怎么办？院系 班号 姓名 学号 成绩一、实验内容矩阵的运算（加法、数乘、乘法、转置、逆）二、预期目标熟悉Mathematica 软件中关于矩阵运算的各种命令．三、常用命令1．矩阵显示命令： 2．求矩阵转置命令： 3．求逆矩阵命令： 4．求矩阵和差命令： 5．求矩阵数乘命令： 6．求矩阵乘命令：四、思考与练习已知矩阵 ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡------------=031948118763812654286174116470561091143A⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡------=503642237253619129113281510551201187851697236421B求：（1） Ａ＇； （2）Ａ-１；（３）Ａ*B ．（1）求Ａ＇的命令： （2）求Ａ-１的命令：Ａ＇= Ａ-１=（3）求Ａ*B 的命令：Ａ*B =（请用矩阵形式表示计算结果）院系 班号 姓名 学号 成绩一、实验内容对矩阵作各种变化,初等变换．二、预期目标1.复习并掌握矩阵初等变换的方法．2.掌握Mathematic 软件中关于矩阵初等变换的相关命令．三、常用命令1．取矩阵元素命令： 2．取矩阵的子矩阵命令： 3．求矩阵维数命令：四、练习内容1．已知矩阵;302 150311101⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=A(1)求A 的行向量组a 1,a 2,a 3, 以及列向量组b 1,b 2,b 3,b 4程序：(2)求A 的一，三，五行，二，三，四列交叉点上的元素做出子矩阵．程序：结果： 2.判断下列向量组是否线性相关(1)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=1211a ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=1302a ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=3123a 程序：结论：(2)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=1121a ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=1112a ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=1353a 程序：结论：实验报告11 行列式运算院系 班号 姓名 学号 成绩一、实验内容行列式的计算.二、预期目标1. 复习矩阵的行列式的求法，矩阵初等变换方法.2. 熟悉Mathematic 软件中关于求一个矩阵的行列式的命令把矩阵进行初等变换的命令以及与其相关的其它命令.三、常用命令1．求矩阵行列式命令：四、练习内容 1.求行列式βααββααββα+++100001000(共10阶)的值计算结果：2.利用克莱姆法则求解下列线性方程组(1)⎪⎪⎪⎩⎪⎪⎪⎨⎧-=-+--=++++=-+-+=+-+--=-+-+3322224343238243214225432154321543215422153321x x x x x x x x x x x x x x x x x x x x x x x x x程序：结果：(2) ⎪⎪⎪⎩⎪⎪⎪⎨⎧=+=++=++=++=+150650650651655454343232121x x x x x x x x x x x x x结果：2.已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=876174114A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=614475914B验证：｜Ａ×Ｂ｜=｜Ａ｜×｜Ｂ｜.程序：实验报告12 求解方程组院系 班号 姓名 学号 成绩一、实验题目求AX=B 的通解．二、实验目的通过本实验,使学生认识到虽然在《线性代数》中求AX=B 的通解比较繁，但在Mathematica 软件中却是比较简单的． 三、常用命令1．矩阵化简命令： 2．解线性方程组命令： 3．求AX=0的基础解系命令：四、练习内容1.求下列矩阵的秩:(1)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=253414312311112A 命令： 结果： (2)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=189513411314311B 命令： 结果：2.解下列线性方程组:(1) ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----512111211121⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡5514321x x x x 程序：结果：(2) ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----1111145212142121⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡3/10324321x x x x结果：实验报告13 特征值、特征向量院系 班号 姓名 学号 成绩一、实验题目计算已知矩阵的特征值和属于每一个特征值的特征向量.二、实验目的1.复习线代中的特征值与特征向量的求法.2.比较Mathematic 软件与普通方法的异同之处.三、常用命令1．求矩阵特征值命令： 2．求矩阵特征向量命令：四、练习内容求出下列矩阵的全部特征值与特征向量:1.⎥⎦⎤⎢⎣⎡-=00a a A ； 程序：结果：2.⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=001010100B ； 结果：3. ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----=1111111111111111C . 结果：实验报告14 离散型随机变量及其相关知识院系班号姓名学号成绩一、实验内容排列、组合的计算，几种离散型随机变量的产生及其相关内容.二、预期目标1.熟练掌握Mathematical软件的基本操作．2.熟悉与排列、组合、离散型随机变量有关的操作命令．3.掌握利用Mathematical软件处理简单的概率问题．三、常用命令1．（双）阶乘运算命令：2．组合数的计算命令：3．排列数的计算命令：4．服从二项分布的随机变量的生成命令：5．服从泊松分布的随机变量命令：6．将离散型随机变量的分布律拟合为函数的命令：四、练习内容1.计算下列结果（1）15！（2）15！！命令：命令：结果：结果：2.计算下列排列组合式的结果（1）P510（2）C510（3）!6!4!2!12⨯⨯命令：命令：命令：结果：结果：结果：3．生成以n=20，p=0.3为参数服从二项分布的随机变量bdist，将其分布律图形显示．程序：4.生成以p=0.4为参数服从几何分布的随机变量bdist，将其分布律图形显示．程序：5.生成以p=0.2为参数服从泊松分布的随机变量bdist，将其分布律图形显示．程序：五、思考与提高1.试分析几种离散型随机变量分布律的最值情况？2.怎样求解离散型随机变量有关的事件概率？实验报告15 连续型随机变量及其相关知识院系班号姓名学号成绩一、实验内容连续型随机变量的产生及其相关内容．二、预期目标1.熟练掌握几种连续型随机变量产生的有关操作命令．2.掌握利用软件对连续型随机变量进行分析的方法．3.掌握利用软件处理简单的概率问题．三、常用命令1．服从均匀分布的随机变量的生成命令：2．服从正态分布的随机变量的生成命令：3．服从t分布的随机变量的生成命令：4．服从χ2分布的随机变量的生成命令：5．服从F分布的随机变量的生成命令：6．求连续型随机变量的概率密度函数的命令：7．求连续型随机变量的分布函数的命令：四、练习内容1.生成以μ=10.05和σ=0.06为参数服从正态分布的连续型随机变量gdist及其概率密度函数、分布函数并图形显示；试求概率P{9.9<gdist<10.17}．程序：2.生成以a=0，b=1为参数服从柯西分布的连续型随机变量gdist及其概率密度函数、分布函数并图形显示．程序：3.生成以n1=4，n2=8为自由度服从F分布的连续型随机变量gdist及其概率密度函数、分布函数并图形显示．程序：4．生成以α=1，β=3为服从威布尔分布的连续型随机变量gdist及其概率密度函数、分布函数并图形显示．程序：五、思考与提高怎样利用软件对随机变量函数的分布进行分析，以及有关事件概率的求解？实验报告16 数字特征院系班号姓名学号成绩一、实验内容随机变量的数字特征及其相关内容．二、预期目标1.熟练掌握随机变量数字特征的有关操作命令．2.掌握利用软件对随机变量的特征函数(母函数)的求解．3.掌握利用软件处理简单的概率问题．三、常用命令1．求随机变量的期望的命令：2．求随机变量的方差的命令：3．求随机变量的标准差的命令：4．求随机变量的函数的方差的命令：5．求数据的协方差的命令：6．求数据的协方差矩阵的命令：7．求两随机变量的相关系数的命令：8．求两数据的相关系数矩阵的命令：四、练习内容1.（1）求以λ为参数服从泊松分布的随机变量的数学期望和方差．（2）求上述随机变量函数（f(x)=x2）的数学期望．（3）求服从参数λ=0.1的指数分布的随机变量的特征函数．程序：结果：2.(1)若样本data={16.5,13.8, 16.6, 15.7, 16.0, 16.4, 15.3},求样本均值、调和均值和中位数．结果：(2)若二维总体的样本data={{1612, 7627}, {1598, 6954},{1804, 8365},{1752, 9469}, {2067, 6410}, {2365, 10327},{1646, 7320}, {1579, 8196}, {1880, 9709}, {1773, 10370},{1712, 7749}, {1932, 6818}, {1820, 9307}, {1900, 6457},{1587, 8309}, {2208, 9559}, {1487, 6255}},求样本均值向量、中位数向量、方差向量和协方差矩阵．程序：结果：实验报告17 估计理论院系班号姓名学号成绩一、实验内容单个和两个总体均值、方差的估计．二、预期目标1.熟练掌握估计理论的相关操作命令．2.熟练掌握利用Mathematical软件对总体均值、方差进行估计．3.掌握利用Mathematical软件处理估计理论相关的实际问题．三、常用命令1．求总体均值的无偏估计的命令：2．求总体方差的无偏估计的命令：3．求总体方差的极大似然估计的命令：4．求单个总体均值的区间估计的命令：5．求两个总体均值之差的区间估计的命令：6．求单个总体方差的区间估计的命令：7．求两个总体方差之比的区间估计的命令：四、练习内容1.若样本data1={4506,4508,4499,4503,4504,4510,4497,4512,4514, 4505,4493,4496,4506,4502,4509,4496}来自正态总体，方差未知(置信度为0.95)：求出总体均值、方差的置信区间．程序：结果：2.若样本data2={4507,4507,4497,4506,4503,4511,4498,4510,4514,4510,4493,4491,4507,4501,4510,4495}来自正态总体，设置信度为0.95：（1）若data1与data2的总体方差都未知，均值之差的置信区间；程序：结果：（2）若data1与data2的总体方差都为40，均值之差的置信区间．程序：结果：3.data1与data2的总体方差之比值的置信区间（置信度为0.95）．程序：结果：实验报告18 假设检验院系班号姓名学号成绩一、实验内容对单个和两个总体均值、方差的假设检验．二、预期目标1.熟练掌握假设检验有关的操作命令．2.熟练掌握利用Mathematical软件对单个总体均值、方差的假设检验．3.掌握利用Mathematical软件对两个总体均值、方差有关的假设检验．三、常用命令1．求单个总体对均值的假设检验的命令：2．求两个总体对均值之差的假设检验的命令：3．求单个总体方差的假设检验的命令：4．求两个总体方差之比值的假设检验的命令：5．求标准正态分布有关概率的命令：6．求t分布有关概率的命令：7．求χ2分布有关概率的命令：8．求F分布有关概率的命令：四、练习内容设有甲、乙两种安眠药,比较其治疗效果．X表示服用甲药后睡眠时间延长时数，Y 表示服用乙药后睡眠时间延长时数，独立观察20个病人，其中10人服用甲药，另10人服用乙药，数据如下表：试就下列两种情况分析这两种药物的疗效有无显著性的差异．（显著性水平为0.05）（1）X与Y的方差相同；（2）X与Y的方差不同．程序：程序：结论：结论：五、思考与提高针对概率论与数理统计中左边、右边假设检验的问题,如何利用软件加以实现？31。

数学实验报告
实验序号：1 日期：2016年 3 月22 日
{ x= x(t) , y= y(t) }的图形命令:
ParametricPlot[{x[t], y[t]} , 要绘图形的参数t的范围, 选择项参数
AspectRatio->Automaic, 表示显示的图形高度与宽度比由
部算法根据函数图形的大小确定; AspectRatio->1, 表示显示的图形高度与宽度比
选项参数名称: Axes
图形是否有坐标轴
该参数的取值为True和None。

该选项参数的默认值为
表示显示的图形有坐标轴；Axes-> None, 表示显示的图形没有坐标轴。

Frame
平面图形是否加框
该参数的取值为True和False。

该选项参数只用于平面图形
False. Frame-> True, 表示显示的图形有框；Frame-> False, 表示显示的图形没有框。

: PlotLabel
是否设置图形名称标记
该参数取值为"字符串"和None, 默认值为None。

PlotLabel-> , PlotLabel->"Figure 1",使显示的图形上标出符号Figure 1
选项参数名称: PlotRange
设置图形的范围PlotRange->Automatic，表示用Mathematica。

数学实验一 Mathematica 基础11. 实验内容① 熟悉软件界面、各个菜单的作用；② 熟练常用运算、函数调用方法、程序运行方法； ③ 熟练使用联机帮助； 2.实验目的熟悉Mathematica 软件的基本操作，熟练使用常用函数。

3.实验要求1） 计算下列表达式的值：2132+，152.5+，100002，π的200位近似值N[]。

362.719950631168807583848837421626835850838234968318861924548520089498529438830221946631919961684036194597899331129423209124271556491349413781117593785932096323957855730046793794526765246551266059895520550086918193311542508608460618104685509074866089624888090489894838009253941633257850621568309473902556912388065225096643874441046759871… N[Pi,200]3.14159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038202） 计算：，tan 0.2，2log 10，200！Sin0.987766 Tan[0.2] 0.20271数学实验一 Mathematica 基础1Log[2,10]L200!7886578673647905035523632139321850622951359776871732632947425332443594499634033429203042840119846239041772121389196388302576427902426371050619266249528299311134628572707633172373969889439224456214516642402540332918641312274282948532775242424075739032403212574055795686602260319041703240623517008587961789222227896237038973747200000000000000000000000000000000000000000000000003） 利用//.计算表达式的值。

目录实验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 无穷大∞数学常数表示精确值。