资源描述:
《《数学实验》上机指导书》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、《数学实验》上机指导书实验题目实验一mathematica软件基本操作(一).实验类型:验证型(二).实验类别:基础实验(三).每组人数:1(四).实验要求:选修(五).实验学时:3个学时(三).实验目的:(1)掌握Mathematica软件的计算器功能;(2)学会使用Mathematica软件求各种类型方程(或方程组)的数值解和符号解;(3)通过本实验深刻理解极限概念;(4)学习并掌握利用Mathematica求极限的基本方法。(5)通过本实验加深理解积分理论中分割、近似、求和、取极限的思想方法;(6)学习并掌握二重积分及线性积分的计算方
2、法;(7)学习常用积分命令。(8)掌握求函数的导函数和偏导数方法;(9)学会使用Mathematica软件进行函数的幂级数展开。(四)【预备知识】(1)方程(或方程组)代数解法的基本理论,函数的零点,方程(或方程组)的解及数值解;(2)本实验所用命令:●用“==”连接两个代数表达式构成一个方程●求方程(组)的代数解:Solve[方程或方程组,变量或变量组]●求方程(组)的数值解:NSolve[方程或方程组,变量或变量组]●从初始值开始搜索方程或方程组的解:FindRoot[方程或方程组,变量或变量组初值]●在界定范围内搜索方程或方程组的解:
3、FindRoot[方程或方程组,变量或变量组范围]14●绘图命令:Plot[表达式,{变量,上限,下限},可选项]●微分方程求解命令:DSolve[微分方程,y[x],x](3)极限、左极限、右极限的概念;(4)本实验所用Mathematica有关命令:●Limit[expr,x->x0]求表达式在时的极限●Limit[expr,x->x0,Direction->1]求左极限●Limit[expr,x->x0,Direction->-1]求右极限(5)定积分的概念、几何意义,二重积分的概念、二重积分化为定积分的过程及其计算方法;(6)本实验
4、所用Mathematica有关命令:●无限积分:Integrate[f,x]●定积分:Integrate[f,{x,上限,下限}](7)函数的导函数、偏导数以及函数的幂级数展开式;(8)本实验所用的Mathematica函数提示:(a)求导数(或偏导数)●D[表达式F,x]求F对于变量x的导数;●D[表达式F,x1,x2,...]按顺序求F关于x1,x2,…的偏导数;●D[表达式F,{x,n}]求F对x的n阶导数。(b)幂级数展开●Series[表达式F,{x,x0,n}]求F关于变量x在x0的n阶泰勒展式。(五).实验内容(1)计算;。(
5、2)对于方程,试用Solve和Nsolve分别对它进行求解,并比较得到的结果,体会代数解即精确解与数值解的差别。(3)先观察函数的图形,然后选择一个初始点求解,并且根据图形确定在某个区间中搜索它的零点。14(4)求方程组的解,然后代入系数和常数项的一组初值,并求解。(5)求微分方程的通解。(6)用Mathematica软件计算下列极限:(1);(2);(3);(4);(5);(6);(7);(8);(9);(10);(11);(12)。(7)求函数的原函数;(8)求;(9)求;(10)求;(11)求。(12)求出被积函数F(x)=的原函数和
6、导函数,并画出被积函数、原函数和导函数的图形,试分辨出哪一条曲线属于哪个函数。14(13)对函数sinx在0点展开10阶和20阶,并以图形方式对比展开的结果和sinx的差别,并分析阶数高的展式对于原来函数的逼近程度是否优于阶数低的展式。(六)【实验操作】(1)学会N[]和expr//N的使用方法,并学会用Precision[]研究误差。In[1]:=N[Exp[3],12]In[2]:=Precision[%]In[3]:=546*54564//NIn[4]:=46545^45676//N(2)学会Solve[]和NSolve[]的使用方法
7、。In[5]:=p=x^4-2x^3-4x^2+3;Solve[p==0,x]In[6]:=NSolve[p==0,x](3)学会Clear[]和FindRoot[]的使用方法In[7]:=Clear[x]In[8]:=f=Sin[x]-Cos[x]In[9]:=Plot[f,{x,-4,4}]In[10]:=FindRoot[f,{x,1}]In[11]:=FindRoot[f,{x,{0,1}}](4)学会用Solve[]求解方程组。In[12]:=Clear[a1,a2,b1,b2,c1,c2]In[13]:=Solve[{a1*x+
8、b1*y==c1,a2*x+b2*y==c2},{x,y}](5)学会DSolve[]的使用方法In[14]:=DSolve[y''[x]+3y'[x]+2y[x]==Exp[x