Julia之坑矩陣運(yùn)算和解ODE

先做一個(gè)簡(jiǎn)單的矩陣計(jì)算，矩陣乘和矩陣除。使用rand 或者randn創(chuàng)造隨機(jī)數(shù)字構(gòu)成的矩陣

A = randn(4,4)
x = rand(4)
b = A*x
x? = A\b 

x? == x
y = A*x? - b
println(y)

得到y = [0.0, 0.0, -2.77556e-17, 0.0]，基本上一致的。

先做一個(gè)矩陣，使用|>pipe 功能直接計(jì)算A + A'，在賦值給A

A = randn(4,4) |> w -> w + w' 
lambda = eigmax(A)
@show det(A - lambda*I)

得到答案為-3.886340427176376e-14，很接近0了

測(cè)試一個(gè)官網(wǎng)的案例，求解洛倫茨方程，具體內(nèi)容不懂，反正很復(fù)雜就是了。方程如下
\frac{dx}{dt} = \sigma(y-x)
\frac{dy}{dt} = x(\rho -z ) - y
\frac{dz}{dt} = xy - \beta z

使用DifferentialEquations 的求解該方程

image.png

成功！

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點(diǎn)，簡(jiǎn)書系信息發(fā)布平臺(tái)，僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容