Благодаря этому ответу на связанный вопрос в Python я написал некоторый код Julia для отображения гистограмм в цикле for.
using PyPlot
γvec = 0.1:0.1:0.9 # γ : fixed paramter
N = 1000 # N : number of samples for each γ
fig, axs = subplots(length(γvec), 1, figsize=(10, 6.18*length(γvec)), constrained_layout=true)
fig.suptitle("Distribution of eigenvalues of ρ at final time", fontsize=20, y=1)
# TL;DR
λ₁vec, λ₂vec = zeros(Complex{Float64}, N), zeros(Complex{Float64}, N);
for (iᵧ, γ) in enumerate(γvec)
# println(γ)
samples = [ myfunction(...) for _ in 1:N ]; # draw samples
for i in 1:N
λ₁₂ = eigvals(samples[i][1]) # calculate (real-valued) eigenvalues of ρ
λ₁vec[i], λ₂vec[i] = maximum(real(λ₁₂)), minimum(real(λ₁₂))
end
axs[iᵧ].hist(λ₁vec) # plot histogram of eigenvalues of ρ
axs[iᵧ].hist(λ₂vec)
axs[iᵧ].set_title("γ = $γ")
axs[iᵧ].set_xlabel("eigenvalues of ρ")
end
gcf() # get current figure
Чтобы создать fig
ure, я добавил синтаксис другого ответа и прочитал пример кода электронной книги на GitHub . Вот фигура, которую я сделал.
Перекрывающиеся заголовок и субтитры, созданные PyPlot Джулии. Что можно сделать, чтобы правильно отобразить заголовок и субтитры?
Отредактировано в ответ на комментарии.
Я использую последнюю версию Юноны и Джулии.
Удалено y=1
, а заголовок идет под первым подзаголовком.