技术文章和通讯

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

玛丽安的分支,MathWorks Russ Minkwitz, MathWorks


常见的测量为研究燃烧化石燃料对气候模式的影响是二氧化碳的水平(有限公司2)在大气中的浓度。一些科学家认为公司的上升趋势2水平可能导致大气温度上升,极地冰盖融化,地球和气候不同的地区从根本上改变。

该数据以来的周期性现象,以及一个向上发展的趋势,是合适的非线性拟合技术,但是工作量将会添加相比拟合线性模型。然而,使用正确的工具和知识的数据,您可以快速确定系数的非线性模型。

本文描述了如何适应可用非线性模型有限公司2使用新发布的数据曲线拟合工具箱。

建模的数据

在这个例子中,我们使用有限公司2收集的数据,在ppm,由美国国家海洋和大气管理局(NOAA)从1979年到1996年在夏威夷的莫纳罗亚山天文台2。捕获的数据下图所示。

modeling_fig1_w.gif
曲线拟合工具中发现的曲线拟合工具箱提供了一个直接的视觉显示的数据。通过使用数据提示功能,我们可以询问具体点。

周期性现象在自然界中可以通过一个正弦信号建模及其高次谐波,也就是说,

modeling_logo_eq1.gif
一个kT / KΦk分别是振幅、周期和相位,k的谐波。正弦条件捕获季节性趋势,但持续的上升趋势需要一个指数项。这些数据,提出模型确定通过检查的记录总世界化石燃料产量从1860年代中期
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平方值(一个衡量拟合优度)从一个不远,暗示全球趋势,但不是当地的趋势。在这种情况下,正弦信号的振幅太小。

改进模型

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

modeling_logo_eq4.gif
这个模型拟合系数d = 349.6α= 0.02185。使用这些值d开始值的正弦模型,我们得到一个新的适合匹配的周期性。
modeling_fig4_w.gif

这个最新型号的幅度大约是3。我们可以验证的阴谋,这是一个合理的值为1(左)。结果面板(下图)也显示了每个参数的置信区间在括号后的系数值。较小的置信区间参数值意味着更多的信心。所有参数的置信区间很小,除了B和d,这可能表明B和d是高度相关的。现在我们可以画出残差仍表现出周期性正弦信号模式,所以另一个正弦信号词是必要的。

现在我们可以画出残差仍表现出周期性正弦信号模式,所以另一个正弦信号词是必要的。


modeling_fig5_w.gif
单一正弦信号模型的残差显示另一个正弦曲线(左)。有两个谐波,适合很好(右)。我们可以用更轻松地创建方程方面,检查每个产生的残差满足直到我们感到满意。

扩展分析

当我们发现了一个令人满意的模型时,我们可以发送适合MATLAB工作区进行进一步分析。我们有一些其他天文台网站检查(例如,萨摩亚、南极和巴罗),所以可以生成一个m文件函数,它可以重建我们的适合。然后我们可以在一个循环中调用该函数自动配合其他数据集,看看在天文台网站生成的模型进行比较。

这个例子中涉及的一些特性,将帮助你适应模型数据。模型是否线性或非线性曲线拟合工具箱允许您快速适应标准模型库或您创建自定义模型。

modeling_fig6_w.gif
m文件中创建的代码复制适合和情节的GUI可以从GUI生成文件菜单。m文件还揭示了工具箱的命令,可以使用在命令行中得到相同的结果。

2002 - 91034 v00出版

引用

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

  2. 卡恩,D。硅藻土,C。纳什,S。数值方法和软件,普伦蒂斯霍尔(1989)

s manbetx 845产品使用

查看相关文章的能力