最小二乘法的编程实现.pdfVIP

百川东到海，何时复西归?少壮不努力，老大徒伤悲。——汉乐府

1、最小二乘法：

1

1）（用AT方法计算逆矩阵）

A

#includestdio.h

#includemath.h

#includestdlib.h

#includestring.h

#includeiostream.h

#defineN200

#definen9

voidGetdata(doublesun[N])//从txt文档中读取数据（小数）

{

chardata;

charsunpot[10]={0000000000};//为防止结果出现‘烫’字

inti=0,j=0;

doubled;

FILE*fp=fopen(新建文本文档.txt,r);

if(!fp)

{

printf(cantopenfile\n);

}

while(!feof(fp))

{

data=fgetc(fp);

if(data!=\n)

{

sunpot[i]=data;

i++;

}

elseif(data==\n)

{

sunpot[i]=\0;//给定结束符

d=atof(sunpot);//将字符串转换成浮点数

sun[j]=d;

j++;

i=0;//将i复位

}

}

}

voidNormal(doublesun[N],doublesun1[N])//将数据进行标准化

{

百川东到海，何时复西归?少壮不努力，老大徒伤悲。——汉乐府

doublemean,temp=0,variance=0;

inti;

for(i=0;iN;i++)

temp=temp+sun[i];

mean=temp/N;

for(i=0;iN;i++)

{

sun1[i]=sun[i]-mean;

}

for(i=0;iN;i++)

{

variance=variance+sun1[i]*sun1[i];

}