最小二乘法直线拟合
3页1、设经验方程是y=F(x),方程中含有一些待定系数an,给出真实值(xi,yi)|i=1,2,n,将这些x,y值代入方程然后作差,可以描述误差:yi-F(xi),为了考虑整体的误差,可以取平方和,之所以要平方是考虑到误差可正可负直接相加可以相互抵消,所以记误差为:e=(yi-F(Xi)2如果经验方程是线性的,形如y=ax+b,就是线性回归。按上面的分析,误差函数为:e=(y-aXj-b)2各偏导为:de/da=2(yy-axi-b)(-xi)=0de/db=-2(y-axb)=O于是得到关于a,b的线性方程组:(2)a+(x)b=yyxi(x)a+nb=y设a=x2,B=x,C=yyxi,D=y,则方程化为:Aa+Bb=CBa+nb=D解出a,b得:a=(Cn-BD)/(An-BB)b=(AD-CB)/(An-BB)这就是我们要进行的算法。C+代码如下:1. #includeiostream2.3.voidMin2Method(double&xyTopX,double&xyTopY,intX,intY,intnCount)4. /5. /功能描述:利用最小二乘法求斜/6./xyTopX-
2、截/7./xyTopY-斜率8.9.10.11.12.13.14.15.16.17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32.33.34.35.36.37.38.39.40.41.42.43.44.45.46./nCount-点数/inti;doubleSumX,SumY,SumXY,SumX2;SumX=0;SumX2=0;for(i=0;inCount;i+)SumX+=Xi;SumX2+=(Xi*Xi);SumY=0;for(i=0;inCount;i+)SumY+=Yi;SumXY=0;for(i=0;inCount;i+)SumXY+=(Xi*Yi);xyTopX=(SumX2*SumY-SumX*SumXY)/(nCount*SumX2-SumX*SumX);xyTopY=(nCount*SumXY-SumX*SumY)/(nCount*SumX2-SumX*SumX);intmain()doublexyTopX=0;doublexyTopY=0;intnCount=10;intx=1,2,3,4,5,6,7,8,9,10;inty=1,2,3,4,5,6,7,8,9,10;Min2Method(xyTopX,xyTopY,x,y,nCount);std:coutxyTopXxyTopXstd:endl;std:coutxyTopYxyTopYstd:endl;system(pause);
《最小二乘法直线拟合》由会员壹****1分享,可在线阅读,更多相关《最小二乘法直线拟合》请在金锄头文库上搜索。
土建监理工作职责职(七篇)
《让书香充盈学生的课外生活》.doc
伏羲班的启示精华本.doc
论文教学评价提高课堂教学质量指挥棒
怎样培养小学生的学习兴趣呢
国内元教育学思考.docx
1万吨无菌冷灌装pet瓶枣汁饮料生产线可研报告书.doc
20XX年四年级班级工作计划.docx
2019年高考物理二轮练习冲刺测试专项7电路
托班班主任工作计划参考范本(4篇)
初一上Unit2Topic2修改单词过关练习.doc
2023年农村清洁工程宣传实施方案.doc
隐藏工程模式帮你挑选iPod
范英男《引人注目的招牌》教学设计.doc
2023年顶岗实习学生总结(4篇).doc
计算机科学与技术学院暨软件学院专业准入实施细则
2023年开学工作会议主持词.docx
叠石水坝工程施工方案正式版(DOC 47页)
六年级日记【推荐】
我眼中的秋天初一作文600字.doc
2024-01-13 21页
2024-01-02 10页
2023-03-23 7页
2023-03-23 7页
2023-09-27 13页
2023-07-20 4页
2023-05-13 24页
2023-11-01 2页
2023-11-19 5页
2023-01-28 16页