c语言拟合
‘壹’ c语言曲线拟合问题
很简单 问题已经转化成 从文件读数据 数据类型转化 赋值到数组 这三步
#include"stdio.h"
#include<stdlib.h>
main()
{
intx[3]={0};//默认文件存了三个数据
FILE*fp1;//定义文件流指针，用于打开读取的文件
chartext[1024];//定义一个字符串数组，用于存储读取的字符
fp1=fopen("d:\forecast1.txt","r");//只读方式打开文件a.txt
inti=0;
while(fgets(text,1024,fp1)!=NULL)//逐行读取fp1所指向文件中的内容到text中
{
inttmp=atoi(text);//输出到屏幕
x[i++]=tmp;
}
fclose(fp1);//关闭文件a.txt，有打开就要有关闭
}
y[i]同理
‘贰’ 如何用c语言编写32阶多项式拟合函数
这与几阶多项式无关，关键是用什么方法拟合。如果只有33个数据点，解一族线性方程就可以了，如果有很多组数据，就用最小而乘法，步骤略微多一些。看看计算方法的书就任意操作了。
‘叁’ 如何用c语言拟合二次方程y=a(x-x0)(x-x0) + b(x-x0) + c,根据多组(x,y)的值计算出参数a,b,c,x0的值
其实很简单的。。就是个搜索穷解问题哎。。
void SOLVE(float x1,float y1,float x2,float y2,float x3,float y3,float x4,float y4)
{
for（a=-10000;a<10000;a++）
for(b=-10000;b<10000;b++)
for(c=-10000;c<10000;c++)
for(x0=-10000;x0<10000;x0++)
{
if(y1=a(x1-x0)(x1-x0) + b(x1-x0) + c)
if(y2=a(x2-x0)(x2-x0) + b(x2-x0) + c)
if(y3=a(x3-x0)(x3-x0) + b(x3-x0) + c)
if(y4=a(x4-x0)(x4-x0) + b(x4-x0) + c)
{
cout<<"a="<<a<<","<<"b="<<b<<","<<"c="<<c<<","<<"x0="<<x0<<endl;
};
};
}
‘肆’ c语言如何对若干离散点进行m阶多项式拟合
m阶多项式拟合，是算法问题，不是计算机语言问题。你先要有你的具体计算策略和方法，数值计算时，用计算机语言写出程序，算出结果。用 c 语言 或 fortran 语言 或别的语言，大同小异，没有原则差别。
你先要定出你的方程形式，例如
y = a1*x^m+a2*x^(m-1)+...am
a1,a2,...,am 是 m 个待定系数
把你的离散点[xi,yi] 代入方程，你可以得到 线性方程式 yi = a1*xi^m+a2*xi^(m-1)+...am
若你有 n 个点，就得到 n 个 方程式。若 n<m 有无穷解，若 n=m 有一解。
通常 离散点的个数 远超过 方程阶数，也就是方程个数超出要解的未知数的个数，n>m，这时，
通常 用 最小二乘法 求解 这个线性方程组。也就是所谓的拟合。
最小二乘法 求解 这个线性方程组 的程序 网上（这里不敢写，写了，你就看不到我的贴子了）可以找到，自己写也不复杂。
曲线拟合文章估计网络文库里也有。
‘伍’ c语言怎么把数据拟合成函数并画线
1、通过一些点拟合出一条直线。
2、参数：pt_input指向传入的点的指针。
3、ptNumbers传入的点数量。
4、k指向拟合直线参数k的指针。
5、b指向拟合直线参数b的指针。
‘陆’ C语言线性曲线拟合
可以直接在线进行拟合，下面是地址（已验证）
http://3.14159.biz/math/Fit.html
建议你直接对数据用Matlab的cftool拟合(可以自由的自定义拟合函数形式)
‘柒’ c语言曲线拟合问题。已经有源代码， 求高手修改符合要求！！！！
以你列举的数据为例，一次运算是哪些数据？
输入拟合多项式的次数:
5
输入给定点的个数n及坐标(x,y):
5
?,?
?,?
?,?
?,?
?,?
你给出了10行数据？能做几次运算？
结果为什么是4个文件？那些数据写第一个文件里去？
我看你的描述，反正得不到这些信息。
‘捌’ 想用C语言编写多项式拟合的程序
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
main()
{
int i,j,m,n=7,poly_n=2;
double x[7]={1,2,3,4,6,7,8},y[7]={2,3,6,7,5,3,2};
double a[3];
void polyfit(int n,double *x,double *y,int poly_n,double a[]);
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]);
getch();
}
/*==================polyfit(n,x,y,poly_n,a)===================*/
/*=======拟合y=a0+a1*x+a2*x^2+……+apoly_n*x^poly_n========*/
/*=====n是数据个数 xy是数据值 poly_n是多项式的项数======*/
/*===返回a0,a1,a2,……a[poly_n]，系数比项数多一（常数项）=====*/
void polyfit(int n,double x[],double y[],int poly_n,double a[])
{
int i,j;
double *tempx,*tempy,*sumxx,*sumxy,*ata;
void gauss_solve(int n,double A[],double x[],double b[]);
tempx=calloc(n,sizeof(double));
sumxx=calloc(poly_n*2+1,sizeof(double));
tempy=calloc(n,sizeof(double));
sumxy=calloc(poly_n+1,sizeof(double));
ata=calloc((poly_n+1)*(poly_n+1),sizeof(double));
for (i=0;i<n;i++)
{
tempx[i]=1;
tempy[i]=y[i];
}
for (i=0;i<2*poly_n+1;i++)
for (sumxx[i]=0,j=0;j<n;j++)
{
sumxx[i]+=tempx[j];
tempx[j]*=x[j];
}
for (i=0;i<poly_n+1;i++)
for (sumxy[i]=0,j=0;j<n;j++)
{
sumxy[i]+=tempy[j];
tempy[j]*=x[j];
}
for (i=0;i<poly_n+1;i++)
for (j=0;j<poly_n+1;j++)
ata[i*(poly_n+1)+j]=sumxx[i+j];
gauss_solve(poly_n+1,ata,a,sumxy);
free(tempx);
free(sumxx);
free(tempy);
free(sumxy);
free(ata);
}
void gauss_solve(int n,double A[],double x[],double b[])
{
int i,j,k,r;
double max;
for (k=0;k<n-1;k++)
{
max=fabs(A[k*n+k]); /*find maxmum*/
r=k;
for (i=k+1;i<n-1;i++)
if (max<fabs(A[i*n+i]))
{
max=fabs(A[i*n+i]);
r=i;
}
if (r!=k)
for (i=0;i<n;i++) /*change array:A[k]&A[r] */
{
max=A[k*n+i];
A[k*n+i]=A[r*n+i];
A[r*n+i]=max;
}
max=b[k]; /*change array:b[k]&b[r] */
b[k]=b[r];
b[r]=max;
for (i=k+1;i<n;i++)
{
for (j=k+1;j<n;j++)
A[i*n+j]-=A[i*n+k]*A[k*n+j]/A[k*n+k];
b[i]-=A[i*n+k]*b[k]/A[k*n+k];
}
}
‘玖’ 已知散点，怎样用C语言求出斜率，线性拟合
拟合直线之后会有一个拟合结果表，在这个表里有直线的方程，斜率值(slop)，截距值(intercept)，以及这些值的误差(error)，和拟合结果的r^2（越接近于1，表明拟合结果越好）。
如果没有看到这个表，那么点击菜单栏上的view
-
results
log，就能看到了。
‘拾’ 急求！！！如何用c语言实现多元线性拟合
Dim a() As String, b() As String
Private Sub Command1_Click()
Dim temp1 As String
Dim i As Integer, k As Integer
temp1 = Text1.Text
a = Split(temp1, ",")
lenolds = Len("->")
For i = 0 To UBound(a)
s = a(i)
j = InStr(s, "->")
Do While j > 0
ReDim Preserve b(i, k)
b(i, k) = Val(Left(s, j + lenolds))
s = Left(s, j - 1) + "->" + Mid(s, j + lenolds + 1)
k = k + 1
j = InStr(s, "->")
Loop
ReDim Preserve b(i, k)
b(i, k) = s
Next i
For i = 0 To UBound(a)
For j = 0 To k
Text2.Text = Text2.Text + b(i, j)
Next j
Next i
End Sub