1樓:馬遠馳
樓上已經給出了完整解法,我用Julia也求解了一遍,樓主可以作為參考。
using
DifferentialEquationsusing
Plots
function
dydt(dy
,y,p
,t)delta
,alpha
,beta
,gamma
,omega=p
dy[1]
=y[2
];dy[2
]=-delta*y
[2]-
alpha*y
[1]-
beta*y
[1]^
3+gamma
*cos
(omega*t
);endy0=
[1.0
,1.0
]tspan=(
0.0,
50.0)p
=(1.0,
1.0,
1.0,
1.0,
1.0)
prob
=ODEProblem
(dydt,y0
,tspan,p
)sol
=solve
(prob
)視覺化
plot
(sol
,vars=(
0,1),
linewidth=3
,title
="Solution of the Duffing Equation"
,xaxis
="Time (t)"
,yaxis=""
,label
="x(t)"
)plot!
(sol
,vars=(
0,2),
linewidth=3
,xaxis
="Time (t)"
,yaxis=""
,label
="x'(t)")相圖
plot
(sol
,vars=(
1,2),
linewidth=3
,title
="phase plot between x(t) and x'(t)"
,xaxis
="x(t)"
,yaxis
="x'(t)"
,label="")
如何用matlab畫加權無向圖?
宇智波帶土 在matlab中有乙個Graph and Network Algorithms模組,該模組可以繪製無向圖和有向圖。下面給乙個簡單的例子 Matlab 無向圖 生成資料,A和鄰接矩陣的形勢相似A magic 10 A dist A index 1 size A 1 names forii ...
如何用mathematica證明下述方程只有唯一實根,且該實根為無理根?
用軟體證明是無理根還是奇怪了點。不過其實你可以直接證明該多項式在Q上面不可約,從而根不可能是有理根。對於這種級數截斷,當你截斷的地方是個質數的時候,可以輕鬆使用Eisenstein判別法判斷。當你截斷的地方不是個質數的時候,就稍微要用到一點點代數知識。利用準素分解就能得到結果,這個東西應該是叫Sch...
如何用Python畫樹?
半個馮博士 先寫乙個樹 class TreeNode def init self key self key key self right None self left None definsert self key if self key key return elif self key ifsel...