Matlab微分

Question

我需要在matlab中編寫for循環來解決使用前向差分方法的導數。使用0.25的步長從0到10導出的函數是10 + 15x + 20x^2。我一直在使用

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);