使用levmar的L-M算法拟合曲线
2015-11-30 18:40
453 查看
本文所要你拟合的曲线公式:
未知参数:Pi
所求参数数目m=4
已知的点数目n=7
初始的参数设定 p[m] = { 1,1,1,1 }
X[i]设定为y: x
= {0,1,2,2,3,4,2};
Data是已知点的x:data[7] = {0,1,2,3,3,4,10};
其他保留为默认。
目标函数:
雅可比函数:
调用:
结果:
未知参数:Pi
所求参数数目m=4
已知的点数目n=7
初始的参数设定 p[m] = { 1,1,1,1 }
X[i]设定为y: x
= {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);
结果:
相关文章推荐
- 大神的指导
- HTML5 基础知识,第 3 部分: HTML5 API 的强大功能
- Android兼容性问题 -- WebP格式图片解码失败
- spark在集群上运行
- Android 百度地图定位指针(系列2)
- 黑马程序员——C语言函数的返回值笔记和练习题
- android studio 编译运行报错
- 借助Babel 6平台使用ES6新特性
- volatile 知识点详解
- Tracking
- 大气科学词典
- wustoj 1056 C语言程序设计教程(第三版)课后习题10.5
- Codeforces 602B Approximating a Constant Range(想法题)
- Java学习笔记(多线程_2)
- vim中常用的命令
- jquery checkbox 点击事件
- Android 图片缓存处理
- WEB打印控件
- wustoj 习题总结
- iOS Reveal 分析他人的UI界面