1樓:有問必答
//最小二乘法
//x[n] y[n] 已知輸入
//n輸入點個數
//a[m] 返回m-1次擬合多項
式的m個係數
//m 擬合多項式的項數,即擬合多項式的最高次為m-1//dt[3] dt[0]返回擬合多項式與各資料點誤差的平方和,dt[1]返回擬合多項式與各資料點誤差的絕對值之和dt[2]返回擬合多項式與各資料點誤差的絕對值的最大值////擬合多項式的輸出
//y(x) = a0 + a1(x-x) + a2(x-x)^2 + …… am(x-x)^m
// 其中x為已知點x的平均值
******************************************/
#include "math.h"
void pir1(x,y,n,a,m,dt)int n,m;
double x,y,a,dt;
c=c/d1; p=p/d1;
a[0]=c*b[0];
if (m>1)
c=c/d2; p=g/d2; q=d2/d1;
d1=d2;
a[1]=c*t[1]; a[0]=c*t[0]+a[0];
}for (j=2; j<=m-1; j++)c=c/d2; p=g/d2; q=d2/d1;
d1=d2;
a[j]=c*s[j]; t[j]=s[j];
for (k=j-1; k>=0; k--)}dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
for (i=0; i<=n-1; i++)return;}
2樓:
先求解一個最小二乘問題,即解一個超定方程組(得到三次曲線係數)再用mfc中的畫圖函式離散的畫出這個三次曲線(注意mfc中的座標是矩陣座標,也就是說原點位於左上角)
最小二乘法三次多項式曲線擬合 演算法 c++ 實現,該怎麼處理
3樓:匿名使用者
/******************************************
//參考 《常用演算法程式集 (c語言描述 第三版)》
//最小二乘法
//x[n] y[n] 已知輸入
//n輸入點個數
//a[m] 返回m-1次擬合多項式的m個係數//m 擬合多項式的項數,即擬合多項式的最高次為m-1//dt[3] dt[0]返回擬合多項式與各資料點誤差的平方和,dt[1]返回擬合多項式與各資料點誤差的絕對值之和dt[2]返回擬合多項式與各資料點誤差的絕對值的最大值////擬合多項式的輸出
//y(x) = a0 + a1(x-x) + a2(x-x)^2 + …… am(x-x)^m
// 其中x為已知點x的平均值
******************************************/
#include "math.h"
void pir1(x,y,n,a,m,dt)int n,m;
double x,y,a,dt;
c=c/d1; p=p/d1;
a[0]=c*b[0];
if (m>1)
c=c/d2; p=g/d2; q=d2/d1;
d1=d2;
a[1]=c*t[1]; a[0]=c*t[0]+a[0];
}for (j=2; j<=m-1; j++)c=c/d2; p=g/d2; q=d2/d1;
d1=d2;
a[j]=c*s[j]; t[j]=s[j];
for (k=j-1; k>=0; k--)}dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
for (i=0; i<=n-1; i++)return;}
求c或c++語言編寫的用最小二乘法進行曲線擬合
4樓:匿名使用者
你的近似解析表示式為y=at+bt^2+ct^2
是不是想寫成為y=at+bt^2+ct^3
但是實際擬合出來的表示式為y=a[3]+a[2]t+a[1]t^2+a[0]t^3會有個常數項的。
簡單的講,所謂擬合是指已知某函式的若干離散函式值,通過調整該函式中若干待定係數f(λ1, λ2,…,λ3), 使得該函式與已知點集的差別(最小二乘意義)最小。如果待定函式是線性,就叫線性擬合或者線性迴歸(主要在統計中),否則叫作非線性擬合或者非線性迴歸。表示式也可以是分段函式,這種情況下叫作樣條擬合。
曲線擬合:
#include
#include
#include
#include
**ooth(double *x,double *y,double *a,int n,int m,double *dt1,double *dt2,double *dt3);
void main()
y = (double *)calloc(n,sizeof(double));
if(y == null)
a = (double *)calloc(n,sizeof(double));
if(a == null)
for(i=1;i<=n;i++)
y[0]=0;
y[1]=1.27;
y[2]=2.16;
y[3]=2.86;
y[4]=3.44;
y[5]=3.87;
y[6]=4.15;
y[7]=4.37;
y[8]=4.51;
y[9]=4.58;
y[10]=4.02;
y[11]=4.64;
/*x[i-1]點對應的y值是擬合已知值*/
呼叫擬合函式*/
for(i=1;i<=m;i++)
printf("a[%d] = %.10f\n",(i-1),a[i-1]);
printf("擬合多項式與資料點偏差的平方和為:\n");
printf("%.10e\n",dt1);
printf("擬合多項式與資料點偏差的絕對值之和為:\n");
printf("%.10e\n",dt2);
printf("擬合多項式與資料點偏差的絕對值最大值為:\n");
printf("%.10e\n",dt3);
free(x); /*釋放儲存空間*/
free(y); /*釋放儲存空間*/
free(a); /*釋放儲存空間*/
} **ooth(double *x,double *y,double *a,int n,int m,double *dt1,double *dt2,double *dt3)//(x,y,a,n,m,dt1,dt2,dt3 )
//double *x; /*實型一維陣列,輸入引數,存放節點的xi值*/
//double *y; /*實型一維陣列,輸入引數,存放節點的yi值*/
//double *a; /*雙精度實型一維陣列,長度為m。返回m一1次擬合多項式的m個係數*/
//int n; /*整型變數,輸入引數,給定資料點的個數*/
//int m; /*整型變數,輸入引數,擬合多項式的項數*/
//double *dt1; /*實型變數,輸出引數,擬合多項式與資料點偏差的平方和*/
//double *dt2; /*實型變數,輸出引數,擬合多項式與資料點偏差的絕對值之和*/
//double *dt3; /*實型變數,輸出引數,擬合多項式與資料點偏差的絕對值最大值*/
t = (double *)calloc(n,sizeof(double));
if(t == null)
b = (double *)calloc(n,sizeof(double));
if(b == null)
z = 0;
for(i=1;i<=n;i++)
z=z+x[i-1]/n; /*z為各個x的平均值*/
b[0]=1;
d1=n;
p=0;
c=0;
for(i=1;i<=n;i++)
c=c/d1;
p=p/d1;
a[0]=c*b[0];
if(m>1)
c=c/d2;
p=g/d2;
q=d2/d1;
d1=d2;
a[1]=c*t[1];
a[0]=c*t[0]+a[0];
} for(j=3;j<=m;j++)
c=c/d2;
p=g/d2;
q=d2/d1;
d1=d2;
a[j-1]=c*s[j-1];
t[j-1]=s[j-1];
for(k=j-1;k>=1;k--) }
*dt1=0;
*dt2=0;
*dt3=0;
for(i=1;i<=n;i++)
/*釋放儲存空間*/
free(s);
free(t);
free(b);
return(1);}
5樓:
#include
void main()
; double y[21] = ;
double midx;
double midy;
midx=0;
midy=0;
///求平均值
for(int i=0;i<21;++i)///求斜率
double linek;
linek=0;
double tempmu;
double tempzi;
tempmu=0;
tempzi=0;
for(int v=0;v<21;++v)if(tempmu==0)
linek=tempzi/tempmu;
////求截距
double lineb;
lineb=midy-linek*midx;
////////////
cout <<"直線方程:"<誤差
double delta;
for(int q=0;q<21;++q) }
最小二乘法資料擬合演算法設計,用c或者c++程式設計出來,急!!!
6樓:明天威威去旅行
#include
#include
void guass(int n,float array[10])for(i=0;i }int main() guass(n,array); return 0;} 最小二乘法一次多項式擬合又稱什麼擬合 7樓:小樂笑了 以及線性迴歸 求用c++實現最小二乘法曲線擬合** 8樓:匿名使用者 你先把最少二乘法的原理搞懂吧,其實最少二乘法不難得。 你可以自己推到,不過這些經典的方法,都有現成公式的了。拿本書找出公式,也就是幾個矩陣的運算,把你的資料帶進去就直接算得出了。 9樓:匿名使用者 用matlab倒是編過 c就沒用過了 最小二乘法擬合多項式問題 10樓:匿名使用者 w=∑(1->n) |yi-(a0+a1*xi+a2*xi^2)/(b0+b1*xi+b2*xi^2+b3*xi^3)|^2 w對各系數求偏導為0 得到7個方程,聯立求解各系數 最小二乘法是一種數學優化技術 它通過最小化誤差的平方和尋找資料的最佳函式匹配。最小二乘法 又稱最小平方法 是一種數學優化技術。它通過最小化誤差的平方和尋找資料的最佳函式匹配。利用最小二乘法可以簡便地求得未知的資料,並使得這些求得的資料與實際資料之間誤差的平方和為最小。最小二乘法還可用於曲線擬合。其他... 給你段程式,你自己照著改吧,需要給初值。function untitled3 clcxi 0.019,0.023,0.027 x data yi 430,380,256 y dataabc0 4,0.5,5 initial valuesabc lsqcurvefit fun,abc0,xi,yi t... 用正交多項式作最小二乘曲線擬合 1 函式語句與形參說明 void spir int n,int m,double x,double y,double a,double dt int n 給定資料點的個數 int m 擬合多項式的項數,即擬合多項式的次數為m 1。要求m n且m 20 double x...最小二乘法是什麼,什麼叫最小二乘法
最小二乘法擬合(y a e bx形式)
以正交多項式為基底,做最小二乘擬合有大神會嗎求