本文所要你拟合的曲线公式:

未知参数:Pi

所求参数数目m=4

已知的点数目n=7

初始的参数设定 p[m] = { 1,1,1,1 }

X[i]设定为y: x[n] = {0,1,2,2,3,4,2};

Data是已知点的x:data[7] = {0,1,2,3,3,4,10};

其他保留为默认。

目标函数:

void ros(double* p, double* x, int m, int n, void* data){LVTYPE* data1 = (LVTYPE*)data;register int i;for (i = 0; i < n; i++){x[i] =(p[0]+p[1] * data1[i]+p[2]*data1[i]*data1[i]+p[3]*data1[i]*data1[i]*data1[i]);}
}

雅可比函数:

void jacros(double* p, double *jac, int m, int n, void* data){LVTYPE* data1 = (LVTYPE*)data;register int i, j;for (i = 0, j = 0; i < n; i++){jac[j++] = 1;jac[j++] = data1[i];jac[j++] = data1[i]*2;jac[j++] = data1[i] * data1[i] * data1[i] * 3;}
}

调用:

int ret = dlevmar_der(ros, jacros, p, x, m, n, 10000, NULL, info, NULL, NULL, (void*)data);

结果:

转载自:http://www.cnblogs.com/jiafenggang/p/5008044.html

使用Levmar的L-M算法拟合曲线相关推荐

  1. Levmar:Levenberg-Marquardt非线性最小二乘算法

    Levmar:Levenberg-Marquardt非线性最小二乘算法 eryar@163.com Abstract. Levmar is GPL native ANSI C implementati ...

  2. 机器学习算法 拟合曲线_制定学习曲线以检测机器学习算法中的错误

    机器学习算法 拟合曲线 机器学习 (Machine Learning) The learning curve is very useful to determine how to improve th ...

  3. 小L的算法课堂——博弈论

    嗨!这里是小L的算法课堂!我是小L 今天我们学习一下博弈论 我们先来讲一个故事: 小L和小Z喜欢玩石头剪刀布. 小L总是打不过小Z赢,于是,小L问小Z:"你是怎么赢我的?" 小Z告 ...

  4. l bfgs算法java代码_优化算法——拟牛顿法之L-BFGS算法

    一.BFGS算法 BFGS算法的校正公式: 利用Sherman-Morrison公式可对上式进行变换,得到 令 ,则得到: 二.BGFS算法存在的问题 在BFGS算法中.每次都要存储近似Hesse矩阵 ...

  5. l bfgs算法java代码_L-BFGS算法介绍

    本文由作者林洋港授权网易云社区发布. 一. L-BFGS是什么 L-BFGS是解无约束非线性规划问题最常用的方法,具有收敛速度快.内存开销少等优点,在机器学习各类算法中常有它的身影.简单的说,L-BF ...

  6. l bfgs算法java代码_理解L-BFGS算法

    理解L-BFGS算法 Mar 30, 2015   #数值优化  #无约束最优化 L-BFGS(Limited-Memory BFGS)是BFGS算法在受限内存时的一种近似算法,而BFGS是数学优化中 ...

  7. l bfgs算法java代码_数值优化:理解L-BFGS算法

    译自<Numerical Optimization: Understanding L-BFGS>,本来只想作为学习CRF的补充材料,读完后发现收获很多,把许多以前零散的知识点都串起来了.对 ...

  8. 数据结构与算法(一)——排序

    虽然之前学过数据结构,但是已时隔四年,大概四月份复习了一遍,但是很多概念也是一知半解,所以重新整理知识点和运行代码的方式来巩固知识. 引言 排序:是计算机程序设计中的一种重要操作,功能是将一个数据元素 ...

  9. Streaming的算法Reservoir Sampling

    转载自   这是一个惊艳了我的算法题 Reservoir Sampling( Reservoir sampling )  这是我在今年求职过程中面试的时候被问到的,因为之前很少接触Streaming的 ...

最新文章

  1. pandas之数据结构
  2. 怎么把分开的pdf放在一起_糖和盐混在一起了要怎么分开?| 趣问万物
  3. 面向对象设计启发规则
  4. java查询SQLServer遇到问题:对象名无效。
  5. C#LeetCode刷题之#747-至少是其他数字两倍的最大数( Largest Number At Least Twice of Others)
  6. hhkb mac设置_我最爱的HHKB
  7. 实验3 OpenGL几何变换
  8. 办公、学习不得不收藏的十个网站,有了它效率至少提升十倍
  9. 一起话单业务量下降问题的排查过程
  10. HTML5+js+css3开心消消乐手机pc端通用源码|H5小游戏
  11. AtCoder Grand Contest 007 题解
  12. android访问服务器405,Http 405 错误
  13. 手机联系人误删怎么恢复?如何找回被删除的电话联系人
  14. wex5 php开发,WeX5开发移动APP(SQLite本地数据优化)
  15. 浙江新青年计算机学校,新青年
  16. 面试官这么表现就有戏?不一定哦!
  17. 网速是什么意思?带宽是什么意思?1M网速的下载速度应是多少?[转.baidu]
  18. C语言 平面向量加法
  19. 三种安装httpd的方法
  20. XDU-OJ-C语言-Test 2

热门文章

  1. Oracle 12C DataGuard部署以及维护
  2. 如何使用 IntelliJ IDEA 2017 配置PHP开发环境 及项目搭建
  3. IntelliJ IDEA Community社区版集成Tomcat or Jetty教程
  4. 在 ASP.NET Core 5.0 中访问 HttpContext
  5. axios的数据拦截(拦截器)
  6. C#算法设计排序篇之07-希尔排序(附带动画演示程序)
  7. mybatis 配置详解
  8. ddos常见攻击报文
  9. 用于Elasticsearch数据可视化和分析的强大工具
  10. linux内核源代码_如何从源代码构建和安装最新的Linux内核