αα ααααααΈ 11 ααααααΆ ααααΆα 1914 ααα ααααΈαααααααΆααααααα Henry Ford ααΆααααα αΆααααα½ααα αααα»αααΆααα New York Times:
βαααα»ααααααΉαααΆ αααα»ααααααααα½αααααΆα ααΎαααΉαα αΆααααααΎαααα·αααααααα’αααα·αααΈα αααα»ααα·αα αΌαα α·ααααα·ααΆαααΏααααααΆααααααΆαααΆααα»ααα ααα»αααααααα»αα ααααααΆααα’αααααΈαααααΆααααααααα»αα ααΆααα·αααΊααΆ ααα Edison αα·ααααα»αααΆαααααΎααΆαααΆα αααΎαααααΆαααΎααααΈαααααΎαααααααα’αααα·αααΈαααααΆαααααααα αα·αααΆαααααααα αα½αααααααΌαααΆααααααΎαα‘αΎαααΆααΆααα·ααααα α αΎαααΎααααα α·ααααααααααΌααα ααΆααααΆαααααααααΊααΆαααααααα αααα αΆβαααααβαααααΆααβααααααβα’αααα·αααΈβαα αΌαβααβαααβαααβαααβααΊβααΆαβαααααΎαβαααβααααααβααααΆαβαααβα’αΆα βααααΎαβααΆαβαααα»αβαααβα ααααΆαβααααΆαβαααβαα·αβα αΆαβααΆα αβαααα αΌαβαααα ααα Edison ααΆαβααΆαααααβαααβααααααβαα½ααααβα αΎαΒ»α
ααα»ααααααΆαα’αααΈαα½ααα»α...
Thomas Edison ααΆαα½α Detroit Electric
ααΆαααααα»αααααααΆααααααΊααΆααΆαααααα‘αΌααΈααααα’ααααααα»ααααααααα»α
ααΎαααΆαααΆαααααααα
α―ααΆ r αααααΆαααααα’ααααΆααααΎαααα
ααααααΈααααΆα αααααΆαααΆααΆαα·ααααααα - αααα»ααααααΆααααααααα αα
αα
αααα·ααααΆααααΈααΆααααααΏαααΉαααααααααΈααααΆα αααααΊα§. ααΆαααααΆααααααΆααα
αα
αααα·ααααΆαα½αααααααΆαααΆααΆαα
αΆααα’αΆαααααααααααΆααααα»αααααΆααααα
αααΎααααααΆαααααΆαααΆααααα½αααααααΆα K αααα»ααα·αααααΆα’αααΈαααααΆαα»ααααααΆαααΌαααΆαααααα
αα
αααα·ααααΆααααΈ α§. αα
αααααααΆααα K αα
αα
αααα·ααααΆααΉααα·αα’αΆα
α
αΆααααααΈααααΆαααΆααααΌαααΆαααααα»ααααααΉαα’αΆα
αααααααααααααΈααααΆααααααΆααΉααα½αα±ααα
αΆααα’αΆααααααααΆααα
αα
αααα·ααααΆαα»αα
ααα ααααΈαααααααΆααααααα αΆααΊααΎααααΈααααααααααΆαααΆααααααα αΆαααΆα ααααααΆαααααΈααΆααααααα»ααΆααααα’αα»ααααΆαα±ααααΎαααααΆααααααΈααΆαα’αα·ααααααααα§ααααΆα ααααααααααα’αααα·αααΈααΉαα’αααααααΎαα
- βααααΆαααΌαααβ ααΊααΆααααΆααααααΆαααααααΆαααααααααααΉαα’αααααααΎααααααααα ααΎαα·αααααααΉαααΆααααΌαααα’αα‘α·α αααα·α (P0=0,5, t=0);
- α’ααααΆααααΎαα ααααααΈααααΆα (r) ααααααα’αααα·αααΈα
αααα»αααααΈαααααΌααα·ααΆαααΆ:
- ααααααα’αααα·αααΈααΉαααααΆααααααΌαααΆαααααα»ααααααααααααΆααααΆαααΈαα
αα ααααΆααααα»α (ICE) α
ααααΈααΈααααΆα (K=1) ααααααααα»ααα·αααΎααααααααα·ααααααααΉαα’αα»ααααΆαα±αααααα
ααααΈααααΆαααααααααΉαα’αααααααΎααααααα
ααΈααααΆααααααΆαααααααααα»αααααα αα·αα§ααααααα·ααααα·αααααΌαααΆααααααααα·α αΆαααΆαα αααα ααααααααΌαααααααα α αΎααα αααα·αααΆαααΈααααΆααααααΆααααααααα’αααα·αααΈαα αααα»αα§ααααΆα αααααααααα
- α₯α‘αΌααααααΎαααααα αααα»α βαααααααΆα’αα·αααααΆαβ (P(t)<0) α αΎααα αααα»ααα»αααΆα ααΎαααΉαααααΎα’α»α αααα·αααΆαααααα ααΉα βααααΆαααΌαααβ αααααΆαααααααααΆααααααΎα (t-t0)α
αααα·αα·ααααα·ααΆααααααααααααΉαα’αααααααΎαααΊααα
ααααΈ
αααα·αα·αααααααααα’αααα·αααΈααααΈ
αααα·αα·αα»αααααΆα 2012 ααΎααΆααααα·αα’αααα·αααΈααΊααααααΆαα α αΎαααΉααα·αααααΌαααΆααααααααα·α αΆαααΆαααα»αααΆααα·ααααΆαααααα
ααΆαααααα ααΎαααΆααα·ααααααααΌα ααΆααααααα
αααααα·ααΈααΎααααΈααααααα Year Zero αα·αα’ααααΆααααΎαααΈααααΆα
import matplotlib.pyplot as plt
import numpy as np
import math
x = np.linspace(2012, 2019, 8)
y1 = np.array([60936407, 63429200, 65708230, 66314155, 69464432, 70694834, 68690468, 64341693]) # ΠΊΠΎΠ»-Π²ΠΎ ΠΏΡΠΎΠΈΠ·Π²Π΅Π΄Π΅Π½Π½ΡΡ
Π»Π΅Π³ΠΊΠΎΠ²ΡΡ
ΠΌΠ°ΡΠΈΠ½
y2 = np.array([52605, 97507, 320713, 550297, 777495, 1227117, 2018247, 1940147]) # ΠΊΠΎΠ»-Π²ΠΎ ΠΏΡΠΎΠΈΠ·Π²Π΅Π΄Π΅Π½Π½ΡΡ
Π»Π΅Π³ΠΊΠΎΠ²ΡΡ
ΡΠ»Π΅ΠΊΡΡΠΎΠΌΠΎΠ±ΠΈΠ»Π΅ΠΉ
y = y2/y1 #Π΄ΠΎΠ»Ρ ΡΠ»Π΅ΠΊΡΡΠΎΠΌΠΎΠ±ΠΈΠ»Π΅ΠΉ Π² ΠΎΠ±ΡΠ΅ΠΌ ΠΏΡΠΎΠΈΠ·Π²ΠΎΠ΄ΡΡΠ²Π΅ Π°Π²ΡΠΎΠΌΠΎΠ±ΠΈΠ»Π΅ΠΉ
ymax=1 #ΠΏΠ΅ΡΠ²ΠΎΠ½Π°ΡΠ°Π»ΡΠ½ΠΎΠ΅ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ΅ ΠΎΡΠΊΠ»ΠΎΠ½Π΅Π½ΠΈΠ΅ ΡΡΠ°ΡΠΈΡΡΠΈΡΠ΅ΡΠΊΠΈΡ
Π΄Π°Π½Π½ΡΡ
ΠΎΡ Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΡΡΠ½ΠΊΡΠΈΠΈ
Gmax=2025 #Π³ΠΎΠ΄ Π΄Π»Ρ Π½Π°ΡΠ°Π»ΠΎ ΠΏΠΎΠΈΡΠΊΠ° "Π½ΡΠ»Π΅Π²ΠΎΠ³ΠΎ Π³ΠΎΠ΄Π°"
rmax=0.35 #Π½Π°ΡΠ°Π»ΡΠ½ΡΠΉ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½Ρ
k=1 #ΠΏΡΠΈΠ½ΡΡ "1" ΠΈΠ· ΠΏΡΠ΅Π΄ΠΏΠΎΡΡΠ»ΠΊΠΈ, ΡΡΠΎ ΡΠ»Π΅ΠΊΡΡΠΎΠΌΠΎΠ±ΠΈΠ»ΠΈ ΠΏΠΎΠ»Π½ΠΎΡΡΡΡ Π·Π°ΠΌΠ΅Π½ΡΡ Π»Π΅Π³ΠΊΠΎΠ²ΡΠ΅ Π°Π²ΡΠΎΠΌΠΎΠ±ΠΈΠ»ΠΈ Ρ ΠΠΠ‘
p0=0.5 # ΠΏΡΠΎΡΠ΅Π½Ρ ΡΡΠ½ΠΊΠ° Π² "Π½ΡΠ»Π΅Π²ΠΎΠΉ Π³ΠΎΠ΄"
for j in range(10): # ΡΠΈΠΊΠ» ΠΏΠ΅ΡΠ΅Π±ΠΎΡΠ° "Π½ΡΠ»Π΅Π²ΡΡ
Π³ΠΎΠ΄ΠΎΠ²"
x0=2025+j
r=0.35
for i in range(10): # ΡΠΈΠΊΠ» ΠΏΠ΅ΡΠ΅Π±ΠΎΡΠ° ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠ° Π² ΠΊΠ°ΠΆΠ΄ΠΎΠΌ "Π½ΡΠ»Π΅Π²ΠΎΠΌ Π³ΠΎΠ΄Ρ"
r=0.25+0.02*i
y4=k*p0*math.e**(r*(x-x0))/(k+p0*(math.e**(r*(x-x0))-1))-y
# print(str(x0).ljust(20), str(r).ljust(20), max(abs(y4)))
if max(abs(y4))<=ymax: # ΠΏΠΎΠΈΡΠΊ ΠΌΠΈΠ½ΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ³ΠΎ ΠΈΠ· ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΡΡ
ΠΎΡΠΊΠ»ΠΎΠ½Π΅Π½ΠΈΠΉ Π²Π½ΡΡΡΠΈ ΠΊΠ°ΠΆΠ΄ΠΎΠ³ΠΎ Π³ΠΎΠ΄Π° ΠΏΡΠΈ ΠΊΠ°ΠΆΠ΄ΠΎΠΌ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠ΅ r
ymax=max(abs(y4))
Gmax=x0
rmax=r
print(str(Gmax).ljust(20), str(rmax).ljust(20), ymax) # Π²ΡΠ²ΠΎΠ΄ "Π½ΡΠ»Π΅Π²ΠΎΠ³ΠΎ Π³ΠΎΠ΄Π°", ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠ° r ΠΈ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ³ΠΎ ΠΈΠ· ΠΎΡΠΊΠ»ΠΎΠ½Π΅Π½ΠΈΠΉ ΠΎΡ ΡΡΠ½ΠΊΡΠΈΠΈ
ααΆαααααααααααααα·ααΈ αααααααΆααααααααααΌαααΆαααααΎαααΎαα
ααααΆαααΌαααααΊ 2028 α
αααα»αααααΎα - 0.37
αααααΆαα’αα·ααααΆαααα·αααααααααα·αα·ααΈααααααα»αααΆαααΊ 0.005255α
ααααΆα αααααα»αααΆααααΆαααααΆα 2012 αα·α 2019 ααΎααα ααΌα αααα
ααααΆα ααα
α»ααααααααΆαα½αααΉαααΆαααααΆαααααα αΌααααααααΆα 2050 ααΎααα
ααΌα
αααα
ααΆααΆααααα αΆαααΈααΆαααΆααααααΆα
α 99% ααααΈααααΆαααΆααααΌααααααΊα§α αα
ααααΆα 2040 ααααααα’αααα·αααΈααΉααααα½αααααααααΆαααααα»αααΆαα½αααΉααααΆαααΈαα
αα ααααΆααααα»αα
αααααα·ααΈααααΆα ααα·α
αα»αααΆα
import matplotlib.pyplot as plt
import numpy as np
import math
x = np.linspace(2012, 2019, 8)
y1 = np.array([60936407, 63429200, 65708230, 66314155, 69464432, 70694834, 68690468, 64341693])
y2 = np.array([52605, 97507, 320713, 550297, 777495, 1227117, 2018247, 1940147])
y = y2/y1
k=1
p0=0.5
x0=2028
r=0.37
y1=k*p0*math.e**(r*(x-x0))/(k+p0*(math.e**(r*(x-x0))-1))
#Π‘ΡΡΠΎΠΈΠΌ Π³ΡΠ°ΡΠΈΠΊ ΡΡΠ½ΠΊΡΠΈΠΈ Π½Π° ΠΎΡΡΠ΅Π·ΠΊΠ΅ ΠΌΠ΅ΠΆΠ΄Ρ 2012 ΠΈ 2019 Π³ΠΎΠ΄Π°ΠΌΠΈ
fig, ax = plt.subplots(figsize=(30, 20), facecolor="#f5f5f5")
plt.grid()
ax.plot(x, y, 'o', color='tab:brown')
ax.plot(x, y1)
#Π‘ΡΡΠΎΠΈΠΌ Π³ΡΠ°ΡΠΈΠΊ ΡΡΠ½ΠΊΡΠΈΠΈ Π½Π° ΠΎΡΡΠ΅Π·ΠΊΠ΅ ΠΌΠ΅ΠΆΠ΄Ρ 2010 ΠΈ 2050 Π³ΠΎΠ΄Π°ΠΌΠΈ
x = np.linspace(2010, 2050)
y2 = [k*p0*math.e**(r*(i-x0))/(k+p0*(math.e**(r*(i-x0))-1)) for i in x]
y3 = 0.99+0*x
fig, ax = plt.subplots(figsize=(30, 20), facecolor="#f5f5f5")
ax.set_xlim([2010, 2050])
ax.set_ylim([0, 1])
plt.grid()
plt.plot(x, y2, x, y3)
ααΆαααααΎα
αααααααΎααΆααααααα·ααααΆααΌα ααααΆααΉαααααα·αααααΆα’αααΈαααααααα·ααααΆαα’αα·αααααααααααααααααΆααααΆαααΈαα αα ααααΆααααα»α αααα»αααΆαααααΆααΆαααααααααΆαααΈααΆαα’αα·ααααααααα§ααααΆα ααααααααααα’αααα·αααΈααΉαα’αααααααΎααααααα’ααααΎαα·αααααααααα·αα·αααααΆαα
αααααααααααα½αααΆααααα αΆαααΆαα ααααΆα 2030 ααΆαααααααΆαααααααααααΉαα’αααααααΎααααααααα ααΎαα·αααααααΉαααΆααααΌαααα’αα‘α·α αααα·α α αΎααα ααααΆα 2040 ααααααααΉαα’αααααααΎααααααΆααααΆαααΈαα αα ααααΆααααα»αααΉαααααΆαααΆααΏαα’ααΈαααΆαα
ααΆααΆααα·αααΆαααααααΆααααΈααααΆα 2030 ααα»ααααα½αα
ααα½αααΉαααΎαααααααααααΆααααααα½αααααΆααα·ααα»αααααΆα 2030 ααα»αααααα½αααααΉαααΉαααΆααΆααα·ααααααΆαααααααα½αααααΉαααααΆαααΆααααααα’αααα·αααΈα
α’ααααΆααααΎααααααΆααααααααα’αααα·αααΈααΊαααααααΆαα’ααααΆααααΎααααααΆαααααααααααααΆααααΆαααΈαα
αα ααααΆααααα»α 4 αα ααααααα αΆαααΆαα
αα
αααα·ααααΆααααΈααααα»αα
αΌααααααα»αααΈαα·αααααααΎαααΆααααααΏα α αΎαααααΆαααΆααααααα½αααααααΆααααααΈαα·ααααα
αΆαααααααααααΎα (αα
ααΈαααααΎαα
αα
αΆαα’αααΈααΌααααααα) .
αααα»αααααααααα»ααααΆαααααΆαααΆααα»α αααα αΆααα Edison αα·αα’αΆα αααααααΆαααΆα αα½αααααααΌαααΆααααααααΆα αααααΊααααααααΆααααααααΆααααααααααΆαα αααααΉαα’αα»ααααΆαα±ααααΆαα ααααααααααΌααααΆαααααΆααΈαααΆααααα
ααΎααααΈαααααΎααααααΆαααααΆααΈααααΆαααααααααααΎααΉααααααΆαααααΆααΈαααααααα₯ααααααααααΆαααααΆαα α αΆαααΆα αααααΌαααααΎααααΎααααααααααΆαα’αααα·αααΈαααααΆαααααΆαααα αααα»αααΈαααα»αααα αα·αααΆαααααααααααΌαα αΆαααα
αα·αααααα»ααααα ααααΎαααααΆααααααααααα’αααα·αααΈααΉαααααΌαααΆααΆαα
PS
ααααα·αααΎ Edison α’αΆα
αααααααΆααααα αΆαααααΆαααααααα±ααααΆααααα "α’αΆαα»ααααα" ααΉααα·αααΆαα
αΆααααααΎαααΌααααΈαα ...
ααΆαααα’αααααααΎααααΆαααααααΆαα
α»ααααααααα»ααααααααα’αΆα
α
αΌααα½ααααα»αααΆααααααααα·αααα
ααΎαα αααααΆαααααα»ααααααααααααΆααΉαααΎαααααααααα’αααα·αααΈ?
-
9,5%αα ααααΆα 2030 α’ααααααααααααΆααΉαααααΌααα ααααΎααααααα’αααα·αααΈ αα·ααααααΆαααααααΆα 18 αααααα
-
20,0%αα ααααΆα 2040 α’ααααααααααααΆααααΆααααΆααααΌααα ααααΎααααααα’αααα·αααΈ38
-
48,4%αα·αααΏαααΆα 2050
-
22,1%ααααααα’αααα·αααΈ ααΉααα·ααααα½αααααααααΆαα α€α’
α’αααααααΎααααΆαα 190 ααΆααααΆαααααααααα α’αααααααΎααααΆαα α‘α€ ααΆααααααΌαααΆαα αΆαααΆααα
ααααα: www.habr.com