燃烧矿物燃料对碳循环影响的常微分方程模型.doc

该【燃烧矿物燃料对碳循环影响的常微分方程模型 】是由【帅气的小哥哥】上传分享，文档一共【6】页，该文档可以免费在线阅读，需要了解更多关于【燃烧矿物燃料对碳循环影响的常微分方程模型 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。燃烧矿物燃料对碳循环影响的常微分方程模型的求解燃烧矿物燃料对大气中二氧化碳有何影响?尽管今天二氧化碳的数量只占大气的大约百万分之350,但其任何增加都有可能对气候产生深远的影响。在LightHouseFoundation维护的网站[1]上有许多相关背景的文章介绍。[2]中提出的模型引入我们的视线,该模型模拟了许多不同形式碳原子的相互转化,包括存在于大气、浅海及深海的三种状态。模型中五个主要的变量都是时间的函数:p,大气中二氧化碳的分压力,浅海中集中溶解的碳,深海中集中的溶解的碳,浅海的碱度,深海的碱度另外在浅海的平衡方程中还涉及到三个量:,浅海中的碳酸氢盐,浅海中的碳酸盐,浅海中气态的二氧化碳的分压力五个主要变量的变化速率由五个常微分方程给出。大气和浅海之间的交换涉及到一个固定的特征传输时间d和源项f(t):下面的方程描述了浅海和深海之间的交换,其中vs和vd分别为两个区域的体积:二氧化碳和浅海中的碳酸盐之间的平衡方程由三个非线性代数方程描述:模型中涉及到的常量数值如下:源项f(t)描述了现代工业燃烧矿物燃料所造成的影响。时间间隔从一千年前到未来几千年后:在t=1000处的初始值:描述了前工业时代的平衡并保持近似不变,其中f〔t〕近似为零。下面的列表大致描述了燃烧矿物燃料,特别是汽油,所释放二氧化碳的源项f〔t〕。其数值在1850年后开始变得显著,在本世纪末到达峰值,然后降低直到燃料消耗完。。下半部的三个图线显示了大气、浅海及深海中碳的变化〔两个碱度值由于其根本保持不变,所以没有绘出〕。开始时,三个区域中的碳维持在平衡状态,所以其数量在1850年之前几乎没有变化。在1850≤t≤2500之间,图的上半部表示了燃烧矿物燃料所产生的碳进入系统,下半局部表示了系统的反响。大气是最先受到影响的,在500年后将增加四倍多,大约有一倍的碳进入浅海并最终转入深海区域内。重新绘制上图。使用pchiptx插入燃料表格,并使用ode23tx在默认的容差下求解此微分方程。比拟三个区域中5000年和1000年时碳的含量。大气中的二氧化碳在何时到达最大值?该方程是适度刚性的,这是由不同时间尺度内多种化学反响的发生造成的。放大图的某些局部,可以发现由于ode23tx采用了小步长所引起的特征锯齿形态。找到这样的一个区域。用MATLAB的常微分方程求解程序,包括ode23,ode45,ode113,ode23s,:odeset〔'RelTol',-6,'AbsTol',-6,'stats','on'〕;哪个方法更适合解决此问题?[1]LighthouseFoundation.[2],NumericalAdventureswithGeochemicalCycles,OxfordUniversityPress,NewYork,:pchiptx是Matlab函数pchip的一个可查看源代码的版本,包含在NCM程序包中。其作用是是实现三次埃米特插值,用于从函数的一组采样点生成该函数的一个近似的连续表达式。其用法可参看“〞刚性问题的定义可参考任何一本计算方法或微分方程数值解的书。也可参看“〞类似ode45这样的以ode开头的函数是Matlab中使用数值方法求解常微分方程的函数,其用法可通过在Matlab的命令窗口中键入helpode45获得,或参看“〞、“〞和“第十五章常微分方程的解法〞