PPy3/Matplotlib - Historia wersji

RobertJB: /* Przykład bardziej zaawansowany */

2021-01-08T09:45:01Z

Przykład bardziej zaawansowany

RobertJB: /* Przykład bardziej zaawansowany */

2021-01-08T09:44:15Z

Przykład bardziej zaawansowany

RobertJB: /* Przykład bardziej zaawansowany */

2021-01-08T09:40:05Z

Przykład bardziej zaawansowany

RobertJB o 12:29, 27 mar 2018

2018-03-27T12:29:35Z

RobertJB: /* Ćwiczenia */

2018-02-15T12:32:35Z

Ćwiczenia

RobertJB: /* Krzywa parametryczna */

2018-02-06T16:30:32Z

Krzywa parametryczna

RobertJB o 16:26, 6 lut 2018

2018-02-06T16:26:23Z

RobertJB: /* Ćwiczenia */

2017-06-30T12:44:02Z

Ćwiczenia

RobertJB: /* Ćwiczenia */

2017-04-21T14:03:00Z

Ćwiczenia

RobertJB: /* Ćwiczenia */

2017-04-21T14:00:18Z

Ćwiczenia

@@ Linia 197: / Linia 197: @@
-n, bins, patches = plt.hist(x, 50, normed=True,
+n, bins, patches = plt.hist(x, 50, density=True,
      facecolor='green', edgecolor='black', alpha=0.75)
 # Tu w jawny sposób odbieramy zwracane przez plt.hist obiekty

@@ Linia 197: / Linia 197: @@
-n, bins, patches = plt.hist(x, 50, normed=True, facecolor='green', alpha=0.75)
+n, bins, patches = plt.hist(x, 50, normed=True,
 # Tu w jawny sposób odbieramy zwracane przez plt.hist obiekty
 # Zmieniamy też:

@@ Linia 187: / Linia 187: @@
 import numpy as np
 import matplotlib.pyplot as plt
-import matplotlib.mlab as mlab
+from scipy.stats import norm
 mi, sigma = 100, 15
@@ Linia 209: / Linia 209: @@
 # zwróconych przez plt.hist w macierzy bins
-y = mlab.normpdf( bincenters, mi, sigma)
+y = norm.pdf( bincenters, mi, sigma)
 # obliczamy wartości w normalnym rozkładzie gęstości prawdopodobieństwa
 # o średniej mi i wariancji sigma**2 dla wartości bincenters

@@ Linia 348: / Linia 348: @@
 ----
-[[PPy3/NumPy|poprzednie]] | [["Programowanie z Pythonem3"|strona główna]] | [[PPy3/Wprowadzenie|od początku ;)]]
+[[PPy3/NumPy|poprzednie]] | [["Programowanie z Pythonem3"|strona główna]] | [[PPy3/TematyDodatkowe]]
 --[[Użytkownik:RobertJB|RobertJB]] ([[Dyskusja użytkownika:RobertJB|dyskusja]]) 14:05, 10 kwi 2017 (CEST)

@@ Linia 329: / Linia 329: @@
 ==Ćwiczenia==
-Spróbuj znaleźć sposób przekształcenia obrazka takiego, jak poniższy portret autora skryptu:
+. Spróbuj znaleźć sposób przekształcenia obrazka takiego, jak poniższy portret autora skryptu:
 [[Grafika:rjb_small.jpg]]
@@ Linia 337: / Linia 337: @@
 [[Grafika:rjb_bw.png]]

@@ Linia 237: / Linia 237: @@
 X = np.cos(Φ)
 Y = np.sin(2*Φ) / 2 # = np.sin(Φ) * np.cos(Φ)
-plt.axes(aspect='equal')
+plt.axes(aspect='equal') # to mi gwarantuje, że okrąg nie wyjdzie „spłaszczony"
 plt.plot(X, Y, 'b')
 YO = np.sin(Φ)

← poprzednia wersja		Wersja z 16:26, 6 lut 2018
Linia 222:		Linia 222:
	[[Grafika:Hist2.gif]]		[[Grafika:Hist2.gif]]

		+	==Krzywa parametryczna==
		+
		+	Rysowanie krzywej danej w postaci parametrycznej niewiele się różni od rysowania wykresu funkcji. Żeby było zabawniej, wykonam jeszcze przy okazji rachunek Monte Carlo oraz jego wizualizację.
		+
		+	Krzywa, jaką tutaj narysujemy to tzw. [https://pl.wikipedia.org/wiki/Lemniskata_Gerona lemniskata Gerona]. Zostanie dodatkowo opisana okręgiem jednostkowym; następnie, wylosuję dużą liczbę punktów o równomiernym rozkładzie leżących w kwadracie opisanym na tym okręgu. Robię tak, ponieważ łatwo mi jest wylosować tablicę wielu liczb z zakresu (0, 1) — jest do tego gotowa funkcja <tt>numpy.random.random</tt>, następnie przeskalować i przesunąć je do odcinka (-1, 1). Wreszcie, traktując je jako pary współrzędnych <tt>(x, y)</tt>, pozostawiam tylko te, które znalazły się we wnętrzu koła jednostkowego (odrzucając pozostałe), i dzielę na te, które znalazły się w obszarze wewnątrz lemniskaty, i te, które znalazły się na zewnątrz.
		+
		+	Intuicyjnie jest dość jasne, że — o ile losowaliśmy punkty rzeczywiście w sposób równomierny — stosunek liczby „trafień" w obszar ograniczony lemniskatą do liczby wszystkich „trafień" w obszar koła będzie w przybliżeniu równy stosunkowi pól powierzchni tych obszarów. I na tym tu poprzestaniemy; to, jak dokładne to może być przybliżenie, czyli jakiej wielkości błędu możemy oczekiwać, to już problem na inne zajęcia...
		+
		+	<source lang=python>
		+	import numpy as np
		+	import matplotlib.pyplot as plt
		+
		+	Φ = np.arange(0, 2np.pi, 2np.pi/600)
		+	X = np.cos(Φ)
		+	Y = np.sin(2Φ) / 2 # = np.sin(Φ) np.cos(Φ)
		+	plt.axes(aspect='equal')
		+	plt.plot(X, Y, 'b')
		+	YO = np.sin(Φ)
		+	plt.plot(X, YO, 'k')
		+	# biorę całą kupę losowych punktów z przedziału (-1, 1)
		+	P = np.random.random((100000, 2)) * 2 - 1 # milion par (x, y)
		+	# wyrzucam te, które są poza kołem x2 + y2 < 1
		+	P = P[(P**2).sum(axis=1) < 1]
		+	# dzielę na te co wewnątrz lemniskaty, i te co na zewnątrz
		+	I_in = np.abs(P[:,1]) < np.abs(P[:,0]) * np.sqrt(1 - P[:,0]**2)
		+	P_in = P[I_in]
		+	P_out = P[~I_in]
		+	# maluję wewnętrzne na niebiesko
		+	plt.plot(P_in[:,0], P_in[:,1], 'b,')
		+	# a zewnętrzne na czerwono
		+	plt.plot(P_out[:,0], P_out[:,1], 'r,')
		+	# powierzchnia to z grubsza stosunek l. punktów wewnętrznych do wszystkich
		+	S = P_in.size / P.size
		+	plt.title('Pole pow. wnętrza lemniskaty to ok. {:.3f} pola pow. koła'.format(S))
		+	plt.show()
		+	</source>
		+
		+	Przy okazji, jak widać w kodzie Pythona (wersja 3) można bezproblemowo używać liter z innych alfabetów, nie tylko łacińskiego. Niestety nie wszystkie systemy operacyjne zapewniają bezproblemową obsługę pełnego zestawu liter wszelkich języków, dlatego na ogół wykorzystywanie liter z egzotycznych systemów pisma nie jest szczególnie zalecane.
		+
		+	W wyniku uruchomienia tego programu powinno pojawić się okno z obrazkiem podobnym do poniższego:
		+
		+	[[File:Lemniskata.png]]

	==Wizualizacja zawartości tablicy dwuwymiarowej==		==Wizualizacja zawartości tablicy dwuwymiarowej==

← poprzednia wersja		Wersja z 12:44, 30 cze 2017
Linia 289:		Linia 289:
	Spróbuj znaleźć sposób przekształcenia obrazka takiego, jak poniższy portret autora skryptu:		Spróbuj znaleźć sposób przekształcenia obrazka takiego, jak poniższy portret autora skryptu:

−	[[Grafika:rjb_small.jpg]] do postaci przypominającej [[Grafika:rjb_bw.png]]	+	[[Grafika:rjb_small.jpg]]
		+
		+	do postaci przypominającej
		+
		+	[[Grafika:rjb_bw.png]]

← poprzednia wersja		Wersja z 14:03, 21 kwi 2017
Linia 289:		Linia 289:
	Spróbuj znaleźć sposób przekształcenia obrazka takiego, jak poniższy portret autora skryptu:		Spróbuj znaleźć sposób przekształcenia obrazka takiego, jak poniższy portret autora skryptu:

−	[[Grafika:rjb_small.jpg]]	+	[[Grafika:rjb_small.jpg]] do postaci przypominającej [[Grafika:rjb_bw.png]]
−
−	do postaci przypominającej
−
−	[[Grafika:rjb_bw.png]]

← poprzednia wersja		Wersja z 14:00, 21 kwi 2017
Linia 286:		Linia 286:

	==Ćwiczenia==		==Ćwiczenia==
		+
		+	Spróbuj znaleźć sposób przekształcenia obrazka takiego, jak poniższy portret autora skryptu:
		+
		+	[[Grafika:rjb_small.jpg]]
		+
		+	do postaci przypominającej
		+
		+	[[Grafika:rjb_bw.png]]
		+
		+

	----		----