数学实验报告.docx
- 文档编号:30202041
- 上传时间:2023-08-07
- 格式:DOCX
- 页数:11
- 大小:76.66KB
数学实验报告.docx
《数学实验报告.docx》由会员分享,可在线阅读,更多相关《数学实验报告.docx(11页珍藏版)》请在冰豆网上搜索。
数学实验报告
集团标准化工作小组#Q8QGGQT-GX8G08Q8-GNQGJ8-MHHGN#
数学实验报告
高等数学数学实验报告
实验人员:
院(系)____学号____姓名___
实验地点:
计算机中心机房
实验一空间曲线与曲面的绘制
一、实验题目:
(实验习题1-2)
利用参数方程作图,做出由下列曲面所围成的立体图形:
(1)
及xOy平面;
(2)
及
二、实验目的和意义
1、利用数学软件Mathematica绘制三维图形来观察空间曲线和空间曲面图形的特点,以加强几何的直观性。
2、学会用Mathematica绘制空间立体图形。
三、程序设计
空间曲面的绘制
作参数方程
所确定的曲面图形的Mathematica命令为:
ParametricPlot3D[{x[u,v],y[u,v],z[u,v]},{u,umin,umax},
{v,vmin,vmax},选项]
(1)
(2)
四、程序运行结果
(1)
(2)
五、结果的讨论和分析
1、通过参数方程的方法做出的图形,可以比较完整的显示出空间中的曲面和立体图形。
2、可以通过mathematica软件作出多重积分的积分区域,使积分能够较直观的被观察。
3、从
(1)中的实验结果可以看出,所围成的立体图形是球面和圆柱面所围成的立体空间。
4、从
(2)中的实验结果可以看出围成的立体图形的上面曲面的方程是
,下底面的方程是z=0,右边的平面是
。
实验一空间曲线与曲面的绘制
一、实验题目:
(实验习题1-3)
观察二次曲面族
的图形。
特别注意确定k的这样一些值,当k经过这些值时,曲面从一种类型变成了另一种类型。
二、实验目的和意义
1.学会利用Mathematica软件绘制三维图形来观察空间曲线和空间曲线图形的特点。
2.学会通过表达式辨别不同类型的曲线。
三、程序设计
这里为了更好地分辨出曲线的类型,我们采用题目中曲线的参数方程来画图,即
输入代码:
ParametricPlot3D
[{r*Cos[t],r*Sin[t],r^2+k*r^2*Cos[t]*Sin[t]},
{t,0,2*Pi},{r,0,1},PlotPoints->30]
式中k选择不同的值:
-4到4的整数带入。
四、程序运行结果
k=4:
k=3:
k=2:
k=1:
k=0:
k=-1:
k=-2:
k=-3:
k=-4:
五、结果的讨论和分析
k取不同值,得到不同的图形。
我们发现,当|k|<2时,曲面为椭圆抛物面;当|k|=2时,曲面为抛物柱面;当|k|>2时,曲面为双曲抛物面。
实验二无穷级数与函数逼近
一、实验题目:
(实验习题2-2)
改变例2中m及
的数值来求函数的幂级数及观察其幂级数逼近函数的情况。
二、实验目的和意义
1.利用Mathematica显示级数部分和的变化趋势。
2.学会如何利用幂级数的部分和对函数进行逼近以及函数值的近似计算。
三、程序设计
若函数
能展开成x-
的幂级数(这里不验证),则根据函数展开为幂级数的展开公式,其展开式为
。
因此首先定义
的n阶导数的函数g(n,
),最后再构成和式即得
的幂级数展开式。
用Mathematica观察幂级数部分和逼近函数的情况。
m=–2,
=2时
输入如下命令:
m=-2;
f[x_]:
=(1+x)^m;
x0=2;
g[n_,x0_]:
=D[f[x],{x,n}]/.x
x0;
s[n_,x_]:
=Sum[
*(x-x0)^k,{k,0,n}];
t=Table[s[n,x],{n,20}];
p1=Plot[Evaluate[t],{x,-1/2,1/2}];
p2=Plot[(1+x)^m,{x,-1/2,1/2},PlotStyle
RGBColor[0,0,1]];
Show[p1,p2]
四、程序运行结果
从输出的图形观察
展开的幂级数的部分和逼近函数
的情况:
五、结果的讨论和分析
从图中可以看到,当n越大时,幂级数越逼近函数。
实验二无穷级数与函数逼近
一、实验题目:
(实验习题2-3)
观察函数
展成的傅里叶级数的部分和逼近
的情况。
二、实验目的和意义
1.利用Mathematica显示级数部分和的变化趋势。
2.学会展示傅里叶级数对周期函数的逼近情况。
三、计算公式
可以展开成傅里叶级数:
,其中
,
四、程序设计
输入代码:
f[x_]:
=Which[-Pi<=x<0,-x,0<=x a[n_]: =Integrate[-x*Cos[n*x],{x,-Pi,0}]/Pi+ Integrate[Cos[n*x],{x,0,Pi}]/Pi; b[n_]: =Integrate[-x*Sin[n*x],{x,-Pi,0}]/Pi+ Integrate[Sin[n*x],{x,0,Pi}]/Pi; s[x_,n_]: =a[0]/2+Sum[a[k]*Cos[k*x]+b[k]*Sin[k*x],{k,1,n}]; g1=Plot[f[x],{x,-2Pi,2Pi},PlotStyle->RGBColor[0,0,1], DisplayFunction->Identity];m=18; For[i=1,i<=m,i+=2, g2=Plot[Evaluate[s[x,i]],{x,-Pi,Pi},DisplayFunction->Identity]; Show[g1,g2,DisplayFunction->$DisplayFunction]] 五、程序运行结果 六、结果的讨论和分析 从图表可以看出,n越大逼近函数的效果越好,还可以注意到傅里叶级数的逼近是整体性的。 实验三最小二乘法 一、实验题目: (实验习题3-2) 一种合金在某种添加剂的不同浓度下进行实验,得到如下数据: 浓度x 抗压强度y 已知函数y与x的关系适合模型: ,试用最小二乘法确定系数a,b,c,并求出拟合曲线。 二、实验目的和意义 1.学会利用最小二乘法求拟合曲线。 2.学会画数据点的散点图及拟合函数的图形,并将两个图画在同一坐标下。 三、计算公式 根据最小二乘法,要求 取最小值,令此函数对各个参数的偏导等于0,解n+1元的方程组便可求得这些参数的最小二乘解。 四、程序设计 输入代码: x=Table[+*i,{i,0,4}]; y={,,,,}; xy=Table[{x[[i]],y[[i]]},{i,1,5}]; q[a_,b_,c_]: =Sum[(a+b*x[[i]]+c*x[[i]]^2-y[[i]])^2,{i,1,5}] NSolve[{D[q[a,b,c],a]==0,D[q[a,b,c],b]==0, D[q[a,b,c],c]==0},{a,b,c}] t1=ListPlot[xy,PlotStyle->PointSize[], DisplayFunction->Identity]; f[x_]: =+*x+*x^2; t2=Plot[f[x],{x,5,35},AxesOrigin->{5,25}, DisplayFunction->Identity]; Show[t1,t2,DisplayFunction->$DisplayFunction] 五、程序运行结果 首先得到a,b,c三个值: {{a->,b->,c->}} 然后得到同一坐标系下的数据点散点图及拟合函数的图形: 六、结果的讨论和分析 观察a,b,c的值以及图像可以发现,二次方项的系数非常小,而所得的图像也非常接近于直线。 实验三最小二乘法 一、实验题目: (实验习题3-3) 在研究化学反应速度时,得到下列数据: 3 6 9 12 15 18 21 24 其中 表示实验中作记录的时间, 表示在相应时刻反应混合物中物质的量,试根据这些数据建立经验公式。 二、实验目的和意义 1.学会利用最小二乘法求拟合曲线。 2.学会由实际经验或相关的学科理论,能够提供拟合函数的可取类型,通过适当的变量代换将拟合函数线性化,建立经验公式。 三、计算公式 在许多场合下,拟合函数不具有线性形式,但是由实际经验或相关的学科理论,能够提供拟合函数的可取类型,而且可以通过适当的变量代换将拟合函数线性化,同样可以建立经验公式。 模型 可以用变量替换 将函数化为线性函数: 。 四、程序设计 输入代码: (1)生成数据并作图观察 t1={3,6,9,12,15,18,21,24}; y1={,,,,,,,}; data1=Transpose[{t1,y1}]; d2=ListPlot[data1,PlotStyle->{RGBColor[0,0,1],PointSize[]}]; (2)确定回归函数的类型 logy=Log[y1]; data2=Transpose[{t1,logy}]; d3=ListPlot[data2,PlotStyle->{RGBColor[0,0,1], PointSize[]}]; (3)对Lny数据进行最小二乘线性拟合 ly=Fit[data2,{1,x},x] y=Exp[ly]//Factor (4)绘图观察回归曲线的拟合效果 g=Plot[y,{x,1,25}, PlotStyle->RGBColor[,,]]; Show[g,d2]; 五、程序运行结果 六、结果的讨论和分析 在实际应用中,可以根据实际背景、理论分析、型值点形态等因素选择适当的拟合曲线。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 实验 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)