直線(xiàn)擬合算法

函數(shù)說(shuō)明：

通過(guò)一些點(diǎn)擬合出一條直線(xiàn)。

參數(shù)：

pt_input
指向傳入的點(diǎn)的指針。

ptNumbers
傳入的點(diǎn)數(shù)量。

k
指向擬合直線(xiàn)參數(shù)k的指針。

b
指向擬合直線(xiàn)參數(shù)b的指針。

返回值：

返回一個(gè)代數(shù)值判斷擬合是否成功。如果成功，則返回 0。如果返回-1，表示點(diǎn)的數(shù)量<2。如果返回-2，表示x坐標(biāo)無(wú)變化。

代碼：

typedef struct tagPOINT_2 {double x;double y; }POINT_2, *PPOINT_2;int LineInfo(POINT_2 *pt_input, int ptNumbers, double * k, double * b) {int nRet = 0;double X_aver = 0;double Y_aver = 0;double A = 0;double B = 0;double _2_XY_av = 0;double _2_XX_av = 0;do{if (ptNumbers < 2){nRet = -1;break;}for (int loop = 0; loop < ptNumbers; loop++){X_aver += pt_input[loop].x;Y_aver += pt_input[loop].y;_2_XY_av += pt_input[loop].x * pt_input[loop].y;_2_XX_av += pt_input[loop].x * pt_input[loop].x;}X_aver /= ptNumbers;Y_aver /= ptNumbers;_2_XY_av /= ptNumbers;_2_XX_av /= ptNumbers;for (int loop = 0; loop < ptNumbers; loop++){A += fabs((pt_input[loop].x - X_aver) * (pt_input[loop].y - Y_aver));B += fabs((pt_input[loop].x - X_aver) * (pt_input[loop].x - X_aver));}if (B < 0.000001){nRet = -2;break;}//方式一//*k = A / B;//*b = Y_aver - *k * X_aver;//方式二*k = (_2_XY_av - X_aver * Y_aver) / (_2_XX_av - X_aver * X_aver);*b = Y_aver - (*k) * X_aver;} while (0);return nRet; }

測(cè)試代碼：

#include <stdio.h> #include <math.h> int main() {POINT_2 pt[10] = {{0.1 , 0.3},{1.2 , 2.0},{2.2 , 4.3},{3.1 , 5.8},{4.0 , 8.5},{5.3 , 10.8},{5.9 , 11.1},{6.8 , 13.9},{8.1 , 16.5},{9.0 , 17.9},};int ret;double k = 0, b = 0;ret = LineInfo(pt, 10, &k, &b);switch (ret){case 0: printf("擬合成功，直線(xiàn)為y = %lf * x + %lf\n", k, b); break;case -1: printf("擬合失敗，傳入的點(diǎn)少于二點(diǎn)\n"); break;case -2: printf("擬合失敗，傳入的點(diǎn)x坐標(biāo)重合\n"); break;default:printf("error\n"); break;}return 0; }

測(cè)試結(jié)果：

總結(jié)

以上是生活随笔為你收集整理的C语言直线拟合函数的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。