用C語(yǔ)言編寫一個(gè)求定積分的程序，用它分別求以下5個(gè)函數(shù)的定積分：

#includestdio.h

創(chuàng)新互聯(lián)建站從2013年成立，是專業(yè)互聯(lián)網(wǎng)技術(shù)服務(wù)公司，擁有項(xiàng)目做網(wǎng)站、成都網(wǎng)站建設(shè)網(wǎng)站策劃，項(xiàng)目實(shí)施與項(xiàng)目整合能力。我們以讓每一個(gè)夢(mèng)想脫穎而出為使命，1280元南縣做網(wǎng)站,已為上家服務(wù),為南縣各地企業(yè)和個(gè)人服務(wù),聯(lián)系電話:18982081108

#includemath.h

float f1(float x)

{

return(1.0+x);

}

float f2(float x)

{

return(2.0*x+3.0);

}

float f3(float x)

{

return(exp(x)+1);

}

float f4(float x)

{

return(pow(1+x,2));

}

float f5(float x)

{

return(pow(x,3));

}

float fsimp(float a,float b,float (*p)(float))

{

float c,s;

c=(a+b)/2;

s=(b-a)/6*(p(a)+4*p(c)+p(b));

return s;

}

int main()

{

float a,b;

printf("請(qǐng)輸入積分下限a的值：");

scanf("%f",a);

printf("請(qǐng)輸入積分上限b的值：");

scanf("%f",b);

printf("%f\n",fsimp(a,b,f1));

printf("%f\n",fsimp(a,b,f2));

printf("%f\n",fsimp(a,b,f3));

printf("%f\n",fsimp(a,b,f4));

printf("%f\n",fsimp(a,b,f5));

}

用C語(yǔ)言求定積分

實(shí)際問(wèn)題描述：

求定積分近似值

程序代碼如下：

#include

#include

void main()

{

int i,n=1000;

float a,b,h,t1,t2,s1,s2,x;

printf("請(qǐng)輸入積分限a,b:");

scanf("%f,%f",a,b);

h=(b-a)/n;

for(s1=0,s2=0,i=1;i=n;i++)

{

x=a+(i-1)*h;

t1=(float)exp(-x*x/2);t2(float)=exp(-(x+h)*(x+h)/2);

s1=s1+t1*h;??????? /*矩形面積累加*/

s2=s2+(t1+t2)*h/2;??????? /*梯形面積累加*/

}

printf("矩形法算得積分值：%f.\n",s1);

printf("梯形法算得積分值：%f.\n",s2);

}

程序運(yùn)行結(jié)果如下：

矩形法算得積分值：0.855821

梯形法算得積分值：0.855624

由上面的比較可知，梯形法的精度要高于矩形法。

c語(yǔ)言 求定積分的通用函數(shù)

對(duì)于一重定積分來(lái)說(shuō)其求解可以使用梯形法進(jìn)行求解，計(jì)算公式如下所示：

其中，f(x)為被積函數(shù)，為橫坐標(biāo)的兩點(diǎn)間的間隔，越小，則計(jì)算出的結(jié)果越精確。

對(duì)于求解此類問(wèn)題可以使用C語(yǔ)言中的回調(diào)函數(shù)編寫通用的計(jì)算函數(shù)，代碼如下：

#include?stdio.h

#include?stdlib.h

#includemath.h

//功能：返回f(x)在積分區(qū)間[a,b]的值

//參數(shù)：FunCallBack?指向用于計(jì)算f(x)的函數(shù)

//??????a??積分區(qū)間的起始值

//??????b??積分區(qū)間的結(jié)束值

//??????dx?橫坐標(biāo)的間隔數(shù)，越小計(jì)算結(jié)果越準(zhǔn)確

double?Calculate(double?(*FunCallBack)(double?x),

double?a,double?b,double?dx)

{

double?doui;

double?total?=?0;????????//保存最后的計(jì)算結(jié)果

for?(doui?=?a;?doui?=?b;?doui?+=?dx)

{

total?+=?FunCallBack(doui)*dx;

}

return?total;

}

double?f2(double?x)

{

return?x*x;

}

double?f(double?x)

{

return?x;

}

double?f3(double?x)

{

return?x*x*x?;

}

int?main()

{

double?total;

total?=?(Calculate(f,?2,?3,?0.000001));

printf("total?=?%lf\n",?total);

total?=?(Calculate(f2,?2,?3,?0.000001));

printf("total?=?%lf\n",?total);

total?=?(Calculate(f3,?2,?3,?0.000001));

printf("total?=?%lf\n",?total);

return?0?;

}

其中，函數(shù)f，f2，f3為自行編寫的關(guān)于x的被積函數(shù)。

運(yùn)行結(jié)果：

total?=?2.500000

total?=?6.333331

total?=?16.249991

C語(yǔ)言求積分

float I_Control(float Input)

{

static float OutData = 0;

OutData+=Input*deltaT；

deltaT+=deltaT；

return OutData;

}

用C語(yǔ)言編寫一個(gè)求定積分的程序

這是辛普森積分法。

給你寫了fun_1( ），fun_2(），請(qǐng)自己添加另外幾個(gè)被積函數(shù)。

調(diào)用方法 t=fsimp(a,b,eps,fun_i);

a,b --上下限，eps -- 迭代精度要求。

#includestdio.h

#includestdlib.h

#include math.h

double fun_1(double x)

{

return 1.0 + x ;

}

double fun_2(double x)

{

return 2.0 * x + 3.0 ;

}

double fsimp(double a,double b,double eps, double (*P)(double))

{

int n,k;

double h,t1,t2,s1,s2,ep,p,x;

n=1; h=b-a;

t1=h*(P(a)+P(b))/2.0;

s1=t1;

ep=eps+1.0;

while (ep=eps)

{

p=0.0;

for (k=0;k=n-1;k++)

{

x=a+(k+0.5)*h;

p=p+P(x);

}

t2=(t1+h*p)/2.0;

s2=(4.0*t2-t1)/3.0;

ep=fabs(s2-s1);

t1=t2; s1=s2; n=n+n; h=h/2.0;

}

return(s2);

}

void main()

{

double a,b,eps,t;

a=0.0; b=3.141592653589793238; eps=0.0000001;

// a definite integral by Simpson Method.

t=fsimp(a,b,eps,fun_1);

printf("%g\n",t);

t=fsimp(a,b,eps,fun_2);

printf("%g\n",t);

// ...

printf("\n Press any key to quit...");

getch();

}

用C語(yǔ)言求積分

基本是這樣的，用梯形發(fā)求定積分，對(duì)應(yīng)于一個(gè)積分式就要有一段程序，不過(guò)你可以改變程序的一小部分來(lái)改變你所要求的積分式。

以c為例：求f(x)=xsinx從1到2的積分

#include math.h

float integral(float(*fun)(float x),float a,float b,int,n)

{float s,h,y;

int i;

s=(fun(a)+fun(b))/2;

h=(b-a)/n; /*積分步長(zhǎng)*/

for(i=1;in;i++)

s=s+fun(a+i*h);

y=s*h;

return y;/*返回積分值*/

}

float f(float x)

{return(x*sinx) /*修改此處可以改變被積函數(shù)*/

}

main()

{float y;

y=integral(f,1.0,2.0,150);/*修改此處可以改變積分上下限和步長(zhǎng)*/

printf("y=%f\n",y);

}

本文標(biāo)題：c語(yǔ)言函數(shù)求積分,C語(yǔ)言積分函數(shù)
URL標(biāo)題：http://www.chinadenli.net/article36/dssoppg.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供網(wǎng)站收錄、靜態(tài)網(wǎng)站、網(wǎng)站內(nèi)鏈、App開發(fā)、手機(jī)網(wǎng)站建設(shè)、Google

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請(qǐng)盡快告知，我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng)，如需處理請(qǐng)聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時(shí)需注明來(lái)源： 創(chuàng)新互聯(lián)