浅析数值分析在数学建模中的应用

网友 分享 时间:

数值分析为数学建模提供了有效的计算方法,能够处理复杂问题,优化模型求解,提高预测精度,广泛应用于工程、经济、物理等领域。下面是阿拉网友收集整理的浅析数值分析在数学建模中的应用优秀范例,欢迎阅读参考,喜欢就支持吧!

一、数值分析在模型建立中的应用

在实际中,许多问题所研究的变量都是离散的形式,所建立的模型也是离散的。例如,对经济进行动态的分析时,一般总是根据一些计划的周期期末的指标值判断某经济计划执行的如何。有些实际问题即可建立连续模型,也可建立离散模型,但在研究中,并不能时时刻刻统计它,而是在某些特定时刻获得统计数据。另一方面,对常见的微分方程、积分方程为了求解,往往需要将连续模型转化成离散模型。将连续模型转化成离散模型,最常用的方法就是建立差分方程。

以非负整数k表示时间,记xk为变量x在时刻k的取值,则称xk=xk+1-xk为xk的一阶差分,称2xk=(xk)=xk+2-2xk+1+xk为xk的二阶差分。类似课求出xk的n阶差分nxk。由k,xk,及xk的差分给出的方程称为差分方程[2]。例如在研究节食与运动模型时,发现人们往往采取节食与运动方式消耗体内存储的脂肪,引起体重下降,达到减肥目的。通常制定减肥计划以周为时间单位比较方便,所以采用差分方程模型进行讨论。记第k周末体重为w(k),第k周吸收热量为c(k),热量转换系数,代谢消耗系数,在不考虑运动情况下体重变化的模型为w(k+1)=w(k)+c(k+1)-w(k)[2],k=0,1,2,,增加运动时只需将改为1+,1由运动的形式和时间决定。

二、数值分析在模型求解中的应用

插值法和拟合法在模型求解中的应用

1.拟合法求解

在数学建模中,我们常常建立了模型,也测量了(或收集了)一些已知数据,但是模型中的某些参数是未知的,此时需要利用已知数据去确定有关参数,这个过程通常通过数据拟合来完成。最小二乘法是数据拟合的基本方法。其基本思想就是:寻找最适合的模型参数,使得由模型给出的计算数据与已知数据的整体误差最小。

假设已建立了数学模型y=f(x,c),其中,c=(c1,c2,,cm)T是模型参数。已有一组已知数据(x1,,y1),(x2,y2),,(xk,,yk),用最小二乘确定参数c,使e(c)=ki=1(yi-f(xi,c))2最小。函数f(x,c)称为数据(xi,,yi)(i=1,2,,k)的最小二乘拟合函数。如果模型函数y=f(x,c)具有足够的可微性,则可用微分方程法解出c。最合适的c应满足必要条件e(c)cj=-2ki=1(yi-f(xi,c))f(xi,c)cj=0,j=1,2,,m。

2.插值法求解

在实际问题中,我们经常会遇到求经验公式的问题,即不知道某函数y=f(x)的具体表达式,只能通过实验测量得到该函数在一些点的函数值,即已知一部分精确的函数值数据(x1,,y1),(x2,y2),,(xk,,yk)。要求一个函数

yi=(xi),i=0,1,,k,(2)

这就是插值问题。函数yi=(xi)称为f(x)的插值函数。xi(i=0,1,,k)称为插值节点,式(2)称为插值条件[2]。多项式插值是最常用的插值方法,在工程计算中样条插值是非常重要的方法。

3.模型求解中的解线性方程组问题

在线性规划模型的求解过程中,常遇到线性方程组求解问题。线性方程组求解是科学计算中用的最多的,很多计算问题都归结为解线性方程组,利用计算机求解线性方程组的方法是直接法和迭代法。直接法基本思想是将线性方程组转化为便于求解的三角线性方程组,再求三角线性方程组,理论上直接在有限步内求得方程的精确解,但由于数值运算有舍入误差,因此实际计算求出的解仍然是近似解,仍需对解进行误差分析。直接法不适用求解n4的线性方程组,因此当n4时,可以采用迭代法进行求解。

迭代法先要构造迭代公式,它与方程求根迭代法相似,可将线性方程组改写成便于迭代的形式。迭代计算公式简单,易于编制计算程序,通常都用于解大型稀疏线性方程组。求解线性方程组的一般设计思想如下,假设建立一个线性规划模型

Ax=b

其中A=a11a12a1na12a22an2an1a12ann,x=x1x2xn,b=b1b2bn,即ARnn,可将A改写为迭代的形式

x=Bx+f

并由此构造迭代法

xk+1=Bxk+f,k=0,1,2,,

其中BRnn,称为迭代矩阵。将A按不同方式分解,就得到不同的迭代矩阵B,也就的带不同的迭代法,例如Jacobi迭代法[5]、高斯-赛德尔迭代法[5]、超松弛迭代法等。

由于计算过程中有舍入误差,为防止误差增大,就要求所使用的迭代法具有稳定性,即迭代收敛,收敛速度越快,误差越小。若x=Bx+f中,1,则认为此迭代法收敛。

4.数值积分在模型求解中的应用 模型求解过程中可能遇到积分求解问题,用求积公式If=bafxdx=Fb-Fa,使定积分计算变得简单,但在实际应用中很多被积函数找不到用解析时表示的原函数,例如10e-x2dx,或者即使找到表达式也极其复杂。另外,当被积函数是列函数,其原函数没有意义,因此又将计算积分归结为积函数值的加权平均值。

假设ax1b,则积分的计算公式[5]为bafxdxb-ani=0ifxi,称其为机械求积公式,其中xi(i=0,1,2,,n)称为求积节点,i与f无关,称为求积系数或权数,机械求积公式是将计算积分归结为计算节点函数值的加权平均,即取ni=0ifxi

得到的。由于这类公式计算极其便捷,是计算机计算积分的主要方法,构造机械求积公式就转化为求参数xi及i的代数问题。

5.数值分析在求解微分方程中的应用

在数学建模中,所建立的模型很多时候是常微分方程或者偏微分方程,这些方程求解析解是很困难的,而且即使能够求得解析解,由于所用数据的误差得到的解也是近似值,所以大部分情况下会采取数值的方法进行求解。

三、误差分析

在数学模型中往往包含了若干参变量,这些量往往是通过观察得到的,因此也带来了误差,这种误差称为观察误差[4]。这些误差是不可避免的,所以我们只能在模型建立和模型求解中避免误差扩大。目前已经提出的误差分析方法有向前误差分析法与向后误差分析,区间分析法,及概率分析,但在实际误差估计中均不可行。不能定量的估计误差,因此在建模过程中更着重误差的定性分析,也就是算法的稳定性分析。

在误差分析中,首先要分清问题是否病态和算法是否稳定,计算时还要尽量避免误差危害。为了防止有效数字的损失,应该注意下面若干原则:一是避免用绝对值小的数作除数;二是避免数值接近相等的两个近似值相减,这样会导致有效数字严重损失;三是注意运算次序,防止大数吃小数,如多个数相加减,应按照绝对值由小到大的次序运算;四是简化步骤,减少算术运算的次数。

四、结论

随着电子计算机的迅速发展、普及以及新型数值软件的不断开发,数值分析的理论和方法无论是在高科技领域还是在传统学科领域,其作用和影响都越来越大,实际上它已成为科学工作者和工程技术人员必备的知识和工具,所以把数值分析的知识正确的应用到数学建模中去不仅是一种趋势,更是用数学的理论解决实际问题的关键

48 3454513
");