Оценка производной с использованием первых принципов
Отношения разностей могут непосредственно использоваться для нахождения не только производной первого порядка, но и производных более высокого порядка. Для начала рассмотрим функцию g и связанное с ней отношение разностей.
In[1]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_28.png)
g[x_] := x^2 Exp[x]
In[2]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_29.png)
dq1[x_] = DifferenceQuotient[g[x], {x, h}]
Out[2]=
![](assets.ru/evaluate-a-derivative-using-first-principles/O_26.png)
Использование предела отношения разностей даёт производную первого порядка.
In[3]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_30.png)
Limit[dq1[x], h -> 0]
Out[3]=
![](assets.ru/evaluate-a-derivative-using-first-principles/O_27.png)
In[4]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_31.png)
Limit[dq1[x], h -> 0];
Simplify[% == g'[x]]
Out[4]=
![](assets.ru/evaluate-a-derivative-using-first-principles/O_28.png)
Производная второго порядка может быть вычислена напрямую от второго отношения разностей без ссылки на производную первого порядка.
In[5]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_32.png)
dq2[x_] = DifferenceQuotient[g[x], {x, 2, h}]
Out[5]=
![](assets.ru/evaluate-a-derivative-using-first-principles/O_29.png)
Предел - это производная второго порядка.
In[6]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_33.png)
Limit[dq2[x], h -> 0]
Out[6]=
![](assets.ru/evaluate-a-derivative-using-first-principles/O_30.png)
In[7]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_34.png)
Limit[dq2[x], h -> 0];
Simplify[% == g''[x]]
Out[7]=
![](assets.ru/evaluate-a-derivative-using-first-principles/O_31.png)
Отношение разностей производной первого порядка является функцией, отличной от отношения разностей второго порядка g, но его предел также является производной второго порядка.
In[8]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_35.png)
dqp[x_] = DifferenceQuotient[g'[x], {x, h}]
Out[8]=
![](assets.ru/evaluate-a-derivative-using-first-principles/O_32.png)
In[9]:=
![Click for copyable input](assets.ru/evaluate-a-derivative-using-first-principles/In_36.png)
Limit[dqp[x], h -> 0]
Out[9]=
![](assets.ru/evaluate-a-derivative-using-first-principles/O_33.png)