一、问题提出

已知一系列的数据点(xi,yi),i=0,1,2,...,n(x_i,y_i),i=0,1,2,...,n(xi,yi),i=0,1,2,...,n，且横坐标等间距(即i≥1i\ge1i≥1且i≤ni\le ni≤n时，恒有xi−xi−1=hx_i-x_{i-1}=hxi−xi−1=h)，如何计算/估计一阶导数f′(xi)f'(x_i)f′(xi)及二阶导数f′′(xi)f''(x_i)f′′(xi)?

二、一阶导数和二阶导数的二阶中心差分格式

如上图，如何利用点B及其左右一个点的坐标(共三个点的坐标)，来估计点B处一阶导数和二阶导数呢？想法非常简单：既然三个点的坐标已知，那么三个点可以插值得到一个一元二次多项式，求得一元二次多项式的系数后，就可以求多项式在B点处的一阶和二阶导数，从而得到一阶导数和二阶导数的二阶中心差分格式。利用matlab的符号计算可以很容易推导得到。
一阶导数的二阶中心差分格式：
f′(x0+h)=f(x0+2h)−f(x0)2h(1)f'(x_0+h)=\frac{f(x_0+2h)-f(x_0)}{2h} \tag 1 f′(x0+h)=2hf(x0+2h)−f(x0)(1)
二阶导数的二阶中心差分格式：
f′′(x0+h)=f(x0)−2f(x0+h)+f(x0+2h)h2(2)f''(x_0+h)=\frac{f(x_0)-2f(x_0+h)+f(x_0+2h)}{h^2} \tag 2 f′′(x0+h)=h2f(x0)−2f(x0+h)+f(x0+2h)(2)

三、一阶导数和二阶导数的四阶中心差分格式

与一阶导数和二阶导数的二阶中心差分推导过程一样，一阶导数和二阶导数的四阶中心差分推导，只需将一元二次多项式插值改为一元四次多项式插值即可。
一阶导数的四阶中心差分格式：
f′(x0+2h)=f(x0)−8f(x0+h)+8f(x0+3h)−f(x0+4h)12h(3)f'(x_0+2h)=\frac{f(x_0)-8f(x_0+h)+8f(x_0+3h)-f(x_0+4h)}{12h} \tag 3 f′(x0+2h)=12hf(x0)−8f(x0+h)+8f(x0+3h)−f(x0+4h)(3)
二阶导数的四阶中心差分格式：
f′′(x0+2h)=−f(x0)−16f(x0+h)+30f(x0+2h)−16f(x0+3h)+f(x0+4h)12h2(4)f''(x_0+2h)=-\frac{f(x_0)-16f(x_0+h)+30f(x_0+2h)-16f(x_0+3h)+f(x_0+4h)}{12h^2} \tag 4 f′′(x0+2h)=−12h2f(x0)−16f(x0+h)+30f(x0+2h)−16f(x0+3h)+f(x0+4h)(4)

四、MATLAB代码

clc
clear
close allsyms x x0 y0 y1 y2 y3 y4 h real%% 一阶导数和二阶导数的二阶中心差分格式
a = [1, x0, x0^2; 1, (x0 + h), (x0 + h)^2; 1, (x0 + 2 * h), (x0 + 2 * h)^2] \ [y0; y1; y2];  %一元二次多项式y(x) = a1 + a2 * x + a3 * x^2的系数
y(x) = a(1) + a(2) * x + a(3) * x^2;
dy(x) = diff(y, 1);
ddy(x) = diff(y, 2);
dy_two_order_central_difference = simplify(dy(x0 + h))
ddy_two_order_central_difference = simplify(ddy(x0 + h))%% 一阶导数和二阶导数的四阶中心差分格式
a = [1, x0, x0^2, x0^3, x0^4; 1, (x0 + h), (x0 + h)^2, (x0 + h)^3, (x0 + h)^4; 1, (x0 + 2 * h), (x0 + 2 * h)^2, (x0 + 2 * h)^3, (x0 + 2 * h)^4; ...1, (x0 + 3 * h), (x0 + 3 * h)^2, (x0 + 3 * h)^3, (x0 + 3 * h)^4; 1, (x0 + 4 * h), (x0 + 4 * h)^2, (x0 + 4 * h)^3, (x0 + 4 * h)^4] \ [y0; y1; y2; y3; y4];  %一元四次多项式y(x) = a1 + a2 * x + a3 * x^2 + a4 * x^3 + a5 * x^4的系数
y(x) = a(1) + a(2) * x + a(3) * x^2 + a(4) * x^3 + a(5) * x^4;
dy(x) = diff(y, 1);
ddy(x) = diff(y, 2);
dy_four_order_central_difference = simplify(dy(x0 + 2 * h))
ddy_four_order_central_difference = simplify(ddy(x0 + 2 * h))%% 验证
n = 50;
x = linspace(0, 2*pi, n);
h = x(2) - x(1);
y = sin(x);
dy = cos(x);
ddy = -sin(x);%一阶导数和二阶导数的二阶中心差分(第一个点和最后一个点无法估算一阶导数和二阶导数)
dy1 = nan * zeros(size(x));
ddy1 = nan * zeros(size(x));
for i = 2 : n - 1dy1(i) = (y(i + 1) - y(i - 1)) / (2.0 * h);ddy1(i) = (y(i - 1) - 2.0 * y(i) + y(i + 1)) / h^2;
end%一阶导数和二阶导数的四阶中心差分(第一二个点和最后两个点无法估算一阶导数和二阶导数)
dy2 = nan * zeros(size(x));
ddy2 = nan * zeros(size(x));
for i = 3 : n - 2dy2(i) = (y(i - 2) - 8.0 * y(i - 1) + 8.0 * y(i + 1) - y(i + 2)) / (12.0 * h);ddy2(i) = -(y(i - 2) - 16.0 * y(i - 1) + 30.0 * y(i) - 16.0 * y(i + 1) + y(i + 2)) / (12.0 * h^2);
endmax_dy1_err = max(abs(dy1(2 : n - 1) - dy(2 : n - 1)));
max_ddy1_err = max(abs(ddy1(2 : n - 1) - ddy(2 : n - 1)));
max_dy2_err = max(abs(dy2(3 : n - 2) - dy(3 : n - 2)));
max_ddy2_err = max(abs(ddy2(3 : n - 2) - ddy(3 : n - 2)));
disp(['一阶导数的二阶和四阶中心差分近似，最大误差分别为:', num2str(max_dy1_err), ',' , num2str(max_dy2_err)])
disp(['二阶导数的二阶和四阶中心差分近似，最大误差分别为:', num2str(max_ddy1_err), ',' , num2str(max_ddy2_err)])

一阶导数和二阶导数的二阶/四阶中心差分格式相关推荐

针对图像的一阶导数和二阶导数
针对图像的一阶导数和二阶导数解释关于图像x,y方向一阶导数和二阶导数的问题.总结如下: 基础知识数字图像中,f(x,y)可表示成一个M*N的二维数字阵列,如下图: 数字阵列: 一阶差分对于图像f ...
关于图形学中图像的一阶导数和二阶导数求法
对于图像求一阶导数和二阶导数可以增强图片或者求图像的边缘. 一阶导数: 我们知道在数学中一阶导数的求法为: 在图形学中由于图片是按像素来离散组成的,所以最小的h取值为1,所以计算后: 所以图像的一阶导 ...
一阶导数和二阶导数的意义
一阶求导是求函数各点的斜率整体就是函数的单调性,而二阶是求函数整体的凹凸性,也就相当于求各点斜率的增减规律,也就是凹凸性拐点不一定是极值点但极值点一定是拐点凹凸性可以纵观函数走势凹凸性的作用 ...
Gamma函数(伽玛函数)的一阶导数、二阶导数公式推导及java程序
本文作者:合肥工业大学管理学院钱洋 email:1563178220@qq.com 内容可能有不到之处,欢迎交流. 未经本人,允许禁止转载. 本文的另外一个地址为:http://www.datal ...
java 伽玛函数_Gamma函数(伽玛函数)的一阶导数、二阶导数公式推导及java程序
Gamma函数(伽玛函数)的一阶导数.二阶导数公式推导及java程序 Gamma函数(伽玛函数)的一阶导数.二阶导数公式推导及java程序本文作者:合肥工业大学管理学院钱洋 email:[ema ...
matlab离散数据差分法,MATLAB中用差分法求解函数 or 离散数据的一阶导数和二阶导数...
再进行编程.以下是我的求解程序: clc;clear all h=0.01; %x属于[a,b] a=-5;b=5; x=a:h:b; n=length(x); %定义y y=sin(0.3*x).* ...
图像的一阶导数和二阶导数MATLAB实现
基于二阶导数的laplace算子 I=imread('11.jpg'); [H,W]=size(I); M=double(I); J=M; for i=2:H-1 for j=2:W-1 J(i,j) ...
一维稳态对流扩散问题，无源项，中心差分格式的python程序
一维稳态对流扩散中心差分格式参考书籍:陶文铨的数值传热学+李人宪的有限体积法基础具体的理论可以参考此文章:有限体积法(5)--对流-扩散方程的离散例题一个长度L为1的规则物体,左边界温度恒为1 ...
由参数方程确定的函数的一阶导数与二阶导数
高阶反函数的导数 {x=u(t)y=v(t)\begin{cases} x=u(t)& \text{}\\ y=v(t)& \text{} \end{cases} {x=u(t)y=v ...
二阶混合差分的几种差分格式
二阶混合差分的几种差分格式对上面两式子对y求偏导,得到: 两式子相减,得到: 将二阶混合项写在左边: ( ∂ u / ∂ y ) i + 1 , j (\partial u/\partial y)_ ...

一阶导数和二阶导数的二阶/四阶中心差分格式

文章目录

一、问题提出

二、一阶导数和二阶导数的二阶中心差分格式

三、一阶导数和二阶导数的四阶中心差分格式

四、MATLAB代码

一阶导数和二阶导数的二阶/四阶中心差分格式相关推荐

最新文章

热门文章