最小二乘法拟合曲线
来源:互联网 发布:返利系统源码 编辑:程序博客网 时间:2024/04/28 06:46
最小二乘法
最小二乘法(又称最小平方法)是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。
利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数据与实际数据之间误差的平方和为最小。
最小二乘法还可用于曲线拟合。
其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达。
人们对由某一变量t 或多个变量t1……tn 构成的相关变量y感兴趣。如弹簧的形变与所用的力相关,一个企业的盈利与其营业额,投资收益和原始资本有关。为了得到这些变量同y之间的关系,便用不相关变量去构建y,使用如下函数模型
- ,
q个相关变量或p个附加的相关变量去拟和。
通常人们将一个可能的、对不相关变量t的构成都无困难的函数类型充作函数模型(如抛物线函数或指数函数)。参数x是为了使所选择的函数模型同观测值y相匹配。(如在测量弹簧形变时,必须将所用的力与弹簧的膨胀系数联系起来)。其目标是合适地选择参数,使函数模型最好的拟合观测值。一般情况下,观测值远多于所选择的参数。
其次的问题是怎样判断不同拟合的质量。高斯和勒让德的方法是,假设测量误差的平均值为0。令每一个测量误差对应一个变量并与其它测量误差不相关(随机无关)。人们假设,在测量误差中绝对不含系统误差,它们应该是纯偶然误差,围绕真值波动。除此之外,测量误差符合正态分布,这保证了偏差值在最后的结果y上忽略不计。
确定拟合的标准应该被重视,并小心选择,较大误差的测量值应被赋予较小的权。并建立如下规则:被选择的参数,应该使算出的函数曲线与观测值之差的平方和最小。用函数表示为:
用欧几里得度量表达为:
最小化问题的精度,依赖于所选择的函数模型。
典型的一类函数模型是线性函数模型。最简单的线性式是y = x0 + x1t,写成行列式,为
直接给出该式的参数解:
- 和
其中,为t值的算术平均值。也可解得如下形式:
[编辑]简单线性模型 y = x0 + x1t 的例子
随机选定10艘战舰,并分析它们的长度与宽度,寻找它们长度与宽度之间的关系。由下面的描点图可以直观地看出,一艘战舰的长度(t)与宽度(y)基本呈线性关系。
以下图表列出了各战舰的数据,随后步骤是采用最小二乘法确定两变量间的线性关系。
仿照上面给出的例子
并得到相应的.
然后确定x1
可以看出,战舰的长度每变化1m,相对应的宽度便要变化16cm。并由下式得到常数项x0:
在这里随机理论不加阐述。可以看出点的拟合非常好,长度和宽度的相关性大约为92%。
[编辑]一般线性情况
若含有更多不相关模型变量t1,...,tq,可如组成线性函数的形式
即线性方程组
通常人们将tij记作数据矩阵 A,参数xj记做参数矢量x,观测值yi记作b,则线性方程组又可写成:
- 即 Ax = b
上述方程运用最小二乘法导出为线性平差计算的形式为:
- 。
[编辑]最小化问题的解
① 一次函数使用polyfit(x,y,1) ②多项式函数使用 polyfit(x,y,n),n为次数 拟合曲线 x=[0.5,1.0,1.5,2.0,2.5,3.0], y=[1.75,2.45,3.81,4.80,7.00,8.60]。 解:MATLAB程序如下: x=[0.5,1.0,1.5,2.0,2.5,3.0]; y=[1.75,2.45,3.81,4.80,7.00,8.60]; p=polyfit(x,y,2) x1=0.5:0.05:3.0; y1=polyval(p,x1); plot(x,y,'*r',x1,y1,'-b') 计算结果为: p =0.5614 0.8287 1.1560 即所得多项式为y=0.5614x^2+0.08287x+1.15560 ③非线性函数使用 lsqcurvefit(fun,x0,x,y)- 最小二乘法拟合曲线
- 最小二乘法拟合二元多次曲线
- 用最小二乘法构造拟合曲线
- 最小二乘法求多项式拟合曲线
- 最小二乘法拟合多项式曲线原理
- Eigen最小二乘法拟合三次曲线
- 最小二乘法拟合曲线:二次函数
- 利用最小二乘法拟合任意次函数曲线
- 最小二乘法拟合曲线:4次函数
- 利用最小二乘法拟合任意次函数曲线(C#)
- 最小二乘法拟合散点组成的椭圆曲线
- 利用最小二乘法拟合任意次函数曲线(C#)
- Python 最小二乘法多项式拟合曲线numpy.polyfit(),numpy.poly1d(),pylab
- 移动最小二乘法(MLS)曲线曲面拟合C++代码实现
- 拟合直线 二次函数曲线 最小二乘法 javascript(p5.js)
- 最小二乘法拟合
- 最小二乘法拟合
- 拟合曲线
- 智能卡的电源电流
- 启用Xdebug使用WinCacheGrind分析脚本执行时间
- 智能卡的激活和去激活序列
- 智能卡的微控制器
- 第十一周实验报告1--学会使用循环控制语句解决实际问题,并强化注重效率的观念
- 最小二乘法拟合曲线
- 智能卡的电气特性
- 智能卡的电源电压
- 智能卡处理器外部时钟
- HashMap漏洞
- 圈存&圈提
- SVN中修改hook实现commit文件之后自动发布的方法
- php数组
- 在VS2005中应用多个其它编译器