matlab - Matlab微分

Question

我需要在matlab中编写一个for循环来使用前向差分法求解导数。导出的函数是从 0 到 10 的 10+15x+20x^2，步长为 0.25。我试过使用

h=.25; 
x=[0:h:10]; 
y = 10+15*x+20*x.^2; 
y(1) = 45; size(x)  
for i=2:47,     
       y(i) = y(i-1) + h*(15+40*x); 
end

Answer 1

我会这样做，作为开始，

h=.25; 
x=[0:h:10]; 
y = 10+15*x+20*x.^2; 
diff(y)./diff(x)

或者，作为替代方案，

syms x;
y = 20.*x.^2 + 15.*x + 10;
dy = diff(y,1);
h=.25; 
xx=[0:h:10];
res = subs(dy,xx);