数值分析C++实现Romberg积分算法

2023-10-18 01:20

本文主要是介绍数值分析C++实现Romberg积分算法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

在这里插入图片描述
算法描述
其算法的伪代码为:在这里插入图片描述
具体的算法:
(1)定义一个二维数组R[n][n]用于存储每次积分的值
(2)首先计算 R[1][1]=(b-a)*(f(a)+f(b))/2;
(3)从j =2开始遍历做以下步骤,直到j=n停止遍历
3.1计算步长h =(b-a)/2^(j-1),计算R[j][1]=
3.2从i =1开始,每次i自增1直到i = 2^(j-2),计算 f(a+(2i-1)h)的累加和sum1。
3.3 计算R[j][1]=R[j-1][1]/2 + h
sum1;
3.4从k =2开始每次k自增1,遍历重复以下操作:
3.4.1 计算R[j][k] = (4^(k-1) *R[j][k-1]-R[j-1][k-1])/(4^(k-1)-1)
(4)打印输出R[j][k]的积分值
源程序代码及运行结果截图

#include<iostream>
#include<vector>
#include<math.h>
#include<iomanip>
#include<string>
using namespace std;
#define MAX  100
float  R[MAX][MAX];
/*
@to do:  获取自变量x对应的函数值
@param x 
*/
float getRombergFuntion(float x) {return log(x);
}
/*
@to do:自定义实现x的y的整数次幂的结果
@param x: x
@param y: y
@return 返回x的y次幂的结果int
*/
int pow(int x,int y)
{int result = 1;for (int  i = 0; i < y; i++){result = result * x;}return result;
}
/*
@to do:获取荣贝格算法的n行打印数据
@param a:积分下限
@param b:积分上限
@param n:要打印的行数
*/
void getRombergValue(float a,float b,int n) 
{//获取初始值R[1][1] = (b-a)*(getRombergFuntion(a)+ getRombergFuntion(b))/2;//步长hfloat h = 0.0f;for (int  i = 2; i <= n; i++){//计算当前步长h = (b - a) / pow(2,i-1);float sum = 0.0f;for (int j = 1; j <= pow(2, i- 2);j++) {sum += getRombergFuntion(a+(2*j -1)*h)*h;}R[i][1] = 0.5*R[i - 1][1] + sum;for (int  k = 2; k <= i; k++){R[i][k] = (pow(4, k - 1)*R[i][k-1]-R[i-1][k-1])/(pow(4,k-1)-1);}}for (int  i = 1; i <= n; i++){for (int j = 1; j <= n; j++){cout << R[i][j] << " ";}cout << endl;}
}int main()
{//测试数据求lnx在区间[1,2]上的积分,并输出四行lnx积分的数据getRombergValue(1,2,4);system("pause");return 0;
}

在这里插入图片描述

这篇关于数值分析C++实现Romberg积分算法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/229252

相关文章

使用Python实现矢量路径的压缩、解压与可视化

《使用Python实现矢量路径的压缩、解压与可视化》在图形设计和Web开发中,矢量路径数据的高效存储与传输至关重要,本文将通过一个Python示例,展示如何将复杂的矢量路径命令序列压缩为JSON格式,... 目录引言核心功能概述1. 路径命令解析2. 路径数据压缩3. 路径数据解压4. 可视化代码实现详解1

PyQt6/PySide6中QTableView类的实现

《PyQt6/PySide6中QTableView类的实现》本文主要介绍了PyQt6/PySide6中QTableView类的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学... 目录1. 基本概念2. 创建 QTableView 实例3. QTableView 的常用属性和方法

PyQt6/PySide6中QTreeView类的实现

《PyQt6/PySide6中QTreeView类的实现》QTreeView是PyQt6或PySide6库中用于显示分层数据的控件,本文主要介绍了PyQt6/PySide6中QTreeView类的实现... 目录1. 基本概念2. 创建 QTreeView 实例3. QTreeView 的常用属性和方法属性

Android使用ImageView.ScaleType实现图片的缩放与裁剪功能

《Android使用ImageView.ScaleType实现图片的缩放与裁剪功能》ImageView是最常用的控件之一,它用于展示各种类型的图片,为了能够根据需求调整图片的显示效果,Android提... 目录什么是 ImageView.ScaleType?FIT_XYFIT_STARTFIT_CENTE

pandas中位数填充空值的实现示例

《pandas中位数填充空值的实现示例》中位数填充是一种简单而有效的方法,用于填充数据集中缺失的值,本文就来介绍一下pandas中位数填充空值的实现,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是中位数填充?为什么选择中位数填充?示例数据结果分析完整代码总结在数据分析和机器学习过程中,处理缺失数

Golang HashMap实现原理解析

《GolangHashMap实现原理解析》HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持高效的插入、查找和删除操作,:本文主要介绍GolangH... 目录HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持

Pandas使用AdaBoost进行分类的实现

《Pandas使用AdaBoost进行分类的实现》Pandas和AdaBoost分类算法,可以高效地进行数据预处理和分类任务,本文主要介绍了Pandas使用AdaBoost进行分类的实现,具有一定的参... 目录什么是 AdaBoost?使用 AdaBoost 的步骤安装必要的库步骤一:数据准备步骤二:模型

Spring Boot中JSON数值溢出问题从报错到优雅解决办法

《SpringBoot中JSON数值溢出问题从报错到优雅解决办法》:本文主要介绍SpringBoot中JSON数值溢出问题从报错到优雅的解决办法,通过修改字段类型为Long、添加全局异常处理和... 目录一、问题背景:为什么我的接口突然报错了?二、为什么会发生这个错误?1. Java 数据类型的“容量”限制

使用Pandas进行均值填充的实现

《使用Pandas进行均值填充的实现》缺失数据(NaN值)是一个常见的问题,我们可以通过多种方法来处理缺失数据,其中一种常用的方法是均值填充,本文主要介绍了使用Pandas进行均值填充的实现,感兴趣的... 目录什么是均值填充?为什么选择均值填充?均值填充的步骤实际代码示例总结在数据分析和处理过程中,缺失数

Java对象转换的实现方式汇总

《Java对象转换的实现方式汇总》:本文主要介绍Java对象转换的多种实现方式,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录Java对象转换的多种实现方式1. 手动映射(Manual Mapping)2. Builder模式3. 工具类辅助映