mathematica 数学实验报告 实验一
数学实验 Mathematic实验一 一元函数的图形

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

13
5、计算的精度
14
6.常用的数学函数
函数 Sqrt[x]
Exp[x]
说
明
函数
Sinh[x] Cosh[x]
说
明
平方根函数 指数函数
双曲正弦函数 双曲余弦函数
Log[x]
Log[a ,x] Sin[x] Cos[x] Tan[x] Cot[x] ArcSin[x]
自然对数函数
以a为底对数函数 正弦函数 余弦函数 正切函数 余切函数 反正弦函数
4
参考教材
由我校编写的 《高等数学实验》
5
实验一 Mathematica软件操作(1)
实验目的 掌握Mathematica的基本操作 掌握Mathematica中的基本运算 掌握常用的数学函数和构造函数的方法
6
实验内容
1、Mathematica的启动
在 Windows 环境下安装好Mathematica。点击开始菜单 中的Mathematica的图标,稍停片刻则显示工作屏幕。
2
数学软件简介
我们在学习数学的过程中,要花费大量的时间 用于推导公式、画图、计算数据。随着电子计算机 技术的不断发展,人们已经开发出一些能够帮助处
理和解决数学问题的软件系统。
常见的通用数学软件包括:Mathematica、
Matlab和Maple,其中Matlab以数值计算见长,
Mathematica和Maple以符号运算、公式推导见长。
相等 不相等 大于 大于或等于 小于 小于等于 都相等 都不相等
24
逻辑运算符 用一个关系式只能表示一个判定条件,要表示 几个判定条件的组合,必须用逻辑运算符将关系表 达式组织在一起,我们称表示判定条件的表达式为 逻辑表达式。常用的逻辑运算符有
mathematica数学实验报告

mathematica数学实验报告本次实验使用Mathematica进行数学建模实验,主要包括以下内容:三角函数、极限和导数、积分和微分方程。
一、三角函数1. 三角函数的绘制使用Mathematica的Plot函数绘制正弦函数和余弦函数的图像。
代码:Plot[{Sin[x], Cos[x]}, {x, -2 Pi, 2 Pi},PlotStyle -> {Blue, Red}, PlotTheme -> "Web"]结果:在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的基本使用

实验一 熟悉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个。
02项目一Mathematica 一元函数微分学

实验一 一元函数微分学实验1 一元函数的图形(基础实验)实验目的 通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Mathematica 作平面曲线图性的方法与技巧.基本命令1. 在平面直角坐标系中作一元函数图形的命令Plot: Plot[f[x],{x,min,max},选项]Plot 有很多选项(Options), 可满足作图时的种种需要, 例如,输入Plot[x^2,{x,-1,1},AspectRatio->1,PlotStyle->RGBColor[1,0,0],PlotPoints->30]则输出2x y =在区间11≤≤-x 上的图形. 其中选项AspectRatio->1使图形的高与宽之比为1. 如 果不输入这个选项, 则命令默认图形的高宽比为黄金分割值. 而选项PlotStyle->RGBColor[1,0,0] 使曲线采用某种颜色. 方括号内的三个数分别取0与1之间. 选项PlotPoints->30令计算机描点作 图时在每个单位长度内取30个点, 增加这个选项会使图形更加精细.Plot 命令也可以在同一个坐标系内作出几个函数的图形, 只要用集合的形式{f1[x],f2[x],…} 代替f[x].2.利用曲线参数方程作出曲线的命令ParametricPlot:ParametricPlot[{g[t],h[t]},{t,min,max},选项]其中)(),(t h y t g x ==是曲线的参数方程. 例如,输入ParametricPlot[{Cos[t],Sin[t]},{t,0,2 Pi},AspectRatio->1]则输出单位圆t y t x sin ,cos ==的图形.3. 利用极坐标方程作图的命令PolarPlot如果想利用曲线的极坐标方程作图, 则要先打开作图软件包. 输入<<Graphics`Graphics`执行以后, 可使用PolarPlot 命令作图. 其基本格式为PolarPlot[r[t],{t,min,max},选项]例如曲线的极坐标方程为,3cos 3t r =要作出它的图形. 输入PolarPlot[3 Cos[3 t], {t,0,2 Pi}]便得到了一条三叶玫瑰线.4. 隐函数作图命令ImplicitPlot这里同样要先打开作图软件包, 输入<<Graphics\ImplicitPlot.m命令ImplicitPlot 的基本格式为ImplicitPlot[隐函数方程, 自变量的范围, 作图选项]例如方程22222)(y x y x -=+确定了y 是x 的隐函数. 为了作出它的图形, 输入ImplicitPlot[(x^2+y^2)^2==x^2-y^2,{x,-1,1}]输出图形是一条双纽线.5. 定义分段函数的命令Which 命令Which 的基本格式为Which[测试条件1, 取值1, 测试条件2, 取值2,…]例如, 输入w[x_]=Which[x<0,-x,x>=0,x^2]虽然输出的形式与输入没有改变, 但已经定义好了分段函数:⎪⎩⎪⎨⎧≥<-=0,0,)(2x x x x x w现在可以对分段函数)(x w 求函数值, 也可作出函数)(x w 的图形.实验举例初等函数的图形例1.1 作出指数函数x e y =和对数函数x y ln =的图形. 输入命令Plot[Exp[x],{x,-2,2}]则输出指数函数x e y =的图形.输入命令Plot[Log[x],{x,0.001,5},PlotRange->{{0,5},{-2.5,2.5}},AspectRatio->1]则输出对数函数x y ln =的图形.注①:x 的, 第二组数{-2.5,2.5}是描述y 的.注②:有时要使图形的x 轴和y 轴的长度单位相等, 需要同时使用PlotRange 和AspectRatio 两个选项. 本例中输出的对数函数的图形的两个坐标轴的长度单位就是相等的. 例1.2 作出函数x y sin =和x y csc =的图形观察其周期性和变化趋势. 为了比较, 我们把它们的图形放在一个坐标系中. 输入命令Plot[{Sin[x],Csc[x]},{x,-2 Pi,2 Pi},PlotRange->{-2 Pi,2 Pi},PlotStyle->{GrayLevel[0],GrayLeve1[0.5]}, AspectRatio->1]注:.例1.3 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势.输入命令Plot[{Tan[x],Cot[x]},{x,-2 Pi,2 Pi},PlotRange->{-2 Pi,2 Pi},PlotStyle->{GrayLeve1[0],GrayLeve1[0.5]},AspectRatio->1]例 1.4 将函数x y ,sin =数的图形间的关系.输入命令p1=Plot[ArcSin[x],{x,-1,1}];p2=Plot[Sin[x],{x,-Pi/2,Pi/2},PlotStyle->GrayLeve1[0.5]]; px=Plot[x,{x,-Pi/2,Pi/2},PlotStyle->Dashing[{0.01}]];Show[p1,p2,px,PlotRange->{{-Pi/2,Pi/2},{-Pi/2,Pi/2}},AspectRatio->1]注 Show[…]命令把称为p1,p2和px 的三个图形叠加在一起显示. 选项PlotStyle->Dashing[{0.01}]使曲线的线型是虚线.例1.5 (教材 例1.1) 给定函数24325555)(x x x x x x f +++++=(a) 画出)(x f 在区间]4,4[-上的图形;(b) 画出区间]4,4[-上)(x f 与)()sin(x f x 的图形. 输入命令f[x_]=(5+x^2+x^3+x^4)/(5+5x+5x^2);g1=Plot[f[x],{x,-4,4},PlotStyle->RGBColor[1,0,0]];则输出)(x f 在区间]4,4[-上的图形.输入命令g2=Plot[Sin[x]f[x],{x,-4,4},PlotStyle->RGBColor[0,1,0]]; Show[g1,g2];则输出区间]4,4[-上)(x f 与)()sin(x f x 的图形.注: Show[…]命令把称为例1.6 在区间]1,1[-画出函数xy 1sin=的图形. 输入命令Plot[Sin[1/x],{x,-1,1}];则输出所求图形,从图中可以看到函数x y 1sin=在0=x 附近来回震荡.二维参数方程作图例1.7 作出以参数方程)20(sin ,cos 2π≤≤==t t y t x 所表示的曲线的图形. 输入命令ParametricPlot[{2 Cos[t],Sin[t]},{t,0,2 Pi},AspectRatio->Automatic]则可以观察到这是一个椭圆.注 在ParametricPlot 命令中选项AspectRatio->Automatic 与选项AspectRatio->1是等效的.例 1.8分别作出星形线)20(s i n 2,c o s233π≤≤==t t y t x 和摆线),sin (2t t x -= )40)(cos 1(2π≤≤-=t t y 的图形.输入命令ParametricPlot[{2 Cos[t]^3,2 Sin[t]^3},{t,0,2 Pi},AspectRatio->Automatic] ParametricPlot[{2*(t-Sin[t]),2*(1-Cos[t])},{t,0,4 Pi},AspectRatio->Automatic]则可以分别得到星形线和摆线的图形.例1.9 画出参数方程⎩⎨⎧==tt t y tt t x 3cos sin )(5cos cos )(的图形:输入命令ParametricPlot[{Cos[5 t]Cos[t],Sin[t]Cos[3t]},{t,0,Pi}, AspectRatio->Automatic];则分别输出所求图形.例1.10 (教材 例1.2) 画出以下参数方程的图形.(1) ⎪⎪⎩⎪⎪⎨⎧+⎪⎭⎫ ⎝⎛-=+⎪⎭⎫⎝⎛-=tt t y tt t x sin 7511sin 5)(cos 7511cos 5)( (2) ⎩⎨⎧-+=-+=t t t t y t t t t x sin )4cos 2sin 1()(cos )4cos 2sin 1()(分别输入以下命令:ParametricPlot[{5Cos[-11/5t]+7Cos[t],5Sin[-11/5t]+7Sin[t]},{t,0,10Pi},AspectRatio->Automatic];ParametricPlot[(1+Sin[t]-2 Cos[4*t])*{Cos[t],Sin[t]},{t,0,2*Pi},AspectRatio->Automatic,Axes->None]; 则分别输出所求图形.例1.11 作出极坐标方程为)cos 1(2t r -=的曲线的图形.曲线用极坐标方程表示时, 容易将其转化为参数方程. 故也可用命令ParametricPlot[…]来作极坐标方程表示的图形.输入命令r[t_]=2*(1-Cos[t]);ParametricPlot[{r[t]*Cos[t],r[t]*Sin[t]},{t,0,2 Pi},AspectRatio->1]可以观察到一条心脏线.极坐标方程作图例1.12 (教材 例1.3) 作出极坐标方程为10/t e r =的对数螺线的图形.输入命令<<Graphics`执行以后再输入PolarPlot[Exp[t/10],{t,0,6 Pi}]则输出为对数螺线的图形.隐函数作图例1.13 (教材 例1.4) 作出由方程xy y x 333=+所确定的隐函数的图形(笛卡儿叶形线). 输入命令<<Graphics\ImplicitPlot.m执行以后再输入ImplicitPlot[x^3+y^3==3x*y,{x,-3,3}]输出为笛卡儿叶形线的图形.分段函数作图例1.14 分别作出取整函数][x y =和函数][x x y -=的图形. 输入命令Plot[Floor[x],{x,-4,4}]可以观察到取整函数][x y =的图形是一条阶梯形曲线.输入命令Plot[x-Floor[x],{x,-4,4}]得到函数][x xy -=的图形, 这是锯齿形曲线(注意: 它是周期为1的周期函数.)例1.15 作出符号函数x y sgn =的图形. 输入命令Plot[Sign[x],{x,-2,2}]就得到符号函数的图形. 点0=x 是它的跳跃间断点.g[x_]: = -1/; x<0; g[x_]: = 0/; x=0; g[x_]: = 1/; x>0; Plot[g[x],{x,-2,2}]便得到上面符号函数的图形. 其中组合符号“/;”的后面给出前面表达式的适用条件例1.16 (教材 例1.5) 作出分段函数⎪⎩⎪⎨⎧>≤=0,,0,cos )(x e x x x h x 的图形.输入命令h[x_]:=Which[x<=0,Cos[x],x>0,Exp[x]] Plot[h[x],{x,-4,4}]则输出所求图形.注:一般分段函数也可在组合符号“/;”的后面来给出前面表达式的适用条件.例1.17 (教材 例1.6) 作出分段函数⎪⎩⎪⎨⎧=≠=0,00,1sin )(2x x xx x f 的图形. 输入命令f[x_]:=x^2Sin[1/x]/;x!=0;f[x_]:=0/; x=0;Plot[f[x],{x,-1,1}];则输出所求图形.函数性质的研究例1.18 研究函数)3(log 3)(35x e x x f x -++=在区间]2,2[-上图形的特征. 输入命令Plot[x^5+3E^x+Log[3,3-x],{x,-2,2}];则输出所求图形. 由图形容易看出, 从左到右, 图形渐渐上升. 因而是增函数.例1.19 判断函数x x x f ππ2cos 2sin )(+=是否为周期函数.任选一个较大的范围, 如取]4,4[-, 在此区间上画出函数)(x f 的图形如图所示.Plot[Sin[2Pi x]+Cos[2Pi x],{x,-4,4}];可以看出函数的图形以某一宽度以单位重复出现.例 1.20 判断函数133)(23+++==x x x x f y 的反函数的存在性. 若存在, 求反函数的表达式, 并画出起图形.先解方程,13323+++=x x x y 求x . 输入命令Solve[y==x^3+3x^2+3x+1,x]; 因此, 所求反函数为.13x y +-= 再输入命令Plot[-1+x^(1/3),{x,-3,3}];则输出反函数在区间]3,3[-内的图形.注:若一个函数满足: 一个y 对应着一个x , 则其反函数一定存在,且在表达式中将y 换成常量求解x , 即将所的表达式中y 换成x , x 换成y 即得到反函数的表达式.作函数图形的动画例1.21 制作函数cx sin 的图形动画, 观察参数c 对函数图形的影响. 输入命令.Do[Plot[Sin[c x],{x,-Pi,Pi},PlotRange->{-1,1}],{c,1,4,1/3}];则输出图形动画.例1.22 (教材 例1.7) 作出函数cx x x f sin )(2+=的图形动画,观察参数c 对函数图形的影响. 输入命令Do[Plot[x^2+Sin[c x],{x,-3,3},PlotRange->{-1,5}],{c,1,5,1/3}];则输出所求动画图形.实验习题1. 把正切函数x tan 和反正切函数x arctan 的图形及其水平渐近线2/,2/ππ=-=y y 和直线 x y =用不同的线型画在同一个坐标系内.2. 作出双曲正切函数x tanh 的图形.3. 输入以下命令Plot[{Sin[x],Sin[2 x],Sin[3 x]},{x,0,2 Pi}, PlotStyle->{RGBColor[1,0,0], RGBColor[0,1,0],RGBColor[0,0,1]}]理解选项的含义.4. 为观察复合函数的情况,分别输入以下命令:Plot[Sqrt[1+x^2],{x,-6,6},PlotStyle->{Dashing[{0.02,0.01}]}] Plot[Sin[Cos[Sin[x]]],{x,-Pi,Pi}]Plot[Sin[Tan[x]]-Tan[Sin[x]]/x^2,{x,-5,5}] Plot[{E^x,ArcTan[x],E^ArcTan[x]},{x,-5,5}]5. 观察函数的叠加, 输入以下命令:a1=Plot[x,{x,-5,5},PlotStyle->{RGBColor[0,1,0]}]a2=Plot[2 Sin[x],{x,-5,5},PlotStyle->{RGBColor[1,1,0]}] a3=Plot[x+2 Sin[x],{x,-5,5},PlotStyle->{RGBColor[1,0,0]}] Show[a1,a2,a3]6. 分别用ParametricPlot 和PolarPlot 两种命令, 作出五叶玫瑰线θ5sin 4=r 的图形.7. 用ImplicitPlot 命令作出椭圆322+=+xy y x 的图形.8. 选择以下命令的一部分输入, 欣赏和研究极坐标作图命令输出的图形.PolarPlot[Cos[t/2],{t,0,4 Pi}] PolarPlot[1-2 Sin[5 t],{t,0,2 Pi}] PolarPlot[Cos[t/4],{t,0,8 Pi}] PolarPlot[t*Cos[t],{t,0,8,Pi}] PolarPlot[t^(-3/2),{t,0,8 Pi}] PolarPlot[2 Cos[3 t],{t,0,Pi}] PolarPlot[1-2 Sin[t],{t,0,2 PI}] PolarPlot[4-3 Cos[t],{t,0,2 Pi}]PolarPlot[Sin[3 t]+Sin[2 t]^2,{t,0,2 Pi}] PolarPlot[3 Sin[2 t],{t,0,2 Pi}] PolarPlot[4 Sin[4 t],{t,0,2 Pi}]PolarPlot[Cos[2 t]+Cos[4 t]^2,{t,0,2 Pi}] PolarPlot[Cos[2 t]+Cos[3 t]^2,{t,0,2 Pi}]PolarPlot[Cos[4 t]+Cos[4 t]^2,{t,0,2 Pi},PlotRange->All]实验2 极限与连续(基础实验)实验目的 通过计算与作图, 从直观上揭示极限的本质,加深对极限概念的理解. 掌握用 Mathematica 画散点图, 以及计算极限的方法. 深入理解函数连续的概念,熟悉几种间断点的图形 特征,理解闭区间上连续函数的几个重要性质.基本命令1.画散点图的命令ListPlot:ListPlot[{{x1,y1},{x2,y2},…{xn,yn}},选项]或者ListPlot[{y1,y2,…yn},选项]前一形式的命令,在坐标平面上绘制点列),(,),,(),,(2211n n y x y x y x 的散点图;后一形式的命令, 默认自变量i x 依次取正整数,,,2,1n 作出点列为),(,),,2(),,1(21n y n y y 的散点图. 命令ListPlot 的选项主要有两个:(1) PlotJoined->True, 要求用折线将散点连接起来; (2) PlotStyle->PointSize[0.02], 表示散点的大小. 2.产生集合或者数表的命令Table:命令Table 产生一个数表或者一个集合. 例如, 输入Table[j^2,{j,1,6}]则产生前6个正整数的平方组成的数表{1,4,9,16,25,36}.3.连加求和的命令Sum:命令Sum 大致相当于求和的数学符号∑. 例如, 输入Sum[1/i,{i,100}]//N执行后得到1001312111++++ 的近似值.与Sum 类似的还有连乘求积的命令Product. 4. 求函数多次自复合的命令Nest: 例如, 输入Nest[Sin,x,3]则输出将正弦函数自己复合3次的函数Sin[Sin[Sin[x]]]5.求极限的命令Limit: 其基本格式为Limit[f[x],x->a]其中f(x)是数列或者函数的表达式, x->a 是自变量的变化趋势. 如果自变量趋向于无穷, 用 x->Infinity.对于单侧极限, 通过命令Limit 的选项Direction 表示自变量的变化方向. 求右极限, 0+→a x 时, 用Limit[f[x],x->a,Direction->-1]; 求左极限, 0-→a x 时, 用Limit[f[x],x->a,Direction->+1]; 求+∞→x 时的极限, 用Limit[f[x],x->Infinity,Direction->+1]; 求-∞→x 时的极限, 用Limit[f[x],x->Infinity,Direction->-1]。
mathematica实验报告

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` 。
数学实验报告

2、画出函数在-3 ( x (源自3 的图形, 且为红色。3、将函数绘制绘制在第一象限内, 区间任选(要求图形高宽比为1)。
4、先画出函数 在-3x25的图形,再显示在平面区域[5,12][5,10]部分的图形以观察局部特征。
5、在同一坐标系中画出三个函数 在-2x2
的图形,并给坐标横轴和纵轴分别标记为x和y。
6、将函数 的图形作在同一坐标系内,并观察直接函数和反函数的图形间的关系(可以选择让图形呈现不同颜色以便观察)。
输入以下命令,观察图形叠加,说明选项意义。
a1=Plot[x,{x,-5,5},PlotStyle-> {RGBColor[0,1,0]}]
a2=Plot[2 Sin[x],{x,-5,5},PlotStyle->{RGBColor[0,0,1]}]
组别
实验小组成员
实验名称
基础实验1:一元函数的图形
成绩
试验序号: 01日期: 2010年月日
实验目的
通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Mathematica作平面曲线图性的方法与技巧。
试验所用版本: Mathematica 5.0
a3=Plot[x+2 Sin[x],{x,-5,5},PlotStyle->{RGBColor[1,0,0]}]
Show[a1,a2,a3]
7、作出分段函数 的图形。(选作)
8、作出隐函数 的图形。(选作)
实验过程记录(含基本步骤,主要程序清单及异常情况记录等)
1.
异常情况记录:
实验结果报告与实验总结:
教师评价
数学软件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= = %
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
数学实验报告
实
验
一
数学与统计学院
信息与计算科学(1)班
郝玉霞
0107
数学实验一
一、实验名:微积分基础
二、实验目的:学习使用Mathematica的一些基本功能来验证或观察得出微积分学的几个基本理论。
三、实验环境:学校机房,工具:计算机,软件:Mathematica。
四、实验的基本理论和方法:利用Mathematica作图来验证高中数学知识与大学数学内容。
五、实验的内容和步骤及结果
内容一、验证定积分
dt
t
s
x
⎰=
1
1
与自然对数
x
b ln=
是相等的。
步骤1、作积分
dt
t
s
x
⎰=
1
1
的图象;
语句:S[x_]:=NIntegrate[1/t,{t,1,x}] Plot[S[x],{x,,10}]
实验结果如下:
图1
dt
t
s
x
⎰=
1
1
的图象
步骤2、作自然对数
x
b ln=
的图象
语句:Plot[Log[x],{x,,10}]实验结果如下:
图2
x
b ln=
的图象
步骤3、在同一坐标系下作以上两函数的图象语句:Plot[{Log[x],S[x]},{x,,10}]
实验结果如下:
图3
dt
t
s
x
⎰=
1
1
和
x
b ln=
的图象
内容二、观察级数与无穷乘积的一些基本规律。
(1)在同一坐标系里作出函数和它的Taylor展开式的前几项构成的多项式函数,,的图象,观察这些多项式函数的图象向的图像逼近的情况。
语句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]}]
实验结果如下:
图4和它的二阶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]}]
实验结果如下:
图5和它的三阶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]}]
实验结果如下:
图6和它的四阶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]}]
实验结果如下:
图7和它的五阶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}]
实验结果如下:
图8 和它的二、三、四、五阶Taylor展开式的图象
(2)分别取n=10,20,100,画出函数在区间[-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]}]
实验结果如下:
图9 n=10时,的图像
语句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]}]
实验结果如下:
图10 n=20时,的图像
语句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]}]
实验结果如下:
图11 n=100时,的图像
(3)分别取5,15,100,,在同一坐标系里作出函数与在区间[-2π,2π]上的图像。
语句1:
p[x_,n_]:=x*Product[1-x^2/(k^2Pi^2),{k,1,n}]
Plot[{Sin[x],p[x,5] },{x,-2Pi,2Pi}]
实验结果如下:
图12 n=5时,与的图像
语句2:
p[x_,n_]:=x*Product[1-x^2/(k^2Pi^2),{k,1,n}]
Plot[{Sin[x],p[x,15] },{x,-2Pi,2Pi}]
实验结果如下:
图13 n=15时,与的图像
语句3:
p[x_,n_]:=x*Product[1-x^2/(k^2Pi^2),{k,1,n}]
Plot[{Sin[x],p[x,100] },{x,-2Pi,2Pi}]
实验结果如下:
图14 n=100时,与的图像
六、实验结果分析
内容一、图1、图2分别作出了定积分
dt
t
s
x
⎰=
1
1
与自然对数
x
b ln=
的图象,
大致看来这两幅图是一样的;由图3在同一坐标系里作出以上两函数的图象,可
以看出这两幅图是完全重合的,由此足以证明:定积分
dt
t
s
x
⎰=
1
1
与自然对数
x
b ln=
是相等的,这与之前我们得出的结论是完全一致的。
内容二、(1)图4、5、6、7分别作出函数和它的二、三、四、五阶Taylor展开式的图象,图8作出了同一坐标系里函数和它的二、三、四阶Taylor展开式的图象,经比较可知,奇数阶的更接近正弦函数;(2)图9、10、11分别作出n=10,20,100时,函数的图像,经观察可知,当n→∞时,这个函数趋向于分段函数;(3)图12、13、14分别作出n=5,15,100时,在同一坐标系里函数与在区间[-2π,2π]上的图像,观察知当n增加时的图像向的图像逼近,且两个函数在x=0处的导数相同,在任何有限的区间上,当n→∞时函数逼近。