2016-12-14 1 views
0

ich das folgende Raster Mal habe definiert:ein Gitter definiert mit np.linspace und Bewertungsfunktion über dieses Gitter

T = 10000 
tmin = 0 
tmax = 10**20 
t = np.linspace(tmin, tmax, T) 

so meine Schrittgrßenregister irgendwo ** 16 um 10. Ich möchte eine Funktion von tmin bis tmax gemäß dieser Schrittgröße auswerten. Wird im Folgenden einfach von tmin bis 10.000 oder über das gesamte Raster ausgewertet?

timedep_PD = np.zeros([t,Y]) 
for t in range(tmin,T): 
    timedep_PD[t]= probdensity_func(x,t,0) 

Wie kann ich beurteilen diese Funktion tminx-tmax in Schrittweiten von (Tmax-Tmin/T), während nach wie vor in Form von timedep_PD[t]

Antwort

1

einfach speichern, um es einfach:

T = 10 
tmin = 0 
tmax = 10**2 
t = np.linspace(tmin, tmax, T) 

timedep_PD = np.zeros((len(t),T)) ##len(t) is T : np.zeros((T,T)) 
dt = ((tmax-tmin)/T) 
k=0 
for t in np.arange(tmin,tmax,dt): 
    timedep_PD[k]= np.sin(t) 
    print timedep_PD[k] 
    k+=1 

Wir können auch np.meshgrid

verwenden