1.数值微分
数值微分是用离散的方法近似地计算函数y=f(x)在某点x=a处的导数值,前差公式:
后差公式:
中心差商:
matlab中可以用diff函数实现前差公式的数值微分:
函数形式:diff(X)
求解X的一阶微分,X是一个向量也可以是一个矩阵,结果:[X(2)-X(1) X(3)-X(2) … X(n)-X(n-1)]
X = [1 1 2 3 5 8 13 21];
Y = diff(X)
Y =
0 1 1 2 3 5 8
diff(X,n)
%求解n阶微分
X = [0 5 15 30 50 75 105];
Y = diff(X,2)
Y =
5 5 5 5 5
ode23函数:
主要形式[t,y]=ode23('fun',t,y0,options)
[t,y]=ode45('fun',t,y0,options)
fun为函数,t=[t0,tf]时,t0、tf为起始与终止值,y0函数的初值
tspan = [0 5];
y0 = 0;
[t,y] = ode23(@(t,y) 2*t, tspan, y0);
%求0至5之间的微分数值,
求解下面的微分方程
y′′