Файл:Population curve of Uruk 2.svg
Исходный файл (SVG-файл, номинально 797 × 522 пкс, размер файла: 49 КБ)
Этот файл находится на Викискладе. Сведения о нём показаны ниже.
Викисклад — централизованное хранилище для свободных файлов, используемых в проектах Викимедиа.
Сообщить об ошибке с файлом |
Краткое описание
ОписаниеPopulation curve of Uruk 2.svg |
English: Estimated population growth of Uruk. |
Дата | |
Источник | Собственная работа |
Автор | Merikanto |
Местоположение камеры | 31° 19′ 20″ с. ш., 45° 38′ 10″ в. д. | Это и другие изображения по их местоположению на OpenStreetMap | 31.322222; 45.636111 |
---|
This plot is based on info from Finnish Wikipedia.
https://fi.wikipedia.org/wiki/Uruk
data
bc;est_pop;pop_lower;pop_upper
5000;0;0;0
4000;5000;5000;5000
3700;11000;8000;14000
3500;14500;14000;15000
3400;20000;20000;2000
3100;45000;40000;50000
2800;80000;80000;80000
2650;65000;50000;80000
2500;50000;50000;50000
1800;75000;75000;75000
Python 3 code
- uruk population Python 3 code
- ver 2
import matplotlib
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
from scipy.interpolate import interp1d
from scipy.optimize import curve_fit
import math as math
- sigmoid
def sigmoid(x, L ,x0, k, b):
y = L / (1 + np.exp(-k*(x-x0)))+b
return (y)
- exp
def func(x, a, c, d):
return a*np.exp(-c*x)+d
df = pd.read_csv(r'uruk3.csv', delimiter=';')
beginx=4200
endx=2800
xnum=20
xnum1=xnum+1
pxlen=100
- bc;est_pop;pop_lower;pop_upper
origx=np.array(df['bc'])
origy=np.array(df['est_pop'])
sekondx=np.linspace(beginx, endx,xnum1)
xdata0=np.linspace(0, xnum,xnum1)
interp1 = interp1d(origx,origy, kind='cubic')
ydata0 = interp1(sekondx)
- print (sekondx)
- rint (xdata0)
- print (ydata0)
xdata1=xdata0
ydata1=ydata0/100000.0
- print(ydata1)
- p0 = [max(ydata1), np.median(xdata1),1,min(ydata1)] # this is an mandatory initial guess
- popt, pcov = curve_fit(sigmoid, xdata1, ydata1,p0, method='dogbox')
popt, pcov = curve_fit(func, xdata1, ydata1, p0=(1, 1e-6, 1))
print ("Popt")
print (popt)
- print (pcov)
x=np.linspace(0, xnum,pxlen)
thirdx=np.linspace(beginx,endx,pxlen)
- y = sigmoid(x, *popt)
y = func(x, *popt)
y2=y*100000
lok1=0
lok2=pxlen-1
rx1=thirdx[lok1]
ry1=y2[lok1]
rx2=thirdx[lok2]
ry2=y2[lok2]
dx1=rx1-rx2 ## years bc
dy1=ry2-ry1
gr1=math.log(ry2/ry1)/dx1
gp1=gr1*100
dout1=math.log(2.0)/gr1
gp1=np.round(gp1,2)
dout1=np.round(dout1,0)
annot1="Population growth rate "+str(gp1)+ "%,\n doubling time "+str(dout1)+ " years"
- print(rx1,rx2, ry1,ry2)
- print(dx1,dy1)
print("Growth rate %", gp1)
- print(gp1)
print(" Doubling time yr",dout1)
- plt.plot(x,y, label='fit')
- plt.plot(xdata1,ydata1)
plt.plot(thirdx,y2, 'g-.', lw=3,label="Theoretical exponent curve")
plt.plot(origx,origy, 'b-', lw=5, label="Data from archaeologists")
plt.xlim(beginx,endx)
plt.xticks(fontsize=15)
plt.yticks(fontsize=15)
plt.xlabel("Years BC" , fontsize=15)
plt.ylabel("Population of Uruk" , fontsize=15 )
plt.title("Population growth of Uruk" , fontsize=18 )
plt.text(4000, 50000, annot1, fontsize=17)
plt.legend()
plt.grid(True)
plt.show()
Лицензирование
- Вы можете свободно:
- делиться произведением – копировать, распространять и передавать данное произведение
- создавать производные – переделывать данное произведение
- При соблюдении следующих условий:
- атрибуция – Вы должны указать авторство, предоставить ссылку на лицензию и указать, внёс ли автор какие-либо изменения. Это можно сделать любым разумным способом, но не создавая впечатление, что лицензиат поддерживает вас или использование вами данного произведения.
- распространение на тех же условиях – Если вы изменяете, преобразуете или создаёте иное произведение на основе данного, то обязаны использовать лицензию исходного произведения или лицензию, совместимую с исходной.
Элементы, изображённые на этом файле
изображённый объект
У этого свойства есть некоторое значение без элемента в
31°19'19.999"N, 45°38'10.000"E
16 июля 2020
История файла
Нажмите на дату/время, чтобы посмотреть файл, который был загружен в тот момент.
Дата/время | Миниатюра | Размеры | Участник | Примечание | |
---|---|---|---|---|---|
текущий | 15:59, 16 июля 2020 | 797 × 522 (49 КБ) | Merikanto | Uploaded own work with UploadWizard |
Использование файла
Нет страниц, использующих этот файл.
Глобальное использование файла
Данный файл используется в следующих вики:
- Использование в fi.wikipedia.org
Метаданные
Файл содержит дополнительные данные, обычно добавляемые цифровыми камерами или сканерами. Если файл после создания редактировался, то некоторые параметры могут не соответствовать текущему изображению.
Ширина | 637.92pt |
---|---|
Высота | 417.6pt |