实验名称：最小二乘拟合学 院：数学学院班级姓名：学号：实验日期 2015 年05月26日自评成绩：97一、实验目的二、实验内容与结果题目1：最小二乘拟合源程序1#include#include using namespace std;float f(int i, int j, int n, float x20)/计算矩阵的函数(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)/计算右端项的函数(float 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为简化平方根法解方程组的函数(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;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 l 10) (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;i+)(coutl it;else(for (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;else(b=f2(r1, n)/f2(r, n);for(i=0;in;i+)(di=r1i+b*di;ri=r1i;int main() /主函数(int n,i,j;cout以最小二乘拟合二次多项式为例endl;float x10,y10,a1010,b10;cout请输入数据的个数：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请输入起始迭代序列l10:endl;float l10,e;for(i=0;i3;i+)(coutlili;coute;f5(3,e,a,b,l);return 0;运行结果1