数值计算 - Richardson外推法求一阶导数（C++实现）

最新推荐文章于 2024-08-02 19:55:53 发布

c橘子

最新推荐文章于 2024-08-02 19:55:53 发布

阅读量4.4k

点赞数 2

CC 4.0 BY-SA版权

分类专栏：【算法】文章标签：数值计算外推法 c++

本文链接：https://blog.youkuaiyun.com/qq_40344307/article/details/90733474

本文介绍了使用Richardson外推法来求解一阶导数的方法，并提供了C++实现代码，适用于数值计算场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

求导只需将 f 中的函数替换即可。

/*本程序使用理查森外推方法就求微分，程序的核心是运用G[i][j] =
 (pow(4, i) * G[i][j-1] - G[i-1][j-1])/(pow(4, i)-1)*/
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
const double le = 1e-10;

////////////////////////////////////////////////////////////////////////////////////
template <typename T> T** Allocation2D(int m, int n){
	T **a;
	a = new T*[m];
	for(int i=0; i<m; i++)
		a[i] = new T[n];
	return a;
}

////////////////////////////////////////////////////////////////////////////////////
double f(double x);
void WaiTui(double t, double h);
void Clear(double **G, int m);
////////////////////////////////////////////////////////////////////////////////////
int main(){
	double t, h;
	cout<<"请输入所求值的导数：";
	cin>>t;
	cout<<"请输入步长的值：";
	cin>>h;
	WaiTui(t, h);
	return 0;
}

////////////////////////////////////////////////////////////////////////////////////
/********************************函数f = sin(x)