《用C语言实现的曲线拟合的最小二乘法》由会员分享，可在线阅读，更多相关《用C语言实现的曲线拟合的最小二乘法（6页珍藏版）》请在金锄头文库上搜索。
1、实验名称：曲线拟合的最小二乘法实验目的了解曲线拟合的最小二乘法实验类型设计型实验环境Windows XP TC实验内容相关知识：已知Ca,b 中函数f(x)的一组实验数据(x,y ) ( i=0,l,m),其中y =f(x )。i iii设9 (x)( j 二 0,1,n;n m)是Ca,b上线性无关函数族。在 = span (x),p (x),p (x)中 j01n找函数f(x)曲线拟合的最小二乘解S*(x) = a 9 (x)，其法方程(组)为: j jj=0 (p ,p )a = d (k = 0,1,n)k j j kj=o其中,x )9 (x )9 (x )i j i k ii=0k
2、=0,l,，n瓦x )f (x )9 (x )三 dii k iki=0特别是，求函数f(x)曲线拟合的线性最小二乘解S*(x) = ax + b的计算公式为:(区x2)(区 y )-(区x )(x y )iiii ib =i=0(m +1)区 x2 - ( x )2iii=0i=0(m +1)区 x yi i-(区 x )( y )ii-(m +1)x )2ix 2-(迟i数据结构：两个一维数组或一个二维数组算法设计：(略)实验用例：已知函数y=f(x)的一张表:x0102030405060708090y6867.166.465.664.661.861.060.860.460试验要求：利用曲
3、线拟合的线性最小二乘法求被逼近函数f(x)在点x=55处的近似值，并画 出实验数据和直线。编写代码：#include#include#includedouble qiuhel(double a102, int p)int i;double y;y=0；for(i=0;i10;i+)y=y+aip;return y;double qiuhe2(double a102, int p)int i;double y=0;for(i=0;i1;i)if(i!=4)ci=n%10+48;n二n/10;cll=O;x1=x;set bkcolor(7);set color(14);setli nes ty
4、le(0,0,3);drawpoly(3,arw);drawpoly(3,arwl);line(120,240,520,240);line(320,40,320,440);x=0;set color(2);setli nes ty le(O,O,l);line(0+320),(in t)(240-nihe(a,0) *q),(90+320),(in t)(240-nihe(a,90) *q); set color(3);outt ex txy(320,30,Y);ou tt ex txy(310,245,0);outt ex txy(525,240,X);outtext xy(xl+330),
5、(240-nihe(a,xl)-10),c);settextsty le(4,0,4);outt ex txy(450,400,Nihe.);for(i=0;i=9;i+)put pixel(ai0+320),(240-ail *q),ll);set color(4);setli nes ty le(l,0,l);line(xl+320),(240-nihe(a,x 1)* q)-80),(xl+320),(240-nihe(a,xl) *q)+120);get ch();closegraph();实验结果（测试用例、实验结果）实验总结与心得通过本次实验，对曲线拟合的最小二乘法有了更深刻的了解!
