技术文章及通讯

大气二氧化碳模型和曲线拟合工具箱

作者:Mary Ann Branch, MathWorks和Russ Minkwitz, MathWorks


研究燃烧化石燃料对天气模式影响的常用测量方法是二氧化碳水平2)大气中的浓度。一些科学家相信CO的上升趋势2这种水平可能导致大气温度上升,极地冰盖融化,以及地球不同地区的气候发生根本变化。

由于这些数据具有周期性现象,以及上升的增长趋势,非线性拟合技术是合适的,但与拟合线性模型相比,可能需要额外的努力。然而,使用正确的工具和数据知识,您可以快速确定非线性模型的系数。

本文描述了如何将一个非线性模型拟合到可用的CO2数据使用新发布的曲线拟合工具箱。

数据建模

在这个例子中,我们使用CO2美国国家海洋和大气管理局(NOAA)从1979年到1996年在夏威夷的莫纳罗亚天文台收集了以百万分率为单位的数据2.捕获的数据如下图所示。

modeling_fig1_w.gif
在曲线拟合工具箱中找到的曲线拟合工具提供了数据的即时可视化显示。通过使用数据提示功能,我们可以询问特定的点。

自然界的循环现象可以用正弦信号及其高次谐波来模拟,即,

modeling_logo_eq1.gif
一个kT / K而且Φk分别是第k次谐波的振幅、周期和相位。正弦项捕获季节性趋势,但持续上升的趋势需要指数项。通过对19世纪60年代中期以来世界化石燃料总产量的记录进行检查,得出了一个建议的数据模型
modeling_logo_eq2.gif

其中t(时间)以年为单位,α被认为约为0.0444/年,其他参数未知(包括所需正弦波的数量)。2

为了构建模型,我们转到曲线拟合工具箱的拟合GUI并创建一个自定义方程,其中包含一个正弦项。

modeling_logo_eq3.gif

由于时间值很大,我们可以选择将时间变量居中并缩放(删除平均值并除以标准偏差)的选项,以改进数值计算的调节。

modeling_fig2_w.gif
在自定义方程面板中,我们输入初始值α = .0444, T = 0.192(定心和缩放后的一年时间),并对其余参数进行随机猜测(如下)。当拟合完成并绘制结果时,我们看到循环模式没有被匹配(左)
modeling_fig3_w.gif
显示的结果表明,r平方值(拟合优度的衡量标准)离1不远,这意味着全球趋势得到满足,但局部趋势没有。在这种情况下,正弦项的振幅太小。

细化模型

我们可以在几个方面改进这个模型。首先,我们可以分别确定指数项的系数。我们使用工具箱内置的指数模型,它会自动从数据中计算出最佳起点,

modeling_logo_eq4.gif
拟合该模型得到的系数为D = 349.6而且α = 0.02185.使用以下值da作为正弦模型的起始值,我们得到了一个新的拟合,很好地匹配了循环性质。
modeling_fig4_w.gif

这个最新模型的振幅大约是3。我们可以在图上验证这是a1(左)的合理值。结果面板(如下)还在系数值后面的括号中显示了每个参数的置信区间。较小的置信区间意味着参数值的可信度更高。除B和d外,所有参数的置信区间都很小,这可能表明B和d高度相关。现在我们可以画出残差仍然显示出周期性正弦模式,所以另一个正弦项是有保证的。

现在我们可以画出残差仍然显示出周期性正弦模式,所以另一个正弦项是有保证的。


modeling_fig5_w.gif
单正弦模型的残差揭示了另一个正弦项(左)。有两个谐波,配合得很好(右)。我们可以很容易地用更多的项创建方程,检查每个结果拟合的残差,直到我们满意为止。

扩展分析

当我们找到一个令人满意的模型时,我们可以将我们的拟合发送到MATLAB工作空间进行进一步分析。我们还有其他几个观测点要检查(例如,萨摩亚,南极和巴罗),因此可以生成一个m文件函数,可以重新创建我们的拟合。然后,我们可以在循环中调用该函数来自动拟合其他数据集,以查看结果模型在不同观测站之间的比较情况。

这个例子涉及到一些可以帮助您将模型与数据匹配的特性。无论模型是线性的还是非线性的,曲线拟合工具箱都使您能够从其库或您创建的自定义模型中快速拟合标准模型。

modeling_fig6_w.gif
可以从GUI File菜单生成用于复制在GUI中创建的拟合和绘图的m文件代码。m文件还显示了可以在命令行上使用的工具箱命令,以获得相同的结果。

发布于2002年- 91034v00

参考文献

  1. GLOBALVIEW-CO2,“合作大气数据集成项目-二氧化碳”(2001)。[公众可在ftp.cmdl.noaa.gov,路径:ccg/co2/GLOBALVIEW]

  2. 卡纳,D.,莫勒,C.,纳什,S.,数值方法和软件,普伦蒂斯霍尔(1989)

s manbetx 845产品使用

查看相关功能的文章