Dessinez une "courbe de seins" dans un graphique 3D (2)

"Dessiner un graphe dans le langage de programmation Julia" J'ai réussi à déplacer l'échantillon pour le moment, mais les données créées par mon propre programme (en langage Julia) J'ai décidé de découvrir comment dessiner un graphique en utilisant.

Cette fois, j'ai mené une expérience pour agrandir les "seins" en me référant à "Calcul de la longueur par équation de surface courbe des seins". J'ai essayé de. Le résultat du calcul obtenu par Julia est affiché sous forme de graphique par Matplotlib. (Veuillez créer un environnement et le déplacer)

Lors de l'enregistrement du code source, l'encodage des caractères doit être utf-8.
Pour lancer l'exemple, double-cliquez sur Oppai2.jl (code source Julia) depuis l'Explorateur.
Augmenter ALPHA dans le code source augmentera les "seins". Augmenter BETA étendra le "mamelon".
Un fichier appelé "Oppai2.png " sera créé dans le répertoire courant.

Remarque: Pour la construction de l'environnement, reportez-vous à "Dessiner des graphes dans le langage de programmation Julia". Ajout des étapes 4, 6 et 11 au 04/07/2014. ** Si vous avez créé l'environnement en consultant l'article ci-dessus avant le 04/07/2014, veuillez effectuer les étapes supplémentaires 4, 6 et 11. **Je vous remercie.

↓ Exemple de $ ALPHA = 1 $. Les "seins" sont petits. Faites glisser le bouton gauche pour faire pivoter le graphique 3D. Vous pouvez effectuer un zoom avant / arrière en faisant glisser le bouton droit.

↓ Exemple de $ ALPHA = 3 $. Les "seins" sont gros. Faites glisser le bouton gauche pour faire pivoter le graphique 3D. Vous pouvez effectuer un zoom avant / arrière en faisant glisser le bouton droit.

`Oppai2.jl`


# -*- coding: utf-8 -*-

using PyPlot

using PyCall
@pyimport matplotlib.cm as cm
@pyimport matplotlib.font_manager as fm
@pyimport matplotlib.pyplot as plt
@pyimport numpy as np

ALPHA = 1.0
BETA  = 1.0
GAMMA = 1.0 / ALPHA

function bust(x, y)
ALPHA * (6 * exp(-((2 / 3 * abs(x) - 1) ^ 2 + (2 / 3 * y) ^ 2) - 1 / 3 * (2 / 3 * y + 0.5) ^ 3) 
+ BETA * 2 / 3 * exp(-2.818 ^ 11 * ((abs(2 / 3 * x) - 1) ^ 2 + (2 / 3 * y) ^ 2) ^ 2) 
- GAMMA * (2 / 3 * x) ^ 4) / 8
end

fig = plt.figure()
ax = fig[:gca](projection="3d")

nx = 201
ny = 201
x = linspace(-3, 3, nx)
y = linspace(-3, 3, ny)
z = zeros(Float64, length(y), length(x))

for iy = 1:length(y)
    for ix = 1:length(x)
        z[iy, ix] = bust(x[ix], y[iy])
    end
end

x, y = np.meshgrid(x, y)
sf = ax[:plot_surface](x,y,z,rstride=4,cstride=4,cmap=cm.coolwarm,linewidth=0.1,antialiased=true)
fig[:colorbar](sf,shrink=0.5,aspect=5)

ax[:set_zlim](-4,3)
ax[:view_init](0,-90)

@windows? (
begin
    #fp = fm.FontProperties(fname="C:\\WINDOWS\\Fonts\\msgothic.ttc")
    fp = fm.FontProperties(fname="C:\\WINDOWS\\Fonts\\msmincho.ttc")
    plt.title("Julia+Échantillon Matplotlib", fontproperties=fp, fontsize=25)
end
: begin
    plt.title("Julia+Échantillon Matplotlib", fontsize=25)
end
)

plt.savefig("Oppai2")
plt.show(block=true)