用C語(yǔ)言如何編寫函數(shù)的求導(dǎo)

求導(dǎo)數(shù)有兩種，一種是表達(dá)式求導(dǎo)，一種是數(shù)值求導(dǎo)。

成都創(chuàng)新互聯(lián)是一家專業(yè)提供臨澤企業(yè)網(wǎng)站建設(shè),專注與網(wǎng)站設(shè)計(jì)、成都做網(wǎng)站、HTML5、小程序制作等業(yè)務(wù)。10年已為臨澤眾多企業(yè)、政府機(jī)構(gòu)等服務(wù)。創(chuàng)新互聯(lián)專業(yè)網(wǎng)站制作公司優(yōu)惠進(jìn)行中。

1.

表達(dá)式求導(dǎo)：需要對(duì)表達(dá)式進(jìn)行詞法分析，然后用常見的求導(dǎo)公式進(jìn)行演算，求得導(dǎo)函數(shù)。在這方面，數(shù)學(xué)軟件matrix，maple做得非常好。如果自己用C進(jìn)行編程，不建議。

2.

數(shù)值求導(dǎo)：利用導(dǎo)數(shù)的定義，用差分計(jì)算，當(dāng)自變量趨于0時(shí)，前后兩次差分收斂到需要精度，計(jì)算結(jié)束。這種方法可以求得某一點(diǎn)的導(dǎo)數(shù)。

例如：

求一階導(dǎo)數(shù)，原函數(shù)

y

=

f(x)，

程序中是float

f(float

x){

...}

dx=0.01; //設(shè) dx 初值

do{

dd1=(f(x0) - f(x0+dx))/dx; //計(jì)算導(dǎo)數(shù)dd1

dx = 0.5 * dx; // 減小步長(zhǎng)

dd2=(f(x0) - f(x0+dx))/dx; //計(jì)算導(dǎo)數(shù)dd2

}while (fabs(dd1-dd2) = 1e-06) //判斷新舊導(dǎo)數(shù)值之差是否滿足精度，滿足則得結(jié)果，不滿足則返回

C語(yǔ)言如何求導(dǎo)函數(shù)

用差分計(jì)算，當(dāng)自變量趨于0時(shí)，前后兩次差分收斂到需要精度，計(jì)算結(jié)束。

例如，一階導(dǎo)數(shù)，寫一個(gè) 函數(shù) y = f(x)：

float f(float x){ ...}

設(shè) dx 初值

計(jì)算 dy

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

導(dǎo)數(shù) 初值

dd1=dy/dx;

Lab:；

dx = 0.5 * dx; // 減小步長(zhǎng)

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

dd2=dy/dx; // 導(dǎo)數(shù) 新值

判斷新舊導(dǎo)數(shù)值之差是否滿足精度，滿足則得結(jié)果，不滿足則返回

if ( fabs(dd1-dd2) 1e-06 ) { 得結(jié)果dd2...}

else { dd1=dd2;goto Lab;};

c語(yǔ)言怎么編求導(dǎo)

//多項(xiàng)式求導(dǎo)數(shù)

intPolyDeri(listnodePolypolyFunc)

{

listnodePoly::iteratoriter;

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

{

if((*iter).ex1)

{

(*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;

}

其中，多項(xiàng)式的定義是listnodePoly，如下：

//多項(xiàng)式節(jié)點(diǎn)結(jié)構(gòu)體定義

typedefstructstuPolynomNode

{

doublecoef;

intex;

}nodePoly;

擴(kuò)展資料

c語(yǔ)言求導(dǎo)數(shù)據(jù)范圍及提示DataSizeHint

#includeiostream

#includecmath

usingnamespacestd;

intmain()

{

intnum=0,i=0;

cinnum;

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

{

if(num%i==0)

break;

}

if(isqrt(num)

coutnum"為素?cái)?shù)"endl;

else

coutnum"不是素?cái)?shù)"endl;

return0;

}