使用Levmar的L-M算法拟合曲线
本文所要你拟合的曲线公式:
未知参数: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算法拟合曲线相关推荐
- Levmar:Levenberg-Marquardt非线性最小二乘算法
Levmar:Levenberg-Marquardt非线性最小二乘算法 eryar@163.com Abstract. Levmar is GPL native ANSI C implementati ...
- 机器学习算法 拟合曲线_制定学习曲线以检测机器学习算法中的错误
机器学习算法 拟合曲线 机器学习 (Machine Learning) The learning curve is very useful to determine how to improve th ...
- 小L的算法课堂——博弈论
嗨!这里是小L的算法课堂!我是小L 今天我们学习一下博弈论 我们先来讲一个故事: 小L和小Z喜欢玩石头剪刀布. 小L总是打不过小Z赢,于是,小L问小Z:"你是怎么赢我的?" 小Z告 ...
- l bfgs算法java代码_优化算法——拟牛顿法之L-BFGS算法
一.BFGS算法 BFGS算法的校正公式: 利用Sherman-Morrison公式可对上式进行变换,得到 令 ,则得到: 二.BGFS算法存在的问题 在BFGS算法中.每次都要存储近似Hesse矩阵 ...
- l bfgs算法java代码_L-BFGS算法介绍
本文由作者林洋港授权网易云社区发布. 一. L-BFGS是什么 L-BFGS是解无约束非线性规划问题最常用的方法,具有收敛速度快.内存开销少等优点,在机器学习各类算法中常有它的身影.简单的说,L-BF ...
- l bfgs算法java代码_理解L-BFGS算法
理解L-BFGS算法 Mar 30, 2015 #数值优化 #无约束最优化 L-BFGS(Limited-Memory BFGS)是BFGS算法在受限内存时的一种近似算法,而BFGS是数学优化中 ...
- l bfgs算法java代码_数值优化:理解L-BFGS算法
译自<Numerical Optimization: Understanding L-BFGS>,本来只想作为学习CRF的补充材料,读完后发现收获很多,把许多以前零散的知识点都串起来了.对 ...
- 数据结构与算法(一)——排序
虽然之前学过数据结构,但是已时隔四年,大概四月份复习了一遍,但是很多概念也是一知半解,所以重新整理知识点和运行代码的方式来巩固知识. 引言 排序:是计算机程序设计中的一种重要操作,功能是将一个数据元素 ...
- Streaming的算法Reservoir Sampling
转载自 这是一个惊艳了我的算法题 Reservoir Sampling( Reservoir sampling ) 这是我在今年求职过程中面试的时候被问到的,因为之前很少接触Streaming的 ...
最新文章
- pandas之数据结构
- 怎么把分开的pdf放在一起_糖和盐混在一起了要怎么分开?| 趣问万物
- 面向对象设计启发规则
- java查询SQLServer遇到问题:对象名无效。
- C#LeetCode刷题之#747-至少是其他数字两倍的最大数( Largest Number At Least Twice of Others)
- hhkb mac设置_我最爱的HHKB
- 实验3 OpenGL几何变换
- 办公、学习不得不收藏的十个网站,有了它效率至少提升十倍
- 一起话单业务量下降问题的排查过程
- HTML5+js+css3开心消消乐手机pc端通用源码|H5小游戏
- AtCoder Grand Contest 007 题解
- android访问服务器405,Http 405 错误
- 手机联系人误删怎么恢复?如何找回被删除的电话联系人
- wex5 php开发,WeX5开发移动APP(SQLite本地数据优化)
- 浙江新青年计算机学校,新青年
- 面试官这么表现就有戏?不一定哦!
- 网速是什么意思?带宽是什么意思?1M网速的下载速度应是多少?[转.baidu]
- C语言 平面向量加法
- 三种安装httpd的方法
- XDU-OJ-C语言-Test 2
热门文章
- Oracle 12C DataGuard部署以及维护
- 如何使用 IntelliJ IDEA 2017 配置PHP开发环境 及项目搭建
- IntelliJ IDEA Community社区版集成Tomcat or Jetty教程
- 在 ASP.NET Core 5.0 中访问 HttpContext
- axios的数据拦截(拦截器)
- C#算法设计排序篇之07-希尔排序(附带动画演示程序)
- mybatis 配置详解
- ddos常见攻击报文
- 用于Elasticsearch数据可视化和分析的强大工具
- linux内核源代码_如何从源代码构建和安装最新的Linux内核