#include iostream
我们提供的服务有:成都网站设计、做网站、微信公众号开发、网站优化、网站认证、富拉尔基ssl等。为1000+企事业单位解决了网站和推广的问题。提供周到的售前咨询和贴心的售后服务,是有科学管理、有技术的富拉尔基网站制作公司
#includemath.h
using namespace std;
double f(double x){//方程
return 2*x*x*x-4*x*x+3*x-6;
}
double bsearch(double l,double r,double p){ //二分求近似值
double ans = (l+r)/2;
while(fabs(f(ans))p){
if(f(ans)*f(l)0){
l = ans;
}else{
r = ans;
}
ans = (l+r)/2;
}
return ans;
}
int main(int argc, char** argv) {
coutbsearch(0,3,0.0001)endl;
return 0;
}
你的意思是说,输入一些数据,让程序拟合出来一个函数。
这个很复杂。比如我举个例子。
多项式拟合函数polyfit之C语言的源码
2010-11-25 10:05:57| 分类: VC++编程 | 标签:double poly_n polyfit int for |字号 订阅
#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;ipoly_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;in;i++)
{
tempx[i]=1;
tempy[i]=y[i];
}
for (i=0;i2*poly_n+1;i++)
for (sumxx[i]=0,j=0;jn;j++)
{
sumxx[i]+=tempx[j];
tempx[j]*=x[j];
}
for (i=0;ipoly_n+1;i++)
for (sumxy[i]=0,j=0;jn;j++)
{
sumxy[i]+=tempy[j];
tempy[j]*=x[j];
}
for (i=0;ipoly_n+1;i++)
for (j=0;jpoly_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;kn-1;k++)
{
max=fabs(A[k*n+k]); /*find maxmum*/
r=k;
for (i=k+1;in-1;i++)
if (maxfabs(A[i*n+i]))
{
max=fabs(A[i*n+i]);
r=i;
}
if (r!=k)
for (i=0;in;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;in;i++)
{
for (j=k+1;jn;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];
}
}
所以这种程序是超级大型的,一般没人能独立写出来。
编程计算器和一些程序能实现,但是源代码都是不公布的。
#include stdio.h
int main()
{
char * file="src.txt";
FILE * fp;
fp=fopen(file,"wt");//打开文件
fwrite("ok",1,2,fp); //写文件
fclose(fp);//关闭文件
return 0;
}
你只要多看看书,深入一点的,才会发现C的功能强大
首先输入要计算什么
比如
sin
cos
...
然后输入要计算的值
接着调用对应的数学函数就可以了
sin
con
tan
cot这些都是有对应数学函数的
最后输出结果。
需要注意的是
C的数学三角函数都是弧度做参数
而不是角度。
提供两个连接给你,自己去下载:
第一个:
第二个:
你所指的软件应该是应用软件,用字符图形作为软件的界面,用C和算法,与数据结构作软件的核心代码,用文件作为数据库(C也可以连接数据库,不过我不会用),做出来的软件界面不好看,但是运行的速度是非常快的,如果你的核心代码写的好,那么这个软件的效率和质量也是非常高的,所以纯C作出来的软件是很优秀的,但是因为开发时间长,开发过程复杂,所以用纯C做应用软件很少
还有C语言有很大的函数库,数学函数只是其函数库中的一部分而以!!