c语言实现最小二乘法（最小二乘法函数）

本篇文章给大家谈谈c语言实现最小二乘法，以及最小二乘法函数对应的知识点，希望对各位有所帮助，不要忘了收藏本站喔。

本文目录一览：

1、单片机c语言的最小二乘法怎么实现，求代码2、计算方法中最小二乘法如何用C语言编程3、用C语言编写最小二乘法线性拟合4、关于最小二乘法的c语言程序5、求一个最小二乘法C语言程序6、c语言最小二乘法 求助

单片机c语言的最小二乘法怎么实现，求代码

#include stdio.h

void main ()

{

int num,i;

float x,y,l,m,n,p,a,b;

i=1;

l=0.0;

m=0.0;

n=0.0;

p=0.0;

printf (“请输入你想计算的x，y的个数：”);

scanf(“%d”,num);

if (num=1)

{

while (i=num);

{

printf(“请输入x的值”);

scanf (“%lf”,x);

printf(“请输入y的值”);

scanf (“%lf”,y);

l+=x;

m+=y;

n+=x*y;

p+=x*x;

i++;

}

a=(num*n-l*m)/(num*p-l*l);

b=(p*m-n*l)/(num*p-l*l);

printf(“最小二乘法所算得的斜率和截距分别为%f和%f\n”,a,b);

}

else printf(“mun”输入有误！);

}

计算方法中最小二乘法如何用C语言编程

#include stdio.h

#include math.h

#define epsilon 1e-6

void nihe1(int n,int m,float sum_x,float sum_y,float sum_xy,float x2);

void nihe2(int n,int m,float sum_x,float sum_y,float sum_xy,float x2,float x2y,float x3,float x4);

int main(){

float x[100]={0.0};

float y[100]={0.0};

int n,i,flag=1;

float sum_y=0.0,sum_x=0,x2=0,sum_xy=0.0,x3=0,x4=0,x2y=0.0;

printf(“请你输入需要测试的数据(先输入x[]，后输入y[])的个数:”);

scanf(“%d”,n);

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

scanf(“%f”,x[i]);}

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

scanf(“%f”,y[i]);}

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

sum_x += x[i];

sum_y += y[i];

sum_xy += x[i]*y[i];

x2 += x[i]*x[i];

x2y += x[i]*x[i]*y[i];

x3 += x[i]*x[i]*x[i];

x4 += x[i]*x[i]*x[i]*x[i];}

printf(“—————请你输入的要拟合的函数——————\n”);

printf(” 1、拟合一次函数\n”);

printf(” 2、拟合二次函数\n”);

scanf(“%d”,flag);

switch(flag){

case 1:

nihe1(n,flag+1,sum_x,sum_y,sum_xy,x2); break;

case 2:

nihe2(n,flag+1,sum_x,sum_y,sum_xy,x2,x2y,x3,x4); break;

default:

printf(“ERROR\n”);}

return 0;}

void nihe1(int n,int m,float sum_x,float sum_y,float sum_xy,float x2){

int i,k,j;

float t,s=0;

float a[2][3] = {{(float)n,sum_x,sum_y},{sum_x,x2,sum_xy}};

n=m;

//if(m == 3)

// a[3][4] = {{n,sum_x,sum_y},{sum_x,x2,x3,sum_xy},{x2,x3,x4,x2y}};

for(k=0;kn-1;k++) {

for(i=k+1;in;i++)

if( abs((int)a[i][k]) abs((int)a[k][k]) )

for(j=k;jn+1;j++) {

t=a[k][j];

a[k][j]=a[i][j];

a[i][j]=t; }

if( abs((int)a[k][k]) epsilon) {

printf(“\nError,主元消去法 cann’t be durable,break at %d!\n”,k+1);

return; }

for(i=k+1;in;i++){

a[i][k]=a[i][k] / a[k][k];

for(j=k+1;jn+1;j++)

a[i][j]=a[i][j]-a[i][k] * a[k][j]; }}

a[n-1][n]=a[n-1][n] / a[n-1][n-1];

for(k=n-2;k=0;k–) {

s=0;

for(j=k+1;jn;j++)

s+=a[k][j]*a[j][n];

a[k][n]=( a[k][n]-s ) / a[k][k]; }

printf(“\n*****The Result*****\n”);

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

printf(” x[%d]=%.4f\n”,i+1,a[i][n]);

printf(“函数为：p(x) = %.4f + (%.4f)*x\n”,a[0][n],a[1][n]);

getchar();}

void nihe2(int n,int m,float sum_x,float sum_y,float sum_xy,float x2,float x2y,float x3,float x4){

int i,k,j;

float t,s=0;

float a[3][4]=

{{(float)n,sum_x,x2,sum_y},{sum_x,x2,x3,sum_xy},{x2,x3,x4,x2y}};

n=m;

for(k=0;kn-1;k++) {

for(i=k+1;in;i++)

if( abs((int)a[i][k]) abs((int)a[k][k]) )

for(j=k;jn+1;j++) {

t=a[k][j];

a[k][j]=a[i][j];

a[i][j]=t; }

if( abs((int)a[k][k]) epsilon) {

printf(“\nError,主元消去法 cann’t be durable,break at %d!\n”,k+1);

return; }

for(i=k+1;in;i++){

a[i][k]=a[i][k] / a[k][k];

for(j=k+1;jn+1;j++)

a[i][j]=a[i][j]-a[i][k] * a[k][j]; } }

a[n-1][n]=a[n-1][n] / a[n-1][n-1];

for(k=n-2;k=0;k–) {

s=0;

for(j=k+1;jn;j++)

s+=a[k][j]*a[j][n];

a[k][n]=( a[k][n]-s ) / a[k][k]; }

printf(“\n*****The Result*****\n”);

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

printf(” x[%d]=%.4f\n”,i+1,a[i][n]);

printf(“函数为：p(x) = %.4f + (%.4f)*x + (%.4f)*x*x\n”,a[0][n],a[1][n],a[2][n]);

getchar();}

用C语言编写最小二乘法线性拟合

首先你要知道最小二乘法的公式，然后使用数组实现。比如定义数组double x[10]代表十个横坐标，double y[10]代表纵坐标，然后套公式用循环计算

关于最小二乘法的c语言程序

已经改正语法的错误，如果有问题

是你的算法有误

#include

int

main

()

{

int

num,i;

float

x,y,l,m,n,p,a,b;

i=1;

l=0.0;

m=0.0;

n=0.0;

p=0.0;

printf

(“请输入你想计算的x，y的个数：”);

scanf(“%d”,num);

while

(i=num)

//去掉分号

{

printf(“请输入x的值：”);

scanf

(“%f”,x);

//加上取地址符

printf(“请输入y的值：”);

scanf

(“%f”,y);

//加上取地址符

l+=x;

m+=y;

n+=x*y;

p+=x*x;

i++;

}

a=(num*n-l*m)/(num*p-l*l);

b=(p*m-n*l)/(num*p-l*l);

printf(“最小二乘法所算得的斜率和截距分别为%f和%f\n”,a,b);

}

求一个最小二乘法C语言程序

#include stdio.h

#include conio.h

#include math.h

#include process.h

#define N 5//N个点

#define T 3 //T次拟合

#define W 1//权函数

#define PRECISION 0.00001

float pow_n(float a,int n)

{

int i;

if(n==0)

return(1);

float res=a;

for(i=1;in;i++)

{

res*=a;

}

return(res);

}

void mutiple(float a[][N],float b[][T+1],float c[][T+1])

{

float res=0;

int i,j,k;

for(i=0;iT+1;i++)

for(j=0;jT+1;j++)

{

res=0;

for(k=0;kN;k++)

{

res+=a[i][k]*b[k][j];

c[i][j]=res;

}

}

}

void matrix_trans(float a[][T+1],float b[][N])

{

int i,j;

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

{

for(j=0;jT+1;j++)

{

b[j][i]=a[i][j];

}

}

}

void init(float x_y[][2],int n)

{

int i;

printf(“请输入%d个已知点：\n”,N);

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

{

printf(“(x%d y%d):”,i,i);

scanf(“%f %f”,x_y[i][0],x_y[i][1]);

}

}

void get_A(float matrix_A[][T+1],float x_y[][2],int n)

{

int i,j;

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

{

for(j=0;jT+1;j++)

{

matrix_A[i][j]=W*pow_n(x_y[i][0],j);

}

}

}

void print_array(float array[][T+1],int n)

{

int i,j;

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

{

for(j=0;jT+1;j++)

{

printf(“%-g”,array[i][j]);

}

printf(“\n”);

}

}

void convert(float argu[][T+2],int n)

{

int i,j,k,p,t;

float rate,temp;

for(i=1;in;i++)

{

for(j=i;jn;j++)

{

if(argu[i-1][i-1]==0)

{

for(p=i;pn;p++)

{

if(argu[p][i-1]!=0)

break;

}

if(p==n)

{

printf(“方程组无解!\n”);

exit(0);

}

for(t=0;tn+1;t++)

{

temp=argu[i-1][t];

argu[i-1][t]=argu[p][t];

argu[p][t]=temp;

}

}

rate=argu[j][i-1]/argu[i-1][i-1];

for(k=i-1;kn+1;k++)

{

argu[j][k]-=argu[i-1][k]*rate;

if(fabs(argu[j][k])=PRECISION)

argu[j][k]=0;

}

}

}

}

void compute(float argu[][T+2],int n,float root[])

{

int i,j;

float temp;

for(i=n-1;i=0;i–)

{

temp=argu[i][n];

for(j=n-1;ji;j–)

{

temp-=argu[i][j]*root[j];

}

root[i]=temp/argu[i][i];

}

}

void get_y(float trans_A[][N],float x_y[][2],float y[],int n)

{

int i,j;

float temp;

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

{

temp=0;

for(j=0;jN;j++)

{

temp+=trans_A[i][j]*x_y[j][1];

}

y[i]=temp;

}

}

void cons_formula(float coef_A[][T+1],float y[],float coef_form[][T+2])

{

int i,j;

for(i=0;iT+1;i++)

{

for(j=0;jT+2;j++)

{

if(j==T+1)

coef_form[i][j]=y[i];

else

coef_form[i][j]=coef_A[i][j];

}

}

}

void print_root(float a[],int n)

{

int i,j;

printf(“%d个点的%d次拟合的多项式系数为:\n”,N,T);

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

{

printf(“a[%d]=%g,”,i+1,a[i]);

}

printf(“\n”);

printf(“拟合曲线方程为:\ny(x)=%g”,a[0]);

for(i=1;in;i++)

{

printf(” + %g”,a[i]);

for(j=0;ji;j++)

{

printf(“*X”);

}

}

printf(“\n”);

}

void process()

{

float x_y[N][2],matrix_A[N][T+1],trans_A[T+1][N],coef_A[T+1][T+1],coef_formu[T+1][T+2],y[T+1],a[T+1];

init(x_y,N);

get_A(matrix_A,x_y,N);

printf(“矩阵A为：\n”);

print_array(matrix_A,N);

matrix_trans(matrix_A,trans_A);

mutiple(trans_A,matrix_A,coef_A);

printf(“法矩阵为：\n”);

print_array(coef_A,T+1);

get_y(trans_A,x_y,y,T+1);

cons_formula(coef_A,y,coef_formu);

convert(coef_formu,T+1);

compute(coef_formu,T+1,a);

print_root(a,T+1);

}

void main()

{

process();

}

]]

/Content

PostDateTime2007-4-19 19:23:57/PostDateTime

/Reply

Reply

PostUserNickName/PostUserNickName

rank一级(初级)/rank

ranknumuser1/ranknum

credit100/credit

ReplyID40389872/ReplyID

TopicID5478010/TopicID

PostUserId1526752/PostUserId

PostUserNamejiangxc2004/PostUserName

Point0/Point

Content

![CDATA[

你可以改一下

不从终端输入,直接在程序中给出参数

请输入5个已知点：

(x0 y0):-2 -0.1

(x1 y1):-1 0.1

(x2 y2):0 0.4

(x3 y3):1 0.9

(x4 y4):2 1.6

矩阵A为：

1 -2 4 -8

1 -1 1 -1

1 0 0 0

1 1 1 1

1 2 4 8

法矩阵为：

5 0 10 0

0 10 0 34

10 0 34 0

0 34 0 130

5个点的3次拟合的多项式系数为:

a[1]=0.408571, a[2]=0.391667, a[3]=0.0857143, a[4]=0.00833333,

拟合曲线方程为:

y(x)=0.408571 + 0.391667*X + 0.0857143*X*X + 0.00833333*X*X*X

]]

/Content

PostDateTime2007-4-19 19:26:11/PostDateTime

/Reply

Reply

PostUserNickName/PostUserNickName

rank一级(初级)/rank

ranknumuser1/ranknum

credit100/credit

ReplyID40390406/ReplyID

TopicID5478010/TopicID

PostUserId1526752/PostUserId

PostUserNamejiangxc2004/PostUserName

Point0/Point

Content

![CDATA[

这样就可以直接调用process()函数了!

二次拟合的话就把宏 T 成2;

拟合点的数目 N 也可以修改!

也可以去到注释的部分进行返回值的调用!

c语言最小二乘法 求助

这个平面上随机点产生我直接用的一维产生方法，不知道你要求的是什么样的分布

#includetime.h

#includestdio.h

#includestdlib.h

////////////////设置区间///////////////////

#define MAX_X (10000)

#define MIN_X (-10000)

#define MAX_Y (10000)

#define MIN_Y (-10000)

#define NODE_NUM (10) //随机点个数

int main()

{

////////////////////生成随机点/////////////////

srand((unsigned int)time(NULL));//设置随机种子

int i;

float interval_x,interval_y;

interval_x = MAX_X-MIN_X;

interval_y = MAX_Y-MIN_Y;

float rand_x[NODE_NUM],rand_y[NODE_NUM];

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

{

rand_x [i] = (rand()/(float)RAND_MAX)*interval_x+MIN_X;

rand_y [i] = (rand()/(float)RAND_MAX)*interval_y+MIN_Y;

printf(“rand_x[%d] = %f, rand_y[%d] = %f\n”,i,rand_x[i],i,rand_y[i]);

}

////////////////////最小二乘法/////////////////

//mean_x，mean_y均值，

//y=ax+b

//sum_xy为x[i]*y[i]的和

//sum_xx为x平方的和

//wiki 最小二乘法给出了两个公式，结果相同 ，选一个就行了

double mean_x,mean_y,a,b,sum_xy,sum_xx;

mean_x = 0;

mean_y = 0;

sum_xx = 0;

sum_xy = 0;

//////////////公式一////////////

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

{

mean_x += rand_x[i];

mean_y += rand_y[i];

sum_xy += rand_x[i]*rand_y[i];

sum_xx += rand_x[i]*rand_x[i];

}

mean_x = mean_x/NODE_NUM;

mean_y = mean_y/NODE_NUM;

a = (sum_xy-NODE_NUM*mean_x*mean_y)/(sum_xx-NODE_NUM*mean_x*mean_x);

b = mean_y-a*mean_x;

//////////////公式二////////////

float fenzi,fenmu,a1,b1;

fenzi = 0;

fenmu = 0;

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

{

fenzi += (rand_x[i]-mean_x)*(rand_y[i]-mean_y);

fenmu += (rand_x[i]-mean_x)*(rand_x[i]-mean_x);

}

a1 = fenzi/fenmu;

b1 = mean_y-a1*mean_x;

printf(“方法一：y = (%f)*x+(%f)\n”,a,b);

printf(“方法二：y = (%f)*x+(%f)\n”,a1,b1);

return 0;

}

关于c语言实现最小二乘法和最小二乘法函数的介绍到此就结束了，不知道你从中找到你需要的信息了吗 ？如果你还想了解更多这方面的信息，记得收藏关注本站。