如何用Matlab 畫達芬方程?

時間 2021-06-03 13:10:37

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...