数学模型与数学建模实验一.docx

实验报告一学院名称：理学院 专业年级： 姓 名： 学 号：课 程：数学模型与数学建模 报告日期：2015年11月10日一、实验题目例1.1.1 人体的体重与身高 利用下表中的一组实际测量的数据，分析人体体重W与身高L的关系。W/kg1217223548546675L/m0.861.081.161.351.551.671.781.85例：1.1.4圆珠体积 利用累次积分、示性函数积分、蒙特卡洛模拟方法计算一个直径为2cm的球状，带有穿过球心，直径为0.2cm的圆柱体的圆珠的体积，并在达到相同精度计算结果下比较这几种方法的计算时间。二、实验目的1、了解MATLAB的基本功能，配合介绍数学建模的最基本的要素，进行误差分析，选择合适的数学模型。2、为了更有效的运用MATLAB指令，简单运用MATLAB编程方法，以四种不同的方式计算圆珠体积，通过比较不同程序的计算用时，通过比较，得到计算用时最少，精度能得到要求的算法。三、问题陈述1、人体的身高与体重 利用下表中的一组实际测量的数据，分析人体体重W与身高L的关系。W/kg1217223548546675L/m0.861.081.161.351.551.671.781.852、圆珠体积的计算 利用累次积分、示性函数积分、蒙特卡洛模拟方法计算一个直径为2cm的球状，带有穿过球心，直径为0.2cm的圆柱体的圆珠的体积，并在达到相同精度计算结果下比较这几种方法的计算时间。四、模型及求解结果1、人体的体重与身高解：根据质量与体积、体积与长度之间的量纲关系，可以建立轮廓模型 关键在于参数的估计。利用实际测量的数据估计参数，有两种方法。估计方法1：用各观测数据的比值 的平均数 来估计估计方法2：用观测数据的均值的比值估计： 两种参数估值不同，估计的好坏可以通过误差平方和的大小进行比较，计算结果通过图像直观地观察，结果如图： （a）估计方法1 （b）估计方法2从数值结果比较，第一种参数估计方法的误差平方和为Q1=342.5765，小于第二种参数估计方法的误差平方和Q2=814.3192，因此结果更好一些，但是它仍然不是使得误差平方和最小的估计，运用拟合模型可以重新 这个问题。2、圆珠体积将求体积的三重积分转为累次积分可以直接得到圆珠的体积 解法1：换成极坐标积分，采用符号函数积分，运行程序在下面的代码中，计算用了61.083812 seconds解法2：用累次积分公式，采用数值函数定积分指令。因为定积分只能计算积分上下限为常数的定积分，虽然可以计算含参量的一重积分，但是再累次积分就只能利用数值梯形积分公式近似计算。通过代码运行可得：计算用时116.902487 seconds解法3：采用示性函数积分方法，集合的示性函数定义为： 因此，很容易采用矩阵的逻辑运算方式定义区域的实行函数，又因为该圆珠可以放在立方体内，由在立方体上的积分就可以得到的体积。计算用时65.264352 seconds.解法4：利用蒙特卡洛积分，在区间内按均匀分布规则书n个点，由落在区域内点的个数s，可以近似计算在中所占的比例，由得体积8，于是可以得到的近似值，计算用时为31.230011 seconds.五、程序代码1、人的升高和体重 W=[12,17,22,35,48,54,66,75]; L=[0.86,1.08,1.16,1.35,1.55,1.67,1.78,1.85]; k1=mean(W./L.^3)%估计方法1.k1 = 13.5892 k2=mean(W)/(mean(L))^3%估计方法2.k2 = 14.5929 plot(L,W,o,L,k1*L.^3,-k,LineWidth,2); figure(2) xlabel(L) ylabel(W) plot(L,W,o,L,k2*L.^3,-k,Linewidth,2); Q1=(W-k1*L.^3)*(W-k1*L.^3)%误差平方和Q1=342.5765Q1 = 342.5765 Q2=(W-k2*L.^3)*(W-k2*L.^3)%误差平方和Q1=814.3192.Q2 = 814.31922、圆珠体积解法1：换成极坐标积分，采用符号函数积分 tic int(int(sym(2*r*sqrt(1-r^2)),0.1,1),0,2*3.14)ans =(15543*11^(1/2))/12500 vpa(ans,6)ans =4.12402 tocElapsed time is 61.083812 seconds.解法2：用累次积分公式 tic f=@(z)quad(@(r)r,0.1,sqrt(1-z.^2)); n=40; x=linspace(0,sqrt(0.99),n); for j=1:n f1(j)=4*3.14*f(x(j));