Julia Bifurcation Kit 延续的记录解决方案

Recording solutions to Julia Bifurcation Kit continuation

提问人:George 提问时间:11/18/2023 更新时间:11/18/2023 访问量:9

问:

我找不到正确的代码行来从 Julia Bifurcation 工具包中获取特定变量的值。我有一个常微分方程系统和一个系统的固定点。当我更改参数时,我可以使用延续函数来获取变量的范数,但我找不到任何方法来获取一个变量的值。我想做的是分别获取所有变量的值并将它们放在文本文件中以便稍后绘制,但我无法让它工作。我尝试使用 and,但是当我打印出分岔点时,我获得的唯一值是变量的范数。printFromSolutionrecordFromSolution

以下是代码的相关部分:

prob = BifurcationProblem(cryptOde, x0, par_sl, (@lens _.I))

opts = ContinuationPar(p_min = -1000., p_max = 2000.,dsmax=0.1,max_steps=100000)

recordedVariables(x, p) = (x1 = x[1], x2 = x[2], nrm = norm(x))

br = continuation(prob, PALC(), opts, bothside = false,printFromSolution=recordedVariables,recordFromSolution=recordedVariables)

当我打印延续曲线的第一个点时,我得到以下结果:

(x = 3119.4774565745524, param = 0.3696058344346101, itnewton = 1, itlinear = 1, ds = 0.014608, n_unstable = 0, n_imag = 0, stable = true, step = 1, eigenvals = ComplexF64[-0.6690887317560925 + 0.0im, -0.7900685522295721 + 0.0im, -0.9999999999999996 + 0.0im, -1.2927168332714745 + 0.0im], eigenvecs = ComplexF64[2.59609717158152e-5 + 0.0im 2.76430039760115e-5 + 0.0im 1.242989057545907e-19 + 0.0im -3.767156349151788e-5 + 0.0im; 9.890022861201471e-6 + 0.0im 6.266810756934409e-6 + 0.0im -1.988920281189839e-22 + 0.0im 8.748783810086575e-6 + 0.0im; 2.5960971715805685e-5 + 0.0im 2.7643003976008426e-5 + 0.0im 0.0 + 0.0im -3.767156349152811e-5 + 0.0im; -0.0008218725707738361 + 0.0im 0.0002776367613707426 + 0.0im 4.854376184199556e-19 + 0.0im -1.270175168102972e-6 + 0.0im; -0.972629467641427 + 0.0im -0.9726279574486716 + 0.0im 0.9999999999988511 + 0.0im 0.9727830062088821 + 0.0im; 0.23229844784378212 + 0.0im 0.23234421809069009 + 0.0im 7.052115215136131e-16 + 0.0im -0.2316755739897396 + 0.0im; 0.005240234288655535 + 0.0im 0.003292223127212381 + 0.0im -2.1416365518527453e-19 + 0.0im 0.004432249430563173 + 0.0im; -0.0010898855800747119 + 0.0im 0.00029980358141065564 + 0.0im 3.4392163450341146e-19 + 0.0im -7.19935385992619e-8 + 0.0im; -0.00015804640699194087 + 0.0im -0.00011594838773710004 + 0.0im 1.5157582958993832e-6 + 0.0im -5.9129236555397425e-5 + 0.0im])

输出的变量是系统的范数。x

我还尝试保存到文件并将文件名传递到延续函数中,但文件没有出现在我的目录中。

朱莉娅

评论


答: 暂无答案