两个长度相同的数组x和y,二者的函数关系不确定,可以通过下列代码求y相对于x的数值微分
diff_array = diff(y)./diff(x)
函数说明
函数名 | diff |
含义 | 求数组中相邻元素的差分 |
语法 | diff(x) |
说明 | 若x是长度为m的向量:[x1, x2, x3, ..., xm] 则diff(x)是长度为m-1的向量:[x(2)-x(1), x(3)-x(2), ..., x(m)-x(m-1)] |
因此,从微分的定义出发,求一个数组相对于另一个数组的数值微分可以转化为求两个数组差分结果的商。
举例
x = 1:1:100;
y = x.^3;
figure;
subplot(3,1,1);
disp(['the length of y is: ', num2str(length(y))])
plot(y);
title('y=x^3');
diff_y_x = diff(y)./d