Fitxer:Plot of various initial mass functions.svg
Salta a la navegació
Salta a la cerca
Mida d'aquesta previsualització PNG del fitxer SVG: 370 × 355 píxels. Altres resolucions: 250 × 240 píxels | 500 × 480 píxels | 800 × 768 píxels | 1.067 × 1.024 píxels | 2.135 × 2.048 píxels.
Fitxer original (fitxer SVG, nominalment 370 × 355 píxels, mida del fitxer: 73 Ko)
| Descripció |
English: Plot of various initial mass functions (IMF). (The vertical axis is actually not ξ(m)Δm, but a scaled version of ξ(m). For m greater than 1 solar mass, it is (m/solar mass)^-2.35.) |
| Data | |
| Font | Treball propi |
| Autor | JohannesBuchner |
Jo, el titular dels drets d'autor d'aquest treball, el public sota la següent llicència:
This file is licensed under the Creative Commons Attribution-Share Alike 4.0 International license.
- Sou lliure de:
- compartir – copiar, distribuir i comunicar públicament l'obra
- adaptar – fer-ne obres derivades
- Amb les condicions següents:
- reconeixement – Heu de donar la informació adequada sobre l'autor, proporcionar un enllaç a la llicència i indicar si s'han realitzat canvis. Podeu fer-ho amb qualsevol mitjà raonable, però de cap manera no suggereixi que l'autor us dóna suport o aprova l'ús que en feu.
- compartir igual – Si modifiqueu, transformeu, o creeu a partir del material, heu de distribuir les vostres contribucions sota una llicència similar o una de compatible amb l'original.
Generating Python Code
import matplotlib.pyplot as plt
import numpy
from numpy import exp, log10 as log
def salpeter55(m):
alpha = 2.35
return m**-alpha
def millerscalo79(m):
return numpy.where(m > 1, salpeter55(m), salpeter55(1))
def chabrier03individual(m):
k = 0.158 * exp(-(-log(0.08))**2/(2 * 0.69**2))
return numpy.where(m <= 1,\
0.158*(1./m) * exp(-(log(m)-log(0.08))**2/(2 * 0.69**2)),\
k*m**-2.3)
def chabrier03system(m):
k = 0.086 * exp(-(-log(0.22))**2/(2 * 0.57**2))
return numpy.where(m <= 1,\
0.086*(1./m) * exp(-(log(m)-log(0.22))**2/(2 * 0.57**2)),\
k*m**-2.3)
def kroupa01(m):
return numpy.where(m<0.08, m**-0.3, numpy.where(m < 0.5, 0.08**-0.3 * (m/0.08)**-1.3, 0.08**-0.3 * (0.5/0.08)**-1.3 * (m/0.5)**-2.3))
plt.figure(figsize=(4,4))
m = numpy.logspace(-2, 2, 400)
for label, imf in zip('Salpeter55 MillerScalo79 Kroupa01 Chabrier03individual Chabrier03system'.split(),\
[salpeter55, millerscalo79, kroupa01, chabrier03individual, chabrier03system]):
plt.plot(m, imf(m)/imf(1), label=label)
plt.gca().set_yscale('log')
plt.gca().set_xscale('log')
plt.xlim(1e-2, 100)
plt.ylim(1e-3, 1e3)
plt.legend(loc='best', prop=dict(size=8))
plt.xlabel('Mass [Solar mass]')
plt.ylabel(r'Mass Function $\xi(m)\Delta m$')
plt.savefig('imf.pdf', bbox_inches='tight')
plt.savefig('imf.svg', bbox_inches='tight')
plt.savefig('imf.png', bbox_inches='tight')
Historial del fitxer
Cliqueu una data/hora per veure el fitxer tal com era aleshores.
| Data/hora | Miniatura | Dimensions | Usuari/a | Comentari | |
|---|---|---|---|---|---|
| actual | 20:14, 4 nov 2015 | | 370 × 355 (73 Ko) | wikimediacommons>JohannesBuchner | Cross-wiki upload from en.wikipedia.org |
Ús del fitxer
La pàgina següent utilitza aquest fitxer:
