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語言編寫最小二乘法線性擬合
首先你要知道最小二乘法的公式，然後使用數組實現。比如定義數組double x[10]代表十個橫坐標，double y[10]代表縱坐標，然後套公式用循環計算
③ 求一個最小二乘法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;i<n;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;i<T+1;i++)
for(j=0;j<T+1;j++)
{
res=0;
for(k=0;k<N;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;i<N;i++)
{
for(j=0;j<T+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;i<n;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;i<N;i++)
{
for(j=0;j<T+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;i<n;i++)
{
for(j=0;j<T+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;i<n;i++)
{
for(j=i;j<n;j++)
{
if(argu[i-1][i-1]==0)
{
for(p=i;p<n;p++)
{
if(argu[p][i-1]!=0)
break;
}
if(p==n)
{
printf("方程組無解!\n");
exit(0);
}
for(t=0;t<n+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;k<n+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;j>i;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;i<n;i++)
{
temp=0;
for(j=0;j<N;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;i<T+1;i++)
{
for(j=0;j<T+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;i<n;i++)
{
printf("a[%d]=%g,",i+1,a[i]);
}
printf("\n");
printf("擬合曲線方程為:\ny(x)=%g",a[0]);
for(i=1;i<n;i++)
{
printf(" + %g",a[i]);
for(j=0;j<i;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>
<PostDateTime>2007-4-19 19:23:57</PostDateTime>
</Reply>
<Reply>
<PostUserNickName></PostUserNickName>
<rank>一級(初級)</rank>
<ranknum>user1</ranknum>
<credit>100</credit>
<ReplyID>40389872</ReplyID>
<TopicID>5478010</TopicID>
<PostUserId>1526752</PostUserId>
<PostUserName>jiangxc2004</PostUserName>
<Point>0</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>
<PostDateTime>2007-4-19 19:26:11</PostDateTime>
</Reply>
<Reply>
<PostUserNickName></PostUserNickName>
<rank>一級(初級)</rank>
<ranknum>user1</ranknum>
<credit>100</credit>
<ReplyID>40390406</ReplyID>
<TopicID>5478010</TopicID>
<PostUserId>1526752</PostUserId>
<PostUserName>jiangxc2004</PostUserName>
<Point>0</Point>
<Content>
<![CDATA[
這樣就可以直接調用process()函數了!
二次擬合的話就把宏 T 成2;
擬合點的數目 N 也可以修改!
也可以去到注釋的部分進行返回值的調用!
④ 關於最小二乘法的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程序演算法
聯系我，我可以把幾種擬合方式的代碼法給你
⑥ 最小二乘法的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;i<n;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;i<T+1;i++)
for(j=0;j<T+1;j++)
{
res=0;
for(k=0;k<N;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;i<N;i++)
{
for(j=0;j<T+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;i<n;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;i<N;i++)
{
for(j=0;j<T+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;i<n;i++)
{
for(j=0;j<T+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;i<n;i++)
{
for(j=i;j<n;j++)
{
if(argu[i-1][i-1]==0)
{
for(p=i;p<n;p++)
{
if(argu[p][i-1]!=0)
break;
}
if(p==n)
{
printf("方程組無解!\n");
exit(0);
}
for(t=0;t<n+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;k<n+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;j>i;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;i<n;i++)
{
temp=0;
for(j=0;j<N;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;i<T+1;i++)
{
for(j=0;j<T+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;i<n;i++)
{
printf("a[%d]=%g,",i+1,a[i]);
}
printf("\n");
printf("擬合曲線方程為:\ny(x)=%g",a[0]);
for(i=1;i<n;i++)
{
printf(" + %g",a[i]);
for(j=0;j<i;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>
<PostDateTime>2007-4-19 19:23:57</PostDateTime>
</Reply>
<Reply>
<PostUserNickName></PostUserNickName>
<rank>一級(初級)</rank>
<ranknum>user1</ranknum>
<credit>100</credit>
<ReplyID>40389872</ReplyID>
<TopicID>5478010</TopicID>
<PostUserId>1526752</PostUserId>
<PostUserName>jiangxc2004</PostUserName>
<Point>0</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>
<PostDateTime>2007-4-19 19:26:11</PostDateTime>
</Reply>
<Reply>
<PostUserNickName></PostUserNickName>
<rank>一級(初級)</rank>
<ranknum>user1</ranknum>
<credit>100</credit>
<ReplyID>40390406</ReplyID>
<TopicID>5478010</TopicID>
<PostUserId>1526752</PostUserId>
<PostUserName>jiangxc2004</PostUserName>
<Point>0</Point>
<Content>
<![CDATA[
這樣就可以直接調用process()函數了!
二次擬合的話就把宏 T 成2;
擬合點的數目 N 也可以修改!
也可以去到注釋的部分進行返回值的調用!
⑦ c語言最小二乘法 求助
這個平面上隨機點產生我直接用的一維產生方法，不知道你要求的是什麼樣的分布
#include<time.h>
#include<stdio.h>
#include<stdlib.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;i<NODE_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;i<NODE_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;i<NODE_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語言進行最小二乘法擬合，這個代碼錯在哪
a=(10.0*sumxy-sumx*sumy)/(10.0*sum-sumx*sumx); /*帶入a=(NΣxy-ΣxΣy)/(NΣx^2(Σx)^2)
NΣx^2(Σx)^2 這邊中間到底有沒有減號？你注釋裡面沒有減號，但是代碼10.0*sum-sumx*sumx中是減法，你看看這邊是否有問題
⑨ 最小二乘法程序（C語言，matlab都可以）
t=20;%t為隨機數個數
x=1:t;
y=rand(1,t);%隨即產生t個隨機數，范圍是0到1，你也可以都乘以某個常數
n=1;%直線擬合n就取1，
p=polyfit(x,y,n);%計算n次多項式系數p
xi=linspace(0,t,100);%linspace用於產生x1,x2之間的N點行矢量
yi=polyval(p,xi);%計算多項式的值
plot(x,y,'.',xi,yi);
xlabel('x'),ylabel('f(x)')
title('最小二乘法擬合')
p%p即為擬合後的直線的一次項系數和常數項
——————————————————————————————————————
注意：隨機產生的點，因為數據方差比較大，擬合成的直線對於某個點來說相差較大。
比如我測試的p=0.00530.4909，那麼，擬合直線就是y=0.0053x+0.4909。
