224 lines
102 KiB
Plaintext
224 lines
102 KiB
Plaintext
|
{
|
||
|
|
"cells": [
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 4,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"import numpy as np\n",
|
||
|
|
"from matplotlib import pyplot as plt\n",
|
||
|
|
"import xarray as xr\n",
|
||
|
|
"from pathlib import Path\n",
|
||
|
|
"\n",
|
||
|
|
"dir_ = Path(globals()[\"_dh\"][0]) # __file__ does not exist in ipynb"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 49,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"def read_wpd_csv(\n",
|
||
|
|
" filepath: Path | str, x_coord: str, ds_coord: str | None = None\n",
|
||
|
|
") -> xr.DataArray:\n",
|
||
|
|
" filepath = Path(filepath)\n",
|
||
|
|
"\n",
|
||
|
|
" header = np.loadtxt(filepath, dtype=str, max_rows=2, delimiter=\",\")\n",
|
||
|
|
" header_values = header[0][::2]\n",
|
||
|
|
" try:\n",
|
||
|
|
" np.array(header_values, dtype=float)\n",
|
||
|
|
" except ValueError:\n",
|
||
|
|
" np.array(header_values)\n",
|
||
|
|
"\n",
|
||
|
|
" values = np.loadtxt(filepath, dtype=float, skiprows=2, delimiter=\",\")\n",
|
||
|
|
"\n",
|
||
|
|
" coords = dict()\n",
|
||
|
|
" dims = [x_coord]\n",
|
||
|
|
" if ds_coord is not None:\n",
|
||
|
|
" dims.append(ds_coord)\n",
|
||
|
|
" coords[ds_coord] = header_values\n",
|
||
|
|
"\n",
|
||
|
|
" for col in range(int(values.shape[1] / 2)):\n",
|
||
|
|
" x = values[:, col * 2]\n",
|
||
|
|
" if col == 0:\n",
|
||
|
|
" coords[x_coord] = x\n",
|
||
|
|
" else:\n",
|
||
|
|
" if np.any(x != coords[x_coord]):\n",
|
||
|
|
" raise ValueError(f\"All webplotdigitizer datasets must share x values\")\n",
|
||
|
|
"\n",
|
||
|
|
" da = xr.DataArray(\n",
|
||
|
|
" values[:, 1::2],\n",
|
||
|
|
" dims=dims,\n",
|
||
|
|
" coords=coords,\n",
|
||
|
|
" attrs=dict(filename=filepath.name),\n",
|
||
|
|
" )\n",
|
||
|
|
" return da"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 78,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"# LED properties\n",
|
||
|
|
"LEDS = {\n",
|
||
|
|
" \"APHHS1005F3C-70MAV\": dict(\n",
|
||
|
|
" emissivity=8e-3, # W/sr\n",
|
||
|
|
" emissivity_current=70e-3, # A\n",
|
||
|
|
" )\n",
|
||
|
|
"}\n",
|
||
|
|
"\n",
|
||
|
|
"# Photodiode properties\n",
|
||
|
|
"PHOTODIODES = {\n",
|
||
|
|
" \"PD15-22B/TR8\": dict(\n",
|
||
|
|
" dark_current=5e-9, # A\n",
|
||
|
|
" light_current=read_wpd_csv(\n",
|
||
|
|
" dir_ / \"data/PD15-22B/reverse_light_current.csv\",\n",
|
||
|
|
" \"irradiance\",\n",
|
||
|
|
" \"vr\",\n",
|
||
|
|
" ).isel(vr=0),\n",
|
||
|
|
" ),\n",
|
||
|
|
" \"TEMD5110X01\": dict(\n",
|
||
|
|
" dark_current=3e-9, # A\n",
|
||
|
|
" light_current=read_wpd_csv(\n",
|
||
|
|
" dir_ / \"data/TEMD5110X01/reverse_light_current.csv\",\n",
|
||
|
|
" \"irradiance\",\n",
|
||
|
|
" \"vr\",\n",
|
||
|
|
" ).isel(vr=0),\n",
|
||
|
|
" ),\n",
|
||
|
|
"}"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 83,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"text/plain": [
|
||
|
|
"Text(0.5, 1.0, 'Photodiode Transfer Function')"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"execution_count": 83,
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "execute_result"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB71klEQVR4nO3dd3hUVf7H8ffMZNJJIISE3mvo0kSp0ptiQ10LoOLqBncVXVcsKFhwrexqlLUA/tZ1F1FRFJWuKKA06RC6KCUJLSF9MnN/f1wzYUgmJJAwk+Tzep55uHPOmXu/M+ROvjn3nnMshmEYiIiIiEiFZ/V1ACIiIiJSNpTYiYiIiFQSSuxEREREKgkldiIiIiKVhBI7ERERkUpCiZ2IiIhIJaHETkRERKSSUGInIiIiUkkosRMRERGpJJTYiVQwc+bMwWKxsH79el+Hcl4HDx7EYrEwZ84cd9nTTz+NxWIp0+M0btyYcePGlek+K5qXXnqJpk2bYrPZ6NSpk6/D8TvffvstFouFb7/91tehiJQrJXYifiI/Yct/BAcH07JlSyZOnEhSUlK5HvvDDz9kxowZ5XqMiuDc/wNvj8aNG/s6VA+LFy/mkUce4corr2T27Nk8//zz5Xq8/OS8qMfMmTPL9djn8+abb3r8ISFS1QT4OgAR8TRt2jSaNGlCdnY2P/zwA2+99RZfffUV27ZtIzQ0tFyO+eGHH7Jt2zYeeOCBctn/2Z544gkeffTRcj/OhejTpw///ve/Pcruvvtuunfvzj333OMuCw8Pv9ShFWv58uVYrVbee+89AgMDL9lx33rrrUKfRY8ePS7Z8Yvy5ptvEh0dXagHt0+fPmRlZV3Sz0fEF5TYifiZYcOG0bVrV8BMKmrWrMmrr77K559/zi233OLj6C5eQEAAAQH++dXTtGlTmjZt6lF277330rRpU2677Tavr8vLy8PlcvksaUhOTiYkJKRMj5+ZmXnePyRuuOEGoqOjy+yY5clqtRIcHOzrMETKnS7Fivi5q666CoADBw54lOfk5DBp0iRq1apFWFgY1157LSkpKYVe/+abb9K2bVuCgoKoW7cu8fHxnD592l3fr18/Fi5cyC+//FLkpcbk5GTuuusuYmNjCQ4OpmPHjrz//vuFjnP69GnGjRtHZGQk1atXZ+zYsR7HyVfUPXZ5eXk888wzNGvWjKCgIBo3bsxjjz1GTk6ORzvDMHj22WepX78+oaGh9O/fn+3btxf5uZ0+fZoHHniABg0aEBQURPPmzfn73/+Oy+Uqsn1J5d83+PLLLzNjxgx3zDt27CA3N5cpU6bQpUsXIiMjCQsLo3fv3qxYscLrPt5++233Prp168a6des82h47dozx48dTv359goKCqFOnDtdccw0HDx4EwGKxMHv2bDIyMtz/f2dfivzggw/o0qULISEhREVFcfPNN/Prr796HKNfv360a9eODRs20KdPH0JDQ3nssccu+jMq6pKoxWLh6aefdj/P/3nYu3cv48aNo3r16kRGRjJ+/HgyMzMLvf6DDz6ge/fuhIaGUqNGDfr06cPixYsB817L7du3891337k/i379+gHe77GbN2+e+/OJjo7mtttu4/Dhwx5txo0bR3h4OIcPH2b06NGEh4dTq1YtHn74YZxO5wV/TiLlwT//bBYRt3379gFQs2ZNj/L777+fGjVq8NRTT3Hw4EFmzJjBxIkTmTt3rrvN008/zdSpUxk4cCD33XcfiYmJvPXWW6xbt45Vq1Zht9t5/PHHSU1N5bfffuO1114DCi41ZmVl0a9fP/bu3cvEiRNp0qQJ8+bNY9y4cZw+fZq//OUvgJlwXXPNNfzwww/ce++9tGnThvnz5zN27NgSvce7776b999/nxtuuIGHHnqIn376ienTp7Nz507mz5/vbjdlyhSeffZZhg8fzvDhw9m4cSODBw8mNzfXY3+ZmZn07duXw4cP88c//pGGDRuyevVqJk+ezNGjR8vkfsLZs2eTnZ3NPffcQ1BQEFFRUaSlpfHuu+9yyy23MGHCBM6cOcN7773HkCFDWLt2baFBDR9++CFnzpzhj3/8IxaLhRdffJHrrruO/fv3Y7fbAbj++uvZvn07999/P40bNyY5OZklS5Zw6NAhGjduzL///W/efvtt1q5dy7vvvgvAFVdcAcBzzz3Hk08+yZgxY7j77rtJSUnh9ddfp0+fPvz8889Ur17dHcuJEycYNmwYN998M7fddhuxsbHn/QxOnjzp8dxms1GjRo0L+jzHjBlDkyZNmD59Ohs3buTdd98lJiaGv//97+42U6dO5emnn+aKK65g2rRpBAYG8tNPP7F8+XIGDx7MjBkzuP/++wkPD+fxxx8HKPZ9zJkzh/Hjx9OtWzemT59OUlIS//jHP1i1alWhz8fpdDJkyBB69OjByy+/zNKlS3nllVdo1qwZ99133wW9Z5FyYYiIX5g9e7YBGEuXLjVSUlKMX3/91fjf//5n1KxZ0wgJCTF+++03j3YDBw40XC6X+/UPPvigYbPZjNOnTxuGYRjJyclGYGCgMXjwYMPpdLrbvfHGGwZgzJo1y102YsQIo1GjRoVimjFjhgEYH3zwgbssNzfX6NmzpxEeHm6kpaUZhmEYn332mQEYL774ortdXl6e0bt3bwMwZs+e7S5/6qmnjLO/ejZt2mQAxt133+1x7IcfftgAjOXLl3u8nxEjRni878cee8wAjLFjx7rLnnnmGSMsLMzYvXu3xz4fffRRw2azGYcOHSr0Xr0JCwvz2PeBAwcMwIiIiDCSk5M92ubl5Rk5OTkeZadOnTJiY2ONO++8s9A+atasaZw8edJd/vnnnxuA8cUXX7hfCxgvvfRSsTGOHTvWCAsL8yg7ePCgYbPZjOeee86jfOvWrUZAQIBHed++fQ3AmDlzZrHHyZf/f3juI/9nKP/9nf3/ng8wnnrqqUL7OvvzMQzDuPbaa42aNWu6n+/Zs8ewWq3Gtdde6/HzbBiGx89D27Ztjb59+xY67ooVKwzAWLFihWEY5s9xTEyM0a5dOyMrK8vd7ssvvzQAY8qUKe6ysWPHGoAxbdo0j3127tzZ6NKlS5GfkYiv6FKsiJ8ZOHAgtWrVokGDBtx8882Eh4czf/586tWr59Hunnvu8bik2bt3b5xOJ7/88gsAS5cuJTc3lwceeACrteBUnzBhAhERESxcuPC8sXz11VfUrl3b494+u93On//8Z9LT0/nuu+/c7QICAjx6Lmw2G/fff3+JjgEwadIkj/KHHnoIwB1n/vu5//77Pd53UQM+5s2bR+/evalRowbHjx93PwYOHIjT6WTlypXnjet8rr/+emrVquVRZrPZ3Pe5uVwuTp48SV5eHl27dmXjxo2F9nHTTTd59HD17t0bgP379wO475v79ttvOXXqVKni+/TTT3G5XIwZM8bjM6hduzYtWrQodHk4KCiI8ePHl+oYn3zyCUuWLHE//vOf/5Tq9We79957PZ737t2bEydOkJaWBsBnn32Gy+ViypQpHj/PwAVNn7N+/XqSk5P505/+5HHv3YgRI2jdunWR50dRMeb/X4n4C12KFfEzCQkJtGzZkoCAAGJjY2nVqlWhX2QADRs29HienyDkJwD5CV6rVq082gUGBtK0aVN3fXF++eUXWrRoUej4bdq08TjGL7/8Qp06dQqNkDz32N6OYbVaad68uUd57dq1qV69uscxAFq0aOHRrlatWoUu/+3Zs4ctW7YUSrzyJScnnzeu82nSpEmR5e+//z6vvPIKu3btwuFwFNv+fP+HQUFB/P3vf+ehhx4iNjaWyy+/nJEjR3LHHXdQu3btYuPbs2cPhmEU+rzy5V/qzVevXr1SD77o06dPmQ2eKO6ziIiIYN++fVitVuLi4srkeN7OD4DWrVvzww8/eJQFBwcX+nmqUaNGqRNukfKmxE7Ez3Tv3t09KrY4NputyHLDMMo6pEuiLCctdrlcDBo0iEceeaTI+pYtW170MUJCQgqVffDBB4wbN47Ro0fz17/+lZiYGGw2G9OnT3ffK3m2kvwfPvDAA4waNYrPPvuMRYsW8eSTTzJ9+nS
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 640x480 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"fig, ax = plt.subplots(tight_layout=True)\n",
|
||
|
|
"for diode, properties in PHOTODIODES.items():\n",
|
||
|
|
" lc: xr.DataArray = properties[\"light_current\"]\n",
|
||
|
|
" lines = ax.loglog(lc.coords[\"irradiance\"], lc, label=diode)\n",
|
||
|
|
"\n",
|
||
|
|
" ax.axhline(properties[\"dark_current\"], color=lines[0].get_color(), linestyle=\"--\")\n",
|
||
|
|
" ax.axvline(\n",
|
||
|
|
" np.max(lc.coords[\"irradiance\"]), color=lines[0].get_color(), linestyle=\"--\"\n",
|
||
|
|
" )\n",
|
||
|
|
"\n",
|
||
|
|
"ax.set_xlabel(\"Irradiance [$W/m^2$]\")\n",
|
||
|
|
"ax.set_ylabel(\"Current [A]\")\n",
|
||
|
|
"ax.grid(True)\n",
|
||
|
|
"ax.legend(loc=\"upper left\")\n",
|
||
|
|
"ax.set_title(\"Photodiode Transfer Function\")"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 84,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"text/plain": [
|
||
|
|
"Text(0.5, 1.0, 'Component comparison')"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"execution_count": 84,
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "execute_result"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAChNElEQVR4nOzdd1gU19cH8O/sLrv03kQRVBREERAbIthQxF5jjAWxJQYrMbYk9l80mihGiS2CLSrR2GJBEQuoWBEbFlAUCyCKgtSF3fv+4cvGlSKwi0s5n+fZ53Fnzt45O1yHw8zcOxxjjIEQQgghhFR7PFUnQAghhBBClIMKO0IIIYSQGoIKO0IIIYSQGoIKO0IIIYSQGoIKO0IIIYSQGoIKO0IIIYSQGoIKO0IIIYSQGoIKO0IIIYSQGoIKO0IIIYSQGoIKO0IIISpjbW2N0aNHqzoNQmoMKuwIqYUePnyIr7/+Gg0bNoS6ujp0dXXh5uaG1atXIycnR9XpVXs7d+5EQECAqtMghNRCHD0rlpDa5ciRIxgyZAhEIhFGjRqF5s2bQywW49y5c/jnn38wevRobNy4UdVpVmu9e/fG7du38fjxY1WnUuXl5eWBx+NBTU1N1akQUiMIVJ0AIeTzSUhIwJdffgkrKyucOnUKderUka3z8/NDfHw8jhw5osIMSW3AGENubi40NDQgEolUnQ4hNQpdiiWkFlm+fDkyMzOxefNmuaKukI2NDaZOnSp7X1BQgMWLF6NRo0YQiUSwtrbG3LlzkZeXJ/c5a2tr9O7dG2fOnEGrVq2goaEBBwcHnDlzBgCwb98+ODg4QF1dHS4uLrh+/brc50ePHg1tbW08evQIXl5e0NLSgoWFBRYtWoSPLypkZWXhu+++g6WlJUQiEWxtbfHrr78WieM4DpMmTcKBAwfQvHlziEQiNGvWDKGhoUW+9/PnzzFmzBiYmZnJ4oKCguRizpw5A47j8Pfff+N///sf6tWrB3V1dXTt2hXx8fGyuE6dOuHIkSN48uQJOI4Dx3GwtrYu+Yfy/3bs2IE2bdpAU1MTBgYG8PDwwIkTJ+Ri/vjjDzRr1gwikQgWFhbw8/PD27dv5WI6deqE5s2b4+bNm+jYsSM0NTVhY2ODvXv3AgDOnj2Ltm3bQkNDA7a2tjh58qTc5xcsWACO43Dv3j188cUX0NXVhZGREaZOnYrc3Fy52ODgYHTp0gWmpqYQiUSwt7fHunXriny3wv5x/PhxWf/YsGGDbN2H99jl5+dj4cKFaNy4MdTV1WFkZIQOHTogLCxMrs1Tp07B3d0dWlpa0NfXR79+/XD37t1iv0t8fDxGjx4NfX196OnpwdfXF9nZ2Z/8mRBSLTFCSK1Rt25d1rBhwzLH+/j4MABs8ODBLDAwkI0aNYoBYP3795eLs7KyYra2tqxOnTpswYIFbNWqVaxu3bpMW1ub7dixg9WvX58tW7aMLVu2jOnp6TEbGxsmkUjktqOurs4aN27MRo4cydauXct69+7NALCffvpJFieVSlmXLl0Yx3Fs3LhxbO3ataxPnz4MAJs2bZpcTgCYo6Mjq1OnDlu8eDELCAhgDRs2ZJqamuzVq1eyuOTkZFavXj1maWnJFi1axNatW8f69u3LALBVq1bJ4k6fPs0AMGdnZ+bi4sJWrVrFFixYwDQ1NVmbNm1kcSdOnGBOTk7M2NiYbd++nW3fvp3t37+/1P28YMECBoC1b9+erVixgq1evZp99dVXbNasWbKY+fPnMwDM09OTrVmzhk2aNInx+XzWunVrJhaLZXEdO3ZkFhYWzNLSkn3//fdszZo1zN7envH5fLZ7925mbm7OFixYwAICAljdunWZnp4ey8jIKLIdBwcH1qdPH7Z27Vo2YsQIBoCNHDlSLu/WrVuz0aNHs1WrVrE1a9aw7t27MwBs7dq1RfqHjY0NMzAwYLNnz2br169np0+flq3z8fGRxc6dO5dxHMfGjx/PNm3axH777Tc2bNgwtmzZMllMWFgYEwgErEmTJmz58uVs4cKFzNjYmBkYGLCEhIQi38XZ2ZkNHDiQ/fHHH2zcuHEMAJs5c2apPxNCqisq7AipJdLT0xkA1q9fvzLFx8TEMABs3LhxcstnzJjBALBTp07JlllZWTEA7MKFC7Jlx48fZwCYhoYGe/LkiWz5hg0bGADZL3bG/isgJ0+eLFsmlUpZr169mFAoZKmpqYwxxg4cOMAAsCVLlsjlNHjwYMZxHIuPj5ctA8CEQqHcshs3bjAAbM2aNbJlY8eOZXXq1JEr9hhj7Msvv2R6enosOzubMfZfYde0aVOWl5cni1u9ejUDwG7duiVb1qtXL2ZlZVXMXi0qLi6O8Xg8NmDAALlit3AfMMbYy5cvmVAoZN27d5eLWbt2LQPAgoKCZMs6duzIALCdO3fKlt27d48BYDwej128eFG2vPBnFBwcLFtWWAz17dtXLpdvv/2WAWA3btyQLSvcNx/y8vIq8sdDYf8IDQ0tEv9xYefo6Mh69epVJO5DTk5OzNTUlL1+/Vq27MaNG4zH47FRo0YV+S5jxoyR+/yAAQOYkZFRqdsgpLqiS7GE1BIZGRkAAB0dnTLFHz16FADg7+8vt/y7774DgCL34tnb28PV1VX2vm3btgCALl26oH79+kWWP3r0qMg2J02aJPt34aVUsVgsu1x49OhR8Pl8TJkypUhOjDEcO3ZMbrmnpycaNWoke9+iRQvo6urKts0Ywz///IM+ffqAMYZXr17JXl5eXkhPT0d0dLRcm76+vhAKhbL37u7uJX6fsjhw4ACkUinmzZsHHk/+kMxxHADg5MmTEIvFmDZtmlzM+PHjoaurW+Rnoa2tjS+//FL23tbWFvr6+mjatKls/wOl/yz8/Pzk3k+ePBnAf/0CADQ0NGT/Tk9Px6tXr9CxY0c8evQI6enpcp9v0KABvLy8StkT7+nr6+POnTuIi4srdn1SUhJiYmIwevRoGBoaypa3aNEC3bp1k8uv0DfffCP33t3dHa9fv5b9nyCkJqHCjpBaQldXFwDw7t27MsU/efIEPB4PNjY2csvNzc2hr6+PJ0+eyC3/sHgDAD09PQCApaVlscvfvHkjt5zH46Fhw4Zyy5o0aQIAstGlT548gYWFRZHitGnTprL1peUEAAYGBrJtp6am4u3bt9i4cSNMTEzkXr6+vgCAly9fltqmgYFBsd+nrB4+fAgejwd7e/sSYwq/l62trdxyoVCIhg0bFvne9erVkxWFhfT09Mr8swCAxo0by71v1KgReDye3Ejf8+fPw9PTU3afm4mJCebOnQsAxRZ2ZbFo0SK8ffsWTZo0gYODA77//nvcvHlTtr6kfQG87wevXr1CVlaW3HJl/8wIqcpoVCwhtYSuri4sLCxw+/btcn3u4wKhJHw+v1zL2WeYaelT25ZKpQCAESNGwMfHp9jYFi1alKvNqqAyfhYf94OHDx+ia9eusLOzw8qVK2FpaQmhUIijR49i1apVsn1b6MOze6Xx8PDAw4cPcfDgQZw4cQJ//vknVq1ahfXr12PcuHFlauNj1eFnRoiyUGFHSC3Su3dvbNy4EVFRUXKXTYtjZWUFqVSKuLg42RkxAEhJScHbt29hZWWl1NykUikePXokO0sHAA8ePAAA2ahSKysrnDx5Eu/evZM7a3fv3j3Z+vIwMTGBjo4OJBIJPD09FfwG/ylrMQy8PxMmlUoRGxsLJyenYmMKv9f9+/flzmqKxWIkJCQoNfdCcXFxcmfZ4uPjIZVKZT+Lf//9F3l5eTh06JDcGbHTp08rvG1DQ0P4+vrC19cXmZmZ8PDwwIIFCzBu3Di5ffGxe/fuwdjYGFpaWgrnQEh1RZdiCalFZs6cCS0tLYwbNw4pKSlF1j98+BCrV68GAPTs2RMAijxBYeXKlQCAXr16KT2/tWvXyv7NGMPatWuhpqaGrl27ynKSSCRycQCwatUqcBwHb2/vcm2Pz+dj0KBB+Oeff4o9k5mamlqBbwFoaWkVuRRZkv79+4PH42HRokVFznIVnlHy9PSEUCjE77//LneWafPmzUhPT6+Un0VgYKDc+zVr1gCAbB8
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 640x480 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"d = np.linspace(0.1, 200, 4000) # downrange distance\n",
|
||
|
|
"area_objective = np.pi * (80e-3 / 2) ** 2\n",
|
||
|
|
"area_spot = np.pi * (3e-3 / 2) ** 2 # approximation\n",
|
||
|
|
"gain_optical = area_objective / area_spot\n",
|
||
|
|
"\n",
|
||
|
|
"# Received power density\n",
|
||
|
|
"power_density = dict()\n",
|
||
|
|
"for led, properties in LEDS.items():\n",
|
||
|
|
" incident_power = properties[\"emissivity\"] * (1 / (d**2)) # W / m^2 at objective\n",
|
||
|
|
" power_density[led] = incident_power * gain_optical\n",
|
||
|
|
"\n",
|
||
|
|
"fig, ax = plt.subplots(tight_layout=True)\n",
|
||
|
|
"for led in LEDS:\n",
|
||
|
|
" lines = ax.loglog(\n",
|
||
|
|
" d,\n",
|
||
|
|
" power_density[led],\n",
|
||
|
|
" label=f\"{led} ({LEDS[led]['emissivity_current']/1e-3:5.1f} mA)\",\n",
|
||
|
|
" )\n",
|
||
|
|
"\n",
|
||
|
|
"ax.set_prop_cycle(None) # reset color cycle to keep colors consistent across plots\n",
|
||
|
|
"for diode, properties in PHOTODIODES.items():\n",
|
||
|
|
" ax.axhline(\n",
|
||
|
|
" np.max(properties[\"light_current\"].coords[\"irradiance\"]),\n",
|
||
|
|
" linestyle=\"--\",\n",
|
||
|
|
" label=f\"{diode} Psat\",\n",
|
||
|
|
" color=ax._get_lines.get_next_color(),\n",
|
||
|
|
" )\n",
|
||
|
|
"\n",
|
||
|
|
"ax.set_xlabel(\"Range [m]\")\n",
|
||
|
|
"ax.set_ylabel(\"RX Power Density [$W/m^2$]\")\n",
|
||
|
|
"ax.grid(True)\n",
|
||
|
|
"ax.legend(loc=\"upper right\")\n",
|
||
|
|
"ax.set_title(\"Component comparison\")"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"metadata": {
|
||
|
|
"kernelspec": {
|
||
|
|
"display_name": ".venv",
|
||
|
|
"language": "python",
|
||
|
|
"name": "python3"
|
||
|
|
},
|
||
|
|
"language_info": {
|
||
|
|
"codemirror_mode": {
|
||
|
|
"name": "ipython",
|
||
|
|
"version": 3
|
||
|
|
},
|
||
|
|
"file_extension": ".py",
|
||
|
|
"mimetype": "text/x-python",
|
||
|
|
"name": "python",
|
||
|
|
"nbconvert_exporter": "python",
|
||
|
|
"pygments_lexer": "ipython3",
|
||
|
|
"version": "3.10.12"
|
||
|
|
}
|
||
|
|
},
|
||
|
|
"nbformat": 4,
|
||
|
|
"nbformat_minor": 2
|
||
|
|
}
|