求曲线的四阶拟合

#include"stdafx.h"
#include <iostream>
#include <vector>
#include <cmath>
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
using namespace std;
void polyfit(int n, double x[], double y[], int poly_n, double a[]);
 int main()
{
int i, n = 11, poly_n = 4;
double x[11] = { 1,2,3,4,6,7,8,9,10,11,25 }, y[11] = { 2,3,6,7,5,3,2,25,12,34,16 };
double a[5];


system("cls");
polyfit(n, x, y, poly_n, a);

for (i = 0; i<poly_n + 1; i++)/*这里是升序排列，Matlab是降序排列*/
printf("a[%d]=%g\n", i, a[i]);
system("pause");
return 0;
}


/*==================polyfit(n,x,y,poly_n,a)===================*/
/*=======拟合y=a0+a1*x+a2*x^2+……+apoly_n*x^poly_n========*/
/*=====n是数据个数 xy是数据值 poly_n是多项式的项数======*/