最小二乘擬合的C++實現(xiàn)(共5頁)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上實驗名稱：最小二乘擬合學(xué) 院：_數(shù)學(xué)學(xué)院_ 班 級 姓 名： 學(xué) 號： 實 驗 日 期 2015 年 05 月 26 日自評成績：97 一、實驗?zāi)康亩?、實驗?nèi)容與結(jié)果題目1：最小二乘擬合源程序1#include#include using namespace std;float f(int i,int j,int n,float x20)/計算矩陣的函數(shù)int k;float s=0;for(k=0;kn;k+)s=s+pow(xk,i+j);return s;float f1(int i,int n,float y20,float x20)/計算右端項的函數(shù)floa

2、t s=0;int k;for(k=0;kn;k+)s=s+pow(xk,i)*yk;return s;float f2(float a10,double n)/f2,f3,f4,f5為簡化平方根法解方程組的函數(shù)float s=0;int i;for(i=0;in;i+)s=s+ai*ai;return s;float f3(float r10,float a1010,int n)float b10,m=0;int i,j;for(i=0;in;i+)float s=0;for(j=0;jn;j+)s=s+rj*aji;bi=s;for(i=0;in;i+)m=m+bi*ri;return m

3、;float *f4(float a1010,float b10,int n)float *r;r=new float10;int i,j;for(i=0;in;i+)float s=0;for(j=0;jn;j+)s=s+aij*bj;ri=s;return r;void f5(int n,float e,float A1010,float B10,float l10)int i,k;float a,b,r10,r110,d10;float *ax;ax=f4(A,l,n);for(i=0;in;i+)ri=Bi-axi;di=ri;if(sqrt(f2(r,n)=e)for(i=0;in;

4、i+)coutlit;elsefor(k=0;kn;k+)a=f2(r,n)/f3(d,A,n);for(i=0;in;i+)li=li+a*di;ax=f4(A,l,n);for(i=0;in;i+)r1i=Bi-axi;if(sqrt(f2(r1,n)=e|k+1=n)for(i=0;in;i+)coutci=lin;break;elseb=f2(r1,n)/f2(r,n);for(i=0;in;i+)di=r1i+b*di;ri=r1i;int main()/主函數(shù)int n,i,j;cout以最小二乘擬合二次多項式為例endl;float x10,y10,a1010,b10;cout請輸入數(shù)據(jù)的個數(shù)：n;for(i=0;in;i+)cout請輸入xi=xi;cout請輸入yi=yi;for(i=0;i3;i+)for(j=0;j3;j+)aij=f(i,j,n,x);for(i=0;i3;i+)bi=f1(i,n,y,x);cout由共軛梯度法求解方程組endl;cout請輸入起始迭代序列l(wèi)10:endl;float