鏈接：https://blog.csdn.net/chen_dsir/article/details/71257631? 博主詳細(xì)介紹了復(fù)合梯形公式，

鏈接：https://wenku.baidu.com/view/b300c53deefdc8d376ee3228.html?詳細(xì)介紹了變步長(zhǎng)的目的，以及變換的方式

下面求解如下問(wèn)題：

1、變步長(zhǎng)：

? ? ? ? 由于多數(shù)情況下，并不知道步長(zhǎng)設(shè)置為多少合適，所以我們可以設(shè)置一個(gè)初始步長(zhǎng)，求解一次，再把步長(zhǎng)變?yōu)樵瓉?lái)的一半，依次接下去，直到達(dá)到我們所需的精度。

2、代碼：

#include<stdio.h> #include<math.h>#define left 1.0 //左端的值 #define right 3.0 //右端的值 #define precision 0.0001 //定義精度 double f(double x) {double m=100/(x*x)*sin(10/x);return m; }int main() {int i;int step=100; //初始化步數(shù)double stepLength=(right-left)/step; //步長(zhǎng)double T=f(left)+f(right); //積分的值 for(i=1;i<=step-1;i++) {double x=left+i*stepLength;T+=2*f(x);}T=stepLength/2*T; //符合梯形公式求解結(jié)果printf("初始步長(zhǎng)得到的積分值:%lf 所分的區(qū)間個(gè)數(shù)：%d\n",T,step);//下面改變步長(zhǎng)，來(lái)細(xì)化精度double T1=T+100;double T2=T;double H; //減小步長(zhǎng)后得到的積分值 while(fabs(T1-T2)>precision/2) {T1=T2;step=step*2;stepLength=(right-left)/step; H=0;for(i=1;i<step;i++){H+=2*f(left+stepLength*i);}H=f(left)+f(right)+H;H=H*stepLength/2; T2=(T1+H)/2; //將上次一的積分結(jié)果和減小步長(zhǎng)后的積分結(jié)果求均值 } printf("初始步長(zhǎng)得到的積分值:%lf 所分的區(qū)間個(gè)數(shù)：%d\n",T2,step);}

3、結(jié)果：

4、驗(yàn)證：

?

總結(jié)

以上是生活随笔為你收集整理的C语言实现变步长求积分算法的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。