当前位置:首页 » 编程语言 » c语言微分求导
扩展阅读
webinf下怎么引入js 2023-08-31 21:54:13
堡垒机怎么打开web 2023-08-31 21:54:11

c语言微分求导

发布时间: 2022-09-08 10:29:15

c语言求导问题

用差分计算,当自变量趋于0时,前后两次差分收敛到需要精度,计算结束。
例如,一阶导数,写一个 函数 y = f(x):

float f(float x)

设 dx 初值
计算 dy
dy = f(x0) - f(x0+dx);
导数 初值
dd1=dy/dx;
Lab:;
dx = 0.5 * dx; // 减小步长
dy = f(x0) - f(x0+dx);
dd2=dy/dx; // 导数 新值
判断新旧导数值之差是否满足精度,满足则得结果,不满足则返回
if ( fabs(dd1-dd2) < 1e-06 )
else ;

② C语言编程求微分,愿意再追加50分.

#include <stdio.h>

main()
{
float x,y,ji_fen,wei_fen;
for(x=0;x<2;x=x+0.01)
{
y=x*x;
ji_fen=ji_fen+y*0.01;
wei_fen=(y-(x-0.01)*(x-0.01))/0.01;
printf("x=%f y=%f ji_fen=%f wei_fen=%f\n",x,y,ji_fen,wei_fen);
}
}

③ 用C语言如何编写函数的求导

求导数有两种,一种是表达式求导,一种是数值求导。
1.
表达式求导:需要对表达式进行词法分析,然后用常见的求导公式进行演算,求得导函数。在这方面,数学软件matrix,maple做得非常好。如果自己用C进行编程,不建议。
2.
数值求导:利用导数的定义,用差分计算,当自变量趋于0时,前后两次差分收敛到需要精度,计算结束。这种方法可以求得某一点的导数。
例如:
求一阶导数,原函数
y
=
f(x),
程序中是float
f(float
x){
...}
dx=0.01; //设 dx 初值
do{
dd1=(f(x0) - f(x0+dx))/dx; //计算导数dd1
dx = 0.5 * dx; // 减小步长
dd2=(f(x0) - f(x0+dx))/dx; //计算导数dd2
}while (fabs(dd1-dd2) >= 1e-06) //判断新旧导数值之差是否满足精度,满足则得结果,不满足则返回

④ 用C语言计算微积分

我给一楼加的注释以及修改:
#include<stdio.h>
#include<math.h>
#define ARRAYBOUND 10001

void main()
{
int i = 0; //辅助变量,最常见那种
int n = 0; //将所求定积分函数曲线在x轴方向,平均分成n等分;n越大,结果越精确;不过限于此算法限制n<ARRAYBOUND,否则溢出.
float x[ARRAYBOUND];//ARRAYBOUND维浮点数组,存放离散的x坐标值
float y[ARRAYBOUND];//ARRAYBOUND维浮点数组,存放每个x坐标对应的函数值;x[i],y[i]满足y[i]=f(x[i]),f是你要求定积分的函数
float x0 = 0.0; //定积分下限
float xn = 0.0; //定积分上限
float h = 0.0; //面积微元宽度
float J = 0.0; //辅助变量

/*f=x^3*/ //这里说明要求定积分的是函数f(x)=x*x*x;(y等于x的立方,x^3是vb的写法)
// printf("input x0,xn,n:");
printf("请分别输入下限(x0),上限(xn),精度(n):");
scanf("%f",&x0);
scanf("%f",&xn);
scanf("%d",&n);

h=(xn-x0)/n;//将函数图形在x方向平分成n份,h是每个面积微元的宽度
x[0]=x0; //将积分下限赋值给x[0]

for(i=0;i<=n && n<ARRAYBOUND;i++)
{
x[i]=x[0]+i*h; //计算n个离散的横坐标值,存入x[]数组
y[i]=(float)pow(x[i],3);//计算n个横坐标对应的函数值,存入y[]数组。在此可以改变要求积分的函数
}

// J=0.0;
for(i=0;i<n;i++)
{
//J=J+y[i]+y[i+1];
J+=y[i];//将所有纵坐标值代数相加,存入J
}
//J=J*h/2.0;
J=J*h;//所有微元面积一次求解,因为∑h*y[i]=h*∑y[i];
printf("\nn=%d \n所求定积分值是: %f\n",n,J);
}

我将//J=J+y[i]+y[i+1]改为J+=y[i];将//J=J*h/2.0;改为J=J*h只是帮助lz理解
其实,这两种表达在理论上是等价的,不过我发现修改后,在n同样大小的情况下,结果的精度有一点点下降,还真不知为什么???

这样的话lz应该能理解了吧,其实一楼的算法还有不少值得改进的地方,希望lz能有所突破!!

⑤ 如何用c语言求函数导数

1、首先要有函数,设置成double类型的参数和返回值。

2、然后根据导数的定义求出导数,参数差值要达到精度极限,这是最关键的一步。

3、假如函数是double fun(doube x),那么导数的输出应该是(fun(x)-fun(x-e))/e,这里e是设置的无穷小的变量。

4、C由于精度有限,因此需要循环反复测试,并判断无穷小e等于0之前,求出上述导数的值。二级导数也是一样,所不同的是要把上述导数公式按定义再一次求导。这是算法,具体的实现自己尝试编程。

C语言的数据长度和精度都有限,因此用C语言编程求的导数并不精确,换句话说C语言编程不适合求导和极限。

(5)c语言微分求导扩展阅读:

举例说明:

一阶导数,写一个函数 y = f(x):

float f(float x){ ...}

设 dx 初值

计算 dy

dy = f(x0) - f(x0+dx);

导数 初值

dd1=dy/dx;

Lab:;

dx = 0.5 * dx; // 减小步长

dy = f(x0) - f(x0+dx);

dd2=dy/dx; // 导数 新值

判断新旧导数值之差是否满足精度,满足则得结果,不满足则返回

if ( fabs(dd1-dd2) < 1e-06 ) { 得结果dd2...}

else { dd1=dd2;goto Lab;}。

⑥ c语言如何做微分

你连题都没有,想要源代码,估计希望不大了,不过有本书里有将如何用C语言解微分方程。
好像是叫计算方法 C语言版。哦,找到了
《计算方法(C语言版)》是作者十多年计算方法研究应用和教学经验的结晶。全书共分9章,主要内容包括算法与误差、非线性方程求根、线性方程组的直接求解和迭代求解、代数插值、数值积分、矩阵特征值与特征向量的计算、常微分方程初值问题的数值解法等。

⑦ 想写个数值微分的C语言程序,用中心差商求一阶导数,初学C语言,哪位大神或者学长学姐能帮忙看看,谢谢

typedefdouble(*FUN_TYPE)(double);
doublederive(doublex0,doubleh,FUN_TYPEfun)
{
return(fun(x0+h)-fun(x0-h))/2.0/h;
}
intmain()
{
doubled=1.15;
doubledd=derive(d,0.01,exp);
cout<<dd<<endl;
}

⑧ c语言怎样用C语言编程来解决数学求“微分

1,一套皆不能初等函数的微分公式;
2,函数和差积商的法则
通过以上两个的明确定义,就可以把一些问题转化成固定的模版上进行计算了。

⑨ 用c语言如何求导

用差分计算,当自变量趋于0时,前后两次差分收敛到需要精度,计算结束。
例如,一阶导数,写一个函数y=f(x):
floatf(floatx){...}
设dx初值
计算dy
dy=f(x0)-f(x0+dx);
导数初值
dd1=dy/dx;
Lab:;
dx=0.5*dx;//减小步长
dy=f(x0)-f(x0+dx);
dd2=dy/dx;//导数新值
判断新旧导数值之差是否满足精度,满足则得结果,不满足则返回
if(fabs(dd1-dd2)<1e-06){得结果dd2...}
else{dd1=dd2;gotoLab;};

⑩ c语言怎么编求导

//多项式求导数

intPolyDeri(list<nodePoly>&polyFunc)

{

list<nodePoly>::iteratoriter;

for(iter=polyFunc.begin();iter!=polyFunc.end();++iter)

{

if((*iter).ex>1)

{

(*iter).coef=((*iter).coef)*((*iter).ex);

(*iter).ex=(*iter).ex-1;

}

elseif(1==(*iter).ex)

{

(*iter).ex=0;

}

elseif(0==(*iter).ex)

{

(*iter).coef=0;

}

}

returnRET_OK;

}

其中,多项式的定义是list<nodePoly>,如下:

//多项式节点结构体定义

typedefstructstuPolynomNode

{

doublecoef;

intex;

}nodePoly;

(10)c语言微分求导扩展阅读

c语言求导数据范围及提示DataSize&Hint

#include<iostream>

#include<cmath>

usingnamespacestd;

intmain()

{

intnum=0,i=0;

cin>>num;

for(i=2;i<=sqrt(num);i++)

{

if(num%i==0)

break;

}

if(i>sqrt(num)

cout<<num<<"为素数"<<endl;

else

cout<<num<<"不是素数"endl;

return0;

}