数学建模讲稿插值拟合方程求根.pptVIP

? 根据微分方程模型得到短时间内喝酒后血液中酒精 浓度与时间的关系为 现通过已有数据来拟合 C_1,C_2,C_3 先建立拟合函数 examp21 ? function f = examp21(c, tdata) ? f=c(1).*(exp(-c(2).*tdata)-exp(-c(3).*tdata)); ? 然后运行以下程序： 3 2 1 ( ) c c y c e e ? ? ? 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 中具体数据，推断得出水泵第一次供 水时间段为［ 9 ， 11 ］，第二次供水时间段为「 20.8 ， 23 ］．