清华数学实验复习试题二(龙格库塔方法_区间估计_约束优化).docVIP

计算方法（数学实验）试题（第2组） 2000.6.22 班级 姓名 学号 说明：（1）1，2题必做，答案直接填在试题纸上； （2）3，4题任选1题，将简要解题过程和结果写在试题纸上； （3）解题程序以网络作业形式提交，文件名用英文字母。 1. 设用数值解法算出 y(1)= 1.1635 ,你用的方法是 龙格-库塔方法 ，调用的 Matlab命令是： ts=0:0.1:2; y0=[1,0]; [x,y]=ode45(@cwf, ts,y0); ode45(@cwf, ts,y0) 算法精度为 4阶 。 %待解常微分方程组函数M文件源程序： function dy=cwf(x,y) dy=[y(2);y(1)*sin(x)]; %应用欧拉方法和龙格-库塔方法求解该常微分方程： ts=:0.1:2; y0=[1,0]; [x,y]=ode45(@cwf, ts,y0); %龙格-库塔方法求数值解 [, y(:,1)] 输出结果： 1.000000000000000 1.163536347595370 注意：ode45/23的步长必须从常微分方程的初值初开始，该命令默认为初值是在步长的起点赋值。 以下为错误程序： ts=-1:0.01:2; y0=[1,0]; [x,y]=ode45(@cwf, ts,y0); [x, y(:,1)] 2． 设总体，未知，现用一容量n=25的样本x对μ作区间估计。若已算出样本均值，样本方差，作估计时你用的随机变量是 ，这个随机变量服从的分布是 t(n-1) ，在显著性水平0.05下μ的的置信区间为 [15.441,17.359]. 若已知样本，对μ作区间估计，调用的 Matlab命令是： [mu, sigma, muci, sigmaci]=normfit(x,alpha) 不可省略其他项！！ 小型火箭初始质量为1200千克，其中包括900千克燃料。火箭竖直向上发射时燃料以15千克/秒的速率燃烧掉，由此产生40000牛顿的恒定推力。当燃料用尽时引擎关闭。设火箭上升的整个过程中，空气阻力与速度平方成正比，比例系数记作k。火箭升空过程的数学模型为 其中为火箭在时刻t的高度，m=1200-15t为火箭在时刻t的质量，T（=30000牛顿）为推力，g (=9.8米/秒2)为重力加速度, t1 (=900/15=60秒)为引擎关闭时刻。 今测得一组数据如下（t~时间（秒），x ~高度（米），v~速度（米/秒））： t 10 11 12 13 14 15 16 17 18 19 20 x 1070 1270 1480 1700 1910 2140 2360 2600 2830 3070 3310 v 190 200 210 216 225 228 231 234 239 240 246 现有两种估计比例系数k的方法： 1．用每一个数据(t,x,v)计算一个k的估计值（共11个），再用它们来估计k。 2．用这组数据拟合一个k。 请你分别用这两种方法给出k的估计值，对方法进行评价，并且回答，能否认为空气阻力系数k=0.5（说明理由）。 1 x=[190 200 210 216 225 228 231 234 239 240 246]; n=length(x); h=1; r(:,1)=(-3*x(:,1)+4*x(:,2)-x(:,3))/2/h; for i=2:n-1 r(:,i)=(x(:,i+1)-x(:,i-1))/2/h; end r(:,n)=(x(:,n-2)-4*x(:,n-1)+3*x(:,n))/2/h; r; t=[10 11 12 13 14 15 16 17 18 19 20]; T=40000;g=9.8; for p=1:11 m=1200-15*t(p); k(p)=(T-m*g-m*r(p))/x(p)^2; end k 4. Inter-Trade公司由中国大陆、菲律宾购买无商标的纺织品，运到香港或台湾地区进行封装和标签后，再运到美国和法国销售。已知两地间的运费如下（美元/吨）： 中国大陆 菲律宾 美国 法国 香港地区 55 72 160 190 台湾地区 67 58 150 210 现Inter-Trade公司从中国大陆和菲律宾分别购得90吨和45吨无标品。假