数值微分的一阶导数和二阶导数公式在具体计算时
题目
数值微分的一阶导数和二阶导数公式在具体计算时
问个数值微分的问题,由于没学过所以想不通,我想如果学过的话应该都遇到过吧.就是对于中心差分格式,f'(x)=[f(x+h)-f(x-h)]/2h,f''(x)=[f(x+h)-2f(x)+f(x-h)]/h^2,但实际在具体已知一系列点(xi,yi),求这系列点的2阶导数值,我可以先求其一阶导数值,然后再用一阶的差分公式求出2阶的导数,也可以直接用2阶导数的差分公式求出,对这两种算法,那个精度更高呢?
答案
答:本题是算是问对人了,
如果你要想深入分析,需要用到函数的泰勒展开.
1) 你说的两种方法都可以用,但是后面的方法精度更高.
f''(x)=[f(x+h)-2f(x)+f(x-h)]/h^2 方法是等效与 f''(x)=[f'(x+h/2)-f'(x-h/2)]/h 是2阶精度
2) 先求其一阶导数值,然后再用一阶的差分公式求出2阶的导数,是1阶精度.
就好比 f'(x)=[f(x+h)-f(x-h)]/2h 是2阶精度,
f'(x)=[f(x+h)-f(x)]/h 是1阶精度.
关键就是在泰勒展开方面
f(x+h)=f(x)+f'(x)*h+b*f''(x)*h^2+c*f'''(x)*h^3 b,c为泰勒展开的系数
f(x-h)=f(x)-f'(x)*h+b*f''(x)*h^2-c*f'''(x)*h^3
可以看出 如果用 [f(x+h)-f(x)]/h =f'(x)+b*f''(x)*h+c*f'''(x)*h^2 后面的是误差项.
如果用 [f(x+h)-f(x=h)]/2h =f'(x) + c*f'''(x)*h^2 可以明显看出第二种方法的误差更小.
同理可以推导二阶导数的精度问题.
如果有不明白的,可以进一步沟通
举一反三
已知函数f(x)=x,g(x)=alnx,a∈R.若曲线y=f(x)与曲线y=g(x)相交,且在交点处有相同的切线,求a的值和该切线方程.
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
最新试题
- 一个无盖的长方体鱼缸,长30厘米,宽20厘米,高20厘米,水深8厘米,放一条鱼,水面高8.5厘米,求鱼的体积
- 已知方程x平方-2x-根号5x+2根号5=0的两个根是直角三角形的两条直角边的长,求直角三角形的斜边长
- 请用孜孜不倦 不耻下问 专心致志 博览群书
- 地球纬度是通过该点的地理纬线与赤道间的角度?
- 正态分布的概率计算
- 如图,在等边△ABC中,D、E分别是BC、AC上的点,且CD=AE,AD与BE相交于点P. (1)求证:∠ABE=∠CAD; (2)若BH⊥AD于点H,求证:PB=2PH.
- 一堆糖果,3个3个数少两个,4个4个数少三个,5个5个数少4个,问这堆糖至少有几颗?
- 当物质的量为amol的氢气和一氧化碳的混合气体完全燃烧时,消耗了0.5amol氧气.则混合气体中.氢气和一氧化碳的物质的量之比为
- 根号72加根号20分之1减根号49分之18
- 在Rt△ABC中,∠C=90°,已知∠A=60°,c=8根号3,解这个直角三角形
热门考点