非自治污染模型的建立与处理.pdf

非自治污染模型 赵钊 5081209401 潘健 5081209393 全球大气中的二氧化碳存量会对全球变暖产生正面影响。这一存量会随着燃 烧化石燃料所带来的二氧化碳的工业排放而增加，但它也会由于海洋和植被的自 然吸收而出现一定数量的下降。在本模型中，我们将会考察全球二氧化碳存量的 动态演变。随着时间的推移，这一存量倾向于随着全球经济活动扩张引发的工业 排放量的增加而增加。另一方面，随着全球变暖问题的恶化，环境污染控制也越 来越严格，这一存量也受到了这个因素的负面影响。 令y 表示二氧化碳的存量，假定存量根据如下方程变化： . y x y （1） 其中，x 表示二氧化碳的工业排放量， ＞0 为参数，它确定了自然环境吸 收二氧化碳的速度。进一步假设，工业排放量随时间的变化如下： . x aebt y （2 ） 其中，a,b和 为正常数。上式中的第一项使排放量随着 t 增加，该项表示 经济活动扩张对二氧化碳工业排放量的影响。该项体现了如下假设：随着污染问 题的恶化，政府对工业二氧化碳的排放量进行了更为严格的控制。 方法1：手工计算 对方程（1）两边微分得到 .. . . y x y . 用方程（2 ）替换x ，重新整理后得到y 的二阶微分方程： .. . y y y ae bt 这是一个可以求解的二阶线性非自治方程，方程的齐次形式为 .. . y y y 0 特征方程为 r2 r  0 特征根为 2   4 r , r  1 2 2 2 因为，已经假定和为正数，所以如果两个根为实根，那么他们都是负的； 如果两个根都是复根，那么实部是负的。齐次形式的解为 r t r t y C e 1 C e 2 h 1 2 为求特解，我们注意到可变项为t 的函数，因此，我们尝试相同形式的特解 y A e bt p 0 为确定系数A0 ，对上式求微分得到 . y bA e bt p 0 再次微分可得 .. 2 bt y b A e p 0 现在将所猜测的特解和它的两个导数代回完全方程。可得 2 bt bt bt bt b A e bA e A e ae 0 0 0 a 求出A0 A0 2   b b   因此，特解为 a bt y p 2 e b  b    微分方程的通解为 a r t r t bt y t C e 1 C e 2  e  1 2 2 b  b    方法2 ：用Maple 求解： 结论： 该解向我们说明了模型中的二氧化碳存量如何随着时间而变化。我们已经确定两 个根都是负数，因此，随着t  ，解中的