From b8208ff23b382502c1a05a241e4753a127296d3e Mon Sep 17 00:00:00 2001 From: Kelvin Ly Date: Tue, 3 May 2022 09:27:25 -0400 Subject: [PATCH] Working transient simulation TODO sweep voltage and measure harmonic distortion --- preamp sim.ipynb | 263 +++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 263 insertions(+) create mode 100644 preamp sim.ipynb diff --git a/preamp sim.ipynb b/preamp sim.ipynb new file mode 100644 index 0000000..45ad246 --- /dev/null +++ b/preamp sim.ipynb @@ -0,0 +1,263 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "id": "9943dcb3", + "metadata": {}, + "outputs": [], + "source": [ + "import matplotlib.pyplot as plt\n", + "\n", + "\n", + "import PySpice.Logging.Logging as Logging\n", + "logger = Logging.setup_logging()\n", + "\n", + "\n", + "from PySpice.Probe.Plot import plot\n", + "from PySpice.Spice.Library import SpiceLibrary\n", + "from PySpice.Spice.Netlist import Circuit\n", + "from PySpice.Unit import *\n", + "\n", + "#spice_library = SpiceLibrary(\"library/\")" + ] + }, + { + "cell_type": "code", + "execution_count": 21, + "id": "cf01c009", + "metadata": {}, + "outputs": [], + "source": [ + "circuit = Circuit(\"JFET preamp\")\n", + "\n", + "v_phantom = circuit.V('phantom_power', 'vdd', circuit.gnd, 48@u_V)\n", + "r_l = circuit.R('1', 'vdd', 1, 6.81@u_kOhm)\n", + "r_r = circuit.R('2', 'vdd', 2, 6.81@u_kOhm)\n", + "r_load = circuit.R('load', 1, 2, 3@u_kOhm)\n", + "\n", + "circuit.model('LSK489A', 'NJF', beta=2.2e-3, betatce=-.5, rd=11, rs=30, lambda_=4.3e-3, vto=-1.13, vtotc=-2.5e-3,\n", + " is_=3e-15,\n", + " isr=0, n=1, xti=0, alpha=30e-6, VK=120, Cgd=3.19e-12, Mj=0.32,\n", + " Pb=0.8, Fc=0.5, Cgs=2.92e-12, Kf=0.0009e-15,\n", + " Af=1, Gdsnoi=2.15, Nlev=3, Mfg='Linear_Systems')\n", + "\n", + "q1 = circuit.JFET('q1', 1, 'g1', 3, model='LSK489A')\n", + "q2 = circuit.JFET('q2', 2, 'g2', 4, model='LSK489A')\n", + "\n", + "rs1 = circuit.R('s1', 3, 5, 100@u_Ohm)\n", + "rs2 = circuit.R('s2', 4, 5, 100@u_Ohm)\n", + "\n", + "rs = circuit.R('s', 5, circuit.gnd, 5@u_kOhm)\n", + "\n", + "# gate biasing\n", + "rg1 = circuit.R('g1', 1, 'g1', 10@u_MOhm)\n", + "rg2 = circuit.R('g2', 2, 'g2', 10@u_MOhm)\n", + "rg3 = circuit.R('g3', circuit.gnd, 'g1', 10@u_MOhm)\n", + "rg4 = circuit.R('g4', circuit.gnd, 'g2', 10@u_MOhm)\n", + "\n", + "cg1 = circuit.C('1', 'g1', 'cin1', 100@u_uF)\n", + "cg2 = circuit.C('2', 'g2', 'cin2', 100@u_uF)\n", + "\n", + "rin1 = circuit.R('in1', 'cin1', 'in1', 0.5@u_MOhm)\n", + "rin2 = circuit.R('in2', 'cin2', 'in2', 0.5@u_MOhm)\n", + "vs = circuit.SinusoidalVoltageSource('input', 'in1', 'in2', amplitude=0.1@u_V, frequency=1e+3)\n", + "rt = circuit.R('tmp', 'in1', circuit.gnd, 5@u_MOhm)" + ] + }, + { + "cell_type": "code", + "execution_count": 22, + "id": "5b07a8e9", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 22, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "simulator = circuit.simulator(temperature=25, nominal_temperature=25)\n", + "\n", + "#analysis = simulator.transient(step_time=1@u_us, end_time=1000@u_us)\n", + "analysis = simulator.operating_point()\n", + "analysis" + ] + }, + { + "cell_type": "code", + "execution_count": 23, + "id": "cd57b8c2", + "metadata": {}, + "outputs": [], + "source": [ + "tran_analysis = simulator.transient(step_time=1@u_us, end_time=1000@u_us)\n" + ] + }, + { + "cell_type": "code", + "execution_count": 24, + "id": "c48c200d", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "WaveForm g1 [17.89245916]@V" + ] + }, + "execution_count": 24, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "analysis['g1']" + ] + }, + { + "cell_type": "code", + "execution_count": 25, + "id": "32f62881", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 25, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD4CAYAAADlwTGnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAuC0lEQVR4nO3deXhU5fn/8fedhWxAQiCBQAghJGELECAsooKsIqKotRUt1t22bkix1farVWttqXWpS4siUHFD61YVcUFUEGQLGvZASMIWQkiA7Hvy/P7Iob80DWQSkpxZ7td1zTUz58wz83kinvuc52xijEEppZTn8bI7gFJKKXtoAVBKKQ+lBUAppTyUFgCllPJQWgCUUspD+dgdoDm6detmoqOj7Y6hlFIuZevWrXnGmLCG012qAERHR5OcnGx3DKWUcikicrCx6ToEpJRSHkoLgFJKeSgtAEop5aG0ACillIfSAqCUUh5KC4BSSnkoLQBKKeWhXOo8ANW2KqprSMspJiu/jGMF5RRXVGOMwctL6BbkR1hnP2LDOhLZJQARsTuuUuocaQHwYLW1hpQj+XyxK4cN6Xnszi6kqqbp+0MEB/gyPCqESQPCmTQgnMguge2QVinV2sSVbgiTlJRk9Ezgc3e8sJzlmw/z1pZDZBeU4+MljOjTheFRIQyLDCEqNJDunf3pHOCDlwjVNYa84gqOF5WTeqyInVkFbMw4SWZeCQCj+4Zy3egopif0wN/X2+beKaUaEpGtxpik/5muBcBzHD5ZyvNfpfH+91lU1xrGx4dx5fCeTOrfneBA32Z/X3puMZ/tPMbbWw5z6GQpYZ38+OWEflw3JkoLgVJORAuABysoq+KZVft4Y9NBRITrRkdx47hoorsFtcr319Ya1u3PY+E36WzIOEH3zn48cMkArkjspfsKlHICWgA8kDGGj7dn89iK3ZworuCaUVHcMzmWiOCANvvNDeknWPBZKtsO5zM6OpTHr0wgrnunNvs9pVTTtAB4mJMllfzm3e18uSeHIb2C+dOVQxgSGdwuv11ba3hn62EWfJpKSWUNv71kADecF42Xl24NKGWHMxUAPQrIDW1IP8G9b//AqZIqHrx0IDed3xfvdlz4enkJ14yKYtKA7tz/3nYe/Xg3X6Ue5+mfJBLWya/dciilzk5PBHMjxhheXJPOdYs3EtjBh/fvGMetF8a068K/vrBOfiy5IYk/XpHAlgMnufyFdew4UmBLFqXU/9IC4CYqqmv49bvbWfBpKjMSIlhx9wUk9GqfIZ+zERHmjO3Du78YhwBXv/gdH6Zk2R1LKYUWALdQUFbF9Ys38+7WI8ydHMcL1w0nyM+5RvcSegXz0d0XMCwyhLlvpfDimnS7Iynl8bQAuLgTxRVcu2gjPxw+xbOzE5k3Nd5pD73s1tGP128dw2XDerLg01T+tHIPrnQQglLuxrlWE1WzHCso56eLN5KVX8biG0YxIf5/7vnsdDr4ePHsNYl0CfRl0doM8ksrWXDVUD1CSCkbaAFwUdkFZVzz0kZOllSy7KbRjInpanckh3l5CY9ePpiQAF+e+2o/3l7C41cM0SKgVDtrcghIRPxFZLOIbBORXSLyqDX9ERHJEpEU6zHjDO3nishOq+299aaHisgqEUmznru0Wq/cXF5xBT9dvIlTJZW8cesYl1r4nyYizJsaz10TY1m++TCPfLxLh4OUameO7AOoACYZY4YBicB0ERlrzXvGGJNoPVY2bCgiCcBtwGhgGDBTROKs2Q8Aq40xccBq671qQkFZFT9bspmj+WUsvWkUw3qH2B2pxUSE+dPi+fn4GF7dcJA/rdxjdySlPEqTBcDUKbbe+loPR1fVBgIbjTGlxphqYA1wpTVvFrDMer0MuMLR0J6qrLKGm1/ZQtrxIl66PolR0aF2RzpnIsIDlwzgxnHRvPxtJou/zbA7klIew6GjgETEW0RSgOPAKmPMJmvWXSKyXUSWnmEIZycwXkS6ikggMAPobc3rbozJBrCew8/w27eLSLKIJOfm5jreMzdTW2v41b9S+P7QKZ6dPdwldvg6SkT4/cxBzBjSg8dX7uGT7dl2R1LKIzhUAIwxNcaYRCASGG0N7SwE+lE3LJQNPNVIuz3AX4BVwGfANqC6OQGNMYuMMUnGmKSwMPdZ6DXXXz5P5dOdx/i/GQOZMSTC7jitzstLePoniYyM6sK8f6WwOfOk3ZGUcnvNOg/AGJMPfANMN8bkWIWhFniZunH+xtosMcaMMMaMB04CadasHBGJALCej7esC+5v+eZDvLQmgzljo7jlgr52x2kz/r7evPyzJCK7BPDz15I5fLLU7khKuTVHjgIKE5EQ63UAMAVIPb3wtlxJ3XBPY+3Dreco4CpguTXrI+AG6/UNwIctyO/2NmWc4MF/72RCfBiPXDbYaU/yai1dgjqw5IZRVNcabn9tK2WVNXZHUsptObIFEAF8LSLbgS3U7QNYATwhIjus6ROBeQAi0lNE6h8R9J6I7AY+Bu40xpyypi8ApopIGjDVeq/qySks5843f6BPaCDPXzccH2/POHG7b7cgnrt2OKnHCvnNe9v18FCl2kiTJ4IZY7YDwxuZfv0ZPn+Uup29p99feIbPnQAmO5zUw1RW13LHG99TWlnNm7eNobN/82/Z6Mom9g/n1xf354nP9jKkV2duH9/P7khKuR3PWKV0QY9/sputB0/xxNVDiffQO2r9ckI/Lh0SwYJPU9mQfsLuOEq5HS0ATujDlCyWbTjIrRf0ZebQnnbHsY2I8MTVQ4nuFsTct37gRHGF3ZGUcitaAJzMoROl/N8HOxnZpwv3XzLA7ji2C/Lz4YVrR5BfVsX8d7ZRW6v7A5RqLVoAnEh1TS1z3/4BEXh2diK+HrLTtymDenbmoUsH8s3eXBav0zOFlWotuoRxIs+tTuOHQ/n86cohRHYJtDuOU5kztg+XJPTgic/28sOhU003UEo1SQuAk9iUcYIXvt7Pj0dGctkwzx33PxMRYcGPhtK9sz/z/7VNzw9QqhVoAXACBWVVzHs7hT5dg3jk8sF2x3FawQG+/PXHQ8nIK+Evn6XaHUcpl6cFwAn8ccVucooq+Ns1iU53L19nM65fN24cF80r3x1g/f48u+Mo5dK0ANjs69TjvLP1CL+YEOPS1/ZvT/dPH0BMWBC/fmcbheVVdsdRymVpAbBRQVkVD7y/nfjuHblnclzTDRQAAR28efonieQUVfCHj3fbHUcpl6UFwEaPrdhNXnElT/54GH4+3nbHcSmJvUO446J+vLv1CF/v1QvJKtUSWgBs8nXqcd61hn6GRobYHccl3TUpltjwjjz4wU5KKpp1mwmlFFoAbFFcUc3vPtihQz/nyM/HmwVXDSErv4ynvthndxylXI4WABs8s2of2QXl/PmqoTr0c46SokOZMzaKV77LZNvhfLvjKOVStAC0s51ZBfxzfSbXjYliZJ/GbqOsmus30wcQ1smP+9/bTlVNrd1xlHIZWgDaUU2t4Xcf7CA0yI/7L9YLvbWWzv6+/GFWAqnHinj5W71WkFKO0gLQjl7feJDtRwp4aOZAggM96wYvbe3iwT2YPrgHz36ZxqETei9hpRyhBaCd5BSW89fP93JhXDcu12v9tImHLx+Et5fw2Cd6boBSjtAC0E4eW7Gbqppa/nhFgtvf2N0uEcEB3D0pjlW7c/TcAKUcoAWgHWzMOMGK7dn88qJ+9OkaZHcct3bLBX2J6RbEox/toqJarxiq1NloAWhj1TW1PPLRLnqFBPCLCXpj87bWwceLRy4fzIETpSz+NtPuOEo5NS0AbWz55kOkHiviwUsH4u+rx/y3h/HxYUwf3IPnv0ojK7/M7jhKOS0tAG3oVEklT63ax3kxXZme0MPuOB7locsGAXWX2lZKNU4LQBt6etU+isqrefjyQbrjt531CgngromxfLrzmN43QKkz0ALQRnYfLeSNTQe5fmwfBvTobHccj3TrhTFEdgngsRW7qak1dsdRyuloAWgDxhge+XgXwQG+zJsSb3ccj+Xv680Dlwwg9VgR7249bHccpZyOFoA28NnOY2zOPMn8af31jF+bXTokghFRITz5xT6K9ZLRSv2XJguAiPiLyGYR2SYiu0TkUWv6IyKSJSIp1mPGGdrPs9rtFJHlIuLfnPauprK6lr98lkpceEdmj+ptdxyPJyI8OHMQuUUVvLQm3e44SjkVR7YAKoBJxphhQCIwXUTGWvOeMcYkWo+VDRuKSC/gHiDJGJMAeAOz633krO1d0ZubDnLgRCm/mzEQH2/dwHIGI6K6cPmwnixam8FRPSxUqf9ocgll6hRbb32tR3P2qPkAASLiAwQCR5ud0kUUllfx7Oo0xvXrykX9w+yOo+r5zfT+GOCvn++1O4pSTsOhVVQR8RaRFOA4sMoYs8madZeIbBeRpSLyPxe3N8ZkAU8Ch4BsoMAY80W9j5y1vfXbt4tIsogk5+bmNqNr7W/hN+mcKq3idzMG6mGfTiaySyC3XtCXD37I0hvHKGVxqAAYY2qMMYlAJDBaRBKAhUA/6oaFsoGnGrazFuqzgL5ATyBIROZYs5tsb/32ImNMkjEmKSzMedeqs/LLWLIukyuH9yKhV7DdcVQjfnlRP7p17MDjn+zBGD0sVKlmDVIbY/KBb4DpxpgcqzDUAi8DoxtpMgXINMbkGmOqgPeBcdZ3OdLeZTxlDS3cd3F/m5OoM+nk78vcKfFsPnBSrxaqFI4dBRQmIiHW6wDqFuqpIhJR72NXAjsbaX4IGCsigVI3JjIZ2GN9lyPtXcLOrAI+SMni5vP70iskwO446ixmj+pNn66BPPHZXj05THk8R7YAIoCvRWQ7sIW6fQArgCdEZIc1fSIwD0BEeorISgBrX8G7wPfADuv3Flnf22h7V2OM4c+f7iEkwJc7JurVPp2dr7cX86f1J/VYER+mZNkdRylbiSuNhSYlJZnk5GS7Y/yXdWl5zFmyid/PHMTNF/S1O45yQG2t4bIX1lFQVsXq+RPw89GrtCr3JiJbjTFJDafrgernwBjDXz9PpVdIAD8dG2V3HOUgLy/h/ukDOHKqjDc3HbI7jlK20QJwDj7flcO2IwXMnRKna5Eu5sK4bozr15UXvtqvl4hQHksLQAvV1Bqe+mIv/cKCuGp4L7vjqGYSEX4zfQAnSipZ/G2G3XGUsoUWgBb6MCWLtOPFzJ/WXy/54KISe4dwSUIPXl6bQV5xhd1xlGp3uuRqgcrqWp75ch8JvTozfbDe6cuV3Xdxf8qra3nhq/12R1Gq3WkBaIG3kw9z+GQZ903rj5eXXvLBlfUL68hPkiJ5Y9NBvX+w8jhaAJqprLKG51enMTo6lAnxzntpCuW4uyfFIQgvfJVmdxSl2pUWgGZatuEAx4sq+PX0/nrBNzfRMySAa0f35p3kIxw6UWp3HKXajRaAZigsr2LhN+lM7B/GqOhQu+OoVnTHxFi8vYTndCtAeRAtAM2weG0GBWVVzJ+mF3xzN907+zNnbB/e//4IGbnFTTdQyg1oAXBQfmklS9cf4NIhEXq5Zzf1y4v64efjzbOrdStAeQYtAA5asi6T4opq7pkcZ3cU1Ua6dfTjhnHRfLTtKPtyiuyOo1Sb0wLggILSKl5Zf4AZQ3rQv0cnu+OoNvTz8TEEdfDhb1/uszuKUm1OC4ADlqzPpEjX/j1Cl6AO3Hx+NCt3HGP30UK74yjVprQANKGgrIp/rs9k+uAeDOjR2e44qh3cckEMnfx9eEa3ApSb0wLQhKXrMikq17V/TxIc6MttF8awancO24/k2x1HqTajBeAsCsqqWLo+k2mDujOop679e5Kbzo8mJNCXp1fpVoByX1oAzuKV9Qd07d9DdfL35efj+/HN3lxSDufbHUepNqEF4AwKy6tYsi6DqYO663H/Hupn5/WhS6Avz+l5AcpNaQE4g2XrD1BYXs1cXfv3WEF+Ptx6YQxfpR7XfQHKLWkBaERReRWL12UyZWC4rv17uJ+d14fgAF+eW633C1DuRwtAI5Z9d4CCsirmTo63O4qyWSd/X265oC9f7slhZ1aB3XGUalVaABoorqhm8bpMJg0IZ0ikrv0ruGFcNJ38fXherxSq3IwWgAaWfXeA/NIqHftX/xEc4MvN5/fl81057MnWs4OV+9ACUE9JRTWLv83gov5hDOsdYncc5URuPr8vnfx0K0C5Fy0A9by64SCndO1fNSI40JcbrWsE7T2mVwpV7kELgKWkopqXv81gQnwYw6O62B1HOaGbz+9LUAdv3QpQbqPJAiAi/iKyWUS2icguEXnUmv6IiGSJSIr1mHGG9vOsdjtFZLmI+FvTQ0VklYikWc+2LnVf23iQkyWVzJ2ia/+qcV2COvCzcdF8siOb/cd1K0C5Pke2ACqAScaYYUAiMF1ExlrznjHGJFqPlQ0bikgv4B4gyRiTAHgDs63ZDwCrjTFxwGrrvS1KK6t5eW0GF8Z1Y4Su/auzuO3CGAJ8vXn+Kz0vQLm+JguAqXP6Jqm+1sM04zd8gAAR8QECgaPW9FnAMuv1MuCKZnxnq3p940FOlFRyr679qyaEBnXg+rF9+HjbUdL13sHKxTm0D0BEvEUkBTgOrDLGbLJm3SUi20VkaWNDOMaYLOBJ4BCQDRQYY76wZnc3xmRbn8sGws/w27eLSLKIJOfm5janbw4pq6xh0doMLojtxsg+oa3+/cr93DY+hg4+XvxdtwKUi3OoABhjaowxiUAkMFpEEoCFQD/qhoWygacatrOKwiygL9ATCBKROc0JaIxZZIxJMsYkhYWFNaepQ97YdJC8Yh37V47r1tGPOWP68O+ULA7kldgdR6kWa9ZRQMaYfOAbYLoxJscqDLXAy8DoRppMATKNMbnGmCrgfWCcNS9HRCIArOfjLetCy5VV1vDimgzG9evKqGhd+1eOu31CDL7eXrzwtW4FKNflyFFAYSISYr0OoG6hnnp64W25EtjZSPNDwFgRCRQRASYDe6x5HwE3WK9vAD5sUQ/OQd3af4Ue96+aLbyTP9eNieKDH7I4dKLU7jhKtYgjWwARwNcish3YQt0+gBXAEyKyw5o+EZgHICI9RWQlgLWv4F3ge2CH9XuLrO9dAEwVkTRgqvW+3ZRX1fDS2gzGxoQyJqZre/60chO/mNAPby9h4RrdClCuyaepDxhjtgPDG5l+/Rk+fxSYUe/9w8DDjXzuBHVbBLZ4c9MhcosqeG72/3RNKYd07+zPNUm9eWvLIe6cGEtkl0C7IynVLB55JnB5VQ0vrklnTN9Qzuuna/+q5X5xUT8AXlyTbnMSpZrPIwvA6xsPcryognun6PX+1bnpFRLA1SN7868tR8guKLM7jlLN4nEFoKSimn98k84Fsd107V+1ijsu6ketMby0JsPuKEo1i8cVgJfWpHOypJL7Lu5vdxTlJnqHBnLViF4s33yI44XldsdRymFN7gR2B8+vTuPfKVlUVNdy5FQZVyT2JFGv969a0Z0TY3nv+ywWrc3gwZmD7I6jlEM8YgugWyc/BvTozOjoUH4xoR9/vHKI3ZGUm+nTNYhZw3ryunVuiVKt5VRJJdcv2dQm96T2iC2Aa0dHce3oKLtjKDd356RYPkjJ4uVvM/jtJQPtjqPcxJJ1mazbn8dDPq2/vu4RWwBKtYd+YR25bGhPXttQd28Jpc5Vfmklr3x3gBkJEcR379Tq368FQKlWdNekWMqqali6LtPuKMoNLF2XSXFFNXdPjm2T79cCoFQriu/eiRkJEbzy3QEKSqvsjqNcWEFpFf9cf4Dpg3swoEfnNvkNLQBKtbK7JsVSXFHN0vW6FaBabun6TIoqqrmnDS9WqQVAqVY2MKIz0wZ1Z+n6TArLdStANV9heRVL12cybVB3BvVsm7V/0AKgVJu4Z3IcReXVLFt/wO4oygW9sv4AReVtu/YPWgCUahMJvYKZPCCcJevrduIp5aii8iqWrMtkysDuJPQKbtPf0gKgVBu5e3Ic+aVVvLbhoN1RlAtZ9t0BCsqq2uVGVVoAlGojib1DmBAfxsvfZlBaqVsBqmnFFdUsXpfJ5AHhDIls27V/0AKgVJu6Z3IsJ0sqeWPjIbujKBew7LsD5JdWtfnY/2laAJRqQyP7hHJ+bFdeWptBeVWN3XGUEyupqGbxtxlc1D+MYe10sUotAEq1sXsmxZFXXMHyzboVoM7s1Q0HOVXaPmP/p2kBUKqNjYnpypi+oby4Jl23AlSjSiqqefnbDCbEhzE8qku7/a4WAKXawdzJceQUVvBO8mG7oygn9PrGugsIzp3Sfmv/oAVAqXZxXr+ujOzThX98k05FtW4FqP+vtLKaRWszuDCuGyPace0ftAAo1S5EhHsmx5FdUM57W7PsjqOcyBsbD3GipJJ723ntH7QAKNVuxsd1Y1jvEP7xzX6qamrtjqOcQFllDS+tTeeC2G6M7BPa7r+vBUCpdiIizJ0cy5FTZXzwg24FKHhj00Hyitt/7P80LQBKtaOJ/cNJ6NWZv3+9n2rdCvBoZZU1vLgmg3H9ujIquv3X/kELgFLtSkS4Z1IcB0+U8tG2o3bHUTZ6c/Mh8oor2vW4/4aaLAAi4i8im0Vkm4jsEpFHremPiEiWiKRYjxmNtO1fb36KiBSKyL2OtlfKHU0d1J0BPTrxwlf7qak1dsdRNiivquHFNemMjQllTExX23I4sgVQAUwyxgwDEoHpIjLWmveMMSbReqxs2NAYs/f0fGAkUAp8UO8jZ22vlDs6fURQRl4Jn+zItjuOssHyzYfILapg7uR4W3M0WQBMnWLrra/1aMlqy2Qg3Rij18ZVHm/64B7EhXfk+dVp1OpWgEc5vfY/um8o5/Wzb+0fHNwHICLeIpICHAdWGWM2WbPuEpHtIrJURJo6g2E2sLzBtCbbi8jtIpIsIsm5ubmOxFXK6Xl5CXdPjiPteDGf7TpmdxzVjt7ecpicwgrutXHs/zSHCoAxpsYaxokERotIArAQ6EfdsFA28NSZ2otIB+By4J16kx1qb4xZZIxJMsYkhYWFORJXKZdw6ZAIYsKCeE63AjxGeVUNC79JZ1R0F9vX/qGZRwEZY/KBb4DpxpgcqzDUAi8Do8/S9BLge2NMTr3vak57pdyOt5dw96RYUo8V6VaAh3hj0yGOFZYzb0o8ImJ3HIeOAgoTkRDrdQAwBUgVkYh6H7sS2HmWr7mWBsM/zWyvlFu6fFgvYsM78vSqfXpEkJsrraxm4Tf7GdevK+Niu9kdB3BsCyAC+FpEtgNbqNsHsAJ4QkR2WNMnAvMARKSniPzniB4RCQSmAu83+N5G2yvlSby9hF9NjWf/8WI+2qZnB7uzZd/VnfU7f5q9R/7U59PUB4wx24HhjUy//gyfPwrMqPe+FPifwa4ztVfK00wf3INBEZ3525dpzBzaE19vPT/T3RSWV/HimnQu6h9myzV/zkT/pSllMy8vYf60eA6eKOXdrUfsjqPawNJ1mRSUVTF/an+7o/wXLQBKOYFJA8JJ7B3Cc6vT9K5hbia/tJIl32Zy8eDuDIkMtjvOf9ECoJQTEBHum9af7IJy3tJ7B7uVRWszKK6sZt5U5xn7P00LgFJO4vzYroyNCeWFr9Mpq9StAHeQV1zBP9cf4LKhPRnQo7Pdcf6HFgClnISIMH9af/KKK3h1wwG746hWsNC6Bagdd/tyhBYApZzIqOhQJsSHsXBNOkXlVXbHUefgWEE5r208yFUjIokJ62h3nEZpAVDKycyfFk9+aRVL1x2wO4o6B89/VXeJDzuv998ULQBKOZmhkSFMG9Sdxd9mcKqk0u44qgUy80p4a8thZo/uTe/QQLvjnJEWAKWc0H0X96ekspq/f73f7iiqBZ78Yi9+Pl7c48Rr/6AFQCmnFN+9Ez8aEcmrGw5y+GSp3XFUM2w7nM8n27O59cIYwjv52x3nrLQAKOWkfjUtHhF4ZtU+u6MoBxljWPBpKqFBHbjtwr52x2mSFgClnFREcAA3nd+XD1Ky2H200O44ygFr0/LYkHGCuyfF0snf1+44TdICoJQT++WEfnT292XBZ6l2R1FNqK01/OXTVHqHBnDdmCi74zhEC4BSTiw40Je7Jsaydl8u6/fn2R1HncXH24+yO7uQ+6b1x8/H2+44DtECoJSTu/68PvQKCWDBp6l660gnVVFdw18/38ugiM5cNrSn3XEcpgVAKSfn7+vN/Gnx7MgqYMWObLvjqEa8sfEQR06Vcf8lA/Dysv9Wj47SAqCUC5iV2IsBPTrx5Od7qayutTuOqie/tJJnV6dxQWw3xsc5x60eHaUFQCkX4O0l/HbGQA6dLNULxTmZv32ZRlF5FQ/OHOgUN3pvDi0ASrmICfFhTIgP49nVaZworrA7jgLSc4t5feNBZo+OcsrLPTdFC4BSLuShmQMprazhaT05zCn86ZM9BPh68ysnvNmLI7QAKOVCYsM7cf3YPizffIg92XpymJ3WpeWxOvU4d06KpVtHP7vjtIgWAKVczL1T4ugc4MtjK3ZjjB4WaoeaWsMfP9lN79AAbjo/2u44LaYFQCkXExLYgXlT4vku/QSrdufYHccjvb3lMKnHivjdJQNd5qSvxmgBUMoF/XRMFHHhHXl85R4qqvX+we2poKyKp1ftZXTfUKYn9LA7zjnRAqCUC/Lx9uKhmYM4eKKUV9YfsDuOR3lm1T5OllTy+5mDXO6wz4a0ACjlosbHhzF5QDjPrU7jWEG53XE8wu6jhby64QBzxvYhoVew3XHOWZMFQET8RWSziGwTkV0i8qg1/RERyRKRFOsxo5G2/evNTxGRQhG515oXKiKrRCTNeu7S6r1Tys09fNlgqmsNj63YbXcUt1dba/j9hzvpEtiB+VP72x2nVTiyBVABTDLGDAMSgekiMtaa94wxJtF6rGzY0Biz9/R8YCRQCnxgzX4AWG2MiQNWW++VUs0Q1TWQuybG8smObNbsy7U7jlt7/4cskg+e4v5LBhAc6PzX+ndEkwXA1Cm23vpaj5YcezYZSDfGHLTezwKWWa+XAVe04DuV8ni3T4ihb7cgHv5wJ+VVukO4LRSUVfHnlXsYERXC1SMi7Y7TahzaByAi3iKSAhwHVhljNlmz7hKR7SKy1IEhnNnA8nrvuxtjsgGs5/DmRVdKAfj5ePOHWYM5cKKUF9ek2x3HLT2zah+nSiv5w6wEl7raZ1McKgDGmBprGCcSGC0iCcBCoB91w0LZwFNnai8iHYDLgXeaG1BEbheRZBFJzs3VTVylGnNhXBgzh0bwj2/SOZBXYncct5JyON+tdvzW16yjgIwx+cA3wHRjTI5VGGqBl4HRZ2l6CfC9Mab+WSs5IhIBYD0fP8NvLjLGJBljksLCwpoTVymP8tDMQXTw9uL3H+3SM4RbSVVNLQ+8t53wTv7cd7F77Pitz5GjgMJEJMR6HQBMAVJPL7wtVwI7z/I11/Lfwz8AHwE3WK9vAD50MLNSqhHdO/szf1o8a/fl8u+ULLvjuIVFazNIPVbEY1ck0NkFbvLeXI5sAUQAX4vIdmALdfsAVgBPiMgOa/pEYB6AiPQUkf8cESQigcBU4P0G37sAmCoiadb8BefcG6U83M/Oi2ZEVAiPfryb3CK9ZPS5SM8t5tnVaVw6JIKpg7rbHadNiCttKiYlJZnk5GS7Yyjl1PYfL2LGs+uYMiicf/x0pN1xXFJtrWH2yxtJzS7ky/kTCO/kb3ekcyIiW40xSQ2n65nASrmZ2PBOzJ0Sx8odx/hU7yHcIm9tOczmzJM8eOkgl1/4n40WAKXc0O3jYxjcszMPfbiL/NJKu+O4lCOnSvnTyj2M69eVHye5zzH/jdECoJQb8vX24omrh5JfWskfPtbLRDiqttZw3zvbMMbwlx8NdfmLvTVFC4BSbmpwz2DumBjL+z9k6VCQg5auz2RjxkkevmwwvUMD7Y7T5rQAKOXG7p4Uy7DIYH77wQ69YmgT9uUU8cTne5kysLvbD/2cpgVAKTfm6+3FM9ckUlFVy33vbKO21nWO+mtPldW1zHs7hU5+Piz40RC3H/o5TQuAUm4uJqwjv79sEOv257F0fabdcZzSM1/uY9fRQv501RCXvcF7S2gBUMoDzB7Vm6mDuvPEZ3vZk11odxynsmZfLgu/SWf2qN5cPNi1b/HYXFoAlPIAIsKCq4YQHOjLnW9+T3FFtd2RnEJOYTm/ejuF/t078fBlg+2O0+60ACjlIbp29OP5a4dzIK+E+9/b7vEXjKupNdyz/AdKK2v4+0+HE9DB2+5I7U4LgFIeZGxMV3598QA+2Z7NK98dsDuOrZ79ch+bMk/yxysSiA3vZHccW2gBUMrD/Hx8DFMGhvP4J3vYevCU3XFs8cWuYzz31X6uHhnJj0Z6xiGfjdECoJSH8fISnvpxIhEh/tzxxlZyCj3r/IC0nCLmvZ3C0Mhg/nhFgt1xbKUFQCkPFBzoy0tzkigqr+a2V5Mpq/SMewkXlFZx26vJBHTw4aXrR+Lv63nj/vVpAVDKQw3q2ZlnZw9nR1aBR5wkVlVTy13Lvycrv4wX54wgIjjA7ki20wKglAebOqg7D0wfwCc7svnb6jS747QZYwy/e38H36bl8ccrEkiKDrU7klPwsTuAUspet4+PYf/xYp5bnUZkSAA/GdXb7kit7pkv03hn6xHumRzHNaOi7I7jNLQAKOXhRITHrxxCTlEFD7y/nc4BPkxPiGi6oYtYvvkQz61O4ydJkcybEmd3HKeiQ0BKKTr4ePHinBEk9g7hnuUprEvLsztSq/gwJYvffbCDCfFhPH6l51zkzVFaAJRSAAR28OGfN44mJiyI219LZnPmSbsjnZMV248y7+0UxvQN5cU5I/H11sVdQ/oXUUr9R3CgL6/eMpqIYH9uWLrZZbcEPtuZzdy3UhjZpwtLbhjlkZd5cIQWAKXUfwnv5M9bt59Hn66B3LxsC1+l5tgdqVne3XqEO9/8gWGRwfzzptEE+emuzjPRAqCU+h9hnfxYfttY4rt35OevbeWDH47YHckhL61J5753tnFeTFdevWUMHXXhf1ZaAJRSjeoS1IE3bh3LiKguzHt7G0+v2ue0VxCtrqnlkY928edPU5k5NIIlNybpwt8BWgCUUmcUHODLa7eM4UcjInludRr3vJVCaaVz3UvgZEklP1u6mVe+O8DN5/fl2dnD8fPRMX9HaIlUSp1VBx8vnvzxUGLCgnjyi7o7iv39uhH072H/JZRTDudz15vfc7yogid/PIyrPfjKni2hWwBKqSaJCHdOjOX1W8aQX1rFrL+v4/WNB227flBVTS1/+3IfP1r4HcbAOz8/Txf+LdBkARARfxHZLCLbRGSXiDxqTX9ERLJEJMV6zDhD+xAReVdEUkVkj4ic15z2SinncX5sN1bOvYCkPqE8+O+dzH55I+m5xe2aYWdWAVcv/I6/fZnGrGE9+fTeCxnWO6RdM7gLaWqnjtSdOhdkjCkWEV9gHTAXmA4UG2OebKL9MuBbY8xiEekABBpj8kXkEUfa15eUlGSSk5Md/bhSqo0YY/hX8mEe/2QP5VW13HRBNHdMiCU40LfNfjOvuIKnV+1j+eZDdA3qwCOXD2bm0J5t9nvuRES2GmOSGk5vch+AqasQp0u8r/VwaLtPRDoD44Ebre+qBCodi6yUclYiwjWjopg4IJwFK1NZtDaD5ZsOceuFMfx0TBRdO/q12m8dLyxn0doMXt90kKoaw03j+nLv1Dg6+7ddsfEUTW4BAIiIN7AViAX+boy531qDvxEoBJKB+caYUw3aJQKLgN3AMOs75hpjShxpb33H7cDtAFFRUSMPHjzYgm4qpdrSnuxCnvx8L6tTj9PBx4tZw3py5YhejI4OxacFl2Coqqll/f483t5ymFW7c6g1hiuG9+LOibH0C+vYBj1wb2faAnCoANT7khDgA+BuIBfIo25r4DEgwhhzc4PPJwEbgfONMZtE5Fmg0BjzkIh0b6p9QzoEpJRzS8sp4pXvDvD+91mUVdXQNagD4+PDGNmnC4m9Q4jqGtjomntReRV7jxWxO7uQjRkn+HZfHkUV1XQJ9OVHIyKZM7YP0d2CbOiRe2iVAmB90cNASf2xexGJBlYYYxIafLYHsNEYE229vxB4wBhzaYPPNdq+IS0ASrmG0spq1uzNZeXOY2xIP0FeccV/5nXy96GTnw++Pl5U1xhOlFRQXlX7n/ndO/sxsX84EweEc1H/MD2mvxW0eB+AiIQBVdaO2wBgCvAXEYkwxmRbH7sS2NmwrTHmmIgcFpH+xpi9wGTqhoNwpL1SyjUFdvDhkiERXDIkAmMMR06VsSOrgCOnSjlyqozSyhqqa2rxEqFrxw507ehHXHhHBkZ0JiLYXy/b3E4cOREsAlhm7QfwAv5ljFkhIq9ZY/wGOAD8HEBEegKLjTGnD+u8G3jDOgIoA7jJmv5EY+2VUu5FROgdGkjv0EC7o6gGmj0EZCcdAlJKqeY70xCQngmslFIeSguAUkp5KC0ASinlobQAKKWUh9ICoJRSHkoLgFJKeSgtAEop5aFc6jwAEckFWno1uG7UXXvIk3hanz2tv+B5fdb+tkwfY0xYw4kuVQDOhYgkN3YihDvztD57Wn/B8/qs/W1dOgSklFIeSguAUkp5KE8qAIvsDmADT+uzp/UXPK/P2t9W5DH7AJRSSv03T9oCUEopVY8WAKWU8lAeUQBEZLqI7BWR/SLygN15WoOI9BaRr0Vkj4jsEpG51vRQEVklImnWc5d6bX5r/Q32isjF9qVvORHxFpEfRGSF9d7d+xsiIu+KSKr13/o8d+6ziMyz/j3vFJHlIuLvbv0VkaUiclxEdtab1uw+ishIEdlhzXtOWnIbNWOMWz8AbyAdiAE6ANuAQXbnaoV+RQAjrNedgH3AIOAJ6u67DPAA8Bfr9SCr735AX+tv4m13P1rQ718Bb1J3D2k8oL/LgFut1x2AEHftM9ALyAQCrPf/Am50t/4C44ERwM5605rdR2AzcB4gwKfAJc3N4glbAKOB/caYDGNMJfAWMMvmTOfMGJNtjPneel0E7KHuf6BZ1C00sJ6vsF7PAt4yxlQYYzKB/dT9bVyGiEQClwKL60125/52pm5hsQTAGFNpjMnHjftM3W1qA0TEBwgEjuJm/TXGrAVONpjcrD6KSATQ2RizwdRVg1frtXGYJxSAXsDheu+PWNPchohEA8OBTUB3Y0w21BUJINz6mDv8Hf4G/AaorTfNnfsbA+QC/7SGvRaLSBBu2mdjTBbwJHAIyAYKjDFf4Kb9baC5fexlvW44vVk8oQA0Ni7mNse+ikhH4D3gXmNM4dk+2sg0l/k7iMhM4LgxZqujTRqZ5jL9tfhQN1Sw0BgzHCihbnjgTFy6z9a49yzqhjp6AkEiMudsTRqZ5jL9ddCZ+tgqffeEAnAE6F3vfSR1m5UuT0R8qVv4v2GMed+anGNtHmI9H7emu/rf4XzgchE5QN0w3iQReR337S/U9eGIMWaT9f5d6gqCu/Z5CpBpjMk1xlQB7wPjcN/+1tfcPh6xXjec3iyeUAC2AHEi0ldEOgCzgY9sznTOrD3+S4A9xpin6836CLjBen0D8GG96bNFxE9E+gJx1O1EcgnGmN8aYyKNMdHU/Tf8yhgzBzftL4Ax5hhwWET6W5MmA7tx3z4fAsaKSKD173sydfu23LW/9TWrj9YwUZGIjLX+Vj+r18Zxdu8Rb6e97jOoO0omHfg/u/O0Up8uoG6TbzuQYj1mAF2B1UCa9Rxar83/WX+DvbTgiAFneQAX8f+PAnLr/gKJQLL13/nfQBd37jPwKJAK7AReo+7oF7fqL7Ccun0cVdStyd/Skj4CSdbfKR14AevKDs156KUglFLKQ3nCEJBSSqlGaAFQSikPpQVAKaU8lBYApZTyUFoAlFLKQ2kBUEopD6UFQCmlPNT/A+f9skKKj/MaAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(tran_analysis['2'])" + ] + }, + { + "cell_type": "code", + "execution_count": 28, + "id": "ae8719b6", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 28, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD4CAYAAAAD6PrjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAtEklEQVR4nO3dd3hUZfr/8fdNeg+B0BIgBBJ6EUIvIgJSdkVFdlVUrIiuvaxY9ru6rquuvWDfgq6KrBWlC6isohA6AVIIJSEhhF7TJs/vjwz8IgZSppwp9+u65po5Z54z535ImE+eU8UYg1JKKQXQyOoClFJKeQ4NBaWUUqdpKCillDpNQ0EppdRpGgpKKaVOC7S6gIZo2rSpSUpKsroMpZTyKqtXr95njIk/VxuvDIWkpCTS09OtLkMppbyKiOysrY1uPlJKKXWahoJSSqnTNBSUUkqdpqGglFLqNA0FpZRSp2koKKWUOk1DQSml1GleeZ6Cco2iIyVs33ec3QdPsu9YKRWVBlulISo0kLiIYBJiw+jYIoqo0CCrS1VKuYiGgh87dKKMxZuL+DazmDW7DlJ4uKROy7VtEs6g9k25oGM8Q1KaEh6sv0ZK+Qr93+xnKisN32UX8/6KnXyfVUxFpaFFdCh928XRu00sKc2iSGgcRnxUCEEBQoAIR0sq2H+8lF0HTrC54Ajr8w/z1foCPlq5i4jgAH7bsxW/79ua89o0trp7SikHiTfeeS0tLc3oZS7qx1Zp+Hztbl5bms2O/SeIjwphYu9ExnZrQY/EGESkXp9XVlFJ+o4DfL52N19vKORkuY1+7eK4c0QKgzs0qffnKaVcT0RWG2PSztlGQ8H3fbO5iKcXbCVn7zG6JUQzdVh7xnRtQXCgc44zOFpSzuz0fN7+fhtFR0oZkBzHYxd3pVOLaKd8vlLKOTQU/NyewyX8ec4mFmYU0T4+gvtHd2RMtxYu+yu+tMLGrJV5vPhNFkdOlnPtwCQeuKgjESG6lVIpT1CXUND/rT7IGMMnq/N5/KvNlNsq+eOYjtw8NJmgANcegRwSGMCUQUlM6NWK5xdlMXPFDpZl7uX5ST1JS4pz6bqVUs6h5yn4mGOlFdzz8Toe+GQDXVtFs+ieYdw2vIPLA6G62PBgnrikG7NuHoCt0vC7t1bw3MJMbJXeNypVyt/oSMGHZBcdZer7q9m5/zj3jkrlDxd0IKCRdTt8+yc3YcHdw3h8TgavLcthw+7DvHJFL2LDgy2rSSl1bjpS8BHfZRVz2es/crSkgo9uHsCdF6ZYGginRIYE8uyknvzt0u6s2LaPi1/7gcw9R60uSyl1FhoKPuD9n3Zyw79XkdA4jC9vH0z/5CZWl/QrV/Vvw6ypAykptzHpzR/5OXe/1SUppWqgoeDFjDG8uiSbP32xieGp8Xxy6yASYsOsLuus+rRtzGe3DaJpVAjX/HMlCzbtsbokpdQZnBIKIjJGRDJFJEdEptfwvojIK/b3N4hIb/v8UBFZKSLrRSRDRB53Rj3+wBjD0wu28vziLC7rncBb1/Qh0gsO/UxsHM6n0wbRtVU0t32wms/X5ltdklKqGodDQUQCgBnAWKALcKWIdDmj2Vggxf6YCrxhn18KjDDG9AR6AWNEZICjNfk6YwyPzcngre9yuXpAG567vCeBbjy6yFGNI4L54Kb+DEhuwn2z1/PF2t1Wl6SUsnPGN0k/IMcYk2uMKQNmARPOaDMBeM9U+QmIFZGW9ulj9jZB9ocet3gOxhiemr+VmSt2cvPQdjwxoRuNPGCHcn2FBwfyjyl96d+uCffOXseX6zQYlPIEzgiFBCCv2nS+fV6d2ohIgIisA/YCi40xP9e0EhGZKiLpIpJeXFzshLK902tLc3j7+1ymDGzLw+M6e/U1hsKCA/jHdWn0axfHPR+vY8GmQqtLUsrvOSMUavpWOvOv/bO2McbYjDG9gESgn4h0q2klxpi3jTFpxpi0+Ph4R+r1Wv/+YfvpfQh//m1Xrw6EU8KDA/nndX3p2TqWO2etY+X2A1aXpJRfc0Yo5AOtq00nAgX1bWOMOQR8C4xxQk0+Z/7GQh77ajOjuzTn7xN7eOUmo7M5tSkpMTaMm2auIqtIz2NQyirOCIVVQIqItBORYOAKYM4ZbeYA19qPQhoAHDbGFIpIvIjEAohIGDAS2OqEmnzKurxD3P3xOnq3ieWVK8/zqp3KdRUXEczMG/oREhTAlH+uZE8db/ijlHIuh79djDEVwO3AQmALMNsYkyEi00Rkmr3ZPCAXyAHeAW6zz28JLBORDVSFy2JjzNeO1uRL8g+e4KaZ6TSLDuGda9MIDQqwuiSXaR0Xzr+v78uRk+Xc8n46JeU2q0tSyu/opbM92PHSCi57/UcKDp/k89sG0aFZlNUlucWijD1MfX81l56XwAu/6+kT+06U8gR1uXS2722H8BHGGP74yQay9x7ljcl9/CYQAEZ3bcG9o1L5fO1u3lmea3U5SvkVDQUP9Y//bWfuxkL+OKYTQ1KaWl2O290xogPjurfg6flb+TZzr9XlKOU3NBQ80E+5+3lq/lbGdG3BLcOSrS7HEiLCc5N6kto8irs/XkfBoZNWl6SUX9BQ8DBFR0q4/cM1tG0SzrOTevj19vTw4EDeuLoP5RWV3PHRWsptlVaXpJTP01DwILZKw92z1nGizMZbV/chKjTI6pIs165pBE9N7MHqnQd5blGm1eUo5fM0FDzIW99vY0Xufh67uCspzf1nx3JtLu7Zisn92/DWd7ks26r7F5RyJQ0FD7Eu7xAvLMpifI+WTOqTaHU5HudPv+lC55bR3Dt7HXuP6IltSrmKhoIHOFZawZ0fraV5dCh/u7S7X+9HOJvQoABeu+o8Tpbb+OOnG/DG82uU8gYaCh7g/77cRP7BE7x0RS9iwnQ/wtm0j4/kobGd+TazmA9+3mV1OUr5JA0Fiy3YVMhna3Zz+4gU+ibFWV2Ox7tmQFuGpjTlyblb2L7vuNXlKOVzNBQsdOB4GY9+sYmuraK5Y0QHq8vxCo0aCc9e3pPgwEbc8/E6KvQwVaWcSkPBQn+ek8Hhk+U8N6knQT545VNXaRETyhOXdGNd3iHe+Hab1eUo5VP0m8giCzYV8tX6Au4YkULnltFWl+N1Lu7Zit/0aMkrS7PJ1vsvKOU0GgoWqL7Z6Nbh7a0ux2s9dnFXIkICefDTDdgq9WgkpZxBQ8ECutnIOZpGhvDn33Zhza5DvLdih9XlKOUT9BvJzZZt3ctX6wu4/QLdbOQMl/RKYHjHeJ5dmEnegRNWl6OU19NQcKOTZTb+9OUmOjSL1M1GTiIiPHlpdwR4+PONelKbUg7SUHCjl5dkk3/wJE9e0o3gQP2nd5aE2DAeHNuJ5dn7+GLdbqvLUcqr6TeTm2TuOcq7y3O5vE8i/ZObWF2Oz7m6f1t6tY7lyblbOXyy3OpylPJaTgkFERkjIpkikiMi02t4X0TkFfv7G0Skt31+axFZJiJbRCRDRO5yRj2eprLS8MjnG4kKDeThcZ2tLscnNWok/PWSbhw4XsoLeoltpRrM4VAQkQBgBjAW6AJcKSJdzmg2FkixP6YCb9jnVwD3GWM6AwOAP9SwrNebnZ5H+s6DPDSuM3ERwVaX47O6JcRw9YC2vP/TTjbtPmx1OUp5JWeMFPoBOcaYXGNMGTALmHBGmwnAe6bKT0CsiLQ0xhQaY9YAGGOOAluABCfU5DH2Hyvlqflb6ZcUp5fEdoP7RnckLiKYR7/YRKWeu6BUvTkjFBKAvGrT+fz6i73WNiKSBJwH/OyEmjzGc4syOV5awV8v7aaXxHaDmLAgHh7XmXV5h/g4Pa/2BZRSv+CMUKjpm+7MP9HO2UZEIoFPgbuNMUdqXInIVBFJF5H04uLiBhfrTpt2H2bWqjyuHZhEqt5JzW0uPS+Bfu3ieHr+Vg4cL7O6HKW8ijNCIR9oXW06ESioaxsRCaIqED4wxnx2tpUYY942xqQZY9Li4+OdULZrGWN4bE4GceHB3DUyxepy/IpI1U7nY6UVvLBYdzorVR/OCIVVQIqItBORYOAKYM4ZbeYA19qPQhoAHDbGFErV9pR/AFuMMS84oRaPMWd9Aek7D3L/RR31xjkWSG0exdX92/Dhz7vYuqfGwadSqgYOh4IxpgK4HVhI1Y7i2caYDBGZJiLT7M3mAblADvAOcJt9/mDgGmCEiKyzP8Y5WpPVTpRV8NS8rXRtFc3v0lrXvoByiXtGpRIdFsQTX2/WM52VqqNAZ3yIMWYeVV/81ee9We21Af5Qw3L/o+b9DV7tzW+3sedICa9edR4BjXyue14jNjyYe0am8uc5GSzeXMTori2sLkkpj6dnNDtZ3oETvPV9Lhf3bKW31/QAk/u3IaVZJE/O20Jphc3qcpTyeBoKTvbU/C00EuGhcZ2sLkUBgQGN+NNvurBz/wn+/cMOq8tRyuNpKDhR+o4DzNu4h2nnt6dlTJjV5Si7YanxXNipGa8uzaH4aKnV5Sjl0TQUnMQYw9/mbaFZVAg3D2tndTnqDI+M70xphY0XFmdZXYpSHk1DwUkWbNrDml2HuHdUKuHBTtl/r5woOT6Syf3bMjs9j5y9ek9npc5GQ8EJyioqeWbBVlKbRzJJD0H1WHeM6EBYUADPLNAT2pQ6Gw0FJ/ho5S527D/BQ2M76yGoHqxJZAi3Dm/P4s1FrNpxwOpylPJIGgoOOlJSzstLshmY3IThHT3/8hv+7obB7WgeHcLf5m3RE9qUqoGGgoPe+m4bB46X8fC4znoVVC8QFhzAvaNSWbvrEAsz9lhdjlIeR0PBAYWHT/Lu8u1M6NWK7okxVpej6mhi70RSmkXyzIJMym2VVpejlEfRUHDA84uyMAbuH93R6lJUPQQGNGL62E5s33ecWSt3WV2OUh5FQ6GBsoqO8umafK4d2JbWceFWl6PqaUSnZvRrF8fLS7I5VlphdTlKeQwNhQZ6flEmEcGB/OGCDlaXohpARHh4XGf2HSvj3eW5VpejlMfQUGiA9XmHWJhRxM1Dk2kcEWx1OaqBerWO5aKuzXl3+XYO6h3alAI0FBrkuUWZxEUEc+NQvZyFt7tvdEeOl1Xw5nfbrC5FKY+goVBPK7btZ3n2Pm49vz2RIXo5C2+X2jyKS3olMHPFDoqOlFhdjlKW01CoB2MMzy3KpHl0CNcMbGt1OcpJ7h6ZQoXN8NrSHKtLUcpyGgr1sCxzL6t3HuTOC1MIDQqwuhzlJG2bRPC7vq2ZtWoXeQdOWF2OUpbSUKijykrDcwuzaBMXrvdd9kF3jkhBRHh5SbbVpShlKQ2FOpq3qZDNhUe4Z1QKQQH6z+ZrWsSEcu2Atny2Jl8vra38mlO+3URkjIhkikiOiEyv4X0RkVfs728Qkd7V3vuniOwVkU3OqMUVKmyVvLAoi9TmkVzcM8HqcpSL3Dq8PWFBAby4WEcLyn85HAoiEgDMAMYCXYArRaTLGc3GAin2x1TgjWrv/RsY42gdrvTZmt3k7jvOfaM76qWxfViTyBBuHNKOuRsL2bT7sNXlKGUJZ4wU+gE5xphcY0wZMAuYcEabCcB7pspPQKyItAQwxnwPeOzF7UsrbLy8JJueiTGM7tLc6nKUi900LJmYsCCeX6Q34lH+yRmhkADkVZvOt8+rb5tzEpGpIpIuIunFxcUNKrQhZqfns/vQSe4b3VEvje0HokODuOX8ZJZlFrNm10Gry1HK7ZwRCjV9U55595K6tDknY8zbxpg0Y0xafLx7bmZTWmHj9WU59GnbmKEpTd2yTmW9KQOTiIsI5qVvdN+C8j/OCIV8oPoxmolAQQPaeJzZ6fkUHi7h7pEpOkrwIxEhgUwdlsz3WcWs3qmjBeVfnBEKq4AUEWknIsHAFcCcM9rMAa61H4U0ADhsjCl0wrpd5tQoIa1tY4Z00FGCv7l2YFv7aCHL6lKUciuHQ8EYUwHcDiwEtgCzjTEZIjJNRKbZm80DcoEc4B3gtlPLi8hHwAqgo4jki8iNjtbkDLNX5dlHCak6SvBD4cGB3DIsmeXZ+1i902OPg1DK6cQbb16elpZm0tPTXfb5pRU2zv/7tyQ2DuO/0wZqKPipE2UVDPv7Mjq1iOY/N/W3uhylHCYiq40xaedqo6fm1uDjVXnsOaKjBH9XNVpoz/9y9rFqh44WlH/QUDhDSbmN15dto29SYwZ3aGJ1OcpiVw9oS9PIEN23oPyGhsIZZqfrKEH9f2HBAUw7P5kfcvazcruOFpTv01CopvooYVB7HSWoKpP7V40WXlysowXl+zQUqtF9CaomYcEB3Dq8PSty9/NT7n6ry1HKpTQU7ErKbbz+bQ79kuJ0lKB+ZXL/NjSL0n0LyvdpKNh9vCqPoiOlevayqlFoUNVo4afcA/y4bZ/V5SjlMhoK/HKUMFBHCeosruzXhubRIbys10RSPkxDAZi1cpeOElStQoMCmHZ+e37efoAV23TfgvJNfh8KJeU23vhuG/3a6ShB1e7KflX7Fl5eovsWlG/y+1D47+p8io6UcteFOkpQtTs1Wvgp94AeiaR8kl+HQllFJW9+u43ebWL1iCNVZ1f1b0N8lO5bUL7Jr0Phi7W72X3oJHfoKEHVw6nRworc/fysowXlY/w2FCpslbz+bQ7dEqIZnuqeO7kp3zH51GhhiY4WlG/x21CYu7GQHftPcPsFOkpQ9RcaFMAtw5L5cZteE0n5Fr8MhcpKw2tLc+jYPIrRXZpbXY7yUqeuiaRHIilf4pehsDBjD9l7j/GHER1o1EhHCaphql9BVe+3oHyFX4bCm9/n0q5pBOO7t7S6FOXlqkYLwXokkvIZfhcK6/IOsT7vENcNSiJARwnKQWHBAafvzpauowXlA5wSCiIyRkQyRSRHRKbX8L6IyCv29zeISO+6Luts7/24g8iQQCb2SXT1qpSfmDygDU0igvVIJOUTHA4FEQkAZgBjgS7AlSLS5YxmY4EU+2Mq8EY9lnWafcdK+XpDIRN7JxAZEuiq1Sg/Ex4cyC3nJ7M8ex+rd+poQbnGsdIK7v/verYVH3PpepwxUugH5Bhjco0xZcAsYMIZbSYA75kqPwGxItKyjss6zayVuyizVXLNwCRXrUL5qasHtKVJRDAv6b4F5SL/+Wknn6zO58jJcpeuxxmhkADkVZvOt8+rS5u6LAuAiEwVkXQRSS8uLm5QoS1jwriyXxs6NIts0PJKnU14cCBTh1WNFtbsOmh1OcrHnCyz8c73uQxNacp5bRq7dF3OCIWa9taaOrapy7JVM4152xiTZoxJi49v2BnIE/sk8tRl3Ru0rFK1uWZgW+Ii9Egk5Xwf/LyT/cfLuOvCFJevyxmhkA+0rjadCBTUsU1dllXKK4QHB3Lz0GS+yypmrY4WlJOUlNt4+/tcBiY3IS0pzuXrc0YorAJSRKSdiAQDVwBzzmgzB7jWfhTSAOCwMaawjssq5TWuHdiWxuFBeiSScprZ6XnsPVrKHRd2cMv6HA4FY0wFcDuwENgCzDbGZIjINBGZZm82D8gFcoB3gNvOtayjNSlllYiQQG4elsy3mcWsyztkdTnKy5VW2Hjj2230TWrMwGT3XN5fjKlxE75HS0tLM+np6VaXoVSNjpVWMPSZpfRqHcu/ru9ndTnKi33w804e+XwT793Qj2FOuJqziKw2xqSdq43fndGslKtFhgRy09BkluloQTmgrKKS15dt47w2sQxNaeq29WooKOUCUwYlERsexCu6b0E10Kdr8tl96KTbbxWsoaCUC0SGVB2JtHTrXtbraEHVU7mtkhnLcujZOpbz3XwTMA0FpVzk2oFtdbSgGuSzNfnkHzzJ3RbcKlhDQSkXiQoN4qYh7ViydS8b8g9ZXY7yEuW2Sl5blkOPxBiGd3T/rYI1FJRyoSmDkogJ09GCqrvP1+4m78BJ7hxhza2CNRSUcqFTo4VvtuxlY/5hq8tRHq7Cvi+hW0I0F3ZuZkkNGgpKudiUwUlEhwbqWc6qVl+sK2Dn/hOWjRJAQ0Epl4sODeKmocl8s6WITbt1tKBqdmqU0KVlNKO6NLesDg0FpdzgOh0tqFp8taGA7fuOc6cFRxxVp6GglBtEhwZx45BkFm/W0YL6NVul4dWlOXRqEcVoC0cJoKGglNtcNziJqNBAPRJJ/crXGwrILT7OXRem0KiRdaME0FBQym1iwoK4cUg7Fm0uIqNARwuqiq3S8MqSbDo2j+Kiri2sLkdDQSl3un5wOx0tqF+Yu7GQbcXHuePCDpaPEkBDQSm3igkL4vrB7ViYUcTmgiNWl6MsVllpeHVJNinNIhnXraXV5QAaCkq53Y2D2xEVoqMFBfM2FZK99xh3eMC+hFM0FJRys5jwIK4fnMSCjD1sKdTRgr+qGiXk0D4+gvHdPWOUABoKSlnihiE6WvB3czcWkll0lDsvTCHAQ0YJoKGglCViw4O5bnAS8zftYeseHS34G1ul4aVvskhtHslve7SyupxfcCgURCRORBaLSLb9ufFZ2o0RkUwRyRGR6dXmTxKRDBGpFJFz3jdUKV9z45B2ROpowS99uW4324qPc8/IVI/Zl3CKoyOF6cASY0wKsMQ+/QsiEgDMAMYCXYArRaSL/e1NwGXA9w7WoZTXiQ0P5rpBSczbuIfMPUetLke5Sbmtkpe+yaZLy2iPOC/hTI6GwgRgpv31TOCSGtr0A3KMMbnGmDJgln05jDFbjDGZDtaglNfS0YL/+XR1PrsOnOC+0Z43SgDHQ6G5MaYQwP5c0wXAE4C8atP59nn1IiJTRSRdRNKLi4sbVKxSnqZxRDBTBrVl3qZCHS34gdIKG68uzaFX61hGdLLmfgm1qTUUROQbEdlUw2NCHddRUxSa+pUJxpi3jTFpxpi0+Hj336JOKVe5aUgy4UEBvLJURwu+bvaqPHYfOsm9o1ItvRLquQTW1sAYM/Js74lIkYi0NMYUikhLYG8NzfKB1tWmE4GCeleqlI+qGi0k8cZ328gqOkpq8yirS1IuUFJu47VlOfRNaszQlKZWl3NWjm4+mgNMsb+eAnxZQ5tVQIqItBORYOAK+3JKKbubhiYTFhSg+xZ82H9+2knRkVLuHdXRY0cJ4HgoPA2MEpFsYJR9GhFpJSLzAIwxFcDtwEJgCzDbGJNhb3epiOQDA4G5IrLQwXqU8kpx9tHC3I2FZBfpvgVfc6Ksgje/28bgDk0Y2L6J1eWck0OhYIzZb4y50BiTYn8+YJ9fYIwZV63dPGNMqjGmvTHmyWrzPzfGJBpjQowxzY0xFzlSj1Le7OZTo4WlOVaXopxs5o872XesjHtHdbS6lFrpGc1KeYhTo4WvNxToNZF8yNGSct76fhvDO8bTp22N5/d6FA0FpTzItGHtiQoJ5NmFevqOr3j7+1wOnSjnPi8YJYCGglIeJSY8iGnD27N0615Wbj9gdTnKQXuPlvDu8u38pkdLuifGWF1OnWgoKOVhrh/UjubRITw9fwvG1PuUHuVBXl2SQ7mtkvtHe8coATQUlPI4YcEB3HVhKmt2HWLx5iKry1ENtGPfcT5auYvf921NUtMIq8upMw0FpTzQ79ISSW4awbMLM7FV6mjBGz2/OIuggEbcdWGK1aXUi4aCUh4oMKAR91/Ukey9x/hsTb7V5ah62rT7MF+tL+CGIUk0iw61upx60VBQykON7daCHokxvLg4i5Jym9XlqHp4ZsFWYsODuOX89laXUm8aCkp5KBHhwTGdKDhcwn9+2ml1OaqOfszZx/LsffxheAeiQ4OsLqfeNBSU8mCDOzRlaEpTZizL4UhJudXlqFoYY3hmwVZaxYRyzcC2VpfTIBoKSnm4P17UiYMnynnz221Wl6JqMW/jHtbnH+buUamEBgVYXU6DaCgo5eG6J8ZwSa9WvPu/7eQfPGF1OeosSsptPL1gC51aRDGxd6LV5TSYhoJSXuCBMZ0Q0MtfeLCZP+4g78BJHh3fhQAPvM1mXWkoKOUFEmLDuHloMl+uK2DtroNWl6POsP9YKa8tzWFEp2YM8eAb6NSFhoJSXmLa8PY0jQzhr3P18hee5qVvsjlRbuPhcZ2sLsVhGgpKeYnIkEDuH53K6p0Hmbdxj9XlKLvsoqN8uHIXk/u3oUMz77+VqoaCUl5kUlprOrWI4ukFWyit0BPaPMHf5m0hPDjA6y5ncTYaCkp5kYBGwqPju5B34CQzf9xhdTl+7/usYpZlFnPHiA40iQyxuhyn0FBQyssMSWnKiE7NeHVJDsVHS60ux2+V2yr569zNtI4LY8qgJKvLcRqHQkFE4kRksYhk259rvNeciIwRkUwRyRGR6dXmPysiW0Vkg4h8LiKxjtSjlL94dHxnSipsPD1/q9Wl+K33Vuwkq+gYj4zrQkigd56oVhNHRwrTgSXGmBRgiX36F0QkAJgBjAW6AFeKSBf724uBbsaYHkAW8JCD9SjlF5LjI7lpaDKfrsln9U69Q5u77T1awkuLsxiWGs9FXZtbXY5TORoKE4CZ9tczgUtqaNMPyDHG5BpjyoBZ9uUwxiwyxlTY2/0EeO9pgEq52e0XdKBlTCj/92WG3nPBzZ6Zn0lJhY3HftsFEe89Ua0mjoZCc2NMIYD9uVkNbRKAvGrT+fZ5Z7oBmH+2FYnIVBFJF5H04uJiB0pWyjdEhATyyPjOZBQc4cOVu6wux2+s3nmAT9fkc9PQZJLjI60ux+lqDQUR+UZENtXwmFDHddQUo7/4s0ZEHgEqgA/O9iHGmLeNMWnGmLT4+Pg6rlop3za+e0sGJjfhuYWZHDheZnU5Ps9WafjTFxm0jAnljhEdrC7HJWoNBWPMSGNMtxoeXwJFItISwP68t4aPyAdaV5tOBApOTYjIFOA3wGSjp2kqVS8iwuMTunK8tIJnF+pOZ1f78OedbC48wiPjOxMeHGh1OS7h6OajOcAU++spwJc1tFkFpIhIOxEJBq6wL4eIjAEeBC42xujlH5VqgNTmUVw3KIlZq/JYo9dFcpl9x0p5blEWg9o3YXz3llaX4zKOhsLTwCgRyQZG2acRkVYiMg/AviP5dmAhsAWYbYzJsC//GhAFLBaRdSLypoP1KOWX7hqZQvOoUB7+bCPltkqry/FJT3y9mZNlNv4yoavP7VyuzqHxjzFmP3BhDfMLgHHVpucB82po55sb5ZRys6jQIP4yoStT31/N29/n8ocL9L+WMy3L3MuX6wq4e2SKT1zf6Fz0jGalfMTori0Y260FLy/JZvu+41aX4zOOl1bw6Oeb6NAskluHt7e6HJfTUFDKhzx2cVdCAhvxyOcb9fLaTvLi4ix2HzrJU5d196kzl89GQ0EpH9I8OpTpYzvx47b9fLpmt9XleL0N+Yf45w/buap/G/omxVldjltoKCjlY67s24a+SY3569zNesE8B5TbKpn+6UaaRobw4Bjvv3lOXWkoKOVjGjUSnrqsOyfKbLoZyQEzluWwufAIf5nQlZiwIKvLcRsNBaV8UIdmUTwwuiOLNhfxmW5GqreN+Yd5bWkOl/RqxZhuvntOQk00FJTyUTcMaUffpMY89lUGBYdOWl2O1ygpt3Hff9fRJDKYxy/uZnU5bqehoJSPCmgkPDepJ7ZKw4OfbtDNSHX04uIssoqO8czEHsSE+89mo1M0FJTyYW2bRPDwuM4sz97Hf37WK6nWJn3HAd5ensuV/dowvGNNF332fRoKSvm4yf3bMDSlKU/O3Ux20VGry/FYR0rKuWf2OhJiw3hkfGery7GMhoJSPk5EeH5STyKCA7n9w7WUlNusLsnjGGN46NONFBwq4eUrziMyxDevgFoXGgpK+YFm0aE8/7ueZBYd5S9fb7a6HI8za1UeczcWct/oVPq0rfFW835DQ0EpPzG8YzNuOT+ZD3/exdwNhVaX4zGyio7y2JwMhnRoyrRhvn9to9poKCjlR+4f3ZFerWOZ/ukG8g7oLUxOlFVw+4driAoN5IXf96RRI9+9JHZdaSgo5UeCAhrx6pXngcAt76/mZJn/7l8wxvDAJxvI2XuMF3/fi2ZRoVaX5BE0FJTyM63jwnnlivPYsueIX5+/8Pb3uczdUMgDF3ViaIre9/0UDQWl/NAFnZpx/+iOzFlfwLvLt1tdjtstzy7mmQVbGd+9JdPOT7a6HI+ioaCUn7pteHvGdW/BU/O3sDy72Opy3GbX/hPc8dFaUppF8ffLe/j0rTUbQkNBKT8lIjx7eU9SmkVx2wdryPKDE9sOnSjjun+vxBh465o+RPjx+Qhn41AoiEiciCwWkWz7c40H+IrIGBHJFJEcEZlebf4TIrJBRNaJyCIRaeVIPUqp+okICeQf16URGhTA9f9aRdGREqtLcpmSchtT31tN/oGTvH1NH5KaRlhdkkdydKQwHVhijEkBltinf0FEAoAZwFigC3CliHSxv/2sMaaHMaYX8DXwfw7Wo5Sqp8TG4fzrur5Vf0X/axVHS8qtLsnpKiurjjRaueMAz/2uJ/2Tm1hdksdyNBQmADPtr2cCl9TQph+QY4zJNcaUAbPsy2GMOVKtXQTgn4dBKGWxbgkxvH51H7KKjnLrf9b41KUwjDE8MXczX60vYPrYTlzcUzdInIujodDcGFMIYH+u6bKCCUBetel8+zwARORJEckDJqMjBaUsc35qPM9M7MH/cvZx2wdrKKuotLokhxljeGZBJv/6YQfXD07ilmF6pFFtag0FEflGRDbV8JhQx3XUtGv/9IjAGPOIMaY18AFw+znqmCoi6SKSXlzsP0dKKOVOl/dJ5MlLu7F0617u+GgN5TbvDoaXl2Tz5nfbmNy/Df/3my56pFEd1BoKxpiRxphuNTy+BIpEpCWA/XlvDR+RD7SuNp0IFNTQ7kNg4jnqeNsYk2aMSYuP1xNNlHKVyf3b8thvu7Awo4i7Z63zyhGDMYZXlmTz0jfZXN4nkScmdNNAqCNHNx/NAabYX08BvqyhzSogRUTaiUgwcIV9OUQkpVq7i4GtDtajlHKC6wa349HxnZm7sZCb3kvnRFmF1SXVWWWl4Ymvt/DC4iwuOy+BZyb20Gsa1YOjofA0MEpEsoFR9mlEpJWIzAMwxlRQtVloIbAFmG2MyTi1vH1T1AZgNHCXg/UopZzkpqHJPDOxO//LLmbyuz9z6ESZ1SXVqtxWyQOfbOCfP2zn+sFJPDepJwEaCPUi3njdk7S0NJOenm51GUr5hQWb9nDnR2tJaBzGO9em0aFZpNUl1ejg8TJu+2ANK3L3c8/IVO68sINuMjqDiKw2xqSdq42e0ayUOqcx3Vrwwc39OXKynEtn/MDSrUVWl/QrWUVHmTDjB1bvPMjzk3py18gUDYQG0lBQStWqb1Icc+4YQpsm4dw4M53nF2V6xJFJxhhmp+dxyYwfOFFmY9YtA5jYJ9HqsryahoJSqk4SYsP4ZNogJvZO5NWlOVz+5gq27ztuWT2HT5Rz+0dr+eMnG+ieEMNXdwymdxv/vpWmM2goKKXqLCw4gOcm9WTGVb3Zse84415ezpvfbXPrYavGGOasL2Dki9+xcNMeHrioIx/ePICWMWFuq8GX6SUClVL1Nr5HS3q3jeVPX2Tw9PytfLI6n0fGd2Z4arxLt+Vv2n2Yp+dv5X85++ieEMM/p/Sle2KMy9bnj/ToI6WUQ5ZsKeLxrzaz68AJerWO5a6RKZyfEu/UcwM25h9mxrIcFmTsISYsiPtGpzK5f1s93LSe6nL0kYaCUsphZRWVfLI6nxnLcth96CRt4sK5sl8bftuzJYmNwxv0mQePl7F4cxEfrNzF+rxDRIUEcuPQdtwwpB3RoUFO7oF/0FBQSrlVWUUl8zcV8sHPu1i5/QAAnVtGc35qPD0SY+jWKoaWsaEEBfxyd6at0lB8tJQte46wPu8QP+XuZ9WOg9gqDe3jI7hmQFsu65OoYeAgDQWllGW27zvO4s17+GbzXtbmHaTcVvVdIwJx4cFEhARSaQzltkr2HSvDVvn/3+/YPIqRnZszqktzeiTG6DkHTqKhoJTyCKUVNjL3HGVzwRH2HClh79FSTpRW0EiEgEZCs+gQWkSH0r5ZJD0SY4nU22S6RF1CQf/llVIuFxIYQI/EWHokxlpdiqqFnqeglFLqNA0FpZRSp2koKKWUOk1DQSml1GkaCkoppU7TUFBKKXWahoJSSqnTNBSUUkqd5pVnNItIMbCzgYs3BfY5sRxvoH32D9pn/+BIn9saY+LP1cArQ8ERIpJe22nevkb77B+0z/7B1X3WzUdKKaVO01BQSil1mj+GwttWF2AB7bN/0D77B5f22e/2KSillDo7fxwpKKWUOgsNBaWUUqd5XSiIyBgRyRSRHBGZXsP7IiKv2N/fICK9a1tWROJEZLGIZNufG1d77yF7+0wRucj1Pfw1d/ZZREaJyGoR2Wh/HuGeXv6qT279OdvfbyMix0Tkftf2rmYW/G73EJEVIpJh/3mHur6Xv+qTO3+3g0Rkpr2vW0TkIff08ld9ckWfJ9l/jpUiknbG59XvO8wY4zUPIADYBiQDwcB6oMsZbcYB8wEBBgA/17Ys8Hdguv31dOAZ++su9nYhQDv78gE+3ufzgFb2192A3b7+c672mZ8C/wXu9/U+U3XXxQ1AT/t0Ez/43b4KmGV/HQ7sAJJ8pM+dgY7At0Batc+q93eYt40U+gE5xphcY0wZMAuYcEabCcB7pspPQKyItKxl2QnATPvrmcAl1ebPMsaUGmO2Azn2z3Ent/bZGLPWGFNgn58BhIpIiIv6djbu/jkjIpcAuVT12Qru7vNoYIMxZj2AMWa/Mcbmor6djbv7bIAIEQkEwoAy4IhrunZWLumzMWaLMSazhvXV+zvM20IhAcirNp1vn1eXNudatrkxphDA/tysHutzNXf3ubqJwFpjTGmDq28Yt/ZZRCKAB4HHnVR/Q7j755wKGBFZKCJrROSPTulF/bi7z58Ax4FCYBfwnDHmgOPdqBdX9dmR9f1CYC0f6GmkhnlnHlN7tjZ1WbYh63M1d/e56gNFugLPUPUXpbu5u8+PAy8aY46J1LS4W7i7z4HAEKAvcAJYIiKrjTFLaivUidzd536ADWgFNAaWi8g3xpjc2gp1Io//DvO2UMgHWlebTgQK6tgm+BzLFolIS2NMoX2Ytrce63M1d/cZEUkEPgeuNcZsc0ov6sfdfe4PXC4ifwdigUoRKTHGvOaMztSRFb/b3xlj9gGIyDygN+DOUHB3n68CFhhjyoG9IvIDkEbVZkN3cVWfHVnfL7lqh4orHlSFWC5VO0xO7Wjpekab8fxyJ83K2pYFnuWXO6b+bn/dlV/upMnF/Tvj3N3nWHu7if7ycz7jcx/Dmh3N7v45NwbWULXDNRD4Bhjv431+EPiX/bMigM1AD1/oc7Vlv+WXO5rr/R1myX96B/9RxwFZVO1Ff8Q+bxowzf5agBn29zee8Q/0q2Xt85tQ9RdStv05rtp7j9jbZwJjfb3PwKNUbXddV+3RzJf7fMZ6H8OCULDod/tqqnasb6KGgPS1PgORVB1dlkFVIDzgQ32+lKpRQSlQBCys9l69vsP0MhdKKaVO87ajj5RSSrmQhoJSSqnTNBSUUkqdpqGglFLqNA0FpZRSp2koKKWUOk1DQSml1Gn/D4HQe80Bgb69AAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(tran_analysis.time, tran_analysis['in1'])" + ] + }, + { + "cell_type": "code", + "execution_count": 30, + "id": "56ee9649", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 30, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAD4CAYAAADcpoD8AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA0i0lEQVR4nO3deXwU9f3H8dcnF4EQIBAC5OJKOMKtATywciriBd5UrVqPWhXqVcWqrW1tq/SyVrFgFUGl3oq2CgqCKCAlCOGQIweBJCAEkCMXuT6/P3b4NWJCAkl2Nruf5+Oxj92dme/M50vCvvOdmZ0RVcUYY4zxliC3CzDGGBNYLHiMMcZ4lQWPMcYYr7LgMcYY41UWPMYYY7wqxO0CfFV0dLR269bN7TKMMaZZWbNmzT5V7XiiZSx4atGtWzfS0tLcLsMYY5oVEdlR1zK2q80YY4xXWfAYY4zxKgseY4wxXmXBY4wxxqsseIwxxniVBY8xxhivalDwiEh7EflERDKc56halhsvIltFJFNEptXVXkSGicg655EuIpOqtVnqrOvY/BhneqKILBGRtSKyXkQmVGtTWW359xvSZ2OMMQ0jDbktgohMBw6o6hNOoESp6oPHLRMMbAPGAXnAamCyqn5dW3sRaQWUqWqFiHQB0oFY5/1S4H5VTTtuO7OAtar6nIikAB+qajdnXqGqtj6ZvqWmpqp9j6dpqCo79hez80Axuw+VcKConMqqKqoU2oSH0KF1C2LbhdO7cxtat7CvmhnTnIjIGlVNPdEyDf1ffSkw0nk9B1gKPHjcMsOATFXNdop6zWn3dW3tVbW4WvtwoD7pqEAb53VbYFe9e2GaXFZBIYs37+HzjH2k5x7kcGlFvdp17dCKM7p34Jxe0YzsHWNBZIwfaOj/4k6quhtAVXcf2+11nDggt9r7PGB4Xe1FZDjwItAVuF5Vq39SzRaRSuBt4HH1DNseAz4WkSlABDC22vLhIpIGVABPqOp7NXVGRG4DbgNITEysR/fNiRwqKefdr/L4139z2brnCAB9Okdy4cBYBsW3pUfH1nRpG06H1mGEBAURJHC4tIL9hUfZsb+YrXuOsD7vIB9t3M3rabmEhwZxXkpnrkpN4OykDoiIyz00xpyKOoNHRBYBnWuY9XA9t1HTp0OdIxhVXQX0E5G+wBwR+UhVS4FrVTVfRCLxBM/1wFxgMvCSqv5ZRM4EXhaR/qpaBSSq6i4R6QF8KiIbVDWrhm3OAmaBZ1dbPftnjrP3SCmzPsvm1VU7KSmvZGB8W359ST/GpnQirl3LE7ZtHxFG+4gwkjtFMjalEwAVlVWszT3I/HX5fJC+m/fTd9Evtg23n9uTC/p3JiTYzpExpjmpM3hUdWxt80Rkj4h0cUYrXYC9NSyWByRUex/P/3aD1dleVTeLSBHQH0hT1Xxn+hERmYdnV95c4GZgvDNvpYiEA9HAXlXd5UzPdo4RDQG+FzymYQqPVvDMp5nMXr6d8soqJg6J48dnd6d/XNsGrTckOIih3doztFt7Hr0ohffW5jNzWTZT/rWW3p0iefSiFEYkRzdSL4wxTa2hfyq+D9zgvL4BmF/DMquBZBHpLiJhwDVOu1rbO8uGOK+7Ar2BHBEJEZFoZ3oocBGw0Wm/ExjjzOuL59hQgYhEiUgLZ3o0cDae40umkagq76fvYsyfl/KPz7KYMKALi+8byV+uGtzg0Dlei5Bgrh6ayKJ7zuXZH55GSXkl172wilvmrCbv2+K6V2CMcV1Dz2rrALwBJOL54L9SVQ+ISCzwT1Wd4Cw3AXgKCAZeVNXf1dH+emAaUA5UAb9R1fdEJAJYBoQ661oE3Kuqlc6ZbM8DrfHsyntAVT8WkbOAmc56goCnVPWFuvpmZ7XVT8GRo0x7ez2Lt+ylf1wbfnNpf05LrPGs+iZxtKKS2ctz+PviDESERy/qy1WpCXb8xxiX1OestgYFjz+z4Knboq/38ODb6zlytIJp4/tww1ndCA5y5wM/90AxD7y1npXZ+xndJ4Y/XTmI9hFhrtRiTCCrT/DYUVlz0iqrlCcXbOGWuWnEtAnn31NG8OMR3V0LHYCE9q149ZbhPHZxCl9k7uPiv39Beu5B1+oxxtTOgseclEMl5dw8ZzXPLc1i8rBE3rvzLHp1inS7LACCgoQbz+7O27efBcCV/1jJ66t3ulyVMeZ4Fjym3nbuL2bSs8v5ImMfj0/szx8uG0CLkGC3y/qeAfFt+feUEQzv0Z4H397AHxduwXYpG+M7LHhMvWzefZjL/7GCA8VlvHrLcK47o6vbJZ1QVEQYs28cyjVDE3h2SRb3vZlOWUWV22UZY2j4lQtMAFidc4Afv7SaiLAQ5v3kTJJ9ZNdaXUKCg/jDZQOIbdeSv3yyjW+LynjuutMJD/W9UZoxgcRGPOaElmfu4/oXVtGxdQve+mnzCZ1jRISpY5L5/aQBLNlawE9eXkNpeaXbZRkT0Cx4TK2+zN7PzXNW061DBG/cfibxUa3cLumU/XB4Ik9ePoBlGQXcOjfNwscYF1nwmBqt2eHZvRbXriWv3DKc6NYt3C6pwa4emsj0ywfyReY+bnt5jR3zMcYlFjzmezbmH+LGF1cTE9mCebee4Rehc8yVqQk8cdkAlm0r4OdvpVNVZWe7GeNtdnKB+Y68b4u56aXVRIaHMO/WM+jUJtztkhrd1UMT2V9UxvQFW4lqFcavLk6xS+wY40UWPOb/HSou58bZqyktr+TVn55FbB23MGjOfnpuT/YXlvHCF9uJadOCO0YmuV2SMQHDgscAnott3vZyGjv2FzHnx8N85moETUVEeHhCXwqOHGX6gq30iI5gfP8ubpdlTECwYzwGVeWhtzewavsB/njFIM7qGRj3tgkKEqZfMZDBCe245/V0NuYfcrskYwKCBY/hxeU5vLM2n7vHJjNxSJzb5XhVeGgws350OlGtQrl1bhp7j5S6XZIxfs+CJ8CtyNrH7z/czHkpnZg6OtntclwRExnO8zekcrC4nJ+8vIajFfYdH2OakgVPAMv7tpi75q2lW4dW/PmqQQS5eFsDt/WLbctfrx7E2p0H+f1/NrtdjjF+zYInQJWWV3L7K2sor6hi1o9SiQwPdbsk143v34Vbz+nOnJU7eD99l9vlGOO3LHgC1O8/3MzG/MP8+apB9OzY2u1yfMYD4/uQ2jWKaW+vJ3PvEbfLMcYvWfAEoAUbdzN35Q5+fHZ3zuvX2e1yfEpocBDP/PA0WoYG89NXvqK4rMLtkozxOxY8ASb3QDEPvLWegfFtmXZBH7fL8Umd24bzt2uGkFlQyG8++NrtcozxOxY8AaS8soqpr62lSuHvk4cQFmI//tqMSI7mJz/oyWurc1mw8Ru3yzHGr9gnTwD5+6eZrN15kD9cNoCuHSLcLsfn3TuuF/3j2vDQO+vZe9i+32NMY7HgCRDrcg/y7JJMLjstjosHxbpdTrMQFhLEU1cPoaS8kvvfWm9XsjamkVjwBIDS8krufWMdnSJb8KuL+7ldTrOSFNOahy9MYdm2AuaszHG7HGP8ggVPAHhywRayC4r445WDaNvSvq9zsq4bnsiYPjH84aMtZBUUul2OMc2eBY+fW5G5j9nLc7jxrG6cnRQYF/9sbCLCHy4fQMvQYB54az2VtsvNmAax4PFjhUcr+Plb6+kRHcGD4+3U6YaIiQznVxensGbHt8xZkeN2OcY0axY8fuxPC7ey61AJf7xyEC3Dgt0up9mbNCSO0X1imL5wCzv2F7ldjjHNlgWPn1qz41vmrMzhhjO7cXrXKLfL8Qsiwu8m9Sc0KIgH7Cw3Y06ZBY8fKquoYtrb6+nSJpz7z+/tdjl+pUvbljxyUV9WbT/Aq//d6XY5xjRLDQoeEWkvIp+ISIbzXOOf1iIyXkS2ikimiEyrq72IDBORdc4jXUQmVWuz1FnXsfkxzvREEVkiImtFZL2ITKjWJlFEPhaRzSLytYh0a0i/fd1zS7PI2FvI45P607qF3d28sV2VmsCIpGimf7TFvlhqzClo6IhnGrBYVZOBxc777xCRYOBZ4AIgBZgsIil1tN8IpKrqYGA8MFNEqn+CXquqg53HXmfaI8AbqjoEuAaYUW35ucAfVbUvMAzYi5/K2HOEZ5ZkcMmgWEb36eR2OX5JRPjtxP4crazit3bvHmNOWkOD51JgjvN6DjCxhmWGAZmqmq2qZcBrTrta26tqsaoeuyxwOFCfnekKtHFetwV2ATghF6KqnzjrLlTV4vp0rrmpqlIeemcDES1C+OXFKXU3MKese3QEd45M4oP0XSzbVuB2OcY0Kw0Nnk6quhvAeY6pYZk4ILfa+zxn2gnbi8hwEdkEbABurxZEALOd3WyPisix22Y+BlwnInnAh8AUZ3ov4KCIvOPshvujMwr7HhG5TUTSRCStoKD5fZi89VUeaTu+5RcT+hLduoXb5fi920f2oEd0BI/O30hpud0u25j6qjN4RGSRiGys4XFpXW2PraKGaXWOYFR1lar2A4YCD4lIuDPrWlUdAJzjPK53pk8GXlLVeGAC8LKIBAEhznL3O+vqAdxYyzZnqWqqqqZ27Nixnt3zDYeKy3nioy2kdo3iitPi3S4nILQICebxif3Zsb+YGUsy3S7HmGajzuBR1bGq2r+Gx3xgj4h0AXCeazp2kgckVHsfj7MbrD7tVXUzUAT0d97nO89HgHl4duUB3Ay84cxbiWcXXbSz/bXOrr4K4D3gtLr63dz86eOtHCwu4zeX9icoqKasN03hrKRoJg2J47nPssjca5fTMaY+Grqr7X3gBuf1DcD8GpZZDSSLSHcRCcNz4P/9E7V3lg1xXncFegM5IhIiItHO9FDgIjwnIgDsBMY48/riCZ4CZ/tRInJsCDMa8Ku7e23MP8Srq3bwozO7kRLbpu4GplH9YkJfWoYG89j7m1C17/YYU5eGBs8TwDgRyQDGOe8RkVgR+RDAGWXcBSwENuM582zTidoDI4B0EVkHvAvcoar7gBbAQhFZD6wD8oHnnTb3AbeKSDrwL+BG9ajEs5ttsYhswLPr71ibZq+qSnl0/kbaR4Rxz7hebpcTkDpGtuDecb34InMfH3+9x+1yjPF5Yn+h1Sw1NVXT0tLcLqNOb6zO5YG31/PnKwdx+el2bMctFZVVTHj6c0rKK/nknnMJD7VLFJnAJCJrVDX1RMvYlQuasUPF5TyxYAtDu0Vx2WlxdTcwTSYkOIhfXdyP3AMl/PPzbLfLMcanWfA0Y39bnMHB4jJ+fUl//ndWuXHL2UnRjO/XmWeXZLH7UInb5Rjjsyx4mqnsgkLmrszh6qEJdkKBD3n4wr5UqvLER1vcLsUYn2XB00z9/sMthIcGc+84uwioL0lo34qf/KAH89ftIi3ngNvlGOOTLHiaoeWZ+1i0eQ93jkqiY6RdocDX/HRkT7q0DeexDzbZrROMqYEFTzNTWaX89t9fEx/VkpvO7uZ2OaYGrcJCeHB8HzbmH2Z+er7b5Rjjcyx4mpk30nLZ8s0RHrqgr52y68MuGRRL/7g2/HHBVruOmzHHseBpRo6UlvPnj7cytFsUEwZ0drsccwJBQcIvJvRl16FSZi/PcbscY3yKBU8zMmNpFvsKy3jkwhQ7fboZOKtnNGP6xDBjSSb7C4+6XY4xPsOCp5nYdbCEF77YzqQhcQxKaOd2OaaeHprQh+LySv7+qV292phjLHiaiacWbQOF+86z67E1J0kxkVw9NIFXvtxBdoFdvdoYsOBpFrbtOcJba/K4/syuxEe1crscc5LuHptMi5Agpi/Y6nYpxvgEC55mYPqCrUSEhXDnqCS3SzGnICYynJ+c25MFm75htX2p1BgLHl+3OucAizbv4faRPWkfEeZ2OeYU3XJOd2IiWzB9wRa7Z48JeBY8Pkyda37FRLawL4s2c63CQpgyJpnVOd+ydGuB2+UY4yoLHh+2aPNe1uz4lrvH9qJVWIjb5ZgGujo1gcT2rZi+cKtdSscENAseH1VRWcX0BVvoER3BVal2gzd/EBYSxL3jerF592E+WL/L7XKMcY0Fj49656t8MvYW8vPzexMSbD8mf3HJoFj6dI7kL59so7yyyu1yjHGFfaL5oNLySv7yyTYGJbRjfH+7NI4/CQoSfn5+b3bsL+aNtFy3yzHGFRY8PuiVL3fwzeFSpo3vY5fG8UOj+8Rwetconl6cQUmZXUDUBB4LHh9TXFbBPz7L4uykDpzZs4Pb5ZgmICI8cH5v9hw+ypyVOW6XY4zXWfD4mJdX7mBfYRn3jLVL4/iz4T06cG6vjjy3NItDJeVul2OMV1nw+JCioxXMXJbNOcnRpHZr73Y5pon9/PzeHCop5/ll2W6XYoxXWfD4kDkrczhQVMY942y0Ewj6x7XlwgFdmL18O98WlbldjjFeY8HjI46UljNrWTYje3fktMQot8sxXnL32GSKyyuZ9bmNekzgsODxEXNW5HCwuNyO7QSY5E6RXDwwljkrcthnN4szAcKCxwccdkY7Y/vG2E3eAtDUMcmUllcyy471mABhweMDXvxiO4dLK7jbRjsBKSmmNRMHxzF3ZQ57j5S6XY4xTc6Cx2WHist54YvtnJfSif5xbd0ux7hkyphkyiuVfyy1UY/xfxY8Lnvhi2yO2Ggn4HWPjuCyIXG8smoHew7bqMf4twYFj4i0F5FPRCTDea7xdCwRGS8iW0UkU0Sm1dVeRIaJyDrnkS4ik6q1Weqs69j8GGd6oogsEZG1IrJeRCY400dVW3adiJSKyMSG9LuxHCwu48XlOVzQvzMpsW3cLse4bMroZKqqlBlLMt0uxZgm1dARzzRgsaomA4ud998hIsHAs8AFQAowWURS6mi/EUhV1cHAeGCmiFS/Ic21qjrYeex1pj0CvKGqQ4BrgBkAqrrk2LLAaKAY+LiB/W4Uz3+eTVGZjXaMR2KHVlyZGs+//pvLroMlbpdjTJNpaPBcCsxxXs8BJtawzDAgU1WzVbUMeM1pV2t7VS1W1QpnejhQn7tmKXBs2NAWqOmGJ1cAH6lqcT3W16QOFJUxe3kOEwZ0oXfnSLfLMT7izlFJKMqzNuoxfqyhwdNJVXcDOM8xNSwTB1S//nueM+2E7UVkuIhsAjYAt1cLIoDZzm6zR+V/l29+DLhORPKAD4EpNdRyDfCv2jojIreJSJqIpBUUNO3tiWcuy6KkvJK7xyQ36XZM8xIf1YqrhybwRlouuQdc//vImCZRZ/CIyCIR2VjD49K62h5bRQ3T6hzBqOoqVe0HDAUeEpFwZ9a1qjoAOMd5XO9Mnwy8pKrxwATgZRH5//6JSBdgALDwBNucpaqpqprasWPHenTt1OwrPMrcFTu4ZFAsyZ1stGO+685RSQhiox7jt+oMHlUdq6r9a3jMB/Y4H+jHPtj31rCKPCCh2vt4/rcbrM72qroZKAL6O+/znecjwDw8u/IAbgbecOatxLOLLrraqq4C3lVV1y8FPPOzLI5WVDLVRjumBl3atmTysATeXJPHzv026jH+p6G72t4HbnBe3wDMr2GZ1UCyiHQXkTA8u7veP1F7Z9kQ53VXoDeQIyIhIhLtTA8FLsJzIgLATmCMM68vnuCpvr9sMifYzeYte4+U8vKXO5g4OI6eHVu7XY7xUXeMSiI4SJix1EY9xv80NHieAMaJSAYwznmPiMSKyIcAzrGZu/Ds4tqM58yzTSdqD4wA0kVkHfAucIeq7gNaAAtFZD2wDsgHnnfa3AfcKiLpeALmRlVVp55ueEZdnzWwvw32j6XZlFcqU2y0Y06gU5twJg9N4K01eeR9a6Me41/E+Ww2x0lNTdW0tLRGXeeew6WcM30JlwyK5U9XDmrUdRv/s/tQCedOX8oVqfH8ftIAt8sxpl5EZI2qpp5oGbtygRc9tzSLyipl6mgb7Zi6dWnbkitT43kzzb7XY/yLBY+X7D5UwrxVO7ny9HgSO7RyuxzTTPx0ZE9U4R+fZbldijGNxoLHS55dkkmVKneOSnK7FNOMxEe14orT43ntv7l8c8iu4Wb8gwWPF+QfLOH11blcNTSBhPY22jEn546RSVSqMnOZjXqMf7Dg8YJnPs1EEBvtmFOS2KEVk4bEMW/VTrtfj/ELFjxNLPdAMW+m5XL10ATi2rV0uxzTTN05Konyyiqet7uUGj9gwdPEnlqUQVCQcMeonm6XYpqx7tERXDo4jle+3Mm+wqNul2NMg1jwNKGMPUd4d20ePzqjK13a2mjHNMydo5Iorajkn59vd7sUYxrEgqcJ/fnjbbQKC+EOO7ZjGkFSTGsuGhjL3JU5HCgqc7scY06ZBU8TWZ1zgAWbvuGWc7rTPiLM7XKMn5gyOoniskpe/MJGPab5suBpAodLy3nw7fXER7Xkth/0cLsc40d6dYpkwoDOvLQih0PFrl9o3ZhTYsHTyHbuL+bc6UvYub+YJy8fSKuwkLobGXMS7hqVTOHRCl5cbqMe0zxZ8DSy+KiWXDCgC/NuPYOzk6LrbmDMSUqJbcN5KZ14cfl2DpfaqMc0rqKjFTT1xaMteBpZUJDw+0kDGNa9vdulGD82dUwyR0ormLM8x+1SjJ958O31TH7+yyYNHwseY5qh/nFtGdMnhn9+sZ3CoxVul2P8RObeI/xnw26GJEYhIk22HQseY5qpqWOSOVRSzpwVOW6XYvzEs0uyCA8J5pYR3Zt0OxY8xjRTgxLaMbJ3R/75eTZFNuoxDbR9XxHz1+Vz3RmJdGjdokm3ZcFjTDM2ZXQy3xaX88qXO9wuxTRzM5ZkEhocxK1e+AqIBY8xzdjpXaM4JzmaWcuyKS6zUY85NbkHinlnbT6ThyUSExne5Nuz4DGmmfvZmGT2F5Xx6pc73S7FNFMzlmYRLMLt53rnYsYWPMY0c6nd2nN2UgdmLsumpKzS7XJMM5N/sIS31uRy1dB4Ordt+tEOWPAY4xemjk5mX+FR5v3XRj3m5Mz8LAtVvDbaAQseY/zC8B4dOKNHe/7xWRal5TbqMfXzzaFSXvtvLlecHk98VCuvbdeCxxg/MXVMMgVHjvKajXpMPc1clkWlKneM9O6tWyx4jPETZ/bowLBu7XnORj2mHvYeKWXeqp1MHBxHYgfvjXbAgscYvyEiTB2TzJ7DR3kzLdftcoyPe35ZNuWVVdw5ynvHdo6x4DHGj5yd1IHTu0YxY2kWRyts1GNqtq/wKC9/uYNLBsXSo2Nrr2/fgscYPyIi/GxMMrsPlfLWmjy3yzE+ataybMoqqrhrdLIr27fgMcbPnJMczeCEdsxYkkVZRZXb5Rgfs6/wKC+v9Ix2kmK8P9oBCx5j/I6I8LOxyeQfLOGdr2zUY77r+WXZlFZUujbagQYGj4i0F5FPRCTDeY6qZbnxIrJVRDJFZFpd7UVkmIiscx7pIjKpWpulzrqOzY9xpieKyBIRWSsi60VkQrU200Vkk4hsFpGnpSlvNGGMDxjZqyMD49vyzJJMyitt1GM89hceZa7Lox1o+IhnGrBYVZOBxc777xCRYOBZ4AIgBZgsIil1tN8IpKrqYGA8MFNEQqqt9lpVHew89jrTHgHeUNUhwDXADGf7ZwFnAwOB/sBQ4NwG9tsYn3bsWE/etyW8+1W+2+UYHzHrc89oZ8po735v53gNDZ5LgTnO6znAxBqWGQZkqmq2qpYBrzntam2vqsWqeuxSu+FAfe7BqkAb53VbYFe16eFAGNACCAX21GN9xjRro/vE0D+uDc8syaTCRj0Bb79zbOfigbEkxUS6WktDg6eTqu4GcJ5jalgmDqj+pYI8Z9oJ24vIcBHZBGwAbq8WRACznd1sj1bbbfYYcJ2I5AEfAlOc9a4ElgC7ncdCVd1cU2dE5DYRSRORtIKCgnr/Ixjji0SEqaOT2XmgmPfW7aq7gfFrz3++nZLySqaOcXe0A/UIHhFZJCIba3hcWlfbY6uoYVqdIxhVXaWq/fDsGntIRI5dNvVaVR0AnOM8rnemTwZeUtV4YALwsogEiUgS0BeIxxN4o0XkB7Vsc5aqpqpqaseOHevZPWN817iUTvTt0oZnbdQT0A4UlTF3ZQ4X+cBoB+oRPKo6VlX71/CYD+wRkS4AzvPeGlaRByRUex/P/3aD1dneGZ0U4Tk+g6rmO89HgHl4duUB3Ay84cxbiWf3WjQwCfhSVQtVtRD4CDijrn4b4w88x3qS2L6viA/W26gnUD3/ebZntOPysZ1jGrqr7X3gBuf1DcD8GpZZDSSLSHcRCcNz4P/9E7V3lg1xXncFegM5IhIiItHO9FDgIjwnIgDsBMY48/riCZ4CZ/q5TttQPCcW1LirzRh/dF5KZ/p0juTvn2ZSWVWfw6XGnxwoKmPuihwuHNCF5E7uj3ag4cHzBDBORDKAcc57RCRWRD4EcI7N3AUsxPOB/4aqbjpRe2AEkC4i64B3gTtUdR+ekwMWish6YB2QDzzvtLkPuFVE0oF/ATeqqgJvAVl4jhWlA+mq+kED+21MsxEUJEwZnUx2QRH/tlFPwPnn59kUl1cydYx739s5nng+m83xUlNTNS0tze0yjGkUVVXK+L8to0ph4d0/IDjIvsoWCL4tKmPEk58yqk8Mz/zwNK9sU0TWqGrqiZaxKxcYEwCOjXoy9xby0cbdbpdjvOSfX/jeaAcseIwJGBMGdKFnxwieXpxhx3oCwMHiMuas2MGEAV3o5SPHdo6x4DEmQAQHCT8b24ttewrtWE8AmLUsm6KyCqa6eE222ljwGBNALhrQhd6dInlqUYZ9r8ePFRw5yuzlOVw8MJbenX1rtAMWPMYElKAg4d7zerF9XxHvrLVruPmr55ZmUVZZxd1jfW+0AxY8xgSc81I6MTC+LX9blGH36/FDuw+V8MqqHVx+WpwrdxetDwseYwKMiHDfeb3JP1jC62m5dTcwzcozn2aiqkzxwWM7x1jwGBOAfpAczdBuUTzzaQal5ZVul2Mayc79xby+OpdrhiaS0L6V2+XUyoLHmAB0bNSz5/BRXvlyh9vlmEbyt8UZBAcJd/nINdlqY8FjTIA6o0cHRiRF89zSLIqOVtTdwPi0zL2FvLs2jx+d2ZVObcLrbuAiCx5jAth95/Vif1EZL63IcbsU00BPLdpGeGgwt5/b0+1S6mTBY0wAG5IYxdi+Mcz8LItDJeVul2NO0de7DvPv9bv58dnd6dC6hdvl1MmCx5gAd8+4XhwureCFz7PdLsWcor98so024SHc+oMebpdSLxY8xgS4frFtuXBAF174YjsHisrcLsecpHW5B1m0eQ+3/aAHbVuGul1OvVjwGGO4Z1wyJeWVPLc00+1SzElQVZ74aDMdIsK48ezubpdTbxY8xhiSYiK57LR45qzcQf7BErfLMfX02bYCvsw+wNQxybRuEeJ2OfVmwWOMATzHegD+/PFWlysx9VFVpTzx0RYS27di8rBEt8s5KRY8xhgA4tq15KazuvHu2ny+3nXY7XJMHean57PlmyPcf35vwkKa10d586rWGNOk7hiZRJvwUJ5csMXtUswJlJZX8qeF2+gf14aLBnRxu5yTZsFjjPl/bVuFcueonny2rYAVmfvcLsfU4pUvPcfipo3vS1CQuF3OSbPgMcZ8x4/O7EZcu5b84aMtVNktsn3O4dJynlmSyTnJ0YxIjna7nFNiwWOM+Y7w0GDuHdeLDfmH+PeG3W6XY44z87MsDhaX8+D4Pm6XcsoseIwx3zNxSBx9Okfyp4Vb7WZxPmT3oRJe+GI7lwyKpX9cW7fLOWUWPMaY7wkOEqZd0IedB4rttgk+ZPqCrVQpPDC+t9ulNIgFjzGmRuf26sjZSR14+tMMDhXbBUTdti73IO+uzeeWEd2Jj/Ldm7zVhwWPMaZGIsLDE1I4XFLOU4u3uV1OQFNVfvPBJqJbt+COUb59k7f6sOAxxtQqJbYNVw9N5OWVO8jcW+h2OQHrg/W7+WrnQX5+fq9mdWmc2ljwGGNO6L7zetEyNJjf/edrt0sJSKXllTz50RZSurThitMT3C6nUVjwGGNOKLp1C6aMSWLJ1gKWbt3rdjkB55+fZ5N/sIRHLupLcDP8smhNLHiMMXW68azudOvQisf/s5nySju92lv2Hi5lxtIszkvpxFk9m+eXRWvSoOARkfYi8omIZDjPUbUsN15EtopIpohMq6u9iAwTkXXOI11EJlVrs9RZ17H5Mc70RBFZIiJrRWS9iEyo1uZJEdnoPK5uSJ+NCURhIUE8fGEKmXsLedVOr/aaJxZsobyyil9M6Ot2KY2qoSOeacBiVU0GFjvvv0NEgoFngQuAFGCyiKTU0X4jkKqqg4HxwEwRqX5E7VpVHew8jo39HwHeUNUhwDXADGf7FwKnAYOB4cDPRaRNA/ttTMAZ2zeGEUnR/HVRht2p1Av+u/0A73yVzy3n9KBbdITb5TSqhgbPpcAc5/UcYGINywwDMlU1W1XLgNecdrW2V9ViVa1wpocD9blglALHAqUtsMt5nQJ8pqoVqloEpOMJM2PMSRARfnlxCkVHK3jyI7t6dVOqqKzil/M3Ets2nCmjm//p08draPB0UtXdAM5zTA3LxAG51d7nOdNO2F5EhovIJmADcHu1IAKY7exme1REjh1tewy4TkTygA+BKc70dOACEWklItHAKKDGU0NE5DYRSRORtIKCgnr+ExgTOHp1iuTmc7rzeloua3YccLscvzVn5Q62fHOEX16cQquw5n/69PHqDB4RWVTt+Ej1x6V1tT22ihqm1TmCUdVVqtoPGAo8JCLhzqxrVXUAcI7zuN6ZPhl4SVXjgQnAyyISpKof4wmiFcC/gJVA9RCrvs1ZqpqqqqkdO3asZ/eMCSxTRycT2zach9/dSIWdaNDo9hwu5a+fbOPcXh05v19nt8tpEnUGj6qOVdX+NTzmA3tEpAuA81zTuZZ5fHeEEc//doPV2V5VNwNFQH/nfb7zfASYh2dXHsDNwBvOvJV4dtFFO+9/5xwPGocnCDPq6rcxpmYRLUL45cUpbPnmCC+tyHG7HL/z+w83U1ZRxWOX9ON/O3T8S0N3tb0P3OC8vgGYX8Myq4FkEekuImF4Dvy/f6L2zrIhzuuuQG8gR0RCnN1liEgocBGeExEAdgJjnHl98QRPgYgEi0gHZ/pAYCDwcQP7bUxAO79fZ0b27shfP9nGN4dK3S7Hb6zM2s/8dbu4/dwedPezEwqqa2jwPAGME5EMYJzzHhGJFZEPAZxjM3cBC4HNeM4823Si9sAIIF1E1gHvAneo6j6gBbBQRNYD64B84HmnzX3ArSKSjmeX2o2qqkAo8LmIfA3MAq477niRMeYkiQi/vqQf5VXKb+2KBo3iaEUlj87fSHxUS3460v9OKKhOPJ/N5nipqamalpbmdhnG+LS/Lcrgr4u28dJNQxnZu6Zzi0x9/eXjrTz9aSazbxrKqGb8bykia1Q19UTL2JULjDGn7Cfn9qBnxwgefncjhUdtR8Kp2vLNYWYszWLi4NhmHTr1ZcFjjDll4aHBTL9iELsOldh3e05RZZXy4NsbaNMylF9e3M/tcrzCgscY0yCnd43iprO68/KXO/gye7/b5TQ7s5dvJz33IL+6OIX2EWFul+MVFjzGmAa7//xeJLZvxbS311NSVul2Oc1G7oFi/vzxNkb3ieGSQbFul+M1FjzGmAZrFRbCE5cPIGd/MX/5ZKvb5TQLVVXK/W+mExwkPD6xv99+Z6cmFjzGmEZxVs9ofjg8kRe+2M6aHd+6XY7Pe3H5dlZtP8AvL04htl1Lt8vxKgseY0yjeeiCPsS2a8k9r6+zs9xOYNueI0xfuJVxKZ248vR4t8vxOgseY0yjiQwP5amrB5P3bTGPvb+p7gYBqKyiirtfW0dkixD+cNmAgNrFdowFjzGmUaV2a8+do5J4a00eH27Y7XY5Pudvi7fx9e7D/OGyAUS3buF2Oa6w4DHGNLqpY5IZlNCOh97ZwO5DJW6X4zOWZ+5jxtIsrkqN5zw/vfJ0fVjwGGMaXWhwEE9dPZjyyirufT2dyiq7NNfeI6X87LV19OzYmscuCYwvitbGgscY0yS6R0fw2CX9WJm9n6cWbXO7HFdVVqlzwkU5z/7wNL+8udvJsOAxxjSZq1ITuCo1nr9/msmSLTXdriswPLskk+WZ+/n1Jf3o3TnS7XJcZ8FjjGlSv7m0P327tOHu19eRe6DY7XK87ouMfTy1aBsTB8dyVWpC3Q0CgAWPMaZJhYcG89y1p1FVpdw57ytKywPnkjo79hdx57yvSI6J5HeTAvPU6ZpY8Bhjmly36Aj+dNUg1ucd4hfvbCAQ7gNWeLSCW+emIQLP/yiViBaBfVynOgseY4xXnN+vM/eO68U7a/OZsTTL7XKaVFWVct8b68jcW8gzk08jsUMrt0vyKRbBxhivmTI6iayCQv64cCs9oiO4YEAXt0tqEk8u2MLCTXt45MK+jEiOdrscn2MjHmOM14gIT14+kNMS23HPG+tYu9P/LiY6e/l2Zi7L5vozunLziO5ul+OTLHiMMV4VHhrMzOtTiYkM56aXVrNtzxG3S2o0H27YzW/+/TXnpXTisUv62ckEtbDgMcZ4XcfIFrxy83DCgoO4/oVVfnGa9fLMfdz9+jqGJLTj6clDCA6y0KmNBY8xxhWJHVrx8s3DKS2v4roXVrH3cKnbJZ2yFVn7uHnOarp3iOCFG4YSHhrsdkk+zYLHGOOa3p0jmX3TUAqOHOWqmSvJP9j8Lij6ZfZ+bn4pjYSoVrx663CiIsLcLsnnWfAYY1x1WmIUL988nP2FZVz1j5Xs3N98drutyNzHj19aTVxUS+bdekbA3ubgZFnwGGNcd3rXKObdegZFZRVcOXNFszjh4IP0Xdww+7/ER7Vk3q3D6RhpoVNfFjzGGJ8wIL4tr912BlUKl89YwbJtBW6XVKvZy7cz9bW1DEmI4s2fnEVMZLjbJTUrFjzGGJ/Rp3Mb3rvzbOKiWnLTS6t5+csdbpf0HaXllUx7ez2//uBrxvXtxNybh9G2VajbZTU7FjzGGJ8S164lb/30LM7t1ZFH39vIva+vo+hohdtlsetgCVfPXMlrq3O5a1QSz113up29dorskjnGGJ/TukUIz/8olb9/msHTizNYl3uQpycPoX9cW1fqeT99F4++t5HKKmXm9adzfgDftrox2IjHGOOTgoOEu8f24tVbPCcdTHx2OdMXbPHqbRX2Hi7lznlfMfVfa+nRMYIPpoyw0GkEDQoeEWkvIp+ISIbzHFXLcuNFZKuIZIrItLrai8gwEVnnPNJFZFK1NkuddR2bH+NM7yoii0VkvbNMfLU2NzjbyBCRGxrSZ2OMd53ZswMLfvYDJg6JY8bSLM5/ahkfbdjdpLdWKKuoYuZnWYz601I+2bSHn5/fmzd/cibdoyOabJuBRBrywxOR6cABVX3CCZQoVX3wuGWCgW3AOCAPWA1MVtWva2svIq2AMlWtEJEuQDoQ67xfCtyvqmnHbedN4N+qOkdERgM3qer1ItIeSANSAQXWAKer6gmvTpiamqppaWknWsQY42UrMvfxy/c3kbm3kIHxbbl7bDIje8UQ1EiXpyktr+T11bnM/CyLXYdKGds3hkcuTKGbBU69icgaVU090TINPcZzKTDSeT0HWAo8eNwyw4BMVc12inrNafd1be1Vtfo3yMLxBEZdUoB7nNdLgPec1+cDn6jqAWf7nwDjgX/VY53GGB9yVlI0C352Du+uzeepRRn8+KU0unZoxfVndOXCgV3o0rblKa138+7DvLUmj/fW5rO/qIzUrlE8ecVAzknu2Mg9MNDw4OmkqrsBVHX3sd1ex4kDcqu9zwOG19VeRIYDLwJdgetVtfppLbNFpBJ4G3hcPcO2dOBy4G/AJCBSRDrUsv24mjojIrcBtwEkJibWo/vGGG8LCQ7iytQELh0cx4JN3zB3RQ6P/2czj/9nM4MS2nFOUjSDEtrRL7YNndqEf+9inRWVVew4UMzm3YdJy/mWz7YVsH1fEaHBwpg+nbjx7G6c0aODS70LDHUGj4gsAmo6mvZwPbdR0xi4zhGMqq4C+olIX2COiHykqqXAtaqaLyKReILnemAucD/wjIjcCCwD8oGKk9m+qs4CZoFnV1tdNRpj3BMWEsQlg2K5ZFAsmXsLWbjpGz7+eg/PfZZFZZXnv29IkBDdugVhIUGIwMHicg6VlP//OsJDgzizRwd+fHY3LhoYa9dZ85I6g0dVx9Y2T0T2iEgXZ7TSBdhbw2J5QEK19/HALud1ne1VdbOIFAH9gTRVzXemHxGReXh25c1V1V3AZU5drYHLVfWQiOTxv915x7a/tK5+G2Oaj6SY1iTFJHHnqCRKyirZtOsQW745wu5DJRQcOUpZRRWVCu1ahtI+IoyE9q3o2yWSpJjWtAix7+J4W0N3tb0P3AA84TzPr2GZ1UCyiHTHMwq5Bvjhido7y+Y6JxN0BXoDOSISArRT1X0iEgpcBCxy2kTjOVGhCngIz246gIXA76udcXeeM98Y44dahgWT2q09qd3au12KqUVDv8fzBDBORDLwnLX2BICIxIrIhwDOsZm78ATAZuANVd10ovbACCBdRNYB7wJ3qOo+oAWwUETWA+vwBNnzTpuRwFYR2QZ0An7nbP8A8Fs8Abga+M2xEw2MMcZ4X4NOp/Zndjq1McacvPqcTm1XLjDGGONVFjzGGGO8yoLHGGOMV1nwGGOM8SoLHmOMMV5lwWOMMcar7HTqWohIAdCQ++5GA/saqZzmItD6HGj9BetzoGhIn7uq6gmvrmrB00REJK2uc9n9TaD1OdD6C9bnQNHUfbZdbcYYY7zKgscYY4xXWfA0nVluF+CCQOtzoPUXrM+Bokn7bMd4jDHGeJWNeIwxxniVBY8xxhivsuCphYiMF5GtIpIpItNqmC8i8rQzf72InFZXWxFpLyKfiEiG8xxVbd5DzvJbReT8pu/h93mzzyIyTkTWiMgG53m0d3r5nf549WfszE8UkUIRub9pe1czF36vB4rIShHZ5Pysw5u+l9/rkzd/r0NFZI7T180i4spNJ5uoz1c6P8cqEUk9bn0n9/mlqvY47gEEA1lADyAMSAdSjltmAvARIMAZwKq62gLTgWnO62nAk87rFGe5FkB3p32wn/d5CBDrvO4P5Ptzf6ut823gTeD+APi9DgHWA4Oc9x0C4Pf6h8BrzutWQA7QzU/63BfP3aCXAqnV1nXSn1824qnZMCBTVbNVtQx4Dbj0uGUuBeaqx5dAOxHpUkfbS4E5zus5wMRq019T1aOquh3IdNbjTV7ts6quVdVdzvRNQLiItGiivtXE2z9jRGQikI2nv27wdp/PA9arajqAqu5X1com6lttvN1nBSJEJARoCZQBh5uma7Vqkj6r6mZV3VrD9k7688uCp2ZxQG6193nOtPosc6K2nVR1N4DzHHMS22tq3u5zdZcDa1X16ClXf/K82l8RiQAeBH7dSPWfCm//jHsBKiILReQrEXmgUXpxcrzd57eAImA3sBP4k6oeaHg3TkpT9bkh2/uOkDpWGKikhmnHn3de2zL1aXsq22tq3u6zZ4Ui/YAn8fx17E3e7u+vgb+qaqFITc29wtt9DgFGAEOBYmCxeG6LvLiuQhuRt/s8DKgEYoEo4HMRWaSq2XUV2oh8/vPLgqdmeUBCtffxwK56LhN2grZ7RKSLqu52hrV7T2J7Tc3bfUZE4oF3gR+palaj9KL+vN3f4cAVIjIdaAdUiUipqj7TGJ2pJzd+rz9T1X0AIvIhcBrgzeDxdp9/CCxQ1XJgr4gsB1Lx7GL1lqbqc0O2911NdYCrOT/wBHI2ngNlxw6w9TtumQv57sG5/9bVFvgj3z0gOd153Y/vHpzLxvsHYb3d53bOcpcHws/4uPU+hjsnF3j7ZxwFfIXnIHsIsAi40M/7/CAw21lXBPA1MNAf+lyt7VK+e3LBSX9+ef0/fHN54DnrYxueMzQedqbdDtzuvBbgWWf+huN+EN9r60zvgOevvQznuX21eQ87y28FLvD3PgOP4NkXvq7aI8Zf+3vcdh/DheBx6ff6OjwnU2ykhhD2tz4DrfGctbgJT+j83I/6PAnP6OYosAdYWG3eSX1+2SVzjDHGeJWd1WaMMcarLHiMMcZ4lQWPMcYYr7LgMcYY41UWPMYYY7zKgscYY4xXWfAYY4zxqv8DeBRhpsyrTHEAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(tran_analysis.time, tran_analysis.branches['vphantom_power'])" + ] + }, + { + "cell_type": "code", + "execution_count": 33, + "id": "e5444251", + "metadata": {}, + "outputs": [], + "source": [ + "vs.amplitude = 1@u_V" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "95fa7e42", + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3 (ipykernel)", + "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.4" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +}