{ "cells": [ { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "from scipy.optimize import curve_fit\n", "import statsmodels.api as sm\n", "from astropy.table import Table\n", "import os,errno\n", "import math \n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "data = pd.read_excel('DATA2.xlsx')\n", "#data = Table.read('measurements.raw.LD.Flow.LD.Occupancy.germany.darmstadt.K41D21.fits', format = 'fits')\n", "#initial_params_Proposed = [6000, 1.5, 4]\n", "#initial_params_Propose = [6000, 1.5, 4, 5]\n", "#initial_params_Prop = [6000, 1.5, 4]\n", "#initial_params_Wang = [6000, 28.2, 0.0069, 0.0028, 0.044]\n", "#initial_params_Edie_Multi_Regime = [6000, 20, 0.05, 1.5, 0.25]\n", "#initial_params_Cheng = [6000, 0.05, 4]\n", "#initial_params_Gaddam = [6000, 1.5, 0.05, 0.6]\n", "#initial_params_Modified_Lee = [6000, 1.5, 10.3, 2.14, 4]\n", "#initial_params_Drake_Two_Regime = [6000, 1.5, 0.25, 5, -15]\n", "#initial_params_van_Aerde = [2000, 0.007, 0.05, 0.5]\n", "\n", "initial_params_Proposed = [100, 160, 4]\n", "#initial_params_Propose = [100, 160, 4,5]\n", "#initial_params_Prop = [100, 160, 4]\n", "initial_params_Wang = [104, 9, 17, 2.1, 0.07]\n", "#initial_params_Edie_Multi_Regime = [104, 65, 30, 175, 20]\n", "initial_params_Cheng = [105, 65, 4]\n", "initial_params_Gaddam = [105, 179, 65, 0.6]\n", "initial_params_Modified_Lee = [105, 179, 10.3, 2.14, 4]\n", "#initial_params_Drake_Two_Regime = [105, 80, 25, 5, -15]\n", "initial_params_van_Aerde = [1000, 0.007, 0.05, 0.5]\n", "initial_params_Kucharski_and_Drabicki = [107.24, 77.94, 8.46, 2.72]\n", "initial_params_MacNicholas = [100, 170, 0.5, 5]\n", "initial_params_Ghandehari_and_Ardekani = [160, 70, 20]\n", "initial_params_Bando = [100, 3, 2]\n", "initial_params_Boardman_and_Lave = [100, 0.3, 0.2]\n", "initial_params_Lee = [115.5, 151.5, 24.87, 2.1]\n", "initial_params_Kerner_and_Konhauser = [105, 40, 0.5, 0.2, 0.8]\n", "initial_params_Del_Castillo_and_Benites_1 = [105, 170, 18]\n", "initial_params_Del_Castillo_and_Benites_2 = [105, 170, 18]\n", "initial_params_Newell_and_Franklin = [100, 160, 2]\n", "initial_params_Papageorgiou = [105, 45, 10]\n", "initial_params_Drake = [105, 45]\n", "initial_params_Underwood = [105, 45]\n", "initial_params_Greenberg = [60, 170]\n", "initial_params_May_and_Keller = [105, 179, 10.3, 2.14]\n", "initial_params_Pipes = [105, 179, 10.3]\n", "initial_params_Drew = [105, 179, 2.14]\n", "initial_params_Greenshields = [105, 179]\n", "\n", "\n", "#initial_params_Yaks = [100, 160, 4] \n", "#data = data[data['OCCUPANCY']> 0] \n", "#data = data[data['FLOW']> 0]\n", "#data = data[data['ERROR_FLAG']<= 2]\n", "\n", "# Transform data into arrays\n", "k_data = np.array(data['DENSITY'])\n", "q_data = np.array(data['FLOW'])\n", "v_data = np.array(data['SPEED'])\n", "\n", "\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", " | FLOW | \n", "DENSITY | \n", "FLOW.1 | \n", "SPEED | \n", "Proposed model | \n", "g(k) | \n", "Proposed model.1 | \n", "
---|---|---|---|---|---|---|---|
0 | \n", "992.0 | \n", "9.719537 | \n", "992.0 | \n", "102.06247 | \n", "103.880665 | \n", "0.246742 | \n", "103.880665 | \n", "
1 | \n", "986.0 | \n", "9.398039 | \n", "986.0 | \n", "104.91551 | \n", "104.047381 | \n", "0.239823 | \n", "104.047381 | \n", "
2 | \n", "976.0 | \n", "9.260186 | \n", "976.0 | \n", "105.39745 | \n", "104.114453 | \n", "0.236833 | \n", "104.114453 | \n", "
3 | \n", "860.0 | \n", "8.268004 | \n", "860.0 | \n", "104.01543 | \n", "104.524115 | \n", "0.214899 | \n", "104.524115 | \n", "
4 | \n", "916.0 | \n", "8.836977 | \n", "916.0 | \n", "103.65536 | \n", "104.304456 | \n", "0.227567 | \n", "104.304456 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
44782 | \n", "940.0 | \n", "9.171828 | \n", "940.0 | \n", "102.48775 | \n", "104.156083 | \n", "0.234909 | \n", "104.156083 | \n", "
44783 | \n", "1088.0 | \n", "10.718139 | \n", "1088.0 | \n", "101.51016 | \n", "103.265398 | \n", "0.267765 | \n", "103.265398 | \n", "
44784 | \n", "1038.0 | \n", "10.199983 | \n", "1038.0 | \n", "101.76487 | \n", "103.603605 | \n", "0.256945 | \n", "103.603605 | \n", "
44785 | \n", "1128.0 | \n", "10.986218 | \n", "1128.0 | \n", "102.67410 | \n", "103.073674 | \n", "0.273291 | \n", "103.073674 | \n", "
44786 | \n", "1000.0 | \n", "9.562135 | \n", "1000.0 | \n", "104.57915 | \n", "103.964111 | \n", "0.243364 | \n", "103.964111 | \n", "
44787 rows × 7 columns
\n", "