非線性最小二乘法C語言代碼

非線性最小二乘法C語言代碼main(){intn,i,flag2;charflag1='y';floatar[50],br[50],x,y,xe,ye,xye,xxe,sx,sy,sxy,sxx,a,b;printf("\n歡迎使用最小二乘法數(shù)據(jù)處理程序\n");printf("\n說明：本程序運(yùn)行結(jié)果保留小數(shù)點(diǎn)后三位\n");for(;flag1=='y'||flag1=='Y';){printf("\n請(qǐng)輸入您要處理的數(shù)據(jù)的組數(shù):");printf("\n*****提示：本程序定義一對(duì)x,y值為一組數(shù)據(jù):");scanf("%d",&n);if(n>50){printf("\n對(duì)不起，本程序暫時(shí)無法處理50組以上的數(shù)據(jù)");continue;}printf("\n請(qǐng)選擇您的數(shù)據(jù)的處理方式:");printf("\n\t1.y與x為一次線性關(guān)系");printf("\n\t2.y與x的2次為線性關(guān)系\n");scanf("%d",&flag2);if(flag2>2||flag2<1){printf("\n對(duì)不起，您的輸入不正確\n");continue;}for(i=0;i<n;i++){printf("\n請(qǐng)輸入第%2d個(gè)x的值\tx%2d=",i+1,i+1);scanf("%f",&ar[i]);printf("\n請(qǐng)輸入對(duì)應(yīng)的y的值:\ty%2d=",i+1);scanf("%f",&br[i]);}if(flag2!=1){for(i=0;i<n;i++)br[i]=br[i]/(ar[i]*ar[i]);}sx=sy=sxx=sxy=0;for(i=0;i<n;i++){sx=sx+ar[i];sy=sy+br[i];sxx=sxx+ar[i]*ar[i];sxy=sxy+ar[i]*br[i];}xe=sx/n;ye=sy/n;xye=sxy/n;xxe=sxx/n;b=(xye-xe*ye)/(xxe-xe*xe);a=ye-b*xe;printf("\n對(duì)您輸入的數(shù)據(jù)的處理已經(jīng)完成，結(jié)果如下:");printf("\n\ta=%8.3f\n\tb=%8.3f\n",a,b);printf("\nb即為擬合直線的斜率，a為截距\n");printf("\n*********************************************************\n");printf("\n是否繼續(xù)使用本程序處理數(shù)據(jù)？(y/n)?");scanf("%c",&flag1);if(flag1=='y'||flag1=='Y')continue;elseif(flag1=='n'||flag1=='N')brea