主題
Search

尤拉-拉格朗日微分方程

下載 Mathematica Notebook下載 Wolfram Notebook

尤拉-拉格朗日微分方程是變分法的基本方程。它指出，如果 J 由 積分 形式 定義

J=intf(t,y,y^.)dt,
(1)

其中

y^.=(dy)/(dt),
(2)

那麼，如果尤拉-拉格朗日微分方程

(partialf)/(partialy)-d/(dt)((partialf)/(partialy^.))=0
(3)

得到滿足，則 J 具有平穩值。

如果時間導數 符號 y^. 被空間導數符號 y_x 替代，則方程變為

(partialf)/(partialy)-d/(dx)(partialf)/(partialy_x)=0.
(4)

尤拉-拉格朗日微分方程的實現形式為EulerEquations[f, u[x], x]，在 Wolfram 語言 包中VariationalMethods` .

在許多物理問題中，f_x （f 關於 x 的偏導數）結果為 0，在這種情況下，對尤拉-拉格朗日微分方程的 manipulation 簡化為大大簡化的和部分積分的形式，稱為 Beltrami 恆等式，

f-y_x(partialf)/(partialy_x)=C.
(5)

對於三個自變數（Arfken 1985, pp. 924-944），方程推廣為

(partialf)/(partialu)-partial/(partialx)(partialf)/(partialu_x)-partial/(partialy)(partialf)/(partialu_y)-partial/(partialz)(partialf)/(partialu_z)=0.
(6)

變分法中的問題通常可以透過求解適當的尤拉-拉格朗日方程來解決。

為了推導尤拉-拉格朗日微分方程，考察

deltaJ=deltaintL(q,q^.,t)dt
(7)
=int((partialL)/(partialq)deltaq+(partialL)/(partialq^.)deltaq^.)dt
(8)
=int[(partialL)/(partialq)deltaq+(partialL)/(partialq^.)(d(deltaq))/(dt)]dt,
(9)

因為 deltaq^.=d(deltaq)/dt。現在，使用 分部積分法 對第二項積分

u=(partialL)/(partialq^.) dv
(10)
=d(deltaq)
(11)
du=d/(dt)((partialL)/(partialq^.))dt v=deltaq,
(12)

因此

int(partialL)/(partialq^.)(d(deltaq))/(dt)dt=int(partialL)/(partialq^.)d(deltaq)=[(partialL)/(partialq^.)deltaq]_(t_1)^(t_2)-int_(t_1)^(t_2)(d/(dt)(partialL)/(partialq^.)dt)deltaq.
(13)

結合 (◇) 和 (◇) 得到

deltaJ=[(partialL)/(partialq^.)deltaq]_(t_1)^(t_2)+int_(t_1)^(t_2)((partialL)/(partialq)-d/(dt)(partialL)/(partialq^.))deltaqdt.
(14)

但是我們只改變路徑，而不是端點，所以 deltaq(t_1)=deltaq(t_2)=0，(14) 變為

deltaJ=int_(t_1)^(t_2)((partialL)/(partialq)-d/(dt)(partialL)/(partialq^.))deltaqdt.
(15)

我們正在尋找平穩值，使得 deltaJ=0。這些必須對任何小的變化 deltaq 消失，由 (15) 得到，

(partialL)/(partialq)-d/(dt)((partialL)/(partialq^.))=0.
(16)

這就是尤拉-拉格朗日微分方程。

J 的變分也可以用引數 kappa 表示為

deltaJ=int[f(x,y+kappav,y^.+kappav^.)-f(x,y,y^.)]dt
(17)
=kappaI_1+1/2kappa^2I_2+1/6kappa^3I_3+1/(24)kappa^4I_4+...,
(18)

其中

v=deltay
(19)
v^.=deltay^.
(20)

並且一階、二階等等變分是

I_1=int(vf_y+v^.f_(y^.))dt
(21)
I_2=int(v^2f_(yy)+2vv^.f_(yy^.)+v^.^2f_(y^.y^.))dt
(22)
I_3=int(v^3f_(yyy)+3v^2v^.f_(yyy^.)+3vv^.^2f_(yy^.y^.)+v^.^3f_(y^.y^.y^.))dt
(23)
I_4=int(v^4f_(yyyy)+4v^3v^.f_(yyyy^.)+6v^2v^.^2f_(yyy^.y^.)+4vv^.^3f_(yy^.y^.y^.)+v^.^4f_(y^.y^.y^.y^.))dt.
(24)

二階變分可以使用以下方式重新表示

d/(dt)(v^2lambda)=v^2lambda^.+2vv^.lambda,
(25)

因此

I_2+[v^2lambda]_2^1=int_1^2[v^2(f_(yy)+lambda^.)+2vv^.(f_(yy^.)+lambda)+v^.^2f_(y^.y^.)]dt.
(26)

但是

[v^2lambda]_2^1=0.
(27)

現在選擇 lambda 使得

f_(y^.y^.)(f_(yy)+lambda^.)=(f_(yy^.)+lambda)^2
(28)

和 z 使得

f_(yy^.)+lambda=-(f_(yy^.))/z(dz)/(dt)
(29)

因此 z 滿足

f_(y^.y^.)z^..+f^._(y^.y^.)z^.-(f_(yy)-f^._(yy^.))z=0.
(30)

然後得出結論

I_2=intf_(y^.y^.)(v^.+(f_(yy^.)+lambda)/(f_(y^.y^.))v)^2dt
(31)
=intf_(y^.y^.)(v^.-v/z(dz)/(dt))^2dt.
(32)

另請參閱

Beltrami 恆等式、最速降線問題、變分法、尤拉-拉格朗日導數、泛函導數、變分

使用 探索

參考文獻

Arfken, G. Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, 1985.Forsyth, A. R. Calculus of Variations. New York: Dover, pp. 17-20 and 29, 1960.Goldstein, H. Classical Mechanics, 2nd ed. Reading, MA: Addison-Wesley, p. 44, 1980.Lanczos, C. The Variational Principles of Mechanics, 4th ed. New York: Dover, pp. 53 and 61, 1986.Morse, P. M. and Feshbach, H. "The Variational Integral and the Euler Equations." §3.1 in Methods of Theoretical Physics, Part I. New York: McGraw-Hill, pp. 276-280, 1953.

在 上被引用

尤拉-拉格朗日微分方程

請引用本文為

Weisstein, Eric W. "尤拉-拉格朗日微分方程。" 來自 Web 資源。 https://mathworld.tw/Euler-LagrangeDifferentialEquation.html

學科分類