插值和拟合概要.ppt

* 供水，上升到35.5英尺时立即停止供水。不考虑水泵因其它原因而中断或需要维修等情况。 (4) 假设从水塔流出的水流速度比水泵供水速度小，以保证不存在断水情况。 (5)假设水的消耗每天相同,每天从水塔中流出的水流速度随时间的变化可以用光滑曲线近似。从统计意义上来看，由于每个用户的用水量远远小于社会群体用水量，以致于个别用户的突然用水基本上不能改变水塔的流水速度。 建立数学模型 (1) 水塔的形状、尺寸都已经清楚地给出，当水泵不供水时，根据水位随时间的变化可以确定从水塔中流出水流的速度 。 * 利用差商确定 ， 表示函数 关于 的函数 表2 观测时间中点与该时间区间平均水流速度 时间中点 区间平均流速 时间中点 区间平均流速 时间中点 区间平均流速 (小时) (×103加仑/小时) (小时) (×103加仑/小时) (小时) (×103加仑/小时) 0 无数据 9.47 无数据 18.49 12.2 0.46 11.2 10.45 无数据 19.50 12.9 1.38 9.7 10.94 无数据 20.40 12.6 2.4 8.6 11.49 15.6 21.43 无数据 3.41 8.1 12.49 16.4 22.49 无数据 4.43 9.3 13.42 15.5 23.42 无数据 5.44 7.2 14.43 13.4 24.43 11.2 6.45 7.9 15.44 13.8 25.32 3.5 7.47 7.4 16.37 12.9 8.45 8.4 17.38 12.2 * (2) 用三次样条曲线拟合流速曲线 利用流速的光滑性假设,采用光滑的三次样条插值近似表示流速曲线,得到如图的插值曲线。 (3) 计算总用水量 通过三次样条插值,可以获得水塔的流速V与时间的关系式V(t),据此, 我们可以求出这个小 镇一天的总的用水量 三次样条插值曲线 * 利用MATLAB软件进行插值计算 一维插值 yi=interp1(x,y,xi,’method’) Method: linear, spline, cubic 高维插值 zi=interp2(x,y,z,xi,yi,’method’) * 加上任何一类边界条件(至少两个)后 一般使用第一、二类边界条件, 即 ------(8) 或 第二类边界条件更常用 * ------(9) * 加以整理后可得 ------(10) ------(11) 未知 * 由条件 由于以上两式相等,得 * ------(12) (12)式称为基本方程组 * 如果问题要求满足第一类(一阶)边界条件: ------(5) 基本方程组(12)化为n-1阶方程组 ------(13) 即 将(13)式化为矩阵形式 * ------(14) 这是一个三对角方程组 如果问题要求满足第二类(二阶自然)边界条件: ------(6) * 由(11)式,可知 ------(15) ----(16) * ------(17) ------(18) 与基本方程组(12)联合,并化为矩阵形式,得n+1阶方程组 -----(19) 与(14)比较 * (19)式与(14)一样,都是三对角方程组,并且都严格对角占优 可以使用追赶法求解