C语言实现变步长求积分算法

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了C语言实现变步长求积分算法，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含1304字，纯文字阅读大概需要2分钟。

内容图文

链接：https://blog.csdn.net/chen_dsir/article/details/71257631  博主详细介绍了复合梯形公式，

链接：https://wenku.baidu.com/view/b300c53deefdc8d376ee3228.html 详细介绍了变步长的目的，以及变换的方式

下面求解如下问题：

1、变步长：

        由于多数情况下，并不知道步长设置为多少合适，所以我们可以设置一个初始步长，求解一次，再把步长变为原来的一半，依次接下去，直到达到我们所需的精度。

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;  //初始化步数
	double stepLength=(right-left)/step; //步长
	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;  //符合梯形公式求解结果
	printf("初始步长得到的积分值:%lf   所分的区间个数：%d\n",T,step);
	//下面改变步长，来细化精度
	double T1=T+100;
	double T2=T;
	double H;      //减小步长后得到的积分值 
	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;    //将上次一的积分结果和减小步长后的积分结果求均值 
		
	}	
	printf("初始步长得到的积分值:%lf   所分的区间个数：%d\n",T2,step);
	
	 
	
}

3、结果：

4、验证：

 

内容总结

以上是互联网集市为您收集整理的C语言实现变步长求积分算法全部内容，希望文章能够帮你解决C语言实现变步长求积分算法所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端

---