数学实验习题
实验1 MATLAB基本特性与基本运算
1． 求解方程ax2?bx?c?0的根。其中
（1）a?1,b?2,c?3 （2）a?1,b??2,c??3 （提示：运用求根公式。结果为（1）
x1,2??1?2i，（2）
x1,2??1,3）
2． 已知圆的半径为15，求圆的周长和面积。
3． 输入例1-6中语句，计算三角形的面积并修改边长值重新计算三角形的面积。 4． 查询表1-4中部分常用函数的功能与用法。
?2A????15． 设
y?x3?1??0?,B??2???2z??2??0?，求矩阵方程AX?2A?B?X的解。
2ln(1?x)x26． 画出
1?x2和在区间[-5,5]上的图形（提示：用 .^ 和 ./ 运算）。 sinx在区间[-5,5]上的图形。
7． 画出f(x)?e8． 设f(x)?e x2sinxcosx?e2cosx2sinxcosx?et2cosxsinx，试在[-5,5]上求出函数的零点及极大、极小值。
9． 求方程? 0(t3?2e?3cost)dt?s?041?x2当s?1、11、21时的根。
??10． 已知
f(x)? 1?dx 0(试证明)，试用不同的积分命令求其近似值(pi=3.14159265358…)。 ?sinax|x|21?e?1/x11．设，试求当
limf(x)x?1存在时a的大小以及极限值。
12．设f(x)?x?sin(x?cosx)，求f(x)在[0,4?]上的极值、拐点。
113．计算积分（1）?xsinxdx；（2）?ln(x?x2?1)dx0。
实验2 MATLAB绘制二维、三维图形
1． 在圆域x2?y2?1上画出上半球面z?1?x2?y2的图形。
x22． 画出椭球面1?y24?z212?1的图形。
?(x?y)223． 在矩形域[-2,2]×[-2,2]区域上画出函数z?xe的图形。
4． 为探测一河床横断面结构，在河面上等距离测量深度如下，试绘图表示。 距 离 测量1 0 10 20 30 40 50 60 70 80 90 0 2.51 8.54 15.60 28.92 20.64 10.53 30.24 15.00 0 测量2 测量3 0 1.81 9.04 16.70 26.99 18.98 11.36 32.56 12.08 0 0 2.02 7.87 14.68 27.87 21.86 12.84 28.43 16.00 0
实验3 MATLAB编程介绍与循环结构 实验4 MATLAB选择结构与应用实验
1、验证“哥德巴赫猜想”，即：任何一个正偶数(n?6)均可表示为两个质数的和。要求编制一个函数程序，输入一个正偶数，返回两个质数的和。
2、编制程序验证一个正整数能否可以表示为多个连续的正整数之和。如：6=1+2+3；15=1+2+3+4+5或=4+5+6或=7+8等等。要求将2到100之间的所有整数给出相应的结果，你能总结出哪些规律。
3、最优截断切割问题（97年数模国赛B题简化） 某行业要进行切割加工，现将一个三条棱长分别 为10、14.5、19cm的长方体，进行6个面的水平和 垂直切割加工处理各切一次，切出一个三条棱长分别 为3、2、4cm的小长方体。若已知切割时，左侧面、 正面、底面的切割厚度分别为6、7、9cm，切割费 用为1元/平方厘米时，求最佳切割方案和切割费用。
4、追击问题
在一边长为1的正方形跑道的四个顶点上各站有1人， 他们同时开始以等速顺时针追逐下一人，在追击过程中， 每个人时刻对准目标，试模拟追击路线。并讨论：
(1) 四个人能否追到一起？
(2) 若能否追到一起，则每个人跑过多少路程？ (3) 追到一起所需要的时间（设速率为1）？ (4) 如果四个人追逐的速度不一样，情况又如何呢？
实验5 开普勒方程近似解与方程求根
1. 某农夫有一个半径10米的圆形牛栏，长满了草.他要将一头牛栓在牛栏边界的栏桩上，但只让牛吃到一半草，问栓牛鼻的绳子应为多长？
2．图5.4所示，为了在海岛I与某城市C之间铺 设一条地下光缆，每千米光缆铺设成本在水下部分是在地下部分是
C2C1， C ，为使得铺设该光缆的总成本最低，光缆 P 的转折点P（海岸线上）应该取在何处？
如果实际测得海岛I与城市C之间水平距离l?30km， I 海岛距海岸线垂直距离
h1?15km,城市距海岸线垂直距离 图5.4 海岛与城市位置
h?10km,C1?3000万元/km,C2?1500万元/km,求P点的坐标（误差?10?3km）.
3．有一艘宽为5m的长方形驳船，欲过某河道的直角 湾，经测量知河道的宽度为10m和12m，如图5.5所示， 试问，要驶过该直角湾，驳船的长度不能超过多少米？ （误差?10?3m） 图5.5 驳船拐弯示意图
4．一个对称的地下油库，内部设计如图5.6所示。 横截面为圆，中心位置上的半径为3m，上下底 上的半径为2m，高为12m，纵截面的两侧是顶 点在中心位置的抛物线。试求：
(1) 油库内油面的深度为h(从底部算起)时， 库内油量的容积V(h)；
(2) 设计测量油库油量的标尺，使当油量 容积V已知时，算出油的深度h，刻出油量大小。 试给出当
V?10m,20m,30m,?333时油的深度. 图5.6 地下油库示意图
5.下面是某报纸2006年3月30日第七版上的一则房产广告： 建筑面积 105.9 m
2
总价 36万
30%首付 10.8万
70%按揭 30年
月还款 1436元
不难算出，买房者向银行总共借了25.2万，30年内共要还51.696万，约为当初借款的两倍。试计算这个案例中贷款年利率是多少？
6. 作为房产公司的代理人，你要迅速准确回答客户各方面的问题。现在有个客户看中了你公司一套建筑面积为120平方米，每平方米单价4200元的房子。他计划首付30%，其余70%用20年按揭贷款（贷款年利率5.5%）。请你提供下列信息：房屋总价格、首付款额、月付还款额。
实验6 Logistic方程求解与混沌
1．Feigenbaum曾对超越函数y??sin(?x)（?为非负实数）进行了分岔与混沌的研究，试利用迭代格式
xk?1??sin(?xk)，做出相应的Feigenbaum图。
2. Henon吸引子是混沌和分形的著名例子，迭代模型为
取初值x0 = 0, y0 = 0, 进行3000次迭代，对于k>1000, 在(xk, yk) 处画一点(注意不要连线)可得所谓Henon引力线图. 试写出程序，画出图形。
实验7 ?的计算与数值积分
1.(排洪量)某河床的横断面如图7.2所示，为了计算最大的排洪量，需要计算它的断面积，试根据图示测量数据（单位：米）用梯形法计算其断面积。
图7.2 河道河床截面图
2.(停产时间)某公司投资2000万元建成一条生产线。投产后，在时刻t 的追加成本和追加收益分别为G(t)= 少？
3.从地面发射一枚火箭，在最初100秒内记录其加速度如下表，试求火箭在100秒时的速度。 t(s) 0 10 20 30 40 50 60 70 80 90 100 (百万元/年), H(t)=
。试确定该生产线在何时停产可获最大利润？最大利润是多
a (m.s2) 30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.67 54.01 57.23 x24.计算椭圆4?y2?1的周长，使结果具有五位有效数字。
5. 已知某直升飞机旋转机翼外形曲线上的部分坐标值如下: x y x y 0.52 5.2880 156.6 36.6000 8.0 13.8400 260.7 31.0000 17.95 20.2000 364.4 20.9000 28.65 24.90000 468.0000 7.8000 50.65 31.1000 507.0 1.5000 104.6 36.5000 520.0 0.2000 且两端点上的一阶导数值分别为1.86548和-0.046115,计算各节点处的一阶和二阶导数值,并计算在[0.52,520.0]内的积分值.
实验8 河流流量估计与数据插值
1. 选择一些函数, 在n个节点上(n不要太大, 取10个左右)分别用拉格朗日﹑分段线性﹑三次样条三种插值方法,计算m个插值点的函数值(m要适中,如50～80).通过数值和图形输出,将三种插值结果与精确值进行比较。适当增加n, 再作比较，由此作初步分析,你能得出什么结论. 下列函数供选择参考:
a. y=cosx, 0≤x≤2π; b. y=(1-x2)1/2,-1≤x≤1; c. y=sin10x, -2≤x≤2; d. y= x exp(-x),-2≤x≤2. 2. 已知f(x)的函数表如下: x y
0.40 0.41075
0.55 0.57815
0.65 0.69675
0.80 0.88811
0.90 1.02652
1.05 1.25382
求四次Lagrange插值多项式,并由此求f (0.596)的近似值.
3.如图所示是瑞士地图,为了算出它的国土面积,首先 对地图作如下测量:以由西向东方向为x轴,由南向北方向为y轴,选作方便的点为原
点,并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段,在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2,这样就得到了如下数据(单位mm ),已知地图的比例是1:2222, 试由测量的值计算瑞士的国土面积,与它的精确值41288平方公里比较 测量数据如下:
x=[7.0, 10.5, 13.0, 17.5, 34, 40.5, 44.5, 48, 56, 61, 68.5, 76.5, 80.5, 91, 96, 101, 104, 106, 111.5, 118, 123.5, 136.5, 142, 146, 150, 157, 158]
y1=[44, 45, 47, 50, 50, 38, 30, 30, 34, 36, 34, 41, 45, 46, 43, 37, 33, 28, 32, 65, 55, 54, 52, 50, 66, 66, 68]
y2=[44, 59, 70, 72, 93, 100, 110, 110, 110, 117, 118, 116, 118, 118, 121, 124, 121, 121, 121, 122, 116, 83, 81, 82, 86, 85, 68]
实验9 人口预测与数据拟合
1．1790年到1980年各年美国人口数的统计数据如下表: 美国人口统计数字(单位: 百万) 年份 1790 统计 3.9 1800 5.3 1810 7.2 1820 9.6 1830 12.9 1840 17.1 1850 23.2 1860 31.4 1870 38.6 1880 50.2 年份 1890 统计 62.0 1900 72.0 1910 92.0 1920 106.5 1930 123.2 1940 131.7 1950 150.7 1960 179.3 1970 204.0 1980 226.5 试根据前100年的数据，分别用Malthus模型和Logistic模型建立美国人口增长的近似曲线（设美国人口总体容纳量为10亿），并预测后100年的人口数，通过与实际数据相比较，对两种预测结果进行分析.
2．由于地球围绕太阳公转，所以气温具有周期性（定性分析）。某地气象台测得当地月平均气温如下表
某地全年气温数据(单位: ℃) 月份 平均气温 1 3.1 2 3.8 3 6.9 4 12.7 5 16.8 6 20.5 7 24.5 8 25.9 9 22 10 16.1 11 10.7 12 5.4 试建立该地全年温度与时间的经验公式。（提示：周期函数可设为傅立叶级数部分和） 3. 某商品的需求量与消费者的平均收入，商品价格的统计数据如下，建立模型并进行检验，预测平均收入为1000，价格为6时的商品需求量。
需求量 收入 价格 100 75 80 70 50 65 90 100 110 60 1000 600 1200 500 300 400 1300 1100 1300 300 5 7 6 6 8 7 5 4 3 9 ，其中V0是电容
4. 用电压U=10伏的电池给电容器充电，t时刻的电压为V(t)=U-(U-V0) 器的初始电压，τ是充电常数．试由下面一组t，V数据确定V0和τ． t (秒) 0.5 1 2 3 4 5 7 9