Mathematica 教學 聯立方程組的係數矩陣

(* 微分方程組 *)
dur=(EA/R-KGA/R*Subscript[d,\[Theta]])*
ur+(EA/R*Subscript[d,\[Theta]]+KGA/R*Subscript[d,\[Theta]])*
ue-KGA*Subscript[d,\[Theta]]*phiz;
due=(EA/R*Subscript[d,\[Theta]]-
KGA/R*Subscript[d,\[Theta]])ue+(EA/R*Subscript[
dd,\[Theta]]+KGA/R)*ue-KGA*phiz;

(* 定義變數 *)
myVars={ur,ue,uz,phir,phie,phiz};

(* 定義微分運算子 *)
myOperator=
Flatten[{Subscript[d,\[Theta]]*#->#'[\[Theta]],
Subscript[dd,\[Theta]]*#->#''[\[Theta]],
(* 所有變數都是\[Theta]的函數 *)
#->#[\[Theta]]&/@myVars}&/@myVars];

test1=Coefficient[#,myVars]&/@{dur,due}
test1//MatrixForm
test2=test1.myVars
test2//MatrixForm
ExpandAll@test2/.myOperator//MatrixForm

Tags: Mathematica 教學, Mathematica 教學 微積分, Teaching

• 編輯文章
• 此文章網址
• 轉寄此文章
• 目前有0 意見