曲線擬合最小二乘法 C 程序

課題八曲線擬合的最小二乘法

實(shí)驗(yàn)?zāi)繕?biāo)：

在某冶煉過程中，通過實(shí)驗(yàn)檢測(cè)得到含碳量與時(shí)間關(guān)系的數(shù)據(jù)如下，試求含碳量y與時(shí)間t

內(nèi)在關(guān)系的擬合曲線。

t05101520253035404550Y01.272.162.863.443.874.154.374.514.584.02最小二乘法：

#include<iostream>#include<math.h>usingnamespacestd;

intArray(double***Arr,intn){

double**p;inti;

p=(double**)malloc(n*sizeof(double*));if(!p)return0;for(i=0;i<n;i++){

p[i]=(double*)malloc(n*sizeof(double));if(!p[i])return0;}*Arr=p;return1;}

voidmain(){

intn,i,j,k;double**A,*B,*S;S=newdouble[3];B=newdouble[3];

cout<<"共2有?DN個(gè)?節(jié)¨2點(diǎn)ì?：êo";cin>>n;

cout<<"請(qǐng)?輸o?入¨?節(jié)¨2點(diǎn)ì?值|ì（ê?§Xi）ê?：êo"<<endl;

double*X;X=newdouble[n];for(i=0;i<n;i++){

cin>>X[i];}

cout<<"請(qǐng)?輸o?入¨?節(jié)¨2點(diǎn)ì?函?￥數(shù)oy值|ì

}

for(i=0;i<n;i++){

B[0]+=Y[i];B[1]+=X[i]*Y[i];}

for(i=0;i<3;i++){

B[i]=0;}}

for(i=0;i<n;i++){

A[0][0]+=1;//cout<<A[0][0];A[0][1]+=X[i];A[0][2]+=X[i]*X[i];A[1][0]=A[0][1];A[1][1]=A[0][2];

A[1][2]+=X[i]*X[i]*X[i];A[2][0]=A[0][2];A[2][1]=A[1][2];

A[2][2]+=X[i]*X[i]*X[i]*X[i];}

else{

for(i=0;i<3;i++){

for(j=0;j<3;j++){

A[i][j]=0;

}

if(!Array(&A,3))

cout<<"內(nèi)¨2存??分¤?配?失o?ì敗?¨1！ê?";

（ê?§Yi）ê?：êo"<<endl;

double*Y;

Y=newdouble[n];for(i=0;i<n;i++){

cin>>Y[i];

曲線擬合最小二乘法 C 程序.doc下載

久久建筑網(wǎng)i5h4u.cn提供大量：建筑圖紙、施工方案、工程書籍、建筑論文、合同表格、標(biāo)準(zhǔn)規(guī)范、CAD圖紙等內(nèi)容。