中心差分法 matlab代码

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

function [y_prime] = center_diff(x, y, h)
% 中心差分法计算导数
% x: 自变量，y: 函数值，h: 步长
n = length(x);
y_prime = zeros(1,n);
for i=2:n-1
y_prime(i) = (y(i+1) - y(i-1)) / (2*h);
end
y_prime(1) = (y(2) - y(1)) / h;
y_prime(n) = (y(n) - y(n-1)) / h;
end

希望对你有所帮助！