Re: [問題] deqn 問題再哪裡呢??

看板Mathematica作者 (養花種魚數月亮賞星星)時間12年前 (2012/09/22 18:53), 編輯推噓4(406)
留言10則, 2人參與, 最新討論串2/2 (看更多)
※ 引述《ntust661 (TOEFL_5!)》之銘言: : http://ppt.cc/ipoL : 我都刪除變數了 : 他一直說我的方程式有錯@@ : why?? : 各位大大救救小弟阿QQ {Pa,Qa,Maa}={1,1,1}; sys={\[Theta]'[s]==4Pi^2*(Maa-Pa-Qa*X[s]), X'[s]==Cos[\[Theta][s]], y'[s]==Sin[\[Theta][s]], M'[s]==Fx[s]*Sin[\[Theta][s]]-Fy[s]*Cos[\[Theta][s]], Fx'[s]==0, Fy'[s]==0}; BC={\[Theta][0]==0,X[0]==0,y[0]==0,Fx[0]==-Pa, Fy[0]==Qa,M[0]==Maa}; ans=NDSolve[ Flatten@{sys,BC},{\[Theta][s],X[s],y[s],M[s],Fx[s], Fy[s]},{s,0,1}][[1]] GraphicsGrid@Partition[Plot[ans[[#,2]],{s,0,1}]&/@Range[6],3] "=" 是 Set 不是 Equal -- 養花種魚數月亮賞星星 http://chungyuandye.twbbs.org -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 218.173.129.126

09/22 19:04, , 1F
老師好._./ 我不知道我錯在哪裡耶QQ...
09/22 19:04, 1F

09/22 19:06, , 2F
一階導數那些要用==,BC那邊最後應該是M[0]==Ma
09/22 19:06, 2F

09/22 19:17, , 3F
我使用你的打法 他最後出現
09/22 19:17, 3F

09/22 19:18, , 4F
{0==4*pi^2*x[s],True,True, ...}
09/22 19:18, 4F

09/22 19:18, , 5F
除了第一個不是True,其他都True!
09/22 19:18, 5F

09/22 19:20, , 6F
好多True,可是沒半個interpolation function QQ
09/22 19:20, 6F

09/22 20:01, , 7F
重開Kernel true 就沒出來了
09/22 20:01, 7F

09/22 20:01, , 8F
然後發現一個重要的東東!!!
09/22 20:01, 8F

09/22 20:04, , 9F
IC = {x'[t],y'[t],z'[t],...} 直接放進去NDSolve
09/22 20:04, 9F

09/22 20:04, , 10F
他竟然不會讀到 X'[t] == ...
09/22 20:04, 10F
文章代碼(AID): #1GNPYGAK (Mathematica)
討論串 (同標題文章)
文章代碼(AID): #1GNPYGAK (Mathematica)