第一章 MATLAB入门 9
9 考虑函数 f(x,y)= y3/9+3x2y+9x2+y2+xy+9 (1)作出f(x,y)在-2 10. 假定某天的气温变化记录如第二章习题5，试用最小二乘方法找出这一天的气温变化规律。考虑下列类型函数, 作图比较效果，并计算均方误差。 (1) 二次函数； (2) 三次函数； (3) 钟形函数f(x)?aeb(t?14)； (4) 函数f(x)?rsin(2?12t??). 11 (化学反应平衡) 一等克分子数一氧化碳(CO)和氧气(O2)的混合物在300K和5bar压力下达到平衡，理论反应方程式为 CO + 0.5 O2 ? CO2 实际反应方程式为 CO + N2 ? x CO + 0.5 (1 +x) O2 + (1 - x) CO2 剩余CO比值x满足化学平衡方程式 Kp?(1?x)1052.?x 0?x?1 x1?xp这里Kp = 3.06, p = 5 bar求x. 12 (月还款额)作为房产公司的代理人，你要迅速准确回答客户各方面的问题。现在有个客户看中了你公司一套建筑面积为180平方米，每平方单价7500元的房子。他计划首付30%，其余70%用20年按揭贷款（贷款年利率5.04%）。请你提供下列信息：房屋总价格、首付款额、月付还款额。如果其中10万元为公积金贷款（贷款年利率4.05%）呢？ 13（栓牛鼻的绳子）农夫老李有一个半径10米的圆形牛栏，里面长满了草，老李要将家里一头牛栓在一根栏桩上，但只让牛吃到一半草，他想让上大学的儿子告诉他，栓牛鼻的绳子应为多长？ ? 14 (弦截法)牛顿迭代法是一种速度很快的迭代方法，但是它需要预先求得导函数。若用差商代替导数，可得下列弦截法 xk?1?xk?xk?xk?1f(xk) f(xk)?f(xk?1) 10 第一章 MATLAB入门 这一迭代法需要两个初值x0, x1，编写一个通用的弦截法计算机程序并用以解习题2。(提示: 函数参数求值用MATLAB函数feval) ? 15 (线性迭代) 迭代过程 x k+1 = g (x k) 的收敛性主要条件是在根的附近满足?g ‘ (x)?<1。从理论上证明线性迭代 x k+1 = a x k + 1 只有两种极限形态：不动点或无穷大。分别就a=0.9, -0.9, 1.1, -1.1 (取x0 =1, 迭代20步)用图形显示迭代过程的不同表现(提示：用subplot将4个子图放在一个图形窗口比较) ? 16 (通道中的细杆) 要运送一根细杆子通过由宽5cm和宽10cm的通道垂直交叉口，在运送过程中必须保持杆子是水平的(如图4.6)，问这根细杆至多可有多长？又通道为园柱形的且细杆不必保持水平，细杆至多可有多长？ 5cm ? ? ? 17 证明当且仅当3 19 (Henon吸引子) 混沌和分形的著名例子，迭代模型为 图4.6 取初值x0 = 0, y0 = 0, 进行3000次迭代，对于k>1000, 在(xk, yk) 处亮一点(注意不要连线)可得 所谓Henon引力线图. 第一章 MATLAB入门 11 习题5 1．某河床的横断面如图5.8所示，为了计算最大的排洪量，需要计算它的断面积，试根据图示测量数据（单位：米）用梯形法计算其断面积。 2．求图5.8各测量点的坡度。 3．作图表示函数z?xe 4. 已知参数方程? ?x2?y3图5.8 ( -1 ?x?lncost?y?cost?tsint, 0 dydy和dxdx的数值解。 x??1 5. 求下列积分的数值解 (1) ?112?310e?x22dx, (2) ?2?0e2xcos3(x)dx ， (3) (6) ??xln(x4)arcsin111sin(x)1?xxdx , , (4), (5)dxdx??200xx2?0d??1?r2sin(?)dr，(7)??(1?x?y2)dydx, D为x2+y2?2x 0D 6 (椭园的周长) 用积分法计算下列椭园的周长 7.(曲面的面积) 求函数z?xe ?x2?y2 x2y2??1 49( -1 12 第一章 MATLAB入门 8 (假奇异积分)试求下列积分, 出现什么问题？分析原因，设法求出正确的解。 9 考虑积分I(k) = I= ?1?1x0.2cos(x)dx ?k?0sin(x)dx=2k, 试分别用trapz（取步长h=0.1或?）, quad 和quadl 求解I(8) 和I(32)。发现什么问题？ 10. (1) 用程序deriv.m求f(x)=x2sin(x2+3x-4)在x=1.3和x=1.5的导数，使精度达到10-3。 (2) 编写用公式(5.21)求函数在某一点二阶导数达到指定精度的算法程序，并用此程序求f(x)=x2sin(x2-x-2)在x=1.4的二阶导数，使精度达到10-3。 11图5.9a和图5.9b中各有两条曲线(粗线为x轴)，辨认每幅图中哪条是f(x)哪条是f (x)的导函数?为什么？ 图5.9b 图5.9a 12 (辛普生积分法)编制一个定步长辛普生法数值积分程序。计算公式为 I?Sn= h(f1+4f2+2f3+4f4+…+2fn-1+4fn+fn+1) 3其中n为偶数，h=(b-a)/n, fi=f(a+(i-1)h). 并取n=5，应用于解习题5(1)。 13 (摩托车)一个重5400kg的摩托车在以速度v=30m/s行驶时突然熄火，设滑行方程为 dv=-8.276 v2 - 2000 dxx为滑行距离，计算要滑行多长距离后, 速度可降至15m/s。 ? 14 一条长凳被牢牢固定在地上，凳面水平。考虑若干块砖在长凳一端叠成阶梯状而尽量向外延伸。一块砖放在长凳右端极端位置是砖的一半在外，但第二块砖若仍放一半(如图5.9)必 5400v 图5.9