中心差分法matlab代码

以下是中心差分法的Matlab代码：

function [dydx] = center_diff(x,y)
%中心差分法求导数
n = length(x);
dydx = zeros(1,n);
for i=2:n-1
dydx(i) = (y(i+1)-y(i-1))/(x(i+1)-x(i-1));
end
dydx(1) = (y(2)-y(1))/(x(2)-x(1));
dydx(n) = (y(n)-y(n-1))/(x(n)-x(n-1));
end

其中，输入参数x和y分别为自变量和因变量，输出参数dydx为y对x的导数。在for循环中，通过中心差分公式计算每个点的导数，首尾两个点使用一阶差分公式计算。