Criando Ruído Branco:

http://en.wikipedia.org/wiki/Colors_of_noise

import numpy as N import wave class SoundFile: def __init__(self, signal): self.file = wave.open('test.wav', 'wb') self.signal = signal self.sr = 44100 def write(self): self.file.setparams((1, 2, self.sr, 44100*4, 'NONE', 'noncompressed')) self.file.writeframes(self.signal) self.file.close() # let's prepare signal duration = 9 # seconds samplerate = 44100 # Hz samples = duration*samplerate ydata = 16384 * N.random.random_sample(samples) signal = N.resize(ydata, (samples,)) ssignal = '' for i in range(len(signal)): ssignal += wave.struct.pack('h',signal[i]) # transform to binary f = SoundFile(ssignal) f.write() |

Algoritmo de Dijkstra para encontrar caminho mínimo em grafo

# Dijkstra's algorithm for shortest paths # David Eppstein, UC Irvine, 4 April 2002 # http://aspn.activestate.com/ASPN/Cookbook/Python/Recipe/117228 from priodict import priorityDictionary def Dijkstra(graph,start,end=None): """ Find shortest paths from the start vertex to all vertices nearer than or equal to the end. The input graph G is assumed to have the following representation: A vertex can be any object that can be used as an index into a dictionary. G is a dictionary, indexed by vertices. For any vertex v, G[v] is itself a dictionary, indexed by the neighbors of v. For any edge v->w, G[v][w] is the length of the edge. This is related to the representation in <http://www.python.org/doc/essays/graphs.html> where Guido van Rossum suggests representing graphs as dictionaries mapping vertices to lists of neighbors, however dictionaries of edges have many advantages over lists: they can store extra information (here, the lengths), they support fast existence tests, and they allow easy modification of the graph by edge insertion and removal. Such modifications are not needed here but are important in other graph algorithms. Since dictionaries obey iterator protocol, a graph represented as described here could be handed without modification to an algorithm using Guido's representation. Of course, G and G[v] need not be Python dict objects; they can be any other object that obeys dict protocol, for instance a wrapper in which vertices are URLs and a call to G[v] loads the web page and finds its links. The output is a pair (D,P) where D[v] is the distance from start to v and P[v] is the predecessor of v along the shortest path from s to v. Dijkstra's algorithm is only guaranteed to work correctly when all edge lengths are positive. This code does not verify this property for all edges (only the edges seen before the end vertex is reached), but will correctly compute shortest paths even for some graphs with negative edges, and will raise an exception if it discovers that a negative edge has caused it to make a mistake. """ final_distances = {} # dictionary of final distances predecessors = {} # dictionary of predecessors estimated_distances = priorityDictionary() # est.dist. of non-final vert. estimated_distances[start] = 0 for vertex in estimated_distances: final_distances[vertex] = estimated_distances[vertex] if vertex == end: break for edge in graph[vertex]: path_distance = final_distances[vertex] + graph[vertex][edge] if edge in final_distances: if path_distance < final_distances[edge]: raise ValueError, \ "Dijkstra: found better path to already-final vertex" elif edge not in estimated_distances or path_distance < estimated_distances[edge]: estimated_distances[edge] = path_distance predecessors[edge] = vertex return (final_distances,predecessors) def shortestPath(graph,start,end): """ Find a single shortest path from the given start vertex to the given end vertex. The input has the same conventions as Dijkstra(). The output is a list of the vertices in order along the shortest path. """ final_distances,predecessors = Dijkstra(graph,start,end) path = [] while 1: path.append(end) if end == start: break end = predecessors[end] path.reverse() return path |

Um decorador simples para fazer cache de dados em funções muito demoradas.

Na primeira chamada de func(), o resultado vai demorar 1 segundo para ser retornado. Já na segunda, será instantâneo.

É claro que só deve ser usado em funções determinísticas.

class Cache: def __init__(self, myfunc): self.elements = {} self.func = myfunc def __call__(self, *args, **kw): key = (args, frozenset(kw.items())) if kw else args try: find = self.elements[key] except KeyError: value = self.func(*args, **kw) self.elements[key] = value return value else: return find def clear(self): self.elements = {} # Test Function from time import sleep @Cache def func(num): sleep(1) #Simulando um calculo demorado return num**num |