数学建模讲稿1.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * 例2: 2019年全国大学生数学建模竞赛C题(酒 后驾车)中给出某人在短时间内喝下两瓶啤酒后,间 隔一定的时间t测量他的血液中酒精含量y(毫克/百 毫升) ,得到数据如表1所示。 根据微分方程模型得到短时间内喝酒后血液中酒精浓度与时间的关系为 现通过已有数据来拟合C_1,C_2,C_3 先建立拟合函数examp21 function f = examp21(c, tdata) f=c(1).*(exp(-c(2).*tdata)-exp(-c(3).*tdata)); 然后运行以下程序： tdata = [ 0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16 ]; ydata = [ 30 68 75 82 82 77 68 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4 ]; c0=[1 1 1]; for i = 1: 50 c = lsqcurvefit(examp21,c0,tdata,ydata); c0 = c End 求解为：c = 114.2472 0.1852 2.0126 （2）lsqnonlin用法：lsqsnonlin ( ′fun′, c0) 。 求含参量非线性函数fun中的参量c,使得各数据点函数值fun的平方和最小。 例如用lsqsnonlin ( ′fun′, c0)命令求解上例 的过程如下: 先建立拟合函数，注意跟上面的不同 function ff=examp22(c,tdata,ydata) tdata = [ 0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16 ]; ydata = [ 30 68 75 82 82 77 68 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4 ]; ff=c(1).*(exp(-c(2).*tdata)-exp(-c(3).*tdata))-ydata; 然后运行： c0=[1 1 1]; for i = 1: 50 c = lsqnonlin(examp22,c0); c0 = c end 求解为：c = 114.2472 0.1852 2.0126 例 估计水塔的水流量 某居民区的民用自来水是由一个圆柱形的水塔提供．水塔高 12.2米，直径17.4米．水塔是由水泵根据水塔内水位高低自动 加水，一般每天水泵工作两次．现在需要了解该居民区用水规律 与水泵的工作功率．按照设计，当水塔的水位降至最低水位，约 8.2米时，水泵自动启动加水；当水位升高到一个最高水位，约 10.8米时，水泵停止工作．可以考虑用用水率（单位时间的用 水量）来反映用水规律，并通过间隔一段时间测量水塔里的水位 来估算用水率．下表是某一天的测量记录数据，测量了28个时 刻，但是由于其中有4个时刻遇到水泵正在向水塔供水，而无水位记录（表1中用符号／／表示）试建立合适的数学模型，推算任意时刻的用水率、一天的总用水量和水泵工作功率． 表1 水位测量记录 时刻(h) 0 0.92 1.84 2.95 3.87 4.98 5.90 水位(cm) 968 948 931 913 898 881 869 时刻(h) 7.01 7.93 8.97 9.98 10.92 10.95 12.03 水位(cm) 852 839 822 // // 1082 1050 时刻(h) 12.95 13.88 14.98 15.90 16.83 17.93 19.04 水位(cm) 1021 994 965 941 918 892 866 时刻(h) 19.96 20.84 20.01 22.96 23.86 24.99 25.91 水位(cm) 843 822 // // 1059 1035 1018 问题分析与数据处理 由问题的要求，关键在于确定用水率函数，即单位时间内 用水体积，记为f(t)，如果能够通过测量数据，产生若干个 时刻的用水率，也就是f(t)在若干个点的函数值，则f(t)的计 算问题就可以转化为插值或拟合问题． 1．假设 1）水塔中水流量是时间的连续光滑函数，与水泵工作 与否无关，并忽略水位高度对水流速度的影响． 2）水泵工作与否完全取决于水塔内水位的高度. 3）水塔为标准圆柱体． 考虑到假设2）结合表1中具体数据，推断得出水泵第一次供 水时间段