當前位置:首頁 » 編程語言 » 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;

}