清华数学实验test2005AB.docVIP

考试课程 数学实验 2005.6.15下午 班级 姓名 学号 得分 [说明] （1）第一、二、三题的答案直接填在试题纸上； （2）第四题将数学模型、简要解题过程和结果写在试题纸上；卷面空间不够时，请写在背面； （3）除非特别说明，所有计算结果小数点后保留4位数字。 （4）考试时间为120分钟。一、（10分）某厂生产A、B两种产品，1千克原料在甲类设备上用12小时可生产3件A，可获净利润64元；在乙类设备上用8小时可生产4件B，可获净利润54元。该厂每天可获得55千克原料，每天总的劳动时间为480小时，且甲类设备每天至多能生产80件A。试为该厂制订生产计划使每天的净利润最大。 以生产A、B产品所用原料的数量x1、x2（千克）作为决策变量，建立的数学规划模型是： min z=-64*x1-54*x2 s.t. 12x1+8x2≤480； x1+ x2≤55,； 0≤3x1≤80；0≤4x2 每天的最大净利润是____3070____元。若要求工人加班以增加劳动时间，则加班费最多 为每小时_____2.5__元。若A获利增加到26元/件，应否改变生产计划？__不改变________ 二、(10分) 已知常微分方程组初值问题 试用数值方法求_______1.73205_________(保留小数点后5位数字)。你用的MATLAB命令是________ode45_________________，其精度为______5级4阶____ (10分) 已知线性代数方程组Ax=b, 其中 , , , 若方程组右端项有小扰动，试根据误差估计式估计___0.0743______（分别表示原问题的解和右端项小扰动后对应的解的变化量）；若取初值，则用高斯-赛德尔迭代法求解Ax=b时，_[1.7160, 0.3926, -0.1306, 0.1381]’_______；对本题而言，此迭代方法是否收敛_____收敛__，原因是______谱半径=0.39681____________。 四、（20分）炮弹射击的目标为一椭圆形区域，在X方向半轴长110m，Y方向半轴长90m.当瞄准目标的中心发射炮弹时，在众多随机因素的影响下，弹着点服从以目标中心为均值的二维正态分布，设弹着点偏差的均方差在X方向和Y方向分别为70m和50m。今测得一组弹着点的横纵坐标如下： X -6.3 -71.6 65.6 -79.2 -49.7 -81.9 74.6 -47.6 -120.8 56.9 Y 28.9 1.6 61.7 -68 -41.3 -30.5 87 17.3 -17.8 1.2 X 100.9 47 9.7 -60.1 -52.7 86 80.6 -42.6 56.4 15.2 Y -12.6 39.1 85 32.7 28.1 -9.3 -4.5 5.1 -32 -9.5 根据这组数据对X方向和Y方向的均值和均方差进行假设检验（设显著性水平为0.05）。 根据这组数据给出随机变量X和Y相关系数的一个点估计。 用蒙特卡罗方法求炮弹落在椭圆形区域内的概率（取10000个数据点；请附程序）。 x=[-6.3 -71.6 65.6 -79.2 -49.7 -81.9 74.6 -47.6 -120.8 56.9 100.9 47 9.7 -60.1 -52.7 86 80.6 -42.6 56.4 15.2]; y=[28.9 1.6 61.7 -68 -41.3 -30.5 87 17.3 -17.8 1.2 -12.6 39.1 85 32.7 28.1 -9.3 -4.5 5.1 -32 -9.5 ]; X`s mean h1=ztest(x,0,70) Y`s mean h2=ztest(y,0,50) X`s sigma h3=chi2test(x,70) Y`s sigma h4=chi2test(y,50) cor=corrcoef(x,y) pause; %均值法 a=1.1;b=0.9;sx=0.7;sy=0.5; n=10000; m=0;z=0; x=unifrnd(0,1.1,1,n); y=unifrnd(0,0.9,1,n); for i=1:n if x(i)^2/a^2+y(i)^2/b^2=1 u=exp(-1/2*(x(i)^2/sx^2+y(i)^2/sy^2))