多项式在建模中的应用.docVIP

多项式在建模中的应用 1 插值公式 在实际问题中，遇到研究的两个变量y与x之间的依赖关系，通过实验或观测可得到这两个变量对应数值：x取值时，y取值，i=0，1，2，…，n．我们希望找一个函数y=f(x)，使得，i=0，1，2，…，n，并且它能尽量准确地反映y与x之间的依赖关系，而计算又比较简单．变量y与x之间的依赖关系是客观存在的，设为．所找的函数f(x)称为的一个插值函数．求插值函数的问题称为插值问题，求插值函数的方法称为插值法． 实际问题中，若我们用某种办法已经知道变量y与x之间的依赖关系是一个很“光滑”的函数，则常常取多项式f(x)确定的函数作为这个函数的插值函数，此时f(x)称为插值多项式．于是，现在的插值问题是：设是数域F上的n+1个不同元素，是F的n+1个元素．我们要在F[x]中找一个次数? n的多项式f(x)，使得，i=0，1，2，…，n．由推论4.5.1，若这个问题的解存在，则它是唯一的．现在我们来证明插值问题的解一定存在． 定理1 设是数域F的n+1个不同的元素,是F的n+1个元素，则在F[x]中存在唯一的次数不超过n的多项式f(x)，使得 ， i=0，1，2，…，n 证 如上所述，若这个插值问题的解存在，则它是唯一的，现在来证存在性． 先看一个特殊情形： 若存在一个次数? n的多项式，使得 ， j=0，1，2，…，n 则 ，…，， ，…，． 这表明是的n个不同的根．由于，所以 (1) 因为，所以由(1)式得 ， (2) 从而得 ． (3) 现在看一般情形．令 ＝， (4) 则degf(x)? n，并且对任意j(j=0，1，2，…，n)有f(cj)=dj ．所以由(4)式给出的f(x)是插值问题的解． ? 公式(4)称为拉格朗日(Lagrange)插值公式． 定理也可以采用待定系数法来证明，设 是插值问题的解，由，i=0，1，2，…，n，可以得到一个含n+1个未知量的由n+1个方程组成的线性方程组，它的系数行列式是Vandermonde行列式．由于是n+1个不同的数，所以这个行列式不等于零．据Cramer法则，该方程组有唯一解，并且解有一个公式表示(这个解也可以用消元法求出)．当然Lagrange插值公式更为方便，因为它既简单又容易记． 上述插值问题的解还可以用牛顿(Newton)插值公式给出： ．(5) 其中系数可以通过把x逐次用代入(5)式求出． 例1 求一个次数不超过3的多项式f(x)，使f(0)=5，f(1)=7，f(? 1)=9，f(? 2)=13 解 用拉格朗日插值公式得 求得插值多项式后，结合分析方法便可解决相关的实际问题 例2 已建立某工厂一个中等水平工人上午上班(8:00开始)的工作效率多项式Q(t)= ? t3+9t2+12t．由之可得这样工人工作效率的高峰时刻在上午11:00． 又如已建立飞机的降落曲线为y=ax3+bx2+cx+d，可指导飞机安全降落． 2 多项式四归(略述) 3 三次样条插值 过去，工程师制图、绘制经过给定点的曲线时，是先用弹性薄木片弯曲使得过给定点，然后画下薄木片弯成的曲线．后来证明，这样得到的曲线就是三次样条插值曲线． 定义 若函数S(x)适合下列条件： 1) S(x)? C (2)[a，b] 2)在每个小区间[xk，xk+1]上，S(x)=Sk(x)是三次多项式； 3)在节点xk上满足S(xk)=f(xk)，k=0，1，…，n， 则S(x)叫做[a，b]上的三次样条插值函数． 这里不介绍三次样条插值的理论陈述(同学们可参考有关数学建模的书)，只简要提一个实例：水塔水量的估计曲线． 美国某州的各用水管理机构要求各社区提供本区每小时多少加仑计的用水率以及每天总的用水量．但许多社区并没有测量流入或流出当地水塔的水量的设备，他们只能代之以每小时测量水塔中的水位，其精度在0.5%以内．更为重要的是，无论何时，只要水塔中的水位降到某一最低水位L时，水泵就启动向水塔重新充水至某一最高水位H，但也无法得到水泵供水量的测量数据．因此，在水泵正在工作时，人们不容易建立水泵中的水位与水泵工作时的用水量之间的关系．水泵每天向水塔充水一次或两次，每次约二小时． 试估计在任何时刻，甚至包括水泵正在工作的时间内，水从水塔流出的流量f(t)，并估计一天的总用水量，表1给出了某小镇一天中的真实数据，其中，从第一次测量开始的时间均以秒为单位，水位测量值为百分之一英尺．例如3316秒后，水塔中水位达到31.10英尺，水塔是一个高40英尺．直径57英尺的直立圆柱．通常当水泵水位降至约27.00英尺时，水泵开始工作，当水位升到35.50英尺时水泵停止工作． 表1 某小镇某天的水塔水位 时间(秒) 水位 (0.01英尺) 时间(秒) 水位 (0.01英尺) 时间(秒) 水位 (0.01英尺) 0 331