实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。曲线拟合（curve fitting）是指选择适当的曲线类型来拟合观测数据，并用拟合的曲线方程分析两变量间的关系。

简介

用连续曲线近似地刻画或比拟平面上离散点组所表示的坐标之间的函数关系的一种数据处理方法。用解析表达式逼近离散数据的一种方法。在科学实验或社会活动中，通过实验或观测得到量x与y的一组数据对(xi，yi)(i=1，2，…m)，其中各xi是彼此不同的。

人们希望用一类与数据的背景材料规律相适应的解析表达式，y=f(x，c)来反映量x与y之间的依赖关系，即在一定意义下"最佳"地逼近或拟合已知数据。f(x，c)常称作拟合模型，式中c=(c1，c2，…cn)是一些待定参数。当c在f中线性出现时，称为线性模型，否则称为非线性模型。

有许多衡量拟合优度的标准，最常用的一种做法是选择参数c使得拟合模型与实际观测值在各点的残差(或离差)ek=yk-f(xk，c)的加权平方和达到最小，此时所求曲线称作在加权最小二乘意义下对数据的拟合曲线。有许多求解拟合曲线的成功方法，对于线性模型一般通 过建立和求解方程组来确定参数，从而求得拟合曲线。

至于非线性模型，则要借助求解非线性方程组或用最优化方法求得所需参数才能得到拟合曲线，有时称之为非线性最小二乘拟合。曲线拟合:贝塞尔曲线与路径转化时的误差。值越大，误差越大;值越小，越精确。

常用函数

指数函数

指数函数(exponential function)的标准式形式为：Y=aebX(12.29)对式（12.29）两边取对数，得lnY=lna+bX

b>0时，Y随X增大而增大；b<0时，Y随X增大而减少。当以lnY和X绘制的散点图呈直线趋势时，可考虑采用指数函数来描述Y与X间的非线性关系，lna和b分别为截距和斜率。

更一般的指数函数：Y=aebX+k，式中k为一常量，往往未知,应用时可试用不同的值。

对数函数

对数函数（lograrithmic function）的标准式形式为：Y=a+blnX，(X>0)，b>0时，Y随X增大而增大，先快后慢；b<0时，Y随X增大而减少，先快后慢，见图12.4(c)、(d)。当以Y和lnX绘制的散点图呈直线趋势时，可考虑采用对数函数描述Y与X之间的非线性关系，式中的b和a分别为斜率和截距。

更一般的对数函数：Y=a+bln(X+k)，式中k为一常量，往往未知。(a)lnY=lna+bX(b)lnY=lna-bX(c)Y=a+blnX(d)Y=a-blnX

幂函数

幂函数（power function）的标准式形式为：Y=aXb(a>0,X>0)，式中b>0时，Y随X增大而增大；b<0时，Y随X增大而减少。对两边取对数，得lnY=lna+blnX所以，当以lnY和lnX绘制的散点图呈直线趋势时，可考虑采用幂函数来描述Y和X间的非线性关系，lna和b分别是截距和斜率。

更一般的幂函数：Y=aXb+k，式中k为一常量，往往未知。

步骤

1、绘制散点图，选择合适的曲线类型，一般根据资料性质结合专业知识便可确定资料的曲线类型，不能确定时，可在方格坐标纸上绘制散点图，根据散点的分布，选择接近的、合适的曲线类型。

2、进行变量变换：Y’=f(Y),X’=g(X)，使变换后的两个变量呈直线关系。

3、按最小二乘法原理求线性方程和方差分析

4、将直线化方程转换为关于原变量X、Y的函数表达式

应用

重力波参数气候特征是确定大气模式中重力波参数化方案的重要条件之一，高垂直分辨率探空资料扰动场是获取重力波参数气候特征的基础数据;目前,获取扰动场的方法较多，但基于不同方法计算的扰动场对重力波参数气候特征影响的研究较少。

基于2014—2017年山西太原气象台高垂直分辨率探空资料，利用2—4阶曲线拟合方法获取下平流层(17—24km高度)温度扰动场、纬向风扰动场和经向风扰动场,经统计发现2阶与3阶曲线拟合方法的扰动场相似程度较高；在此选取相似度较高的2阶、3阶曲线拟合方法的扰动场分别计算大气重力波参数，并对大气重力波参数间的气候差异特征进行研究。

结果表明：(1)不同阶曲线拟合方法扰动场的变化振幅及随高度变化趋势存在差异，且扰动场间的相关较弱;(2)2阶、3阶曲线拟合方法扰动场得到的重力波参数大小、年内变化趋势及在不同区间范围内占有率均存在差异，且相关较弱;(3)1—12月，相对3阶曲线拟合方法的扰动场，基于2阶曲线拟合方法的扰动场得到的重力波群速、水平波长、垂直波长、周期、固有相速均较大，而重力波能量上传百分比在某些月份较大。

因此,不同阶曲线拟合方法扰动场间存在差异,会导致计算得到的大气重力波参数气候特征存在差异,最终对研制大气模式中的大气重力波参数化方案产生影响。