Matlab与数学实验梁锦桐2007160036

1、bu深圳大学考试答题纸(以论文、报告等形式考核专用)二 九 二 十 学年度第 一 学期课程编号2316000801课程名称Matlab与数学实验主讲教师仇高新评分学 号2007160036姓名梁锦桐专业年级07微电一班教师评语：题目：Matlab与数学实验期末作业1. 用数值积分公式计算 (结果保留小数点后8位)：(1) 取积分步长, 用梯形公式计算S= 6.24764132 。>> format long;>> x=0:pi/2:2*pi;>> trapz(x,(1-(0.15.2)*(sin(x).2).0.5)(2) 要求相对误差为10-6, 用Sim

2、pson公式S= 6.24769188 ，Matlab命令是_ quad('(1-(0.15.2)*(sin(x).2).0.5',0,2*pi)_.>> format long;>>quad('(1-(0.15.2)*(sin(x).2).0.5',0,2*pi)2. 设用数值解法算出 y(1)= 1.1635 ,你用的方法是 Runga-Kutta 方法 ，调用的 Matlab命令是 x,t=ode45('verderpol',0:0.1:1.2,1,0) ，算法精度为 4阶 。解：建立函数文件verderpol.mf

3、unction xprime = verderpol(x,t) xprime = t(2); t(1)*sin(x);运行命令：x,t=ode45('verderpol',0:0.1:1.2,1,0)3. 设 用数值解法算出y(1)= 0.2713 （精确到4位小数）, 你用的方法是 Runga-Kutta 方法 ，调用的 Matlab命令是 x,t=ode45('verderpol2',0:0.1:1.2,1,0) ，算法精度为 45阶 。解：建立函数文件 verderpol2.mfunction xprime = verderpol2(x,t) xprime

4、 = t(2);t(2)*sin(x)-t(1)*exp(x);运行命令：x,t=ode45('verderpol2',0:0.1:1.2,1,0)然后在Matlab中用命令：x,t=ode45('verderpol2',0:0.1:1.2,1,0)4. 用电压V=14伏的电池给电容器充电，电容器上t时刻的电压满足：，其中是电容器的初始电压，是充电常数。试用下列数据确定和。t（秒）0.30.51.02.04.07.0v(t)5.68736.14347.16338.862611.032812.6962 你用的方法是 线性最小二乘法 ，结果是= 5.0001 ，=

5、3.6165 。5. 小型火箭初始质量为900千克，其中包括600千克燃料。火箭竖直向上发射时燃料以15千克/秒的速率燃烧掉，由此产生30000牛顿的恒定推力。当燃料用尽时引擎关闭。设火箭上升的整个过程中，空气阻力与速度平方成正比，比例系数为0.4（千克/米）。重力加速度取9.8米/秒2.A. 建立火箭升空过程的数学模型（微分方程）；B. 求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点的时间和高度。解：A.火箭升空应分为两个过程：1.有燃料产生向上推力的过程：2.燃料用完，引擎关闭的运动过程：对于第一个过程：持续的时间为：t0=60015=40s向上加速过程火箭的质量为：M=900-

6、15*t (1) t40sF空气阻力与速度平方成正比，即 f=k*v2，方向是竖直向下根据牛顿第二定律得在加速过程有：F-f-M*g=M*a即F-k*v2-M*g=M*a (2) 又有： v=dydt (3 ) 和 a=d2ydt2 (4)联立（1）、（2）、（3）和（4）得：MgF-k*dydt2-900-15*t*g=900-15*t*d2ydt2f代入数值得：向上加速过程30000-0.4*dydt2-900-15*t*9.8=900-15*t*d2ydt2即:d2ydt2=30000-0.4*dydt2-900-15*t*9.8/900-15*t , 0 t40s初始条件为：y(0)=

7、0 ; y(0)=0对于第二个过程：t>40如右图所示：火箭只受到重力和阻力的作用由于燃料已经用完,则剩余质量为：m=900-600=300kg对物体进行受力分析得：-f-mg=ma (5)即：-k*dydt2-m*g=m*d2ydt2 (6)mg代入数值得 ：f-0.4*dydt2-300*9.8=300*d2ydt2 t>40 （7）即： d2ydt2=-0.4*dydt2/300-9.8 t>40 （8）初始条件由第一个过程的终值给出。 B.在Matlab中计算得，引擎关闭瞬间火箭的高度8323米，速度259米/秒，引擎关闭前瞬间加速度0.7709米/秒2，引擎关闭后瞬

8、间加速度为，99.2291米/秒2；到达最高点的时间51秒，高度9192米。M文件为，answer5.m6. 冰淇淋的下部为椎体，上部为半球。设它由锥和球面围成，用蒙特卡罗方法计算它的体积。解：方程是一个以（0，0，1）为球心，半径为1，所求的椎体和球面围成的体积包含在球体里面，设计一个与球相切的正方体，边长为2，正方体的体积为8. 在正方体内随机的投点，则点落在圆锥和球面围城的体积里面的概率应该等于圆锥和球面围城的体积和正方体的体积之比。假说总投点数为n,落在圆锥和球面围城的体积内的点数为M，则应有： VV正=mn =>V=mnV正在matlab里面编写随机点进行试验求出所求的体积为：

9、3.1315 ，程序如下：answe6.mn=100000;a=0;b=1;m=0;for i=1:n x=rand(1)*2-1;y=rand(1)*2-1;z=rand(1)*2; if(x2+y2)0.5<=z && x2+y2+(z-1)2<=1) m=m+1; endendfprintf('计算出来的体积为：%fn',8*m/n);7. 容器盛满水后，低端直径为的小孔开启。根据水力学知识，当水面高度为时，水从小孔中流出的速度（为重力加速度，0.6为孔口收缩系数）。若容器为倒圆锥形，现测得容器高和上底面直径均为1.2m，小孔直径为3cm，问水

10、从小孔中流完需要多少时间？2分钟时水面高度是多少？解：容器中总共的水的体积为：Vtot,经过t时间后，假设从容器中流出的水的体积为Vout,剩余水的体积为：Vleft ,高度为h; ; 则有：Vout=Vtot-Vleft; (1)而：Vout=*sdt （2） ; Vleft=13h*S=13h*(h2)2=112*h3; （3）将（2）、（3）代入（1）得*sdt=Vtot-112*h3; （4）对式（4）两边求导得：*(d02)2=-112*3*h2dhdt （5）即：0.6*2*g*h*(d02)2=-14*h2dhdt （6）代入数值进一步化简得：dhdt=-5.4*10-4*2*9.8*h/h2 (7