对给定的试验数据点(xi,yi)(i=1,2,……,n),可以构造m次多项式

数据拟合的最简单的做法就是使误差p(xi)-yi的平方和最小

当前任务就是求一个P(x)使得

从几何意义上讲就是寻求给与定点(xi,yi)距离的平方和最小的曲线y=p(x),函数p(x)称为拟合函数或者是最小二乘解,求拟合函数p(x)的方法称为曲线拟合的最小二乘法

由极值条件得到

即求得

通过简单运算可以得出系数是下面线性方程组的解。

在matlab中编程实现多项式曲线拟合函数为:LeastSquareDetail

功能:求已知数据点的多项式曲线拟合插值法多项式

调用格式:A=LeastSquareDetail(x,y,m,x0)

其中:x为已知数据点的x坐标向量

y为已知数据点的y坐标向量

m为拟合多项式的次数

A为拟合多项式的系数向量

多项式曲线拟合的matlab实现

LeastSquareDetail<span style="font-family: Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255);">.m</span>
<pre name="code" class="plain">function f = LeastSquareDetail(x,f,m,x0)
%x代表输入的x的值
%y代表输入的x对应的y值
%m代表拟合的次数
%x0为待求的点的x坐标format long g
a=length(x);
b=length(f);
if(a~=b)disp('xy的维数不相等!!');
end%求得矩阵A
A=zeros(m+1);
for i=1:m+1for j=i:m+1if(i==1&&j==1)A(i,j)=a;elsefor k=1:aA(i,j)=A(i,j)+x(k)^(i+j-2);endendA(j,i)=A(i,j);end
end%求得等式右边的矩阵
B=zeros(1,m+1);
for i=1:m+1if(i==1)for k=1:aB(i)=B(i)+f(k);endelsefor k=1:aB(i)=B(i)+x(k)^(i-1)*f(k);endend
end
B=B';
%求解系数
X=A\B;%求得到的方程
syms t;
c=length(X);
f=0;
for i=1:cf=f+X(i)*t^(i-1);f=collect(f);f=vpa(f,6);
endif nargin==4f=subs(f,'t',x0);
end
LeastSquareDetail<span style="font-family: Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255);">Insert.m</span>
<pre name="code" class="plain">% x=[36.9 46.7 63.7 77.8 84.0 87.5];
% y=[181 197 235 270 283 292];
% LeastSquareDetail(x,y,1)x=0:2*pi;
y=sin(x);
xx=0:0.1:2*pi;
yy=LeastSquareDetail(x,y,3,xx);
plot(x,y,'b',xx,yy,'r')

多项式曲线拟合最小二乘法相关推荐

  1. 最小二乘法多项式曲线拟合原理与实现--转

    原文地址:http://blog.csdn.net/jairuschan/article/details/7517773/ 概念 最小二乘法多项式曲线拟合,根据给定的m个点,并不要求这条曲线精确地经过 ...

  2. 最小二乘法多项式曲线拟合数学原理及其C++实现

    目录 0 前言 1 最小二乘法概述 2 最小二乘法求解多项式曲线系数向量的数学推导 2.1 代数法 2.2 矩阵法 3 代码实现 4 总结 参考 0 前言 自动驾驶开发中经常涉及到多项式曲线拟合,本文 ...

  3. 多项式曲线拟合之最小二乘法推导

    1.多项式曲线拟合之最小二乘法 1.1 问题来源 1801年,意大利天文学家朱赛普·皮亚齐发现了第一颗小行星谷神星.经过40天的跟踪观测后,由于谷神星运行至太阳背后,使得皮亚齐失去了谷神星的位置.随后 ...

  4. 最小二乘法多项式曲线拟合原理与实现

    概念 最小二乘法多项式曲线拟合,根据给定的m个点,并不要求这条曲线精确地经过这些点,而是曲线y=f(x)的近似曲线y= φ(x). 原理 [原理部分由个人根据互联网上的资料进行总结,希望对大家能有用] ...

  5. 最小二乘法多项式曲线拟合原理与实现(错误地方已经修改底层补充自己写的java实现)

    目录(?) [-] 概念 原理 运行前提 代码 运行效果 概念 最小二乘法多项式曲线拟合,根据给定的m个点,并不要求这条曲线精确地经过这些点,而是曲线y=f(x)的近似曲线y= φ(x). 原理 [原 ...

  6. c++ 三次多项式拟合_最小二乘法多项式曲线拟合数学原理及其C++实现

    本文使用 Zhihu On VSCode 创作并发布 0 前言 自动驾驶开发中经常涉及到多项式曲线拟合,本文详细描述了使用最小二乘法进行多项式曲线拟合的数学原理,通过样本集构造范德蒙德矩阵,将一元 N ...

  7. 基于自动分段最小二乘法对股票的多项式曲线拟合

    基于自动分段最小二乘法对股票的多项式曲线拟合 摘 要 针对传统的分段最小二乘法确定分段步长时经验成分较多的不足,提出一种通过比较拟合优度,自动确定相对最优的步长.通过实际数据的验证,验证了此方法的拟合 ...

  8. 最小二乘法多项式曲线拟合及其python实现

    最小二乘法多项式曲线拟合及其python实现 多项式曲线拟合问题描述 最小二乘法 针对overfitting,加入正则项 python实现 运行结果 多项式曲线拟合问题描述 问题描述:给定一些数据点, ...

  9. 趋势预测方法(一) 多项式拟合(最小二乘法)_函数拟合

    多项式拟合(最小二乘法) a基本原理: b拟合函数原理: c方法优缺点: 优点: 思想简单,实现容易.建模迅速,对于小数据量.简单的关系很有效. 解决回归问题,拥有很好的解释性. 是很多非线性模型的基 ...

最新文章

  1. JDK1.8源码分析:线程安全的CopyOnWriteArrayList与CopyOnWriteArraySet
  2. android application常见错误
  3. 添加MSDN搜索的正确方法
  4. python删除为空值的列_python如何删除列为空的行
  5. SAP UI5 jQuery.sap.formatMessage
  6. Java随机生成长宽的矩形_java – 将正方形或矩形分解为大量随机大小的正方形或矩形...
  7. android studio 制作表格_红爆网络的旅游排行榜视频制作其实很简单,赶快动手尝试...
  8. tcp报文格式_34.TCP取样器
  9. android文件体系,Android文件体系-基础
  10. 逾期怎么处理_信用卡3万逾期三个月,催收说已经起诉生效,立案处理了,该怎么办?...
  11. 工作开不开心,都在钱上
  12. [.NET] 在Windows系统中搭建基于.NET的iPhone应用程序虚机开发环境
  13. 编译redis-5.0.9报错zmalloc.h:50:31: 致命错误:jemalloc/jemalloc.h:没有那个文件或目录问题解决
  14. CentOS 7系统上制作Clonezilla(再生龙)启动U盘并克隆双系统
  15. 职称专业技术工作小结 计算机,申报计算机专业中级职称任职来专业技术工作总结(范文).doc...
  16. 立创eda学习笔记二十九:原理图转PCB
  17. 千万不要死于无知——心理状态
  18. 【Algorithm】Prim
  19. 熬之滴水穿石:JAVA的世界(2)
  20. 完整的DJ表演工具包djay Pro 2 Mac,DJ打碟神器

热门文章

  1. mysql 密码1045_mysql登录1045错误时 修改登录密码
  2. java thread_Java(多线程Thread)
  3. java servlet接口开发_servlet三种实现方式之一实现servlet接口
  4. 写代码?程序猿?你不能不懂的八大排序算法的Python实现
  5. CF--思维练习--CodeForces - 220C Little Elephant and Shifts (STL模拟)
  6. 图论--Dijkstra算法总结
  7. linux 头文件汇总
  8. 一位强人写的电子工程师起码技术要求
  9. uCOS-II任务建立示例
  10. 关于keil环境的 三个红点(备忘)