C语言版的线性回归函数总结.docVIP

C语言版的线性回归分析函数 分类： C/C++ 2007-08-03 23:39 13840人阅读 评论(31) 收藏 举报 语言c数学计算delphisystem ??????? DOS下的程序及代码，看来目前也没什么用了，想打个包刻在光碟上，却发现有些代码现在可能还能起作用，其中就有计算一元回归和多元回归的代码，一看代码文件时间，居然是1993年的，于是稍作整理，存放在这，分析虽不十分完整，但一般应用是没问题的，最起码，可提供给那些刚学C的学生们参考。 ? //?求线性回归方程：Y?=?a?+?bx//?dada[rows*2]数组：X,?Y；rows：数据行数；a,?b：返回回归系数//?SquarePoor[4]：返回方差分析指标:?回归平方和，剩余平方和，回归平方差，剩余平方差//?返回值：0求解成功，-1错误int?LinearRegression(double?*data,?int?rows,?double?*a,?double?*b,?double?*SquarePoor){????int?m;????double?*p,?Lxx?=?0.0,?Lxy?=?0.0,?xa?=?0.0,?ya?=?0.0;????if?(data?==?0?||?a?==?0?||?b?==?0?||?rows??1)????????return?-1;????for?(p?=?data,?m?=?0;?m??rows;?m?++)????{????????xa?+=?*p?++;????????ya?+=?*p?++;????}????xa?/=?rows;?????????????????????????????????????//?X平均值????ya?/=?rows;?????????????????????????????????????//?Y平均值????for?(p?=?data,?m?=?0;?m??rows;?m?++,?p?+=?2)????{????????Lxx?+=?((*p?-?xa)?*?(*p?-?xa));?????????????//?Lxx?=?Sum((X?-?Xa)平方)????????Lxy?+=?((*p?-?xa)?*?(*(p?+?1)?-?ya));???????//?Lxy?=?Sum((X?-?Xa)(Y?-?Ya))????}????*b?=?Lxy?/?Lxx;?????????????????????????????????//?b?=?Lxy?/?Lxx????*a?=?ya?-?*b?*?xa;??????????????????????????????//?a?=?Ya?-?b*Xa????if?(SquarePoor?==?0)????????return?0;????//?方差分析????SquarePoor[0]?=?SquarePoor[1]?=?0.0;????for?(p?=?data,?m?=?0;?m??rows;?m?++,?p?++)????{????????Lxy?=?*a?+?*b?*?*p?++;????????SquarePoor[0]?+=?((Lxy?-?ya)?*?(Lxy?-?ya));?//?U(回归平方和)????????SquarePoor[1]?+=?((*p?-?Lxy)?*?(*p?-?Lxy));?//?Q(剩余平方和)????}????SquarePoor[2]?=?SquarePoor[0];??????????????????//?回归方差????SquarePoor[3]?=?SquarePoor[1]?/?(rows?-?2);?????//?剩余方差????return?0;} ? 为了理解代码，把几个与代码有关的公式写在下面（回归理论和公式推导就免了，网上搜索到处是，下面的公式图片也是网上搜的，有些公式图形网上没找到或者不合适，可参见后面多元回归中的公式）： 1、回归方程式： 2、回归系数：??? ???其中： ??????? ? ??????? ?? 3、回归平方和： 4、剩余平方和： 实例计算： double?data1[12][2]?=?{//????X??????Y????{187.1,?25.4},????{179.5,?22.8},????{157.0,?20.6},????{197.0,?21.8},????{239.4,?32.4},????{217.8,?24.4},????{227.1,?29.3},????{233.4,?27.9},????{242.0,?27.8},???