求導(dǎo)數(shù)有兩種,一種是表達式求導(dǎo),一種是數(shù)值求導(dǎo)。
創(chuàng)新互聯(lián)專注于興文企業(yè)網(wǎng)站建設(shè),響應(yīng)式網(wǎng)站建設(shè),商城網(wǎng)站開發(fā)。興文網(wǎng)站建設(shè)公司,為興文等地區(qū)提供建站服務(wù)。全流程按需設(shè)計網(wǎng)站,專業(yè)設(shè)計,全程項目跟蹤,創(chuàng)新互聯(lián)專業(yè)和態(tài)度為您提供的服務(wù)
表達式求導(dǎo):需要對表達式進行詞法分析,然后用常見的求導(dǎo)公式進行演算,求得導(dǎo)函數(shù)。在這方面,數(shù)學(xué)軟件matrix,maple做得非常好。如果自己用C進行編程,不建議。
數(shù)值求導(dǎo):利用導(dǎo)數(shù)的定義,用差分計算,當(dāng)自變量趨于0時,前后兩次差分收斂到需要精度,計算結(jié)束。這種方法可以求得某一點的導(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;????//計算導(dǎo)數(shù)dd1
dx?=?0.5?*?dx;??//?減小步長
dd2=(f(x0)?-?f(x0+dx))/dx;????//計算導(dǎo)數(shù)dd2
}while?(fabs(dd1-dd2)?=?1e-06)?//判斷新舊導(dǎo)數(shù)值之差是否滿足精度,滿足則得結(jié)果,不滿足則返回
1、首先要有函數(shù),設(shè)置成double類型的參數(shù)和返回值。
2、然后根據(jù)導(dǎo)數(shù)的定義求出導(dǎo)數(shù),參數(shù)差值要達到精度極限,這是最關(guān)鍵的一步。
3、假如函數(shù)是double fun(doube x),那么導(dǎo)數(shù)的輸出應(yīng)該是(fun(x)-fun(x-e))/e,這里e是設(shè)置的無窮小的變量。
4、C由于精度有限,因此需要循環(huán)反復(fù)測試,并判斷無窮小e等于0之前,求出上述導(dǎo)數(shù)的值。二級導(dǎo)數(shù)也是一樣,所不同的是要把上述導(dǎo)數(shù)公式按定義再一次求導(dǎo)。這是算法,具體的實現(xiàn)自己嘗試編程。
C語言的數(shù)據(jù)長度和精度都有限,因此用C語言編程求的導(dǎo)數(shù)并不精確,換句話說C語言編程不適合求導(dǎo)和極限。
擴展資料:
舉例說明:
一階導(dǎo)數(shù),寫一個函數(shù) y = f(x):
float f(float x){ ...}
設(shè) dx 初值
計算 dy
dy = f(x0) - f(x0+dx);
導(dǎo)數(shù) 初值
dd1=dy/dx;
Lab:;
dx = 0.5 * dx; ?// 減小步長
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;}。
//多項式求導(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;
}
其中,多項式的定義是listnodePoly,如下:
//多項式節(jié)點結(jié)構(gòu)體定義
typedefstructstuPolynomNode
{
doublecoef;
intex;
}nodePoly;
擴展資料
c語言求導(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"為素數(shù)"endl;
else
coutnum"不是素數(shù)"endl;
return0;
}
用差分計算,當(dāng)自變量趨于0時,前后兩次差分收斂到需要精度,計算結(jié)束。
例如,一階導(dǎo)數(shù),寫一個函數(shù)y=f(x):
floatf(floatx){...}
設(shè)dx初值
計算dy
dy=f(x0)-f(x0+dx);
導(dǎo)數(shù)初值
dd1=dy/dx;
Lab:;
dx=0.5*dx;//減小步長
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;gotoLab;};