logo

어웨이

어웨이

줄리아에서 미분방정식 풀이 조기 종료하는 법


줄리아에서 암시적 미분방정식 풀이하는 법

DAE 미분 대수 방정식 https://kr.mathworks.com/help/matlab/math/solve-differential-algebraic-equations-daes.html

import DifferentialEquations as DE
import Sundials
import DiffEqBase
# function f2(out, du, u, p, t)
#     out[1] = 2u[1] + 0.15u[1]^2 - 0.005u[1]^3 - 0.1u[1]*u[2] - 0.1u[1]*du[1] - du[1]
#     out[2] = -u[2] - 0.1u[2]^2 + 0.1u[1]*u[2] - 0.15u[2]*u[3] - 0.015u[1]*u[2]*u[3] + 0.01u[1]*u[2]^2 - 0.1u[1]*du[2] - 0.1u[2]*du[2] - 0.01u[1]*u[2]*du[2] - du[2]
#     out[3] = -0.7u[3] + 0.065u[2]*u[3] - 0.05u[2]*du[3] - du[3]
# end
# prob = DE.DAEProblem(f2, collect(traj0[1, 5:end]), collect(traj0[1, 2:4]), (0, 200), differential_vars = [true, true, true])
# sol = DE.solve(prob, Sundials.IDA())
# plot(eachrow(stack(sol.u))...)

줄리아 좌표축 틱의 위치를 바꾸는 방법

using Plots

heatmap(rand(10, 10), xmirror = true, yflip = true)

줄리아 좌표축 틱에 다른 내용 쓰는 방법 xticks = (1:10, 'a':'j')

using Plots

plot(randn(10), xticks = (1:10, 'a':'j'))
plot(randn(10), xticks = (2:2:10, 'A':'E'))

줄리아 패키지 설치 속도 확인하는 법

https://discourse.julialang.org/t/any-way-to-speed-up-loading-large-precompiled-packages/98762/11

@time_import


줄리아에서 데이터 읽어들일 때 깨지는 인코딩 문제 해결

https://github.com/JuliaStrings/StringEncodings.jl