数学实验大作业
2016-2017学年第一学期《数学实验-大作业》
专业班级测绘工程1502
姓名闻小玖
学号1501060212
开课系室计算数学系
完成日期2016年11月6日
一、(基本题目，每题10分，共8题)用mathematica 求解下面问题，写出程序和结果。 1. 求隐函数()y y x =的导数
arctan
ln y
x =解：代码：
In[1]=
Out[1]=
解题过程：
结果： =
2. 求2,sin x t y t ==在(2,sin 2)处的切线和法线方程。 解：代码：
In[1]=
所以切线方程为（
y-Sin
）=Cos[
]/(2
)*(x-2)
因为切线斜率与法线斜率乘积为-1 输入：-1/Cos[()/(2
)
得-2
Sec[()]
所以法线方程为（y-Sin
）=-2
Sec[(
)]*(x-2)
解题过程：
2
2
2
22
2
2
2
2
2
结果：切线方程为（
y-Sin ）=Cos[]/(2)*(x-2)
法线方程（y-Sin
）=-2
Sec[(
)]*(x-2)
3.
求
?
-π
3)cos(2
dx x e x 的近似值.
In[2]= Out[2]= 解题过程：
结果： 4. 计算x
x x ln lim
20+
→
解：代码
ln[4]:=Limit[x^2*lnx, x -> 0, Direction -> -1] Out[4]=
解题过程：
结果：0
5. 作出柱面2221x y +=和圆柱面2
2
1x z +=相交的图形 解：代码
In[9]=
2222
2
2
解题过程：
6. 求()61222,2244+--+=y x y x y x f 的极值. 解：代码
In[50]=Clear[f];
f[x_,y_]=x^4+2*y^4-2*x^2-12*y^2+6; fx=D[f[x,y],x] fy=D[f[x,y],y] -4 x+4 x3 -24 y+8 y3
(-4+12 x2) (-24+24 y2)
data={x,y,fxx,disc,f[x,y]}/.critpts;
TableForm[data,TableHeadings →{None,{"x","y","fxx","disc","f"}}] Out[52]=
所以：
当x=-1,y=-Sqrt[3]时判别式judge=384>0,A=fxx=8>0,函数有极小值-13; 当x=-1,y=Sqrt[3]时判别式judge=384>0,A=fxx=8>0,函数有极小值-13; 当x=0,y=0时判别式judge=96>0,A=fxx=-4<0,函数有极大值6;
当x=1,y=-Sqrt[3]时判别式judge=384>0,A=fxx=8>0,函数有极小值-13
;
当x=1,y=-Sqrt[3]时判别式judge=384>0,A=fxx=8>0,函数有极小值-13;
当x=-1,y=0或x=0,y=-Sqrt[3]或x=0,y=Sqrt[3]或x=1,y=0时判别式judge<0,函数无极值。
解题过程：
结果：
当x=-1,y=-Sqrt[3]时判别式judge=384>0,A=fxx=8>0,函数有极小值-13; 当x=-1,y=Sqrt[3]时判别式judge=384>0,A=fxx=8>0,函数有极小值-13; 当x=0,y=0时判别式judge=96>0,A=fxx=-4<0,函数有极大值6;
当x=1,y=-Sqrt[3]时判别式judge=384>0,A=fxx=8>0,函数有极小值-13; 当x=1,y=-Sqrt[3]时判别式judge=384>0,A=fxx=8>0,函数有极小值-13;
当x=-1,y=0或x=0,y=-Sqrt[3]或x=0,y=Sqrt[3]或x=1,y=0时判别式judge<0,函数无极值。
7. 求曲面积分,22zdxdy y x ??∑
其中∑为球面2222a z y x =++的下半部分的外侧.
解：代码
In[3]= Clear[A, fa, ds]; A = {0, 0, z*x^2*y^2}; fa = {x, y, z}/a; ds = a^2*Sin[u];
Integrate[(A.fa /. {x -> a*Sin[u]*Cos[v], y -> a*Sin[u]*Sin[v], z -> a*Cos[u]})*ds // Simplify, {u, Pi/2, Pi}, {v, 0, 2*Pi}] Out[7]=
解题过程：
结果：
8. 数列n x 的递推公式为： 1x =
n x =
（1） 试求出数列充分多的项，观察它的极限是否存在； （2） 解：代码
1__1.41421356237 2__1.54221082541 3__1.57598084511 4__1.58453826525 5__1.58668487076 6__1.58722197637 7__1.58735628117 8__1.58738985915 9__1.58739825376 10__1.58740035241 11__1.58740087708 12__1.58740100825 13__1.58740104104 14__1.58740104924 15__1.58740105129 16__1.58740105180 17__1.58740105193 18__1.58740105196 19__1.58740105197 20__1.58740105197 21__1.58740105197 22__1.58740105197 23__1.58740105197 24__1.58740105197
26__1.58740105197 27__1.58740105197 28__1.58740105197 29__1.58740105197 30__1.58740105197 31__1.58740105197 32__1.58740105197 33__1.58740105197 34__1.58740105197 35__1.58740105197 36__1.58740105197 37__1.58740105197 38__1.58740105197 39__1.58740105197 40__1.58740105197 41__1.58740105197 42__1.58740105197 43__1.58740105197 44__1.58740105197 45__1.58740105197 46__1.58740105197 47__1.58740105197 48__1.58740105197 49__1.58740105197 50__1.58740105197 51__1.58740105197 52__1.58740105197 53__1.58740105197 54__1.58740105197 55__1.58740105197 56__1.58740105197 57__1.58740105197 58__1.58740105197 59__1.58740105197 60__1.58740105197 61__1.58740105197 62__1.58740105197 63__1.58740105197 64__1.58740105197 65__1.58740105197 66__1.58740105197 67__1.58740105197 68__1.58740105197
70__1.58740105197
71__1.58740105197
72__1.58740105197
73__1.58740105197
74__1.58740105197
75__1.58740105197
76__1.58740105197
77__1.58740105197
78__1.58740105197
79__1.58740105197
80__1.58740105197
解题过程：
结果：故存在极限1.587
（3）用Mathematica画出它的散点图；
（4）通过求解方程，求出此数列的极限。
因为存在极限，所以当x->+∞时，有x= Sqrt[2 + Sqrt[x]]
故：{{m→1.83118}}
二、(综合题，本题20分) 画出由半个椭球面和若干个小球组成的图形，其中小球的半径远小于椭球的短半径。
（下面两种情况任选其一）
1情况一：半个椭球面固定，小球以某种能够平衡的方式放在椭球面里面，要求球面和球面，球面和椭球面相切，小球的个数大于等于3。
2情况二：半个椭球面固定，小球在椭球面里面运动，小球的个数大于等于1，做出动画图。解：
一、问题背景与提出
半个椭球面固定，小球以某种能够平衡的方式放在椭球面里面，要求球面和球面，球面和椭球面相切，小球的个数大于等于3，画出由半个椭球面和若干个小球组成的图形
二、实验目的
画出题目要求组成的图形
三、实验原理与数学模型
通过数学建模的方法，寻找一种可行的方法使椭球与球能满足问题要求，并绘制图像。根据题意，相切时，再切点两曲面法向量方向相同；三个球相切有且只有一种情况，可以考虑先固定球体，然后在寻找半椭球，使半椭球与球相切；椭球与球都具有对称性。
四、实验内容（要点）
1.实验要求绘制出图形，并没对球与半椭球的相关几何性质做要求，再加之三球相切的位置关系有且只有一种，所以可以考虑固定三个球体的位置，然后再寻找一个与其相切的椭球。
2.为了简化椭球的表示，可以选用有一个截面是圆的椭球，椭球的参数可以减小，且根据经验，半椭球的中轴线应该位于三个球的中心，此时只要有一个球与椭球相切就可以保证其他的球与椭球也相切。
3.为了使计算更为直观，可以将椭球与球都投影到平面上，转化成椭圆与圆只有一个交点的问题。
五、实验过程记录（含基本步骤、主要程序清单及异常情况记录等）
1.取3个半径为1的球：
b1=ParametricPlot3D[{Sin[φ]Cos[θ],Sin[φ]Sin[θ]+Q,Cos[φ]},{φ,0,π},{θ,0,2π},PlotRa nge→{{-2.1,2.1},{-2.5,2.5},{-2.1,2.1}}];
b2=ParametricPlot3D[{Sin[φ]Cos[θ]+1,Sin[φ]Sin[θ]-Q/2,Cos[φ]},{φ,0,π},{θ,0,2π}];
b3=ParametricPlot3D[{Sin[φ]Cos[θ]-1,Sin[φ]Sin[θ]-Q/2,Cos[φ]},{φ,0,π},{θ,0,2π}];
Show[b1,b2,b3]
2.这时，做出三个球在ZOY 平面上的投影
Rightsightc1=ParametricPlot[{Cos[θ]+Q,Sin[θ]},{θ,0,2π},PlotRange →{{-3,3},{-1.2,1.2}}];
Rightsightc2=ParametricPlot[{Cos[θ]-Q/2,Sin[θ]},{θ,0,2π}]; Rightsightc3=ParametricPlot[{Cos[θ]-Q/2,Sin[θ]},{θ,0,2π}]; Show[Rightsightc1,Rightsightc2,Rightsightc3]
2. 为了简便分析，可以选取一个平行于XOY 平面方向截面为圆、有一根对
称轴与Z 轴重合的椭球进行研究。其标准方程为：，现只要确定参数p ，似的半椭球与三个球分别有且只有一个交点。 3. 为了直观，将所有的立体研究对象都投影到ZOY 平面上，且为了讨论符
合习惯，将X 轴代替投影后的Y 轴，将Y 轴代替投影后的Z 轴，建立新的XOY 平面直角坐标系。
el=ParametricPlot[{3Cos[θ],Sqrt[3]Sin[θ]},{θ,0,2π}]; Show[el,Rightsightc1,Rightsightc2,Rightsightc3]
321
1.0
0.5321
1.5
1.00.5
中的参数只是为了绘图而设置）
4.根据对称性，只要考虑一个圆与椭圆相切即可保证其他圆于椭圆相切。
在图中记最右侧的圆为
Rightsightc1：,
其直角坐标方程为Rightsightc1：。
椭圆方程为,
直角坐标系中的方程为：。
取，得圆上一点,过该点切线斜率为。
椭圆若与其相切，则椭圆要求过该点，且过该点切线斜率为，即：
用mathematica求解：
Clear[a,b,p];
Solve[{7b 3a(1+2p),49/(12a)+(1/2+p)2/b 1},{a,b}]
{{a→(7 (4+p))/6,b→1/2 (4+9 p+2 p2)}}
（程序中为了表达方便，用a，b代替了其平方，最后的解其实为
）。
5.绘制平面图可有如图效果：
El=ParametricPlot[{Sqrt[a]*Cos[θ],Sqrt[b]*Sin[θ]+p},{θ,0,2π}];
REC=ParametricPlot[{Cos[θ]+Sqrt[4/3],Sin[θ]},{θ,0,2π},PlotRange→{{-3,3},{-
1.2,1.2}}];
po=ListPlot[{7*Sqrt[3]/6,-0.5}];
Show[El,REC,po]
21
1
6.取p=0.5，可以得，带入椭球方程：,转化为参数方程：
7.再绘制图像：
El3D=ParametricPlot3D[{Sqrt[a]*Sin[φ]Cos[θ],Sqrt[a]*Sin[φ]Sin[θ],Sqrt[b]*Cos[
φ]+p},{θ,0,2π},{φ,0.5π,π},PlotRange→{{-2.7,2.7},{-2.7,2.7},{-1.9,2}}];
Show[El3D,b1,b2,b3]
六、实验结果报告与实验总结
1.实验结果：根据分析与计算，我们完成了实验目的，绘制出可行的结果。
2.实验总结：此次实验大体分为四个步骤：
绘制相切的三个球；
a)将球投影到平面上；
b)再平面内寻找合适的椭圆；
c)将平面图形转化为立体图。
3.先发现了固定球体，而不是固定椭球是突破口；其次，将立体问题适当
的经过投影转换为平面问题可以很好的化简；然后，问题没有规定具体
的方程参数，合理的给未知参数赋值可以使方程求解更加直观。
七、参考文献
1.中国石油大学（华东）数学实验初稿（2016-2017学期）
2.网页
https://www.360docs.net/doc/f719018013.html,/link?url=RzhG4JHvl9ChxXqGYEEsoSCDf AqqCWGGBmj6bnMRuxwetOeKv7X2tQ0r6SHeP7U_9ZQRpc88eBMgXa--64 IaKWvyZZLrssghglrB1bSsJ8e
3.费祥历，亓建，《高等数学（上册）》《高等数学（下册）》，高等教育出版
社（2015年8月）
最新版北京科技大学第三次数学实验报告
《数学实验》报告 实验名称Matlab三维曲面绘图 学院东凌经济管理学院 专业班级 姓名 学号 2016年3月
一、【实验目的】 1.了解并掌握Matlab三维曲面绘图； 2.进一步掌握绘图程序格式和意义； 3.初步掌握meshgrid, mesh, surf, colordef, colormap, light等使用。 二、【实验任务】 79-7 79-9 三、【实验程序】 79-7 t1=-3:0.1:3; [x1,y1]=meshgrid(t1); z1=x1.^2+y1.^2;
subplot(1,2,1);colordef white;light('position',[20,20,5]);colormap(pin k); mesh(x1,y1,z1),title('x^2+3.*y^2'); subplot(1,2,2);colordef white;light('position',[20,20,5]);colormap(pin k); surf(x1,y1,z1),title('x^2+3.*y^2') 79-9 t=-2:0.1:2; [x,y]=meshgrid(t); z1=5-x.^2-y.^2; subplot(1,3,1),mesh(x,y,z1),title('抛物面') z2=3*ones(size(x)); subplot(1,3,2),mesh(x,y,z2),title('平面') r0=abs(z1-z2)<=0.2; zz=r0.*z2;yy=r0.*y;xx=r0.*x; subplot(1,3,3),plot3(xx,yy,zz,'x'),title('交线') 四、【实验结果】 79-1
数学实验上机汇总未完成
数学实验上机作业整理∈hyd 实验一 1. 计算球体体积（半径r=5） r=5;v=(4/3)*pi*r^3 v =523.5988 2.设矩阵1234567891023416A ?? ? = ? ??? （1）提取A 的第2列赋值给B; A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];B=A(:,2) B = 2 7 3 （2）提取A 的第2行前3个数给C ； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];C=A(2,[1,2,3]) C = 6 7 8 （3）提取A 第1,3行和2, 4列相交位置元素构成子矩阵D ； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];D=A([1,3],[2,4]) D = 2 4 3 1 （4）构造矩阵E 使得E 的结构为：132213C E D C ???? ?= ? ?? A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];E=[D [C;C]] E = 2 4 6 7 8 3 1 6 7 8
（5）把A 中间的8换为0； A(2,3)=0;A A = 1 2 3 4 5 6 7 0 9 10 2 3 4 1 6 （6）去掉A 的第2行； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6]; A(2,:)=[] A = 1 2 3 4 5 2 3 4 1 6 3.写出完成下列操作的命令 （1） 建立10阶单位矩阵A; A=eye(10) （2）建立5×6的随机矩阵A ，其元素为[100，200]范围内的随机数； A=rand(5,6)*100+100 （3）将A 对角线元素加30 A+eye(5,6)*30 4.（选做题）设有分块矩阵333223E R A O S ????? =? ??? ，其中E,R,O,S 分别为单位矩阵、随机矩阵、零矩阵和对角矩阵，试通过数值计算验证2 2 E R RS A O S +?? =? ??? 。 S=[1 1;1 1]; E=eye(3);R=rand(3,2); O=zeros(2,3); [E R;O S]^2 [E R+R*S;O S^2] 实验二 1.设矩阵1215346562A -?? ? = ? ?-?? （1）求A 的秩、A 的每个元素3次方; A=[1 2 -1;5 34 6;-5 6 2];
大学数学实验
大学数学实验 项目一 矩阵运算与方程组求解 实验1 行列式与矩阵 实验目的 掌握矩阵的输入方法. 掌握利用Mathematica (4.0以上版本) 对矩阵进行转置、加、减、数乘、相乘、乘方等运算, 并能求矩阵的逆矩阵和计算方阵的行列式. 基本命令 在Mathematica 中, 向量和矩阵是以表的形式给出的. 1. 表在形式上是用花括号括起来的若干表达式, 表达式之间用逗号隔开. 如输入 {2,4,8,16} {x,x+1,y,Sqrt[2]} 则输入了两个向量. 2. 表的生成函数 (1) 最简单的数值表生成函数Range, 其命令格式如下: Range[正整数n]—生成表{1,2,3,4,…,n }; Range[m, n]—生成表{m ,…,n }; Range[m, n, dx]—生成表{m ,…,n }, 步长为d x . (2) 通用表的生成函数Table. 例如,输入命令 Table[n^3,{n,1,20,2}] 则输出 {1,27,125,343,729,1331,2197,3375,4913,6859} 输入 Table[x*y,{x,3},{y,3}] 则输出 {{1,2,3},{2,4,6},{3,6,9}} 3. 表作为向量和矩阵 一层表在线性代数中表示向量, 二层表表示矩阵. 例如,矩阵 ??? ? ??5432 可以用数表{{2,3},{4,5}}表示. 输入 A={{2,3},{4,5}} 则输出 {{2,3},{4,5}} 命令MatrixForm[A]把矩阵A 显示成通常的矩阵形式. 例如, 输入命令: MatrixForm[A] 则输出 ??? ? ??5432 但要注意, 一般地, MatrixForm[A]代表的矩阵A 不能参与运算. 输入 B={1,3,5,7} 输出为 {1,3,5,7} 输入 MatrixForm[B] 输出为
MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案
“”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;
>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =
0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1
6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5
8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =
数学实验作业
练习2﹒1 画出下列常见曲线的图形（其中a=1，b=2,c=3）。 1. 立方抛物线y = 解： x=-4:0.1:4; y=x.^(1/3); plot(x,y) -4 -3-2-101234 0.20.40.60.811.21.4 1.6 2.高斯曲线2 x y e -= 解： fplot('exp(-x^2)',[-4,4])
-4 -3 -2 -1 1 2 3 4 00.10.20.30.40.50.60.70.80.9 1 3、笛卡儿曲线23 3 2 2 33,(3)11at at x y x y axy t t = = +=++ 解：ezplot('x^3+y^3-3*x*y',[-4,4])
-4 -3-2-1 01234 -4-3-2-10123 4x y x 3+y 3-3 x y = 0 或：t=-4:0.1:4; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)
-1.5 -1-0.500.51 1.5 00.5 1 1.5 2 2.5 3 4、蔓叶线233 2 2 2 ,()11at at x x y y t t a x = = = ++- 解：t=-4:0.1:4; x=t.^2./(1+t.^2); y=t.^3,/(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)
00.10.20.30.40.50.60.70.80.91 -4 -3-2-10123 4 或： ezplot('y .^2-x.^3/(1-x)',[-4,4])
清华大学数学实验报告4
清华大学数学实验报告4
————————————————————————————————作者: ————————————————————————————————日期： ?
电13 苗键强2011010645
一、实验目的 1.掌握用 MATLAB 软件求解非线性方程和方程组的基本用法, 并对结果作初步分析； 2.练习用非线性方程和方程组建立实际问题的模型并进行求解。 二、实验内容 题目1 【问题描述】 (Q1）小张夫妇以按揭方式贷款买了1套价值２0万元的房子，首付了5万元，每月还款1000元，15年还清。问贷款利率是多少？ (Q2)某人欲贷款50 万元购房,他咨询了两家银行,第一家银行 开出的条件是每月还45０0元,15 年还清；第二家银行开出的条件是每年还４5０00 元,2０年还清。从利率方面看,哪家银行较优惠(简单假设:年利率=月利率×12）? 【分析与解】 假设初始贷款金额为x０,贷款利率为ｐ,每月还款金额为x,第i 个月还完当月贷款后所欠银行的金额为x i,(i=1,2，3,.．....，n）。由题意可知： x1=x0(1+p)?x x2=x0(1+p)2?x(1+p)?x x3=x0(1+p)3?x(1+p)2?x(1+p)?x ……
x n=x0(1+p)n?x(1+p)n?1???x(1+p)?x =x0(1+p)n?x (1+p)n?1 p =0 因而有： x0(1+p)n=x (1+p)n?1 p (1) 则可以根据上述方程描述的函数关系求解相应的变量。 （Ｑ1) 根据公式(１)，可以得到以下方程： 150p(1+p)180?(1+p)180+1=0 设 f(p)=150p(1+p)180?(1+p)180+1，通过计算机程序绘制ｆ（p）的图像以判断解ｐ的大致区间,在Ｍatlab中编程如下： ｆorｉ = 1:25 t = ０．0001＊i; p(i) = t； f(i) =１50*ｔ＊(１+t).^180-(１+t).^180+1; ｅｎd; ｐｌｏｔ(ｐ,ｆ),hｏld ｏn,grid on; 运行以上代码得到如下图像：
c++大作业学生实验报告
学生实验报告 实验课名称： C++程序设计 实验项目名称：综合大作业——学生成绩管理系统专业名称：电子信息工程 班级： 学号： 学生： 同组成员： 教师：
2011 年 6 月 23 日 题目：学生成绩管理系统 一、实验目的： （1）对C++语法、基础知识进行综合的复习。 （2）对C++语法、基础知识和编程技巧进行综合运用，编写具有一定综合应用价值的稍大一些的程序。培养学生分析和解决实际问题的能力，增强学生的自信心，提高学生学习专业课程的兴趣。 （3）熟悉掌握C++的语法和面向对象程序设计方法。 （4）培养学生的逻辑思维能力，编程能力和程序调试能力以及工程项目分析和管理能力。 二、设计任务与要求： （1）只能使用/C++语言，源程序要有适当的注释，使程序容易阅读。 （2）至少采用文本菜单界面（如果能采用图形菜单界面更好）。 （3）要求划分功能模块，各个功能分别使用函数来完成。 三、系统需求分析： 1.需求分析: 为了解决学生成绩管理过程中的一些简单问题，方便对学生成绩的管理 （录入，输出，查找，增加，删除，修改。） 系统功能分析: (1):学生成绩的基本信息：学号、、性别、C++成绩、数学成绩、英语成绩、 总分。 (2):具有录入信息、输出信息、查找信息、增加信息、删除信息、修改信息、 排序等功能。 2.系统功能模块（要求介绍各功能） （1）录入信息(Input)：录入学生的信息。 （2）输出信息(Print)：输出新录入的学生信息。 （3）查找信息(Find)：查找已录入的学生信息。 （4）增加信息(Add)：增加学生信息。 （5）删除信息(Remove)：在查找到所要删除的学生成绩信息后进行删除并输出删除后其余信息。 （6）修改信息(Modify)：在查到所要修改的学生信息后重新输入新的学生信息从而进行修改，然后输出修改后的所有信息。 （7）排序(Sort)：按照学生学号进行排序。 3.模块功能框架图
北京科技大学参数检测实验报告全
北京科技大学参数检测实验报告全
实验六工业热电偶的校验 摘要：本实验重在了解热电偶的工作原理并通过对热电偶进行校正验证镍铬热电偶的准确性并了解补偿导线的使用方法。 关键词：热电偶校正标准被校补偿导线 1 引言 （1）实验目的 1．了解热电偶的工作原理、构造及使用方法。了解热电势与热端温度的关系。了解对热电偶进行校正的原因及校正方法，能独立地进行校正实验和绘制校正曲线。 2．了解冷端温度对测量的影响及补偿导线的使用方法。 3．通过测量热电势掌握携带式直流电位差计的使用方法。 （2）实验设备 1．铂铑－铂热电偶（标准热电偶）１支 ２．镍铬－镍硅热电偶（被校正热电偶）１支 3．热电偶卧式检定炉（附温度控制器）１台 4．携带式直流电位差计 1台 5．酒精温度计 1支 6．广口保温瓶 1个 7．热浴杯及酒精灯各1个 2 内容 1．了解直流电位差计各旋钮、开关及检流计的作用，掌握直流电位差计的使用方法。 2．热电偶校正 （１）实验开始，给检定炉供电，炉温给定值为400oC。当炉温稳定后，用电位差计分别测量标准热电偶和被校正热电偶的热电势，每个校正点的测量不得少于四次。数据记录于表6－１。 （２）依次校正600oC、 800oC、 1000oC各点。 （３）将测量电势求取平均值并转换成温度，计算误差，根据表6－３判断被热电偶是否合格。绘制校验曲线。 3．热电偶冷端温度对测温的影响及补偿导线的使用方法。 （１）1000oC校正点作完后，保持炉温不变。测量热浴杯中的水温，然后用电位差计分别测量镍铬－镍硅热电偶未加补偿导线和加补偿导线的热电势。数据记录于表6－２中。 （２）用酒精灯加热热浴杯，当水温依次为30oC、 40oC、 50oC时，用电位差计分别测量镍铬－镍硅热电偶未加补偿导线和加补偿导线的热电势。数据记录于表6－２中。 （３）用铂铑－铂热电偶测量炉温，检查实验过程中炉温是否稳定，分析若炉
数学建模作业——实验1
数学建模作业——实验1 学院：软件学院 姓名： 学号： 班级：软件工程2015级 GCT班 邮箱： 电话： 日期：2016年5月10日
基本实验 1.椅子放平问题 依照1.2.1节中的“椅子问题”的方法，将假设中的“四腿长相同并且四脚连线呈正方形”，改为“四腿长相同并且四脚连线呈长方形”，其余假设不变，问椅子还能放平吗？如果能，请证明；如果不能，请举出相应的例子。 答：能放平，证明如下： 如上图，以椅子的中心点建立坐标，O为原点，A、B、C、D为椅子四脚的初始位置，通过旋转椅子到A’、B’、C’、D’，旋转的角度为α，记A、B两脚，C、D两脚距离地面的距离为f(α)和g(α)，由于椅子的四脚在任何位置至少有3脚着地，且f(α)、g(α)是α的连续函数，则f(α)和g(α)至少有一个的值为0，即f(α)g(α)=0，f(α)≥0，g(α)≥0，若f(0)＞0，g(0)=0，
则一定存在α’∈（0，π），使得 f(α’)=g(α’)=0 令α=π（即椅子旋转180°,AB 边与CD 边互换），则 f(π)=0，g(π)＞0 定义h(α)= f(α)-g(α)，得到 h(0)=f(0)-g(0)＞0 h(π)=f(π)-g(π) ＜0 根据连续函数的零点定理，则存在α’∈（ 0，π），使得 h(α’)= f(α’)-g(α’)=0 结合条件f(α’)g(α’)=0,从而得到 f(α’)=g(α’)=0，即四脚着地，椅子放平。 2. 过河问题 依照1.2.2节中的“商人安全过河”的方法，完成下面的智力游戏：人带着猫、鸡、米过河，船除需要人划之外，至多能载猫、鸡、米之一，而当人不在场时，猫要吃鸡、鸡要吃米，试设计一个安全过河的方案，并使渡河的次数尽量的少。 答： 用i =1,2,3,4分别代表人，猫，鸡，米。1=i x 在此岸，0 =i x 在对岸， ()4321,,,x x x x s =此岸状态，()43211,1,1,1x x x x D ----=对岸状态。安全状态集合为 :
MATLAB实验练习题(计算机) 南邮 MATLAB 数学实验大作业答案
“MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6)
1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408
3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans =
数学实验报告-6
《数学实验》报告 实验名称常微分方程的求解 学院材料科学与工程 专业班级材料1209 姓名曾雪淇 学号 41230265 2014年 5月
一、【实验目的】 掌握常微分方程求解和曲线拟合的方法，通过MATLAB求解一阶甚至是二阶以上的高阶微分方程。 二、【实验任务】 P168习题24，习题27 三、【实验程序】 习题24：dsolve('Dy=x*sin(x)/cos(y)','x') 习题27：function xdot=exf(t,x) u=1-2*t; xdot=[0,1;1,-t]*x+[0 1]'*u; clf; t0=0; tf=pi; x0t=[0.1;0.2]; [t,x]=ode23('exf',[t0,tf],x0t) y=x(:,1); Dy=x(:,2); plot(t,y,'-',t,Dy,'o') 四、【实验结果】 习题24：ans = -asin(-sin(x)+x*cos(x)-C1) 习题27： t = 0.014545454545455 0.087272727272727 0.201440113885487 0.325875614772746 2
0.462108154525786 0.612058884594697 0.777820950596408 0.962141414226468 1.148168188604642 1.276725612086219 1.405283035567796 1.518837016595503 1.670603286779598 1.860122410374634 2.089084425249819 2.356884067351406 2.654570124097287 2.968729389456267 3.141592653589793 x = 0.100000000000000 0.200000000000000 0.103024424647132 0.215787876799993 0.121418223032493 0.288273863806750 0.159807571438023 0.379808018692957 0.211637169341158 0.447918********* 0.275587792496926 0.484712850141869 0.348540604264411 0.481263088285519 3
大学数学实验心得体会
大学数学实验心得体会 [模版仅供参考，切勿通篇使用] 大学数学实验心得体会（一） 数学，在整个人类生命进程中至关重要，从小学到中学，再到大学，乃至更高层次的科学研究都离不开数学，随着时代的发展，人们越来越重视数学知识的应用，对数学课程提出了更高层次的要求，于是便诞生了数学实验。 学期最初，大学数学实验对于我们来说既熟悉又陌生，在我们的记忆中，我们做过物理实验、化学实验、生物实验，故然我们以为数学实验与它们一样，当我们在网上搜索有关数学实验的信息时，我们才知道，大学数学实验作为一门新兴的数学课程在近十年来取得了迅速的发展。数学实验以计算机技术和数学软件为载体，将数学建模的思想和方法融入其中，现在已经成为一种潮流。 当我们怀着好奇的心情走进屈静国老师的数学实验课堂时，我们才渐渐懂得，数学实验是一门有关计算机软件的课程，就像c语言一样，需要编辑运行程序，从而进行数学运算，它不需要自己来运算，就像计算器一样，只要我们自己记下重要程序语句，输入运行程序，便可得到运行结果，大大降低了我们的运算量，
给我们生活带来许多便捷，在大一时，我学过c语言，由于这样的基础，让我能够更快的学会并应用此软件。 时间飞逝，转眼间，我们就要结课了，这学期我们学习了mathematics的基础，微积分实验，线性代数实验，概率论与数理统计实验，数值计算方法及实验。通过这学期的学习，我也积累了些自己的学习方法和心得。首先，我们要在平时上课牢记那些mathematics语言和公式，那些东西就想单词和公式一样，只需要背诵;然后，我们要看几遍书，并多看一下例题;最后，我们要多应用mathematics软件去练习。正所谓熟能生巧，我坚信，只要我们能够做到这三步，我们就能很好的掌握这门课程。 通过学习使用数学软件，数学实验建模，使我们能够从实际问题出发，认真分析研究，建立简单数学模型，然后借助先进的计算机技术，最终找出解决实际问题的一种或多种方案，从而提高了我们的数学思维能力，为我们参加数学竞赛和数学建模打下了坚实的基础，同时也为我们进一步深造和参加工作打下一定的实践基础! 大学数学实验心得体会（二） 在此期间我充分利用研修活动时间学习，感到既有辛苦，又有收获。既有付出，又有新所得。这次远程研修让我有幸与专家和各地的数学精英们交流，面对每次探讨的主题，大家畅所欲言，
北京科技大学数学实验第五次讲解学习
北京科技大学数学实 验第五次
精品资料 《数学实验》报告 实验名称 Matlab拟合与插值 2013年12月
一、【实验目的】 1.学习Matlab的一些基础知识，主要多项式及其相关计算等； 2.熟悉Matlab中多项式的拟合，编写一些相关的Matlab命令等； 3.熟悉Matlab中多项式的插值，并编写一些相关的Matlab命令等； 4.完成相关的练习题。 二、【实验任务】 1.在钢线碳含量对于电阻的效应的研究中，得到以下数据.分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形，计算当x=0.45时的电阻值. 碳含量 0.10 0.30 0.40 0.55 0.70 0.80 0.95 x 电阻y 15 18 19 21 22.6 23.8 26 2.在某种添加剂的不同浓度之下对铝合金进行抗拉强度试验，得到数据如下，现分别使用不同的插值方法，对其中没有测量的浓度进行推测，并估算出浓度X=18及26时的抗压强度Y的值. 浓度X 10 15 20 25 30 抗压强度Y 25.2 29.8 31.2 31.7 29.4 3.用不同方法对在(-3,3)上的二维插值效果进行比较.
三、【实验程序】 1.在钢线碳含量对于电阻的效应的研究中，得到以下数据.分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形，计算当x=0.45时的电阻值. M文件 clc; clf; x=[0.1 0.3 0.4 0.55 0.7 0.8 0.95]; y=[15 18 19 21 22.6 23.8 26]; p1=polyfit(x,y,1); p3=polyfit(x,y,3); p5=polyfit(x,y,5); x1=0.1:0.05:1; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); plot(x,y,'rp',x1,y1,'b-',x1,y3,'g-.',x1,y5,'m--'); legend('拟合点','一次拟合','三次拟合','五次拟合'); disp('以下为当x=0.45时的电阻值：') disp('一阶拟合函数值'),g1=polyval(p1,0.45) disp('三阶拟合函数值'),g3=polyval(p3,0.45) disp('五阶拟合函数值'),g5=polyval(p5,0.45)
北理工数学实验作业
一． 1. 1/e 2. 3 3.1 4.e3 5. ∞ 6. 0 7.∞ 8.0 9.1/2 10.0 11.e2c12.不存在13. 1/12 Matlab实验过程： 1.1/exp(1) syms n; f=(1-1/n)^n; limit(f,n,inf) ans = 1/exp(1) 2.3 syms n; f=(n^3+3^n)^(1/n); limit(f,n,inf) ans = 3 3. 1 syms n; f=(1+sin(2*n))/(1-cos(4*n)); limit(f,n,pi/4) ans = 1 4.e^3 syms x; f=(1+cos(x))^(3*sec(x)); limit(f,x,pi/2) ans = exp(3) 5.inf syms x; f=(x^2)*exp(1/(x^2));
limit(f,x,0) ans = Inf 6.0 syms x; f=(x^2-2*x+1)/(x^3-x); limit(f,x,1) ans = 7.inf syms x; f=((2/pi)*atan(x))^x; limit(f,x,+inf) ans = Inf 8.0 syms x y; f=(1-cos(x^2+y^2))/((x^2+y^2)*exp(x^2+y^2)); limit(limit(f,x,0),y,0) ans = 9.1/2 syms x; f=(1-cos(x))/(x*sin(x)); limit(f,x,0) ans = 1/2 10.0 syms x;
f=atan(x)/(2*x); limit(f,x,inf) ans = 11.exp(2*c) syms c; f=sym('((x+c)/(x-c))^x'); limit(f,'x',inf) ans = exp(2*c) 12.极限不存在 syms x; f=cos(1/x); limit(f,x,0) ans = limit(cos(1/x), x = 0) 13.1/12 syms x; f=1/(x*log(x)^2)-1/(x-1)^2; limit(f,x,1) ans = 1/12 二．观察函数logbx,当b=1/2,1/3,1/4和b=2,3,4时函数的变化特点，总结logbx的图形特点。
matlab与数学实验大作业
《数学实验与MATLAB》 ——综合实验报告 实验名称：不同温度下PDLC薄膜的通透性 与驱动电压的具体关系式的研究学院：计算机与通信工程学院 专业班级： 姓名： 学号： 同组同学： 2014年 6月10日
一、问题引入 聚合物分散液晶(PDLC)是将低分子液晶与预聚物Kuer UV65胶相混合，在一定条件下经聚合反应，形成微米级的液晶微滴均匀地分散在高分子网络中，再利用液晶分子的介电各向异性获得具有电光响应特性的材料，它主要工作在散射态和透明态之间并具有一定的灰度。聚合物分散液晶膜是将液晶和聚合物结合得到的一种综合性能优异的膜材料。该膜材料能够通过驱动电压来控制其通透性，可以用来制作PDLC型液晶显示器等，具有较大的应用范围。已知PDLC薄膜在相同光强度及驱动电压下，不用的温度对应于不同的通透性，不同温度下的阀值电压也不相同。为了尽量得到不同通透性的PDLC薄膜，有必要进行温度对PDLC薄膜的特性的影响的研究。现有不同温度下PDLC 薄膜透过率与驱动电压的一系列数据，试得出不同温度下PDLC薄膜通透性与驱动电压的具体关系式，使得可以迅速得出在不同温度下一定通透性对应的驱动电压。 二、问题分析 想要得到不同温度下PDLC薄膜通透性与驱动电压的具体关系式可以运用MATLAB多项式农合找出最佳函数式，而运用MATLAB多项式插值可以得出在不同温度下一定通透性所对应的驱动电压。 三、实验数据 选择10、20、30摄氏度三个不同温度，其他条件一致。
（1）、10摄氏度 实验程序： x=2:2:40; y=[5.2,5.4,5.8,6.4,7.2,8.2,9.4,10.8,12.2,14.0,16.6,22.0, 30.4,39.8,51.3,55.0,57.5,58.8,59.6,60.2]; p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); disp('三次拟合函数'),f3=poly2str(p3,'x') disp('五次拟合函数'),f5=poly2str(p5,'x') disp('七次拟合函数'),f7=poly2str(p7,'x') x1=0:1:40; y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,'rp',x1,y3,'--',x1,y5,'k-.',x1,y7); legend('拟合点','三次拟合','五次拟合','七次拟合') 实验结果：
数学实验报告-2
《数学实验》报告 实验名称 MATLAB绘图 学院材料科学与工程 专业班级材料1209 姓名曾雪淇 学号 41230265 2014年 5月
学会用MATLAB绘制二维曲线、三维曲线，掌握gtext, legend, title,xlabel,ylabel,zlabel,axis 等指令用法,并学会图形的标注。二、【实验任务】 P79 习题1，习题3，习题5 三、【实验程序】 习题一： x=0:pi/10:4*pi; y1=exp(x./3).*sin(3*x); y2=exp(x./3); y3=-exp(x./3); plot(x,y1,'b*',x,y2,'r-.',x,y3,'r-.') 习题二： x1=-pi:pi/10:pi; y1=x1.*cos(x1); x2=pi:pi/10:4*pi; y2=x2.*tan(1./x2).*sin(x2).^3; x3=1:0.1:8; y3=exp(1./x3).*sin(x3); subplot(1,3,1);plot(x1,y1,'r*')，grid on,title(‘y1= x1*cosx1’) subplot(1,3,2) ;plot(x2,y2,’b-‘),grid on,title (‘y2=x2*tan(1/x2)*sinx2^3’) subplot(1,3,3);plot(x3,y3,'g+'),grid on,title (‘y3=exp(1/x3)*sinx3’) gtext(‘y1=x1cos(x1)’),gtext(‘y2=x2tan(1/x2)sin(x2)^3’), gtext(‘y3=exp(1/x3)sin(x3)’) legend(‘y1= x1*cos(x1)’, ‘y2=x2tan(1/x2)sin(x2^)3’ ‘y3=exp(1/x3)sin(x3)’) xlabel(‘x轴’),ylabel(‘y轴’),axis xy 习题三： t=0:pi/10:20*pi; x=t.*cos(pi/6.*t); y=t.*sin(pi/6.*t); z=2*t; plot3(x,y,z,'r*'),grid on title(‘圆锥螺线的图像’) xlabel(‘x轴’),ylabel(‘y轴’),zlabel(‘z轴’)
数学实验作业汇总终审稿)
数学实验作业汇总 文稿归稿存档编号：[KKUY-KKIO69-OTM243-OLUI129-G00I-FDQS58-
（1）产生一个5阶魔方矩阵M：M=magic(5) （2）将矩阵M的第3行4列元素赋值给变量t：t=M(3,4) （3）将由矩阵M第2，3，4行第2,5列构成的子矩阵赋给变N：N=M(2:4,2:3:5) （4）将由矩阵M的前3行赋给变量N：N=M(1:3,:) （5）将由矩阵M的后3列赋给变量N：N=M(:,end:-1:end-2) （6）提取M的主对角线元素，并以这些对角线元素构成对角矩阵N：N=diag(diag(M))或N=tril(triu(M)) (7)随机产生1000个100以内的整数赋值给变量t：t=round(rand(1,1000)*100) （8）随机产生100*5个100以内的实数赋值给变量M：M=rand(100,5)*100 （1）删除矩阵M的第7个元素M(7)=[] （2）将含有12个元素的向量t转换成3*4的矩阵：reshape(t,3,4) （3）产生和M同样大小的单位矩阵：eye(size(M)) （4）寻找向量t中非零元素的下标：find(t) （5）逆序显示向量t中的元素：t(end:-1:1) （6）显示向量t偶数位置上的元素：t(2:2:end) （7）利用find函数，将向量t中小于10的整数置为0：t(find(t<10&rem(t,1)==0))=0（8）不用find函数，将向量t中小于10的整数置为0：t(t<10&rem(t,1)==0)=0 （9）将向量t中的0元素用机器0（realmin）来代替：t(find(t=0))=realmin （10）将矩阵M中小于10的整数置为0：M(find(M<10)&rem(M,1)==0)=0
东南大学高等数学数学实验报告上
高等数学数学实验报告实验人员：院（系） ___________学号_________姓名____________ 实验地点：计算机中心机房 实验一 一、实验题目： 根据上面的题目，通过作图，观察重要极限：lim（1+1/n）n=e 二、实验目的和意义 方法的理论意义和实用价值。 利用数形结合的方法观察数列的极限，可以从点图上看出数列的收敛性，以及近似地观察出数列的收敛值；通过编程可以输出数列的任意多项值，以此来得到数列的收敛性。通过此实验对数列极限概念的理解形象化、具体化。 三、计算公式（1+1/n）n 四、程序设计 五、程序运行结果 六、结果的讨论和分析 当n足够大时，所画出的点逐渐接近于直线，即点数越大，精确度越高。对于不同解题方法最后均能获得相同结果，因此需要择优，从众多方法中尽可能选择简单的一种。程序编写需要有扎实的理论基础，因此在上机调试前要仔细审查细节，对程序进行尽可能的简化、改进与完善。 实验二 一、实验题目 制作函数y=sin cx的图形动画，并观察参数c对函数图形的影响。 二、实验目的和意义 本实验的目的是让同学熟悉数学软件Mathematica所具有的良好的作图功能，并通过函数图形来认识函数，运用函数的图形来观察和分析函数的有关性态，建立数形结合的思想。 三、计算公式：y=sin cx 四、程序设计 五、程序运行结果
六、结果的讨论和分析 c 的不同导致函数的区间大小不同。 实验三 一、实验题目 观察函数f(x)=cos x 的各阶泰勒展开式的图形。 二、实验目的和意义 利用Mathematica 计算函数)(x f 的各阶泰勒多项式，并通过绘制曲线图形，来进一步掌握泰勒展开与函数逼近的思想。 三、计算公式 四、程序设计 五、程序运行结果 六、结果的讨论和分析 函数的泰勒多项式对于函数的近似程度随着阶数的提高而提高，但是对于任一确定次数的多项式，它只在展开点附近的一个局部范围内才有较好的近似精确度。 实验四 一、实验题目 计算定积分的黎曼和 二、实验目的和意义 在现实生活中许多实际问题遇到的定积分，被积函数往往不能用算是给出，而通过图像或表格给出；或虽然给出，但是要计算他的原函数却很困难，甚至原函数非初等函数。本实验目的，就是为了解决这些问题，进行定积分近似计算。 三、计算公式 四、程序设计 五、程序运行结果 六、结果的讨论和分析 本实验求的近似值由给出的n 的值的不同而不同。给出的n 值越大，得到的结果越接近准确的
李萨如图模拟(Matlab大作业)
《数学实验》报告 实验名称李萨如图模拟（Matlab大作业） 2011年11月8日
一、【实验目的】 运用数学知识与MATLAB相结合，运用数学方法，建立数学模型，用MATLAB软件辅助求解模型，解决实际问题。 二、【实验任务】 一个质点沿 X轴和 Y轴的分运动都是简谐运动,分运动的表达式分别为: x=Acos ( w1t+beta ) , y=Acos(w2t+beta ) 。如果二者的频率有简单的整数比, 则相互垂直的简谐运动合成的运动将具有封闭的稳定的运动轨迹, 这种图称为李萨如图。 1，用matlab分别画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的图像（未合成）2，用matlab画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像 3，用matlab画出x轴方向和y轴方向传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。（李萨如图） 三、【实验分析及求解】 1，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，y =Acos ( w1t+beta ) 分别画出两个波的传播图像。 2，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，y =Acos ( w1t+beta )， 用matlab画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。
3，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，画出x轴方向和y 轴方向传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。（李萨如图）。
数学实验作业一
数学实验作业一 对以下问题,编写M文件: (1)用起泡法对10个数由小到大排序. 即将相邻两个数比较,将小的调到前头. 解： 代码如下： zuoye1 clear all;clc; a=[7 2 1 0 9 4 5 -3 8 6]; n=length(a); for ii=1:n-1 if a(ii+1)>=a(ii) t1=a(ii); a(ii)=a(ii+1); a(ii+1)=t1; end for jj=1:n-1 if a(jj+1)>=a(jj) t2=a(jj); a(jj)=a(jj+1); a(jj+1)=t2; end end end a 运行结果显示如下： a = 9 8 7 6 5 4 2 1 0 -3
(2)有一个 矩阵,编程求出其最大值及其所处的位置. 解： 代码如下：zuoye2.m clear; clc; a=[1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6] max=-1; flage1=0; flage2=0 for i=1:4 for j=1:5 if (a(i,j)>max) t=max; max=a(i ,j); a(i,j)=t; flage1=i; flage2=j ； end end end max flage1 flage2 运行结果显示如下： a = 1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6 flage2 = max = 45′
9 flage1 = 2 flage2 = 5 结果： （3）编程求∑=20 1 !n n 。 解： 代码如下：zuoye3.m clear; clc; sum=0; for i=2:11 sum=sum+gamma(i); end sum
MTLB实验练习题计算机南邮MATLAB数学实验大作业答案
“M A T L A B ”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) a = 1.10447+1.05983*i -1.00450+1.06095*i -. -1.00450-1.06095*i
1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2 sin cos0 x x x -=所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6)
0.7022 3、求解下列各题： 1）3 0sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans = (-32)*exp(x)*sin(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> sym x; >> vpa((int(exp(x^2),x,0,1/2)),17)