δ函数与微分方程的建模.ppt

* 第 六章 δ函数与微分方程的建模 微分方程建模是数学建模的重要方法，因为许多实际问题的数学描述将导致求解微分方程的定解问题。把形形色色的实际问题化成微分方程的定解问题，大体上可以按以下几步： 1. 根据实际要求确定要研究的量（自变量、未知函数、必要的参数等）并确定坐标系。 2. 找出这些量所满足的基本规律（物理的、几何的、化学的或生物学的等等）。 3. 运用这些规律列出方程和定解条件。 列方程中常见的方法有： 1. 按规律直接列方程（如牛顿第二定律、放射性物质的放射性规律等）。 2. 微元分析法与任意区域上取积分的方法。 3. 模拟近似法（在一定的假设下，给出实际现象所满足的规律，然后利用适当的数学方法列出微分方程）。 在实际的微分方程建模过程中，往往是上述方法的综合应用。 第一节 传染病的流行与宣传运动的关系的数学描述 开展预防传染病（例如艾滋病等）的流行的宣传运动对防止传染病的蔓延起多大作用？这个宣传运动要持续多长时间？要具有多大的强度？ 这些问题的讨论将导致含有间断的非齐次项或非齐次项含有δ函数的一阶线性常微分方程。 我们从最简单的情形——不开展宣传运动的情形开始。 最简单的模型——不开展宣传运动的情形 设总人数为N是不变的，t时刻得病人数为x(t)，它传染给正常人的传染率为r。 显然，从t到t+Δt时间内平均传染率为 t时刻的传染率为 从而微分方程为 求解得 令 得 这表明：最终每个人都要染上疾病。 我们得到最简单的数学模型： 可分离变量 持续宣传的作用 为了预防传染病的流行，宣传有如在人们耳边敲起警钟“要当心！”因而是很重要的。 假设：宣传的开展将使得传染上疾病的人数x(t)减少，减少的速度与总人数N成正比，这个比例常数称为宣传强度。 若从t=t00开始，开展一场持续的宣传运动，宣传强度为a，则所得的数学模型为 其中 即Heaviside函数。 （单位阶梯函数） 解法1 分段求解 当 时， 当 时， 解一阶线性非齐次方程的初值问题： 得 两段表达式合起来，得 解法2 用拉氏变换法求解 的Laplace变换 的Fourier变换 令 的Laplace变换为 ，对方程两边求拉氏变换并利用初值得 求逆变换得 令 得 这表明：持续的宣传是起作用的，最终会使发病率减少（0ar）。 进一步地，如果我们从一系列时刻t=t1,t2,…,tm开始，分别进行强度为a1,a2,…,am的宣传运动，则得到如下的数学模型： 同样可用分段求解法求解上述初值问题，但计算比较繁琐。若用拉氏变换法求解，则是十分方便的。 令 得 解得 顺时宣传情形的数学模型 如果宣传运动是短暂进行的，这在日常生活中是常见的，例如仅仅是听一个报告，或街头散发传单等等，即在t=t1,t2,…,tm等m个时刻进行m次宣传，宣传强度分别为a1,a2,…,am，则数学模型变成 这里我们用 来数学的近似描述时间极为短暂的宣传。这时，模型中含有顺时点源项。 用拉氏变换求解这类方程是方便的。 求解得 令 得 这表明：短暂的宣传（即使是多次的短暂宣传）是不起作用的，最终还是所有的人都传染上疾病。 思考：实际情况，即使短暂的宣传也往往会在人们的头脑中留下一定时间延续的印象，应该怎样描述呢？请大家自己作为练习来建模。 第二节 附录——δ函数及其性质 一、δ函数的物理模型 对于连续分布的量，一般用密度函数来描述。现设有一条无穷长的杆，沿杆建立了一维坐标系，点的坐标为x，杆的分布密度为ρ(x)。(-∞,x]段上杆的质量为m(x)，则它们的关系是 对于集中分布的量，如何来描述它的密度？ 仍以无穷长的杆为例。设无穷长杆的质量为1，质量集中于点x0，则(-∞,x]段杆的质量为 若仍形式地套用连续变化时的ρ(x)来描述它的密度，则应该有 且 特别有