Re: [問題] 如何求解方程式"符號運算"與"簡化"
看板Mathematica作者chungyuandye (養花種魚數月亮賞星星)時間12年前 (2012/12/23 10:28)推噓1(1推 0噓 0→)留言1則, 1人參與討論串4/4 (看更多)
※ 引述《jurian0101 (Hysterisis)》之銘言：
: 直接貼到Mathematica讀
: line1=L*D[i[L][t],t]==V[g][t]-r[L]*i[i][t]-v[0][t];
: line2=i[i][t]==(L*D[i[L][t],t])/R[S]+i[L][t];
: line3=v[0][t]==c*r[c]*D[v[c][t],t]+v[c];
: line4=L*D[i[L][t],t]==-(((r[c]*R[S])*(c*D[v[c][t],t]))/(r[L]+R[S]))+
: (R[S]*V[g][t])/(r[L]+R[S])-(v[c]*R[S])/(r[L]+R[S])+
: (i[L][t]*((-r[L])*R[S]))/(r[L]+R[S]);
: line5=c D[v[c][t],t]==R[L]/(r[c]+R[L])(1/R[S]*L*D[i[L][t],t]+i[L][t])
: -v[c]/(r[c]+R[L])
(* 由方程組直接求解 *)
myans1=L*Derivative[1][i[L]][t]/.
Solve[{line1,line2,line3,line5},{i[L]'[t],i[i][t],v[0][t],
v[c]'[t]}][[1]]
(* Eq 2, 3 代入 Eq 1, 再將 Eq 5代入 Eq 1，最後用Solve做移項整理 *)
myans2=L*Derivative[1][i[L]][t]/.
Solve[(line1/.({line2,line3}/.Equal->Rule))/.(line5/.
Equal->Rule),Derivative[1][i[L]][t]][[1]]
(* 係數 *)
mycoef=Table[FullSimplify@Coefficient[z,#]&/@{i[L][t],v[c],
V[g][t]},{z,{myans1,myans2}}]
SameQ@mycoef
--
養花種魚數月亮賞星星
http://chungyuandye.twbbs.org
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.131.237
推
12/23 13:15, , 1F
12/23 13:15, 1F
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 4 之 4 篇)：
Mathematica 近期熱門文章
PTT數位生活區 即時熱門文章
