{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# ARMマイコンでつくるダイレクト・サンプリングSDRを試す\n", "トラ技2021/01から2021/05で連載されていた「ARMマイコンでつくるダイレクト・サンプリングSDR」(以下ARM-SDRと略す)を実際につくって試しみます。\n", "\n", "記事を参考に作ったARM-SDRです。\n", "\n", "\n", "\n", "## トラ技2021/01のおさらい\n", "\n", "ブレッドボードの配線をトラ技2021/01の図3から引用します。\n", "\n", "\n", "\n", "### NUCLIEO-F446REのアプリ\n", "NUCLIEO-F446REのアプリは、mbedのサイトで公開されています。「SDR」で検索すると「SDR_AM_Rx_CIC_FIR」のプロジェクトが見つかりました。\n", "\n", "富山の放送局に合わせてmain.cppのF_C_配列を以下のように修正しました。\n", "```C++\n", "// 各放送局の搬送波の周波数,単位: Hz\n", "const float F_C_[] = { 648.0e3f, // NHK 第1\n", " 738.0e3f, // 北日本放送\n", " 954.0e3f, // TBS ラジオ\n", " 1035.0e3f, // NHK 第2 \n", " 1134.0e3f, // 文化放送\n", " 1242.0e3f, // ニッポン放送\n", " 1422.0e3f}; // ラジオ日本\n", "```\n", "\n", "### ARMでダイレクト・サンプリングができた理由\n", "コラム1にも書いてありますが、NUCLIEO-F446REに搭載されているSTM32F446REは、フラッシュ・メモリのプログラムをノーウェイトでアクセスできることが大きいです。他のCPUで試す場合には、CPUクロックの他にフラッシュ・メモリのアクセス時間にも注意必要です。\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## RFフロントエンドについて\n", "ここでの目的は、コラム2にあるAM受信機用アンプとバンドパスフィルタの総合周波数特性を求めることです。\n", "\n", "RFフロントエンドの回路をトラ技2021/01の図2から引用します(ハイパス・フィルター、バンドパス・フィルターの値は筆者が付記))。\n", "\n", "\n", "\n", "コラム2では、周波数特性を以下のように求めたそうです。\n", "- OPアンプの周波数1MHzのオープン・ループ・ゲインが30dB, ゲインの減衰率が20dB/dec\n", "- バンドパスフィルタの周波数特性は、4端子回路のF行列の縦続接続の関係を使用\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "OPアンプの増幅は、R1=1kΩ、R3=4.7kΩの非反転増幅とすると、\n", "$$\n", "G = \\frac{V_{out}}{V_{in}} = \\frac{R_1 + R_3}{R_1} = 5.7倍 \\simeq 15dB\n", "$$\n", "\n", "NJU77701のデータシートの利得/位相 対 周波数 特性例を引用します。\n", "この図から20dB/decの減衰率が求まります。\n", "\n", "" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 4端子行列(F行列)\n", "できるだけ覚えることを少なくして、特性を求めることを考えてみましょう。\n", "\n", "抵抗R、コンデンサーC、コイルLのインピーダンスは、$ Z_R = R, Z_C = \\frac{1}{sC}, Z_L = sL $となります。\n", "\n", "以下のような$ Z_1, Z_2 $の直列と並列のインピーダンスの回路について考えます(トラ技2013/8月号付録p26から引用)。\n", "\n", "\n", "\n", "この場合のF行列は、直列、並列のインピーダンスを表す2つのマトリックスの積で表すことができます。\n", "\n", "$$\n", "\\begin{bmatrix}\n", " V_1 \\\\\n", " i_1\n", "\\end{bmatrix} = \n", "\\underbrace{\n", " \\begin{bmatrix}\n", " 1 & Z_i \\\\\n", " 0 & 1\n", " \\end{bmatrix}\n", "}_{直列インピーダンス}\n", "\\underbrace{\n", " \\begin{bmatrix}\n", " 1 & 0 \\\\\n", " 1/Z_2 & 1\n", " \\end{bmatrix}\n", "}_{並列インピーダンス}\n", "\\begin{bmatrix}\n", " V_2 \\\\\n", " i_2\n", "\\end{bmatrix}\n", "$$\n", "\n", "出力端子を開放したとき、$i_2 = 0 $なので、これからゲインGは、以下のようになります。\n", "\n", "$$\n", "\\begin{bmatrix}\n", " V_1 \\\\\n", " i_1\n", "\\end{bmatrix} = \n", "\\begin{bmatrix}\n", " 1 + Z_1/Z_2 & Z_1 \\\\\n", " 1/Z_2 & 1\n", "\\end{bmatrix}\n", "\\begin{bmatrix}\n", " V_2 \\\\\n", " 0\n", "\\end{bmatrix}\n", "$$\n", "\n", "$$\n", " G = \\frac{V_2}{V_1} = \\frac{Z_2}{Z_1 + Z_2}\n", "$$" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 4端子の特性\n", "MATLAB互換のライブラリPython-Controlを使って4端子の伝達関数を求め、特性曲線を出力してみましょう。\n", "\n", "最初にPython-Controlをインストールします。\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "!pip -q install control" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "次に、$ Z_1, Z_2 $にバンドパス・フィルターの左のローパス・フィルター(LPF)の特性を計算します。\n", "\n", "$$\n", "Z_1 = R + sL\n", "$$\n", "$$\n", "Z_2 = \\frac{1}{sC}\n", "$$\n", "から$ G_L $は、次のようになります。\n", "$$\n", "G_L = \\frac{1}{LC s^2 + RC s + 1}\n", "$$" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", " 1\n", "-----------------------------\n", "7.26e-15 s^2 + 8.91e-08 s + 1\n", "\n" ] } ], "source": [ "from control import matlab\n", "\n", "R = 270.0\n", "L = 22.0e-6\n", "C = 330.0e-12\n", "# tfの引数のnumrator, denominatorは、sの次元の高い順に指定係数を配列で指定する\n", "GL = matlab.tf([1], [L*C, R*C, 1])\n", "print(GL)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd3xUVfr48c8zk15pIUAoIfQqEEQRqTZsa++6awNddXXXde37c91dv+q66xa7ImAF27oqilgggPQiIl0ggAGkSiC9zPP7YyYYQ8okU5M879drXpk599xzn3BITu49TVQVY4wxxhuOUAdgjDGm8bBGwxhjjNes0TDGGOM1azSMMcZ4zRoNY4wxXrNGwxhjjNciQh1AoLVp00bT09MDeo38/Hzi4+NDUo4359SVp7bjNR2rLt3btGAI9zrxJl+g6iVUdeLPa9e3nHCvk/rE6G81XXfFihX7VTXlmAOq2qRfmZmZGmhz5swJWTnenFNXntqO13SsunRv04Ih3OvEm3yBqpdQ1Yk/r13fcsK9Try5dqDUdF1guVbzO7XRPZ4SkfEislFENovIvaGOxxhjmpNG9XhKRJzAM8BpQA6wTEQ+VNV1oY2scVNVSstdlLuUcpdSdvSriwOFLr4/WHD0eGm5+9iWQ+XEbztIWbk7X5lLWbW3jKI1uylzqSddWZtTSs7i7Z5zXUfLLytXyl0uSj3X+vkxF+UuiItykhgTQWJMJIkxEbSMi6J9ixjSWsTSJiEap0NC/U9nTLPTqBoNYBiwWVW3AojIdOA8wO+Nxn++/I69R4pQBZcCKKq4XxXvcX/+4YdiPtyzyvNZcR09pnhO/emcas6vyLf/QBFTs5dWOuZe4kWrnO+qVO6Phwr519oFR3/hl7uUctWffc4vLCLiqy8od7mqzeNSYNbMmv8x5s6pPn3xomPTVq48Nm3NmhqLjnQKEQ4HEQ4hwik4HQ4inYJDhPySMo4UlVHuOnapmwiH0LFlLL3bJdGnfRK92ycyuFML2ibF1Px9GGN8JhW/mBoDEbkYGK+qN3o+XwOcoKq3Vck3EZgIkJqamjl9+vR6X+uRxYX8kO9CBEAQgYq/awU86e73qi4cDsfR9Mr58Hw+5r07zp+V5XKVE+FwHs1X3bV+Vr6AlpcTGRmBQ8Ap4Kh4AQ4HOBBc5aVER0X+lAdwOORonvLSEmKio46e6xQ5WlZZSTFxsdE4RIioOO6AkqIi4uNicXo+OwWKiwpJiIvF6XCf7xQoKiwgMSEOpwgOwV2GoyLWuu8UVJUSFxSWKkdK4UChi4NFysEi5Yd8F98fcbG3wNOIAqlxQq9WTjLiSxnWMZ64SN/uRvLy8khISAjIOXXlq+14TceqS6+a1pDvyV/8de36lhPudVKfGP2tpuuOHTt2haoOPeZAdR0d4foCLgEmVfp8DfBUbedYR3jT79zLLy7VFdsP6gtzN+v1U5Zq/4c+1S73zNBu932sV720WKcuyNY9uYUNKts6wv3LOsLrlxYM9e0Ib2yPp3KATpU+dwR2hSgWEybioiIY0rklQzq3ZOKobpS7lMkfzGZ/TAc+X7eHhz5cy8MfrWVkjxQuHJLG6X3bERvlDHXYxjRKja3RWAb0EJGuwE7gcuDK0IZkwo3TIfRo6WTCmD7cd2YfNu89wv++3sX7X+/kjumrSIiO4BeDOvDL4V3o3S4p1OEa06g0qkZDVctE5DZgFuAEJqvq2hCHZcJc97aJ3HVGL+48rSdLsg/yzorveXdFDm8u2cGwrq345fAunNGvHZHORjcC3Ziga1SNBoCqfgJ8Euo4TOPjcAjDu7VmeLfW/PHsvry9/HteX7Kd2978mraJ0fzqpHSuPqELyXGRoQ7VmLBlf1qZZqllfBQ3je5G1l1jmXztUHq1S+SJWRsZ/tiXPPzRWnJ+LAh1iMaEJZ/uNESkLTAC6AAUAmtw97i7/BCbMQHndAjjeqcyrncq63cf5qV5W3lt0XZeXbSdswe0JzO+PNQhGhNWGtRoiMhY4F6gFfA1sBeIAc4HuonIu8A/VPWwvwI1JtD6tE/iycsGcdcZvZiyIJtpS7/nw+IyPtu7mFvHdGd4t9aIF3NLjGnKGnqncRYwQVV3VD0gIhHAObiX+njPh9iMCYkOLWJ54Oy+/OaUHvz1zTlk7c7jyklLGNSpBbeN7c4pfdpa42GarQb1aajqH6prMDzHylT1f6pqDYZp1JJiIjkrI4p5d4/lr+f3Z39eMTe+upwz/z2fD7/ZVe3yJsY0dQ3uCBeR0SIy0PP+UhF5WkR+JyLR/gvPmNCLiXRy9YldyLprDE9eehxlLuX2aV9zyj+yeGvZDkrKrAvPNB8N7dN4BhgIRIvIJiAB+BQ4CZgMXOW3CI0JExFOBxcO6cj5g9L4bN0PPDNnC/e89y3/+uI7Jo7K4PLjO9tMc9PkNbRPY6yq9hWRGNwzs9uqarmIvACs9l94xoQfh0MY3789Z/Rrx7zv9vPMnM08/NE6np69metP7krXMntsZZquhjYaRQCqWiQi21W13PNZRaTUb9EZE8ZEhNE9UxjdM4Vl2w7y9OzNPDFrI7ERsM61ketGpNM6wZ7WmqaloY1GWxG5E/cq3RXv8Xw+dk9ZY5q449Nb8cr1w1izM5c/vb2IZ7I28/JX2VwxrDMTRnWlfXJsqEM0xi8a2mi8BCRW8x5gkk8RGdOI9U9L5rbBMXTsm8mzWVt4ZdE2Xlu8jQsHd+Sm0RlkpIRmHwtj/KVBjYaqPuzvQIxpSrq3TeTJSwfxu1N78tL8rby17HveXvE9Z/Zvxy1jutM/LTnUIRrTIA0dPfWf2o6r6u0NC8eYpqVTqzj+fF5/fjOuB1MWZPPaou188u0PjOzRhpNalDNa1SYKmkalofM0VnheMcAQ4DvPaxBgi/UYU0VKYjR3j+/NgvvGcff4XqzffZjHlxVx4XML+XzdHlw2UdA0Eg19PPUKgIhci3v4bann8/PAZ36LzpgmJikmklvGdOf6EV15ZNps5uwuZsKry+mVmsjNYzI4d2AHImxfDxPGfP3f2YGfd4IneNKMMbWIiXRySudIsu4awz8vOw5F+d1b3zDm71m8tmgbRaV2w27Ck6+bMD0GfC0iczyfRwN/8rFMY5qNCKeDCwZ35Lzj0pi9YS/PZm3mjx+s5d9fbmZMBxdDTiwlKcY2hTLhw6dGQ1WniMhM4ARP0r2q+oPvYRnTvDgcwql9UzmlT1uWZB/k2awtvLtpH7Menc01w7tw/cldaWMTBU0YaOjoqXRV3QbgaSQ+qHJcgDRVzfE5QmOaERHhxIzWnJjRmqkffMmyvJY8N3cLL3+VzWXHd2LCyAw6tYoLdZimGWvoncYTIuLA3VisAPbhHknVHRgLnAI8BFijYUwDpSc7ufa8IWzdl8cLc7cybekO3liyg18c14EhcbayrgmNho6eukRE+uJezfZ6oD1QAKwHPgEeUdWihpQtIk8A5wIlwBbgOlU9JCLpnvI3erIuVtWbG3INYxqTjJQEHr94IL89rQeT5mczbekO3i8p58v9S7lpVDdOzGhlcz1M0DS4T0NV1wEP+DGWCp8D96lqmYg8DtwH3OM5tkVVBwXgmsaEvfbJsfzxnL7cPq4Hf542h7k7c7nipcUc1zGZm0Z344x+7UIdomkGwm5AuKp+pqplno+LgY6hjMeYcJMcF8m53aL46p5x/N8FA8gtLOWWN1Yy7h9ZzN5RasN1TUCJavjORBWRj4C3VPV1z+OptcAm4DDwoKrOr+G8icBEgNTU1Mzp06cHNM68vDwSEnxfiK4h5XhzTl15ajte07Hq0r1NC4ZwrxNv8nlbLy5VVuwp55PsUrJzXSRGwWldIhnXKZKEKDmmnFDViT+vXd9ygl0ndaU3hp+VsWPHrlDVocccUNWgv4AvgDXVvM6rlOcB4H1+atiigdae95nA90BSXdfKzMzUQJszZ07IyvHmnLry1Ha8pmPVpXubFgzhXife5KtvvbhcLn3uvS/02slLtMs9M7TPH2fqwx+u1Xc/+bJB8QVCqOolVHVSU3pj+FkBlms1v1N9mqfhGVp7FZChqn8Wkc5AO1VdWtt5qnpqHeX+CjgHOMUTPKpaDBR73q8QkS1AT2C5L9+DMU2FiNC7lZObLxzGhh8O8+K8rby6aBsuVb46vIqJozLo0z4p1GGaRs7XPo1ngeHAFZ7PR4BnfClQRMbj7vj+haoWVEpPERGn530G0APY6su1jGmqerdL4slLBzH37rGc1jmCWWt/4Mx/z+dXk5ey/kB5xd28MfXma6Nxgqreyk/bv/4IRPlY5tO417P6XERWeRZBBBgFrBaRb4B3gZtV9aCP1zKmSUtrEcsVfaJZdO8p/OGMXqzdlcvjy4o4/5kFfPLtbsptdV1TT76uPVXq+etfwX03APg060hVu9eQ/h7wni9lG9NcJcdFcuvY7txwclcemz6buXvKuOWNlXRpHceEkRlcnNmRmEhnqMM0jYCvdxr/wd1Z3VZEHgG+Av7P56iMMQERE+lkTKdIvrhzNM9fPYQWcVE8+L81jHhsNk99+R2HCkpCHaIJc74uWPiGiKzAvWyIAOer6nq/RGaMCRinQxjfvz1n9GvHkuyDvDB3C//4fBPPzd3CZcd34oaTu9Kxpa1xZY7V0AULW1X6uBeYVvmY9TUY0zhUXiCxYsTVa4u28+qi7Zw9oD0TR2XYfubmZxp6p7ECdz+GAJ2BHz3vWwA7gK5+ic4YEzQVI67uOr0XUxdu480lO/jwm12c1K01E0ZlMKZniq1xZRrWp6GqXVU1A5gFnKuqbVS1Ne65Ff/1Z4DGmODq0CKW+8/qw8L7xnH/Wb3Zui+f66Ys44x/zeOd5d9TXGbLlDRnvnaEH6+qn1R8UNWZuHfvM8Y0ckkxkUwc1Y15d4/lyUuPwyHCH95dzcjH5/Bc1hZyC0tDHaIJAV+H3O4XkQeB13E/rroaOOBzVMaYsBEV4eDCIR25YHAa87/bz4vztvL4pxt4evZ3XD6sM9eNSLdO82bE10bjCtybLb3v+TyPn2aHG2OaEBFhVM8URvVMYe2uXCbNz+aVhduYunAb5wxsz4SR1mneHPg65PYgcIefYjHGNBL9OiTzz8sG8YczejH5K/fGUB+s2sWI7q2ZMDKD0dZp3mT5umDhHDyzwStT1XG+lGuMaRw6tIjlwXP68ptTejB96Q4mL8jm2inL6JWayIRRGfziuA5ERYTdtj3GB74+nrqr0vsY4CKgrIa8xpgmKjk2kptGd+O6EV356JtdvDR/K3e98w1PzNrAdSO6csWwziTHRoY6TOMHvj6eWlElaYGIzPWlTGNM4xUV4eCizI5cOCSNed/t56V5W3ls5gae+tLdad7X6dPSdCYM+Pp4qvLMcAfuzZFso2JjmjkRYXTPFEb3TGHNzlxemr+VqQu3oarMz/2aG63TvNHy9fFU5ZnhZUA2cIOvQRljmo7+acn8+/LB3D2+N3+eNo/P1+3hf55O84mjujGqRxvrNG9EfG00+qhqUeUEEYn2sUxjTBNUsbfH364dwbSlO5iyIJtfTV5K73aJTBiZwbnWad4o+FpDC6tJW+RjmcaYJiw5NpKbR3dj/t3j+Pslx6EKv3/nG0b+bTbPz93C4SKbaR7OGrrKbTsgDYgVkcG4H08BJAE2NdQYU6eoCAcXZ3bkoiFpzN20jxc9neZPz97M5cd34vqTu9KhRWyowzRVNPTx1BnAtUBH4MlK6UeA+32MyRjTjIgIY3q1ZUyvtqzZmcuL87YyxTPT/NzjOnDjyK7062Cd5uGiQY2Gqr4CvCIiF3m2YTXGGJ/1T0vmP1cM5u7xvZiyYBvTl+7g/a93cnL3NkwclcFI6zQPuYY+nrpaVV8H0kXkzqrHVfXJak7ztuw/AROAfZ6k+ytW0hWR+3CPzioHblfVWQ29jjEmfHVsGccfz+nL7eN68Kan0/yXnk7ziaMySHIdsxCFCZKGPp6K93xN8FcgVfxTVf9eOUFE+gKXA/2ADsAXItJTVW1xf2OaqOS4SH49phvXn5zOh6vcM83vfPsbWkYLNzu3cMUJnUmKsZnmwdTQx1MveL4+7N9wanUeMF1Vi4FsEdkMDMNGaxnT5EVHOLlkaCcuzuxI1qZ9/O2DFTw6cwNPzd7MFcM6cd0I6zQPFl9nhKfgfpSUXrksVb3et7C4TUR+CSwHfq+qP+IerbW4Up4cT5oxppkQEcb2aosMi6V198G8OH8rkxdsY8oC6zQPFlFt+LNBEVkIzMc9M/zoY6K6OsdF5AuqX27kAdwNw37cM83/ArRX1etF5BlgkacvBRF5GfikumuJyERgIkBqamrm9OnTG/DdeS8vL4+EBN+f1DWkHG/OqStPbcdrOlZdurdpwRDudeJNvkDVS6jqxJ/XrlzOvgIXn20vZV5OGcXl0K+1g/HpkfRv4zzaaR7udVKfGP2tpuuOHTt2haoOPeaAqjb4Bazy5Xwvyk8H1nje3wfcV+nYLGB4XWVkZmZqoM2ZMydk5XhzTl15ajte07Hq0r1NC4ZwrxNv8gWqXkJVJ/68dnXlHMov0WfmfKfH//Vz7XLPDD3jn3P1neXfa3FpedjXiTfXDpSargss12p+p/o6I3yGiJzlYxk/IyLtK328AFjjef8hcLmIRItIV6AHsNSf1zbGNF7JcZHcMqY78+8ZyxMXD8Slyl2emeYfby2xPc39xNe1p+4A7heRYqAU98xwVdUkH8r8m4gMwv14ahtwE+5C14rI28A63Isj3qo2csoYU0XlTvO5m/bx0vytvLPpAJ88+qXtae4Hvu6nkeivQCqVeU0txx4BHvH3NY0xTU/lmeavfPglXxe0YqpnpvnZA9x7mg/oaJ3m9eXr6Kkh1STnAttV1XbwM8aEhS5JTn71C/fy7FMWZDNt6fd8+M0uhme0ZuIo957mDofNNPeGr4+nngWGAN96Pg8AvgFai8jNqvqZj+UbY4zfdGgRywNnu/c0f2vp90xekM11U5fRo20CE0Zm0NJmmtfJ147wbcBgVc1U1UxgEO6O61OBv/lYtjHGBERSTCQTRmUw7+6x/OuyQUQ4Hdz93mp+n1XIM3M2c6igJNQhhi1fG43eqrq24oOqrsPdiGz1sVxjjAm4SKeD8wen8cntJ/PGjSfQJcnBE7M2MvzR2fzpw7XsOFAQ6hDDjq+PpzaKyHNAxey5y4BNnt37bHybMaZREBFGdG/D74fG0K73ECbNz+aNJdt5ddE2zuzfnhtHdmVw55ahDjMs+NpoXAvcAvwW93Dbr4C7cDcYY30s2xhjgq53uyT+fslx/OGMXkxduI3XF2/n4293Myy9FcNblTHKpc2609zXIbeFwD88r6ryfCnbGGNCKTUphnvG9+bWsd15e9n3vPxVNv/eVsxHO+Zy48gMLhySRkykM9RhBp1PfRoi0kNE3hWRdSKyteLlr+CMMSbUEqIjuP7krsz9wxh+fVw08dER3P/+t4x4bDb//uI7DuY3r05zXx9PTQEeAv6J+3HUdfy0X7gxxjQZEU4HJ7SP4O7LR7Ak+yAvzdvKP7/YxLNZm7lkaEf6R7pCHWJQ+NpoxKrqlyIiqrod+JOIzMfdkBhjTJMjIpyY0ZoTM1qzee8RJs3P5u1lObxR7mLOweVMHJVBZpdWoQ4zYHwdclskIg7gOxG5TUQuANr6IS5jjAl73dsm8thFA/nq3rGc0y2SJdkHuei5RVz47AI+XbOb8iY4WdDXRuO3QBxwO5AJXAP8ytegjDGmMWmbGMNFPaJYeO84/nxeP/bnlXDz6ysZ948svtxRSmFJ01lb1dfRU8s8b/Nw92cYY0yzFRcVwS+Hp3PVCV34bO0PvDBvK6+tK2DGY19yzYlduGZ4OimJ0aEO0ycNajRE5MPajqvqLxoWjjHGNH5Oh3DmgPaM79+OSf+bzbIjyTw1ZzPPz9vKRUPSGBDVeDvNG3qnMRz4HpgGLMFGTBljzDFEhB4tnUy4YChb9+Xx8lfZvLsih2llLmYfWMaEkRkM69rq6La0jUFD+zTaAfcD/YF/A6cB+1V1rqrO9VdwxhjTVGSkJPDIBQNYeO84zu8eycodh7jsxcWc/8wCZqze1Wg6zRvUaKhquap+qqq/Ak4ENgNZIvIbv0ZnjDFNTOuEaM7v7u40f+SC/hwuKuO2N7/mnvmFTFmQTX5xeG9F1ODRU569ui8EXgduBf4D/NdfgRljTFMWE+nkqhO68OWdo3nxmkxaxQgPf7SO4Y9+yd8+3cDew0WhDrFaDe0IfwX3o6mZwMOqusavURljTDPhcAin92tH1L5YkjKOY9L8rTw/dwsvzd/K+YPSmDAqg56pft9Zu8Ea2hF+DZAP9ARur9SJI4CqapIfYjPGmGZlSOeWPHtVJtsP5DP5q2zeXp7DOytyGNMrhYkjMxjerXXIO80b1Gioqq+TAmskIm8BvTwfWwCHVHWQiKQD64GNnmOLVfXmQMVhjDGh0qV1PA+f15/fntqTN5ZsZ+rC7Vw5aQn9OiQxcVQGZw1oT6QzYL+Ga+Xr2lN+p6qXVbwXkX8AuZUOb1HVQcGPyhhjgq9lfBS3jevBjSMz+GDVTl6an80d01fx+MwNXH9yVy47vhOJMZFBjSk0TZUXxH0PdinuuSDGGNNsxUQ6uez4znz221FMvnYonVvH8deP13PSo7N59JP17M4tDFosohqeY4NFZBTwpKoO9XxOB9YCm4DDwIOqOr+GcycCEwFSU1Mzp0+fXl02v8nLyyMhISEk5XhzTl15ajte07Hq0r1NC4ZwrxNv8gWqXkJVJ/68dn3LCfc6qU+MFbJzy/k0u5Rle8oR4IT2EYxPj6BzUv02hqrpumPHjl1R8fv3Z1Q16C/gC2BNNa/zKuV5Dvh9pc/RQGvP+0zcM9KT6rpWZmamBtqcOXNCVo4359SVp7bjNR2rLt3btGAI9zrxJl+g6iVUdeLPa9e3nHCvE2+uXZMdB/L14Q/Xap8/ztQu98zQqyct1rkb96rL5fLq/JquCyzXan6nhqRPQ1VPre24iEQAF+JuHCrOKQaKPe9XiMgW3KO3lgcwVGOMCWudWsXx/87tyx2n9ODNpTuYsiCbX05eSu92ibx6/TDaJsX49Xph1xHucSqwQVVzKhJEJAU4qKrlIpIB9ABsa1ljjAGS4yL59Zhu3HByVz78ZhezN+wJyIq64dpoXM6xHeCjgD+LSBlQDtysqgeDHpkxxoSxqAgHF2d25OLMjgEpPywbDVW9tpq094D3gh+NMcaYCmE75NYYY0z4Cdsht/4iIvuA7QG+TDI/n4QYzHK8OaeuPLUdr+lYdenVpbUB9tcRXyCEe514ky9Q9RKqOqkulmCVE+51AuH3s9JFVVOOSa1uSJW96j2E+MVQlePNOXXlqe14TceqS68hrdphe829TkJZL6Gqk1DWS7jXSSjrpb7/lvZ4yj8+CmE53pxTV57ajtd0rLp0f/07+EO414k3+axe/FeO1UnN6hVLk388ZUJLRJZrdbNKTchYnYSnxlIvdqdhAu3FUAdgjmF1Ep4aRb3YnYYxxhiv2Z2GMcYYr1mjYYwxxmvWaBhjjPGaNRrGGGO8FpZrT5mmSUQcwF+AJNwTmV4JcUgGEJGRwFW4fx/0VdWTQhxSsycinYGncc8Q36Sqj4U4pKPsTsP4REQmi8heEVlTJX28iGwUkc0icq8n+TwgDSgFcqqWZfynPvWiqvNV9WZgBmANeYDU82elJ/Cxql4P9A16sLWwRsP4aiowvnKCiDiBZ4Azcf+Hv0JE+gK9gEWqeifw6yDH2dxMxft6qXAlx25JYPxnKt7XydfA5SIyG5gT5DhrZY2G8YmqzgOq7msyDNisqltVtQSYjvsuIwf40ZOnPHhRNj/1rJeKxyG5qno4uJE2H/Wsk+uAh1R1HHB2cCOtnTUaJhDScO/hXiHHk/Zf4AwReQqYF4rAmrma6gXgBmBK0CMyNdXJp8DtIvI8sC0EcdXIOsJNIEg1aaqqBbh/OZnQqLZeAFT1oSDHYtxq+llZA1wc7GC8YXcaJhBygE6VPncEdoUoFvMTq5fw0+jqxBoNEwjLgB4i0lVEonDv+f5hiGMyVi/hqNHViTUaxiciMg1YBPQSkRwRuUFVy4DbgFnAeuBtVV0byjibG6uX8NNU6qTRrXIrIuOBfwNOYFI4TXoxxpimrlE1Gp4xzZuA03A/C1wGXKGq62o6p02bNpqenh7QuPLz84mPjw9JOd6cU1ee2o7XdKy6dG/TgiHc68SbfIGql1DViT+vXd9ywr1O6hOjv9V03RUrVuzXxr5HODAcmFXp833AfbWdk5mZqYE2Z86ckJXjzTl15anteE3Hqkv3Ni0Ywr1OvMkXqHoJVZ3489r1LSfc68SbawdKTdelhj3LG9udxsXAeFW90fP5GuAEVb2tSr6JwESA1NTUzOnTpwc0rry8PBISEkJSjjfn1JWntuM1Hasu3du0YAj3OvEmX6DqJVR14s9r17eccK+T+sTobzVdd+zYsSu0uu1nq2tJwvUFXIK7H6Pi8zXAU7Wd09A7jT2HCzW3sERdLledecP9r1r76ym45YT7X7V2p9HwfM3pZ4Ua7jQa2+S+oI1pvmbSUjbuOYLTISTHRtIiNpLkOPfXFnFRJMdGutPjItm9qwzdsJdW8VG0SYymTUIU0RHOQIRljDEh1dgajaNjmoGduMc0XxmIC/3mlO78kFvEoYJSDhWWcKiglNzCUvbnlbB5Xx6HCko5UlR2NP+Lq5f97PykmAjaJEaTkhB99GuKp0FpmxhDhxaxtG8RQ1JMZCDCN8aYgGhUjYaqlolIxZhmJzBZAzSm+ZyBHerMU1bu4khRGbOyvqLngMEcyCthf14x+48Usy+v2PO+hPW7DjPvSDFHisuOKSMxOoL2LWJonxyLFBTzbfl3dGgRS1rLWNJbx9M2MRqHo7qVBowxJvgaVaMBoKqfAJ+EOg6ACKeDlvFRtIt3MKRzyzrzF5WWsz+vmD2Hi9h1qIjduYXsOlTErkOF7MotZPu+MrJyNv3snJhIB+mt40lvHU+XNnF0bR1Pept4MtrEk5IYHahvzRhjqtXoGo3GLCbSSceWcXRsGUdml2OPZ2VlcYvlzYkAAB5jSURBVOKIkezOLSLnxwK2HShg2/58th/I57u9R5i9YS8l5a6j+VvERZIaXc7s3DX0TE2kd7tEeqQmkhxrj7yMMYFhjUaYiYl00rVNPF3bxDOyx8+PlbuUXYcK2XYgny1789i4J49lm3L478qd5FV69JXWIpaBHZMZ2LEFx3VKprCs8QyrNsaEt0Y1T6M+RORc4Ny0tLQJr7/+ekCvFeo5AfHx8RwoUnKOuMjJc7HjsIvsXBf7Ct11Kyjt4x10TXbSs5WDPq2cpMQKIlLndZva2PNglBPucwJsnkbD8zWnn5UmMU+jIa/mPCP8QF6xztmwR3/70iy9fspSHfLnz7TLPTO0yz0z9IRHvtA7pq3UaUu26/QZX9a7/MY69jwY5YT7nACbp9HwfDZPo/HN0zD10Co+ijG92sLuKMaMOR5VZcu+fBZvPcDirQdYsOUA/1vlnuby3Lo5jOudyrjebRnWtRVREbYAsjHmWNZoNCMiQve2CXRvm8DVJ3ZBVdm6P58pnyxipyue15dsZ/KCbOKjnJzcow1nDWhPtPWHGGMqsUajGRMRuqUkcGqXSMaMGUZhSTkLt+znyw17mb1+L7PW7iHKAaftXcm5AzswplcKMZE2092Y5swaDXNUbJSTU/qkckqfVFznKcu3/8jzM5exeMsBPl69m6SYCC4c0pErhnUOdajGmBCxRsNUy+EQhnVtRUHfaE4eOYqFWw7w7ooc3lyyg6kLt9G9hYMDiTmcc1x7W2fLmGbEhtz6QbgP7/TnMMIjJcpXO8uYs6OYvYVCi2jh9C4RjOkUias4P+yHEQajnHAf3mlDbhuez4bcYkNu/SHch3cGYhjh7Nmzde7GvXrlS4u0yz0ztP//+1R//fwsPZBXXO/4AiHc68SbfDbk1n/lhHudeHPtQKnvkFsbV2kaREQY1TOFN248kY9uO5nRvVKYmV3K6L/N4akvvyO/msUZjTGNn/VpGJ8N6JjM01cOYXjSbOb9mMQ/Pt/EK4u2ccepPUlroo8/jWmurNEwfpOW4OCFc4aycseP/O3TDfzxf2vokuSgZbcfGezFKsDGmPBX6+MpEekoIneJyAciskxE5onIsyJytogE9NGW57oqIm08n9NFpFBEVnlezwfy+qbhhnRuybQJJ/L0lYPJLVYufG4h97//LYeLSkMdmjHGRzXeaYjIFCANmAE8DuwFYoCewHjgARG5V1Xn+TsoEekEnAbsqHJoi6oO8vf1jP+JCOcM7IBz70ZWFKUyZeE25m7cxxOXDOSkbm1CHZ4xpoFqezz1D1VdU036GuC/IhIFBGqW1z+Bu4EPAlS+CZLYCOHBc/py9sD2/P7tb7jypSVcNyKde8b3ttnlxjRCYTdPQ0R+AZyiqneIyDZgqKruF5F0YC2wCTgMPKiq82soYyIwESA1NTVz+vTpAY053OcEhMvY8+Jy5Z2NJXyxo4z0JAe3DoomJS4wTznDvU68yWfzNPxXTrjXSX1i9De/z9MAvgVWV3nNx3030Lqu82so8wvcdyxVX+cBS4BkT75tQBvP++iK6wGZwPdAUl3Xsnka4Tf2fNaa3dr/oU914J9m6ez1e2qNraHCvU68yWfzNPxXTrjXiTfXDpRALI0+EygH3vR8vtzz9TAwFTjXizKqNlSnVpcuIgOArsA3ng2COgIrRWSYqv4AFHvOXyEiW3D3ryyv7/VNaJ3erx0z2iVy8+sruW7qMu49szc3jco4uimUMSZ8edNojFDVEZU+fysiC1R1hIhc7c9gVPVboG3F5yqPp1KAg6paLiIZQA9gqz+vb4KnS+t43r/lJH7/zjc8NnMD2w8U8Jfz+hHhtPmmxoQzbxqNBBE5QVWXAIjIMKDiAVgwp/2OAv4sImW473xuVtWDQby+8bOYSCdPXT6Yzq3ieC5rCzsPFfLsVUNIiLbpQ8aEK29+Om8EJotIRUNxBLhRROKBRwMWGaCq6ZXevwe8F8jrmeBzOIR7xvemc6s4HvzfGn758hKmXj+MpJjIUIdmjKlGnc8CVHWZqg4ABgGDVXWgqi5V1XxVfTvwIZrm4IphnXnmysF8uzOXq15awo/5JaEOyRhTjTobDRFJFZGXgemqekhE+orIDUGIzTQz4/u358VrhrJxzxGueGkx+/OKQx2SMaaKOudpiMhMYArwgKoeJyIRwNeeu4+wZftpeJ8n3MaerztQzr9WFNEu3sE9w2KIj6z/qKpwrxNv8tk8Df+VE+51Up8Y/S0Q8zSWeb5+XSltVV3nhcvL5mk0zrHnWRv3avf7P9YLn12g+cWl9T4/3OvEm3w2T8N/5YR7nXhz7UAJxH4a+SLSGlAAETkRyG1oq2aMN0b3TOE/lw/m6x0/ctNrKyguKw91SMYYvOjTAO4EPgS6icgC4FXgNwGNyhjgzAHtefyigcz/bj93v7u64i7XGBNCdQ65VdWVIjIa6AUIsFFVbY1rExSXDO3E3iPFPDFrI51bxfH703uFOiRjmrXalka/sIZDPUUEVf1vgGIy5mduGdONHQcKeGr2Zjq1iuPSoZ1CHZIxzVZtdxoVa0q1BU4CZns+jwWyAGs0TFCICH+9oD+7cgu5/7/fktYilhHdbU8OY0Khxj4NVb1OVa/D3QHeV1UvUtWLgH5Bi84Yj0ing2euGkJGSjy3vLGS7QfyQx2SMc2SN/M01qhq/0qfHcDqymnhyOZpeJ+nMY0931vg4uFFhbSMFh48MZaYiOrncIR7nXiTz+Zp+K+ccK+T+sTob4GYp/E0MAu4FvgV7qXSn6rrvHB52TyNpjf2fN6mvdr13hl682vL1eVyVZsn3OvEm3w2T8N/5YR7nXhz7UDx+zwNVb0NeB44Dvf6Uy+qqg25NSEzskcK953Zh5lrfuDZrC2hDseYZqW20VPiaW1Q1feB92vLY0ww3TiyK2t25fL3zzbSp30i43qnhjokY5qF2u405ojIb0Skc+VEEYkSkXEi8grux1V+JSLHicgiEflWRD4SkaRKx+4Tkc0islFEzvD3tU3jISI8duFA+rZP4o5pq9i6Ly/UIRnTLNTWaIzHvdnRNBHZJSLrRCQb+A64Avinqk4NQEyTgHvVvSDi+8AfAESkL+6tZvt5YntWRJwBuL5pJGKjnLxwTSYRTuGm11aQVxzMPcGMaZ5qG3JbpKrPqnur1y7AKbj30+iiqhNUdVWAYuoFzPO8/xy4yPP+PNzLsxerajawGRgWoBhMI9GxZRzPXDmELfvyuOvtb2ypEWMCrM4ht8EmIguBx1X1AxG5E3hYVRNF5Glgsaq+7sn3MjBTVd+tpoyJwESA1NTUzOnTpwc05nAf3tkchhF+ml3K9I0lXNwjknO6RYV9nXiTz4bc+q+ccK+T+sTob34fchuIF/AFsKaa13lAb+AzYAXwEHDAc84zwNWVyngZuKiua9mQ2+YxjNDlculv3lyp6ffO0Dkb9oR9nXiTz4bc+q+ccK8Tb64dKPUdcuvNHuF+p6qn1pHldAAR6Qmc7UnLASovOtQR2OX/6ExjJCI8ftFAvtubx+3TvubBYbbHuDGB4M3S6IhIFxE51fM+VkQSAxWQiLT1fHUAD+KeIwLu5dkvF5FoEekK9ACWBioO0/jERjl54epMRIT/rCyioMQ6xo3xN2/2CJ8AvAu84EnqCPwvgDFdISKbgA247ySmAKjqWuBtYB3wKXCrqtrOPOZnOreO46krBrMzT20PDmMCwJs7jVuBEcBhAFX9DvfKtwGhqv9W1Z6e171a6adeVR9R1W6q2ktVZwYqBtO4jeqZwsU9I5mxejcvzd8a6nCMaVK86dMoVtUSEffCcCISgWfrV2PC1VldI8mPbs1jMzfQt30yJ/ewpdSN8Qdv7jTmisj9QKyInAa8A3wU2LCM8Y2I8MTFx9G9bQK3TVvJ9wcLQh2SMU2CN0ujO4AbcI9oEtwr3k7SMH9YbEuje5+nKY8935PvXkq9TayDB06MIdpZ/VLqdZUTiHNsnkbwygn3OqlPjP4W0HkaQCtgYH3OCfXL5mnY2PPZG/Zo+r0z9KZXl2tZefVLqXtTjr/PsXkawSsn3OvEm2sHit+XRheRLBFJEpFWwCpgiog86WvrZkywjO3VlgfP7suna3/grx+vC3U4xjRq3nSEJ6vqYRG5EZiiqg+JyOpAB2aMP91wcld2HSrk5a+ySWsRy40jM0IdkjGNkjcd4REi0h64FJgR4HiMCZgHzurDWQPa8deP1zNjtS0mYExDeNNo/Bl35/dmVV0mIhm4l0c3plFxOIQnLx3E8ekt+d1bq5izYW+oQzKm0fFmu9d3VHWgqt7i+bxVVS+q6zxjwlFMpJOXrz2e3u2SuOn1FXz13f5Qh2RMo+JNR3iMiNwqIs+KyOSKVzCCMyYQkmIiefX6YWS0iefGV5exZOuBUIdkTKPhzTyNd3CvA3Ul7kdVVwHrVfWOwIfXcDZPw/s8zWnseWWHi5VHlxZyoEi5Y3AM/docuxGkzdPwL5unUb+0YPD7PA3ga8/X1Z6vkcDsus4Ll5fN07Cx57XZc7hQz/jnXO1x/yc689vdDS6nIefYPI3glRPudeLNtQPF7/M0gFLP10Mi0h9IBtIb0KAZE3baJsbw1sTh9EtL4tY3V/LuipxQh2RMWPOm0XhRRFoCf8S9p8U64G8BjcqYIEqOi+T1G05geEZr7nrnG578bCMuV1ivkmNMyHgzemqSqv6oqnNVNUNV26rq83WdVxsRuURE1oqIS0SGVkofJiKrPK9vROSCSseyRGRjpeMBW57dND/x0RFMvvZ4Lhvaif/M3sxt01ZSWGLbtRhTVZ0zwkUkGrgI9yOpo/lV9c8+XHcNcCE/bexUOX2oqpZ5JhR+IyIfqWrFFmxXqepyH65rTI2iIhw8dtEAurdN4P9mrmfHwYX8spsr1GEZE1a8WUbkAyAXWAEU++Oiqroe3MtXV0mvvH51DLZvhwkyEWHCqAwyUuL53VureGhPGQmddnPWgPahDs2YsOBNo9FRVccHPBIPETkBmAx0Aa6pdJcB7sUSy4H3gL96eviN8btT+qTy8e0j+dULc7nljZVcc2IX7jurN3FR3vzIGNN0eTNP40XgKVX9tl4Fi3wBtKvm0AOq+oEnTxZwV3WPnESkD/AKMEpVi0QkTVV3ikgi7kbjdVV9tYZrTwQmAqSmpmZOnz69PqHXm83TqF9aMPjruocO5zFzVySztpWREitc3z+aPq2Pnc/RkGvbPI3glRPudVKfGP3Nb/M0gG+B1bhHS5UCGz2fv8UzZ8PXF5CFuw+jpuNzqjsOXAs87c01bJ6GjT33RzmLtuzXUX+brV3umaEPvv+tHi4s8fnaNk8jeOWEe514c+1Aqe88jdrutc9pWLvVcCLSFfhe3R3hXYBewDbPvuQtVHW/iER6Yvsi2PGZ5uvEjNbMvGMkf5+1iSkLs5m5Zjd3n9GbizM74nDUbzdAYxqz2obc7gEuAP4AjAd2qur2ipcvFxWRC0QkBxgOfCwiszyHTsY9YmoV8D5wi6ruB6KBWZ59PFYBO4GXfInBmPqKi4rg/53blw9uHUHnVnHc/d5qzn92Acu2HQx1aMYETW13Gq/gfiw1HzgT6Av4Zb0pVX0fd6NQNf014LVq0vOBTH9c2xhfDezYgvd+fRIfrNrFozPXc8nzixjZow13ntaTwZ1bhjo8YwKqtkajr6oOABCRl4GlwQnJmPAnIpw/OI0z+rXjtcXbeH7uVi54diHjerdleHI5o1WPGVJuTFNQW6NRseYUnj6GIIRjTOMSG+Vk4qhuXHlCF15ZuI2X5m9ldkEpM3Yu4IaRGZzZvx2RTm9W6zGmcahxyK1nPkR+xUcgFijwvFdVTQpKhA1kS6N7n6c5DSMMdDnFZcrs7HyydjvYU6C0ihFGdYzg5LQI2sT+vPGwIbfBK8eG3NbM70ujN/aXDbm1YYTBLmfOnDlaXu7SL9b9oFdPWqzp987Q9Htn6NWTFusHq3ZqYUmZV2XbkFv/lWNDbmvmzyG3xpgGcjiEU/qkckqfVHJ+LODdFTm8szyH26d9TXyUk3F9UukiZZxYWk5MZO2TBY0JJ9ZoGBNgHVvG8dtTe3L7uB4s2nqAGat3MWvtHj7KL2Hyus8Z27st43q1ZVTPFFISo0MdrjG1skbDmCBxOIQR3dswonsb/nKeixfen0OOoy2fr9vDx6t3AzAgLZnRPVMY0yuFMtvTw4QhazSMCYEIp4N+bZzcOmYAj5zfn3W7D5O1cS9ZG/fxbNZmnp6zmSgnHL91McPSWzOsaysGd25hj7JMyFmjYUyIORxC/7Rk+qclc9u4HuQWlLJwy37++9VqduaX8q8vN6EKkU6hX4dkjuuYTFReKWl7jpCRkoDTljExQWSNhjFhJjkukjMHtCf2wEbGjBlJbmEpK7YfZEn2Qb7ecYh3VuRQUFLOS9/OIz7KSb+0ZAakJdMrNZG8Q+UcX1xGfLT9aJvAqHNp9MbK5ml4n6c5jT0PRjmBnhPgUmXLvnz2lkazNdfFtlwX3x9xUVJpk8E2sUJagoOOCQ5aRZTQpVUsbeMdJEZCfn6+zdNoYL7m9LNi8zQCKNRzAnzN05zGngejnFDMCSgrd+m2/Xn6j+mf61NfbtLb3lyppz85V7vd97F2uWfG0Vf/hz7V0Y98ore+sUKf+HSDvr1shy7cvF+378/XotIyr+L2hc3TqF9aMNg8DWOaIadD6NI6niGpEYwZ0+NoekmZi3c/zSK1Wz+2HShg+4F8Vn6Xw7c7c5m55gfKq4zQSkmMpkNyDB1axNI+OZYOLWJonxxLalI0bRKiSUmMtkdfzZzVvjFNWFSEgw4JDsb0ST2alpW1nzFjxlBa7iLnx0J2HSpk56FCdh8qYtehQnblFrJpzxHmbtpHQUn5MWXGRjpJSYymTUKU56u7MWmdEE3LuEhaxkWRHBtJi7hIWsRFER/ltMUbmxBrNIxppiKdDrq2iadrm/hqj6squYWl7DpUxL68YvYfKf7Z131Hisnen8/S7IP8WFBabRnu6wjJsVG0iIvEUVrI69uX0yIuksSYCBKjI0iIiSA+OoKE6AgSYyJIiI48+r4iPSrCFn0MFyFpNETkEuBPQB9gmHr2CPfsyjcJGOKJ7VVVfdRzLBOYinvhxE+AOzzP3YwxASAitIiLokVcVJ15S8tdHMwv4VBBKYcKSjhUWEpuQSk/et4fKiglt7CE7J0F7DxUyLpduRwpLiOvuAxvfoqjIhwkRkcQG+UkNtJJWVEhz25cREykk9hIB7GRTmKjnJ7PzmM+b/2hjPL1e4h0OoiK8LycDqIjHD9LKyhVikrLiXI6bEfGGoTqTmMNcCHwQpX0S4BoVR0gInHAOhGZpqrbgOeAicBi3I3GeGBm8EI2xtQk0ukgNSmG1KSYWvNlZWUxZszIo59VlYKScvKKyzhSVEa+pyE5UuT+mldUSn5JOUeKyjhSVEphaTlFpeXk/FCIQyC3sJS9h8spLC2nsKT86PHS8mpaolXLvftmvvwUgAiHHG1MIp3uRqa8pIgWq+bhdAgRTgcRDnG/dwiHDxXx6rZlRz87HUKk08G+vcV8emD10fQIp4NdO0tYWrTh6GenQ9i+rYTNzq04Pec6xP1yOqj0XhBx92E5RRBPmtPhbuSdlfOIMDS9ld/n8YSk0VDV9UB1zzkViPfsCR4LlACHRaQ9kKSqizznvQqcjzUaxjRqIkJ8tPsxVGo9NltwNz7DazxeWu6iqNTTiJS4mLtwMQMHDaGk3EVpmYvichclZe5XacX7chfrNnxH565djx6rSC8td1Fc5iJn1w+0aBVHuUspc6nnqztfQZmy70ixJ91FWbk7T35BOd8d2Xv0c7lLKSktQ3dsPbZx27i+gf+S1dvwl/E4Hf5dRSCk8zREJAu4q8rjqdeAU4A44Heq+qKIDAUeU9VTPflGAveo6jk1lDsR910JqampmdOnTw/o9xHucwJs7Hlwywn3OQG2n0bD8/m7TlQVBcoVDh/JJyYuHlVwKbjQn96r+y/qo+8VXLjn7Pz0nmPy927lwFHHIISwmacBfIH7MVTV13mV8mQBQyt9HgG8AUQCbYGNQAZwPPBFpXwjgY+8icPmadjY82CXE+5zAmw/jYbna04/KwR7noZ67grq6UrgU1UtBfaKyAJgKDAf6FgpX0dgl+9RGmOMqY9wG8e2AxgnbvHAicAGVd0NHBGRE8XdEfJL4INQBmqMMc1RSPo0ROQC4CkgBTgErFLVM0QkAZgC9MW9F/kUVX3Cc85QfhpyOxP4jXoRvIjsA7YH4vuoJBnIDVE53pxTV57ajtd0rLr06tLaAPvriC8Qwr1OvMkXqHoJVZ1UF0uwygn3OoHw+1npoqopx6RW98zKXvXuv3kxVOV4c05deWo7XtOx6tJrSKv2uWhzr5NQ1kuo6iSU9RLudRLKeqnvv2W4PZ5qrD4KYTnenFNXntqO13SsunR//Tv4Q7jXiTf5rF78V47VSc3qFUuTXRrdhAcRWa7VDdszIWN1Ep4aS73YnYYJtBdDHYA5htVJeGoU9WJ3GsYYY7xmdxrGGGO8Zo2GMcYYr1mjYYwxxmu2CZMJGhFxAH8BknCPSX8lxCEZji4AehXu3wd9VfWkEIfU7IlIZ+Bp3JP9NqnqYyEO6Si70zA+EZHJIrJXRNZUSR8vIhtFZLOI3OtJPg9IA0qBnGDH2pzUp15Udb6q3gzMAKwhD5B6/qz0BD5W1etxr5ARNqzRML6aintDrKNExAk8A5yJ+z/8FSLSF+gFLFLVO4FfBznO5mYq3tdLhSuBacEKsBmaivd18jVwuYjMBuYEOc5aWaNhfKKq84CDVZKHAZtVdauqlgDTcd9l5AA/evKUBy/K5qee9VLxOCRXVQ8HN9Lmo551ch3wkKqOA84ObqS1s0bDBEIa8H2lzzmetP8CZ4jIU8C8UATWzNVULwA34F4s1ARXTXXyKXC7iDwPbAtBXDWyjnATCNVtFaaqWoD7l5MJjWrrBUBVHwpyLMatpp+VNcDFwQ7GG3anYQIhB+hU6bNtmhUerF7CT6OrE2s0TCAsA3qISFcRiQIuBz4McUzG6iUcNbo6sUbD+EREpgGLgF4ikiMiN6hqGXAbMAtYD7ytqmtDGWdzY/USfppKndiChcYYY7xmdxrGGGO8Zo2GMcYYr1mjYYwxxmvWaBhjjPGaNRrGGGO8Zo2GMcYYr1mjYZoNESkXkVWVXumhjslfRGSwiEzyvL9WRJ6ucjxLRIbWcv50EekR6DhN42drT5nmpFBVB1V3QEQE97wlV5Bj8pf7gb/6cP5zwN3ABP+EY5oqu9MwzZaIpIvIehF5FlgJdBKRP4jIMhFZLSIPV8r7gGejnC9EZJqI3OVJP/oXvIi0EZFtnvdOEXmiUlk3edLHeM55V0Q2iMgbngYLETleRBaKyDcislREEkVkvogMqhTHAhEZWOX7SAQGquo3XnzPv6h0p7VRRLI9h+YDp4qI/SFpamX/QUxzEisiqzzvs4Hf4d4Y6jpVvUVETgd64N7jQIAPRWQUkI97TaDBuH9mVgIr6rjWDbj3pzheRKKBBSLymefYYKAf7oXpFgAjRGQp8BZwmaouE5EkoBCYBFwL/FZEegLRqrq6yrWGAmuqpF0mIidX+twdQFU/xLO2kYi8Dcz1pLtEZDNwnBffm2nGrNEwzcnPHk95+jS2q+piT9LpntfXns8JuBuRROB9z9LuiIg3C8qdDgwUkYrlrZM9ZZUAS1U1x1PWKiAdyAV2q+oygIrNkETkHeCPIvIH4Hrcu79V1R7YVyXtLVW9rdL3mlX5oIjcjfvf45lKyXuBDlijYWphjYZp7vIrvRfgUVV9oXIGEfktnn0nqlHGT495Y6qU9RtVnVWlrDFAcaWkctw/h1LdNVS1QEQ+x72b26W47yqqKqxy7VqJyCnAJcCoKodiPGUZUyPr0zDmJ7OA60UkAUBE0kSkLe5dBi8QkVhP/8G5lc7ZBmR63l9cpaxfi0ikp6yeIhJfy7U3AB1E5HhP/sRK/QuTgP8Ay1S16nah4F4dtbs336CIdAGeBS5V1aoNRE8grFdYNaFndxrGeKjqZyLSB1jk6ZvOA65W1ZUi8hawCtiOu9O4wt+Bt0XkGmB2pfRJuB87rfR0dO8Dzq/l2iUichnwlIjE4v6L/1QgT1VXiMhhatiOVVU3iEiyiCSq6pE6vs1rgdbA+57vcZeqniUiqbgfV+2u43zTzNnS6MbUk4j8Cfcv878H6XodgCygd01DgkXkd8ARVZ3UwGv8Djisqi83OFDTLNjjKWPCmIj8ElgCPFDHHJLn+HlfSX0dAl7x4XzTTNidhjHGGK/ZnYYxxhivWaNhjDHGa9ZoGGOM8Zo1GsYYY7xmjYYxxhivWaNhjDHGa/8fkV3nQbfPRRcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "from matplotlib import pyplot as plt\n", "\n", "matlab.bode(GL, Hz=True)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "同様にハイパス・フィルター(HPF)の特性を計算します。\n", "\n", "$$\n", "Z_1 = \\frac{1}{sC} + R \n", "$$\n", "$$\n", "Z_2 = sL\n", "$$\n", "から$ G_L $は、次のようになります。\n", "$$\n", "G_L = \\frac{LC s^2}{LC s^2 + RC s + 1}\n", "$$" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", " 1.551e-13 s^2\n", "-----------------------------\n", "1.551e-13 s^2 + 3.3e-07 s + 1\n", "\n" ] } ], "source": [ "R = 1000.0\n", "L = 470.0e-6\n", "C = 330.0e-12\n", "# tfの引数のnumrator, denominatorは、sの次元の高い順に指定係数を配列で指定する\n", "GH = matlab.tf([L*C, 0, 0], [L*C, R*C, 1])\n", "print(GH)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deXxU1fn48c8zM9kXdsK+KLgAghhE6wquuIKirVZtXVq/trV2+WrVqrVq/Vp/Vlut2mrdaxWtuyiiKBFElEXZ950AEpYA2SeZeX5/zAWGMEkmsyd53q/XNXPPOffcJzmSk3vvOeeKqmKMMcaEw5XsAIwxxrQc1mkYY4wJm3UaxhhjwmadhjHGmLBZp2GMMSZs1mkYY4wJmyfZAcRb586dtV+/fnE9R0VFBTk5OUmpJ5xjmirTWH5DeaHSw01LhFRvk3DKxatdktUmsTx3c+tJ9TZpToyxFuq8nTt3ZvLkyZNVdcxBB6hqq94KCws13qZOnZq0esI5pqkyjeU3lBcqPdy0REj1NgmnXLzaJVltEstzN7eeVG+TcM4dLw2dF5ijIX6n2u0pY4wxYWv1t6eMSRafX1m7vZzV2ypYt72CHRVedlV6Ka+pQxW2bavmteK5ZKW7ycvwkJeZRl6mhw7Z6XTOS2f9bh+H7aqic24G6R77+86khhbXaYjIGOBRwA08o6p/TnJIxuyzcWclk9bW8sLaWcxZV0p5Td2+vMw0F+2y0sjLTMMlUFHhZ7eWU+n1UVZdS3lNHf56q/rcM/MzANplpdE5N50ueRl0zs2gelcNC3wr6ZSbTqecdDrmZNApN51yr+L3Ky6XJPLbNm1Ii+o0RMQNPAGcCRQDs0XkPVVdktzITFtW5fXx9rebmDB7AwuKdwMwsGsVFx7dg2P6dGBg11z6d8khPzPtgOOKiooYNerUffuqSoXXR2mFl+3lNRR9NZeCvoexvbxm/1bmZfHmPWzdVceUDStCxuMpmkSHnEBn0iE7HW9FNRO3zScvM3A1U1JcS0nOxn37eZkecjI8ZKW7yUoLbBkel3U8JqQW1WkAI4FVqroGQEQmAGMB6zRMwu0or+Ff09fy6qwN7K6q5cju+dx+zhF0rFjPpeee2nQF9YgIuRkecjM89O6Yze41HkYd1ydk2aKiIk446RRKK73sKPeyo6KGnRVevpq3hI7deztpXkorvGyt8LNl1Xb2VNftu/KZsHxBk/FkprnITveQleYmM821r1PJTHPjcQlpbhdpHhfpbldgv97nNLeLdLfgce//vGpDLVtmbcAtgssluATcLsElgc3twvnqpLkkUFYIfHaOWb3LR8fiXbhkf8cmAoKwNyl4f1O5n1UlZXtLOnmBn7kEHV9S6Wfjzsp9+wSVEYHSaj9b91Q7tQT+s/cce2qUHeU1B/0c93iVnRXeA9LKnLTgbrncq5Q65aSB/lpoMCOk/EwP0lBlEWppnUZPYGPQfjFwXJJiMW1UWXUt/5q+lmenr6Gq1seYId24+oT+HNuvAyJCUdHGpiuJgXSPi4L8TAryM/eltdu1klGjjjigXOCKZhQQeM7y0adFDBtxHHuq6iirrqWsuo4Kbx3VtT6qvD6qav1U1fqC9gNbtfO5vKaOOp9S6/M7m9b7uj89pCULY/MD+GpG88p/MS28ctOmNp5f9GnDeVOnhE7/7JPI06Kw7L4xZKa5Y1qnaAtaGl1ELgXOVtWfOPtXASNV9Zf1yl0PXA9QUFBQOGHChLjGVV5eTm5ublLqCeeYpso0lt9QXqj0cNMSIR5t4ldlWnEdb6zwUl4LIwrcXDwwnR65rgaPiSbGeLVLotpEVfEp+PxQp1Dnh/KKCrKys/ErqIJ/78befQ3ksT/vgM9OmcqqatIzM53zBMrs/QxB+85/qqqryczM3LevQWX2xgpQVV1DRkbG/vSgOhSorqkhMyPjwHM6X6vrHbtXTc3B6TU1NaSHSMvIyDgwsOAYQyc3mA5wRh8P7iZuMzb0/8Po0aPnquqIg0+YAnMpwt2A7wGTg/ZvB25v7Bibp2Fjz2NRz8LiXTr28S+0760T9dJ/fqkLNu6K+tw2TyNx9aR6m4Rz7nhp7jyNlnZ7ajYwUET6A5uAy4AfJjck05rV+JR73l/Mi1+uo2NOOo98fxgXDe8Z8/vExrQULarTUNU6EbkRmExgyO1zqro4yWGZVurbDaXcPaOK7yrXcdXxfbn57MNpl5XW9IHGtGJRdRoi0hU4EegBVAGLCFzS+GMQW0iq+iHwYbzqN8Zb5+fvn63kiamraJ8hvPKT4zhhQOdkh2VMSoio0xCR0cBtQEfgW6AEyATGAYeKyBvAw6q6J1aBGpMIK7aW8dvX57Fo0x7GH9OL0zvstA7DmCCRXmmcC/xUVTfUzxARD3A+gQl4b0YRmzEJ4/Mrz32xloc+Xk5ehod/XlnImCHdKCoqSnZoxqSUiDoNVb2lkbw64J2IIzImwTburOR//zufWWt3cuagAh64+Cg65x48dNIYE8UzDRE5FShV1QUi8n3gFGA18KSqHjwt0pgUo6r8d04x97y/GBHhoUuGcklhLxsZZUwjIn2m8QQwFMgQkRVALvARcALwHHBFzCI0Jg5Kyqr5/VsLmbK0hOMP6chfLh1Grw7ZyQ7LmJQX6ZXGaFUdJCKZBOZLdFVVn4g8BTS9qI0xSTRp4RZ+//ZCKrw+7jp/ENec0M8W5zMmTJF2GtUAqlotIutV1efsq4jUxiw6Y2Jod1Ut97y3mLe+3cRRPdvxyPeHMbAgL9lhGdOiRNppdBWR3xJYW3HvZ5z9LjGJzJgY+mLldm55Yz4lZTX86vSB3HjaANLc9mIjY5or0k7jX0BeiM8Az0QVkTExVOmt48+TlvHSzPUc0iWHt352AsN6t092WMa0WJEOub0n1oEYE2tz15fyv6/PY92OSq49sT+3nH04WemxXSbamLYm0tFTjzWWr6o3RRaOMdGrqfPx3+VeJk3+ku7tsnj1p8fzvUM7JTssY1qFSG9PzXW+nggMAl5z9i8NyjMm4ZZs3sNvX5/Hsu9q+cGI3tx5/pHkZdoig8bESqS3p14EEJGrCQy/rXX2/wl8HLPojAlTnc/PU9PW8LcpK2ifnc6vj8ng15cMTXZYxrQ60S6N3oPAQ/Cdzn6uk2ZMwqzeVs7/vj6feRt3cd7Q7vxp7BDmz/4y2WEZ0ypF22n8GfhWRPa+VPdU4I9R1mlMWPx+5cWZ63jwo2VkeNw8dvlwLhxmf7MYE09RdRqq+ryITAKOc5JuU9Xvog/LmMat217BrW8u4Ou1Oxl1eBceHD+UgvzMZIdlTKsX6eipfqq6DsDpJN6tly9AT1UtjjpCY4L4/MrzM9byl4+Xk+Zy8eD4o/j+iN62yKAxCRLplcZDIuIi0FnMBbYReAnTAGA0cDpwN2CdhomZVSXl/O6N+XyzYRenHdGV/7voKLq1s6sLYxIp0tFTl4rIIAKr2V4LdAcqgaUEXsV6v6pWxyxK06bV+fw8PX0Nf5uykqw0N3/9wTDGHd3Tri6MSYKIn2mo6hLgjhjGYsxBln23h1v+u4CFm3YzZnA37h03mK55dnVhTLJEO3rKmLjw1vl5smgVT0xdRX5mGk/88BjOG9o92WEZ0+ZZp2FSzsLi3dzyxnyWfVfGhcN68McLB9MxJz3ZYRljSMFOQ0QeAi4AvAReH3uNqu5y8m4HrgN8wE2qOjlpgZqYq6718dinK3lq2ho65aTz9FWFnDW4W7LDMsYEieqFAhJwpYj8wdnvIyIjo4zpE2CIqg4FVgC3O3UPAi4DBgNjgCdFxJYsbSW+XL2dcx6dzpNFqxl3dE8++c2p1mEYk4KivdJ4EvADpwH3AmXAm8CxkVaoqsFrV30FXOJ8HgtMUNUaYK2IrAJGAjMjPZdJvl2VXv7vw6W8PqeYPh2z+fd1Izl5oL3Hy5hUJaoa+cEi36jqMSLyraoOd9Lmq+qwmAQn8j7wmqq+LCKPA1+p6stO3rPAJFV9I8Rx1wPXAxQUFBROmDAhFuE0qLy8nNzc3KTUE84xTZVpLL+hvFDp4aYBqCqzvvPxn6U1lNfCmH5pjB2QRoY7NsNoU71NwikXr3aJ1c8mEslql1Rvk+bEGGsNnXf06NFzVXXEQRmqGvEGfA24gW+c/S7At2EcNwVYFGIbG1TmDuBt9ndsTwBXBuU/C4xv6lyFhYUab1OnTk1aPeEc01SZxvIbyguVHm5acWmlXvP8LO1760Q9/7HpurB4V6PxRSLV2ySccvFql1j9bCKRrHZJ9TYJ59zx0tB5gTka4ndqtLenHnN+sXcVkfsJ3Eq6s6mDVPWMxvJF5MfA+cDpTvAQmF3eO6hYL2BzJEGb5PD5lRe/XMdfPl6OKtx53pFcfUI/PPaubmNajGgXLPyPiMwlsGyIAONUdWk0dYrIGOBW4FRVrQzKeg94RUQeIbD8+kBgVjTnMomzdMsebntrIfM37uLUw7rwp3FD6N0xO9lhGWOaKdIFCzsG7ZYArwbnqerOg48K2+NABvCJs0zEV6p6g6ouFpHXgSVAHfALVfVFcR6TAF6f8uBHy/jXtDW0y0rj0cuO5sJhPWwJEGNaqGhe96oEri76AKXO5/bABqB/pAGp6oBG8u4H7o+0bpNYM1Zt584ZVZRUruaSwl7cce6RdLBJesa0aJEuWNgf9r3e9T1V/dDZPwdo9HmFaf1KK7zc/+FS3phbTNds4ZWfHMcJAzonOyxjTAxE+yD8WFW9Ye+Oqk4SkfuirNO0UKrKm3OLuf/DpeyuquVnow5leNoW6zCMaUWi7TS2i8idwMsEblddCeyIOirT4qwqKefB2dUs2zmf4X3ac/+4oxjUI5+iInuRozGtSbSdxuUEXrb0trM/zUkzbUR1rY8npq7in5+vJk2U/7voKC47tjculz3oNqY1inbI7U7gVzGKxbQw01Zs4653F7F+RyUXDe/JqPaljD2uT7LDMsbEUVSdhohMJXBb6gCqelo09ZrUVrKnmifnVTPro1kc0jln34PuoqKiZIdmjImzaG9P3Rz0ORMYT2AOhWmFfH7lP1+v56GPllNV6+M3ZxzGDaMOIcNjiw0b01ZEe3tqbr2kGSLyeTR1mtS0aNNu7nh7IfOLd3PSgM6c372Cy84YmOywjDEJFu3tqeCZ4S6gELCXILQiVXXKve8v4YUv19IxJ2PfjO7PP7e/DYxpi6K9PRU8M7wOWEvgzXqmhVNVPlr0Hb+fXsUu71quOK4Pt5x9BO2y0pIdmjEmiaLtNI5U1ergBBHJiLJOk2Qbd1Zy93uL+WxZCb3zXDx73fc4pk+HZIdljEkB0XYaXwLH1EubGSLNtAC1Pj/PTF/Lo5+uwCXCnecdSf/a9dZhGGP2iXSV225ATyBLRIYTuD0FkA/Yetct0Ox1O7nj7YWs2FrO2YMLuPuCwfRon0VR0YZkh2aMSSGRXmmcDVxN4EVIjwSllwG/jzImk0ClFV6eW1TDtI9m0rN9Fs/8aARnDCpIdljGmBQV6Sq3LwIvish4VX0zxjGZBFBV3vxmE//34VJ2VdbxP6ccwq/OGEh2erR3LI0xrVmkt6euVNWXgX4i8tv6+ar6SIjDTIpY/l0Zd72ziFnrdnJMn/ZcdLSbq849MtlhGWNagEj/rMxxvubGKhATf+U1dTw6ZQXPzVhHfqaHB8cfxaWFvZk2zeZcGGPCE+ntqaecr/fENhwTD6rKBwu2cN/EJXy3p5rLR/bmd2cfYW/RM8Y0W7QzwrsAPwX6BdelqtdGF5aJlTXbyvnLnGoW7/iGwT3yefLKY2wIrTEmYtE+9XwXmA5MAXzRh2NipbrWx5NTV/HPz9fgFj9/vGAQVx7fF4/blezQjDEtWLSdRraq3hqTSEzMfLZsK3e/t5iNO6sYd3QPRnUoZdyJ/ZMdljGmFYj2z86JInJuTCIxUSsureT6l+Zw7QtzyPC4eeWnx/G3y4bTPsOuLowxsRHtlcavgN+LSA1QS2BmuKpqfrSBicjNwENAF1Xd7qTdTmBBRB9wk6pOjvY8rYG3zs8zX6zhsU9XIgi3jjmC607qT7rHOgtjTGxF+z6NvFgFEkxEegNnAhuC0gYBlwGDgR7AFBE5TFXb9LOUpTt83PfoNFZvq+DswQX84YLB9GyfleywjDGtVLSjp0ItTLgbWK+q0bzB76/A7wg8aN9rLDBBVWuAtSKyChhJYIHENqdkTzX3f7iUd+dV06djNs9ffSyjj+ia7LCMMa2cqB70iu/wDxb5isCKtgudpKOA+UAn4AZV/TiCOi8ETlfVX4nIOmCEqm4XkceBr5yZ6IjIs8AkVX0jRB3XA9cDFBQUFE6YMKH531wzlJeXk5sb/TzHcOrx+ZXPNtTx1iovtT44s5dy0RE5pLulwWOaqrex/IbyQqWHm5YIiWyTSI9JVrskq01iee7m1pPqbdKcGGOtofOOHj16rqqOOChDVSPegAnA4KD9QcDzwCHAvEaOmwIsCrGNBb4G2jnl1gGdnc9PAFcG1fEsML6pGAsLCzXepk6dmpB65q7fqef8bZr2vXWiXvXs17pmW3lY526qTGP5DeWFSg83LRES1SbRHJOsdklWm8Ty3M2tJ9XbJJxzx0tD5wXmaIjfqdE+CD9CVRcHdUBLRGS4qq4RafgvX1U9I1S6iBwF9AfmO8f3Ar4RkZFAMdA7qHgvYHOU8bcIpRVeHvxoGRNmb6RbfiZPXnEM5wzphoiwPtnBGWPalGg7jeUi8g8CVxwAPwBWOG/vq21uZaq6ENh3Y77e7an3gFdE5BECD8IHArOijD+l+f3K63M28uBHy9hTXcf1pxzCTacPJDfDVqI1xiRHtL99rgZ+DvyawHDbL4CbCXQYo6Os+wCqulhEXgeWEHgf+S+0FY+cWrx5N3e+s4hvN+xiZL+O3DduCId3i8tgNWOMCVu0Q26rgIedrb7yaOp26u9Xb/9+4P5o601le6preeTjFbw0cx0dc9J5+NJhXHxMTxq73WeMMYkS7ZDbgcADBB6AZ+5NV9VDooyrzVFVZm6u45aHP2d7eQ1XHteXm886nHbZackOzRhj9on29tTzwN0E5lWMBq5h//vCTZhWlZRx1zuLmbmmhmG92vHsj0cwtFf7ZIdljDEHibbTyFLVT0VEVHU98EcRmU6gIzFNqPTW8ffPVvHM9DVkpbn50aB07r7yRNwu63eNMakp2k6jWkRcwEoRuRHYRNDoJxOaqvLxkq3c+/4SNu2q4pLCXtx2zhEsmjPTOgxjTEqLttP4NZAN3ATcB5wG/DjaoFqzDTsq+eP7i/lsWQmHF+Tx3xu+x7H9OiY7LGOMCUu0o6dmOx/LCTzPMA2orvXx9LQ1PDF1FR6XcOd5R/LjE/qRZi9FMsa0IBF1Gs5Euwap6oWRhdM6TVuxjT+8u4h1Oyo5b2h37jpvEN3aZTZ9oDHGpJhIrzS+B2wEXiWwVpTdiA9hy+4q/jRxKR8s3EL/zjn8+7qRnDywS7LDMsaYiEXaaXQj8L6Ly4EfAh8ArwavQ9WW1fr8vDBjHX+dsgKfX/nfMw/j+lMPIcPjTnZoxhgTlYg6DWf5jo+Aj5x1pi4HikTkXlX9eywDbGlmrd3JXe8sYvnWMk47oiv3XDiY3h2zkx2WMcbERMQPwp3O4jwCHUY/4DHgrdiE1fJsL6/hgQ+X8eY3xfRsn8XTVxVy5qACW/7DGNOqRPog/EVgCDAJuEdVF8U0qhYk8FKkWm4qKqKq1sfPRx3KjacNIDvdVqI1xrQ+kf5muwqoAA4Dbgr6a1oAVdX8GMSW8uZv3MVd7y5iQbGXEw7txL1jhzCga3LeiGaMMYkQ6TONNj25YHdlLQ99vIz/fL2BzrkZ3DA0g1svP85uRRljWj27h9IMqsqb32zigQ+XUlrp5eoT+vGbMw/jm69mWIdhjGkTrNMI0/LvyrjrnUXMWreTY/q056XrRjK4R7tkh2WMMQllnUYTymvqeHTKCp6bsY78TA8Pjj+KSwt747KFBY0xbZB1Gg1QVT5c+B33TlzM1j01XD6yN787+wg65KQnOzRjjEka6zRCqPP5ue7FOXy+YhuDuufzjysLOaZPh2SHZYwxSWedRgget4vDu+Ux+vAuXHl8Xzy2Eq0xxgDWaTTo9+cemewQjDEm5dif0MYYY8JmnYYxxpiwiaomO4a4EpFtwPo4n6YdsDtJ9YRzTFNlGstvKC9Ueqi0zsD2JuKLh1Rvk3DKxatdktUmoWJJVD2p3iaQWv9WtgOo6piDSquqbVFuwNPJqiecY5oq01h+Q3mh0htIm2Ntklrtkqw2SWa7pHqbJLNdmvuztNtTsfF+EusJ55imyjSW31BeqPRY/RxiIdXbJJxy1i6xq8fapGHNiqXV354yySUic1R1RLLjMPtZm6SmltIudqVh4u3pZAdgDmJtkppaRLvYlYYxxpiw2ZWGMcaYsFmnYYwxJmzWaRhjjAmbdRrGGGPCZgsWmoQSkVHAfcBiYIKqFiU1IIOIuAi0ST6BCWYvJjmkNk9ETgauIPA7epCqnpDkkPaxKw0TNRF5TkRKRGRRvfQxIrJcRFaJyG1OsgLlQCZQnOhY24pmtslYoCdQi7VJ3DSnTVR1uqreAEwEUqoTtyG3JmoicgqBjuAlVR3ipLmBFcCZBH4RzQYuB5apql9ECoBHVPWKJIXdqjWzTS4ESlX1KRF5Q1UvSVLYrVpz2kRVlzj5rwM/UdU9yYn6YHalYaKmqtOAnfWSRwKrVHWNqnqBCcBYVfU7+aVARgLDbFOa0yYEflmVOmV8iYuybWlmmyAifYDdqdRhgD3TMPHTE9gYtF8MHCciFwNnA+2Bx5MRWBsWsk2AR4G/O/fRpyUjsDasoTYBuA54PuERNcE6DRMvEiJNVfUt4K1EB2OAhtukksAvKJN4IdsEQFXvTnAsYbHbUyZeioHeQfu9gM1JisUEWJuknhbXJtZpmHiZDQwUkf4ikg5cBryX5JjaOmuT1NPi2sQ6DRM1EXkVmAkcLiLFInKdqtYBNwKTgaXA66q6OJlxtiXWJqmntbSJDbk1xhgTtqQ8CBeRS4E/AkcCI1V1jpOeBjwDHOPE9pKqPuDkFQIvAFnAh8CvNIwer3PnztqvX7/YfxNBKioqyMnJSUo94RzTVJnG8hvKC5UebloipHqbhFMuXu2SrDaJ5bmbW0+qt0lzYoy1UOft3LkzkydPnqyp8o5wAp3F4UARMCIo/YcElpYAyAbWAf2c/VnA9wiMNpgEnBPOuQoLCzXepk6dmrR6wjmmqTKN5TeUFyo93LRESPU2CadcvNolWW0Sy3M3t55Ub5Nwzh0vDZ2XBt5ZnpRnGqq6VFWXh8oCckTEQ+CKwgvsEZHuQL6qznS+mZeAcYmL2BhjDKTeg/A3gApgC7AB+Iuq7iQwASZ4TZxiJy1utpfXsKO8hrLqWmrqfHuvhIwxpk2L24NwEZkCdAuRdYeqvuuUKQJu1v3PNE4Efg5cDXQApgPnAJ2AB1T1DKfcycDvVPWCBs59PXA9QEFBQeGECROaHf/vv6hkc/mBPxuPC9Jce78KHuezS/1keNwHpQfKCmnO53S3kOGGDLeQXu9rhhvqaqrokJd9QJrbFWruz37l5eXk5uZGVaax/IbyQqWHm5YIsTpvJPWEe0yy2iVZbRLLcze3nlRvk+bEGGsNnXf06NFzVXVE/fS4PQjf+wu+mX4IfKSqtUCJiMwARhDoPHoFlWt0AoyqPo3zkvYRI0boqFGjmh3I7R02s6O8Bq/Pj7cusNUEffbW+fflbd66jfz2HagJyiv3+amp8e0/ts5PVW0djffRAlQdkJLmFrLS3GSne8jLDGy5mWnkZXrIz/RQWuJl8MCe5GZ6yHPSA3mBzx1y0pk78wsa+xkUFRU1mN9QXqj0cNMSIVbnjaSecI9pqly82iVZbRLLcze3nlRvk+bEGGvNPW+qLSOyAThNRF4m8CD8eOBvqrpFRMpE5Hjga+BHwN/jGcgFw3qEXTbwQz+uyXKqGug8vD6qan1Uen1UO1+ran3M/mY+hx5+RGDf2Sprna/eOsqqA9vuSi/FOyvZU13H7spaPlq3otHzugU6fjmFjtnptM9Oo2NOOh1y0umYHfi6dXMdLC+hID+TrnkZdMhOx9XEFY4xpm1K1pDbiwj80u8CfCAi81T1bOAJAgt0LSLwZ/fzqrrAOexn7B9yO8nZWhQRITPNTWaamw4h8nWzh1HDe4XIaVhRUREnnHQK5TV1lFXXUlZdx569X6tq2VVZy7ylq8jr3JWdFV52VdaysqSc0govpZVe/M6Vz9MLZu+rM80tdM3LpEteBm5vNZ/tXkTXvAy65mfSs30WvTpkUee3ZzzGtEVJ6TRU9W3g7RDp5cClDRwzBxgS59BapHSPi46edDrmpIfML/JvYNSooQel+/3KnupaPpr6BQMGH01JWQ1b91SzdU8NJWXVlOypYW2pn1XzNrO7qvaAYwXo/vWn9OqQTa8OWfTqmE3F1lqy1uzg0K65dMpJR8SuVoxpbVLt9pRJIJdLaJ+dTrccFyP6dQxZZu/9zupaHyV7ati0q4ri0kq++HYp7nadKC6t4uu1O3ln3ib8Cs8u+gqA9tlpHNoll+y6GpbLag7tksuArrn06Zhtt76MacGs0zBhyUxz06dTNn06ZQOd6FK+mlGjjt6X763z887HRRQMOIrVJeWs3hbYFnznY/qkZfvK5WZ4OLJ7HoN7tGNQj3wG98hnYNc80j2pNvrbGBNKq117SkQuAC7o2bPnT19++eW4nivVh3cmexihZOTwXYWf4nI/G/b4Wb/Hz8YyPzXOO+I8An3zXQzs4GZgh8DX/PTorkZSvU3CKWdDbmNXT6q3SXNijLXmDrlNyjIiidxsGZHUXBqhzufXVSVl+u68TXr/B0t0/JMzdODvP9S+t07UvrdO1NEPTdVb35ivE+dv1l0V3kbjb07ciagn1ZessGVEIi9ny4io3Z4yyeF2CYd2yeXQLrlc6Axvrq71sWjTbuasL2XOup18sHALE2ZvxCUwrHd7ThnYhdOO6MrQXu3sIbsxSWKdhkkZmWluRvTrGHgof+qh1Pn8zC/execrtjNtxTYe+2wlj366kp7tszhnSGbRKhUAABw/SURBVDfOOaobw3t3sAfrxiSQdRomZXncLgr7dqSwb0d+e+ZhlFZ4+WxZCZMWbeGlmet55ou1dG+XyfhjevH9Eb2dh/TGmHiyTsO0GB1y0hlf2Ivxhb0oq67l06UlvDNvE08WreLxqav43iGduGxkb849qnuyQzWm1bJOw7RIeZlpjBvek3HDe7JldxVvzi3m9TnF/GrCPP48aRmndPNzzPG15GemJTtUY1oVGxxvWrzu7bK48bSBFN08iuevPpZ+nXJ4bbmXEx74jAc/WsbuytqmKzHGhMXmacRAqs8JaItjz5d8V07Rdx5mf+cjywPnHpLGmX3TyHA376G5zdOILZun0by0RLB5GjZPo9llWvPY88Wbdus1z8/SvrdO1JH3f6Lvz9+kfr+/2fVEcu5oy9k8jdjVk+ptEs6546VFvO7VmEQZ1COf564+ltf/53t0zs3gxle+5cfPz2b9jopkh2ZMi9RopyEivUTkZhF5V0Rmi8g0EXlSRM4TEetwTIsxsn9H3v3Fidx9wSC+WV/KWX+dxr9nrrPX+BrTTA3+4heR54HnAC/wIHA5gVexTgHGAF+IyCmRnFREHhKRZSKyQETeFpH2TvqZIjJXRBY6X08LOqZIRJaLyDxn6xrJuU3b5XG7uObE/kz57akcf0gn7np3Mde9OIdtZTXJDs2YFqOxq4WHVfUsVX1MVb9U1VWqukhV31LVXwKjaOSVq034BBiiqkOBFcDtTvp24AJVPQr4MfDvesddoapHO1tJhOc2bVy3dpm8cM2x3HPhYL5YtZ3zHpvONxtKkx2WMS1Cg52Gqi5q7EBV9arqqkhOqqofq2qds/sVzvu/VfVbVd3bES0GMkUkI5JzGNMYEeHHJ/TjvRtPJDPNzWVPfcV/52xMdljGpLwmn0s4t4oW1Numi8hfRaRTDGK4ltCvbh0PfKuqwfcOnnduTd0ltmKdiYEjuuXz7i9O5Nj+HbjljQU8MGmpPecwphFNztMQkf8H+IBXnKTLnK97gJNU9YIGjpsCdAuRdYeqvuuUuQMYAVysQYGIyGDgPeAsVV3tpPVU1U0ikge8Cbysqi81cO7rgesBCgoKCidMmNDo9xgtm6fRvLREaO55fX7lP0u9fLaxjpN7erh6cDpul9g8jRizeRrNS0uEmM/TAGY0lAYsbOr4Rur9MTATyK6X3ovAc44TGzn2auDxcM5j8zRs7Hm4/H6/Pvzxcu1760S9/qXZWlPrs3kaMWbzNJqXlgjxmKeRKyLH7d0RkZHA3m6pLvQhjRORMcCtwIWqWhmU3h74ALhdVWcEpXtEpLPzOQ04H2j0mYsxzSUi/PbMw/jD+YOYvHgrv3ltHj6/3aoyJlg4Cxb+BHhORPZ2FGXAT0QkB3ggwvM+DmQAnziPJr5S1RuAG4EBwF0icpdT9iygApjsdBhuAsN+/xXhuY1p1LUn9cevyp8+WErpDg+jR6m9s8MYR5OdhqrOBo4SkXYEnoHsCsp+PZKTquqABtL/BPypgcMKIzmXMZH4ycmHUOn18cgnK/jTB0v5wwWDkh2SMSkhnNFTBSLyLDBBVXeJyCARuS4BsRmTVL88bQBn9vXw3Iy1/HvmumSHY0xKCOeZxgvAZKCHs78C+HW8AjImVYgIlx+RzhlHduXu9xYzdZnNJzUmnCG3s1X1WBH5VlWHO2nzVPXohEQYIVsaPfwybWkYYST1eDJzeGBWNVsr/PzxhCy65TT+t1aqD++0IbeRl2tL/1aiGXJbBHQCvnH2jwc+b+q4VNlsyK0NI4xFPcWllXr0PZP17L9+rpU1dTE5tw25TVw9qd4m4Zw7XuIx5Pa3BCbaHSoiM4CXgF9G1qcZ0zL1bJ/F3y4bzvKtZdz5ziKbNW7arHBGT30jIqcChwMCLFdVe3+maXNOPawLN502kEc/Xcmx/Tpw2cg+yQ7JmIRrsNMQkYsbyDpMRFDVt+IUkzEp66bTBzJ3fSn3vL+E4w7pRP/OOckOyZiEauz21AXOdh3wLHCFsz0DXBn/0IxJPW6X8JdLh5HucfGb1+ZR5/MnOyRjEqqxpdGvUdVrAAUGqep4VR0PDE5YdMakoG7tMrn/oiHM27iLJ6auTnY4xiRUOA/C+6nqlqD9rcBhcYrHmBbh/KE9GHd0Dx77bCXzN+5q+gBjWolw5mk8DgwEXiVw1XEZsEoDb+9LWTZPI/wybWnseSzrqahV7ppRRaYb7jkxizRnfapUnxNg8zQiL9eW/q1EPE/D6VQuAv7qbBeFc0yqbDZPw8aex7Oez5Zt1b63TtSHP17e7HPbPI3E1ZPqbRLOueOlufM0Ghs9Jc6BqOrbwNuNlTGmLRp9eFcuGt6TJ6eu4pwh3Tiye36yQzImrhp7pjFVRH4pIgcMRheRdBE5TUReJPAiJWPatLvOH0S7rDRue3OBvX/DtHqNdRpjCLzm9VUR2SwiS0RkLbASuBz4q6q+EMlJReQhEVnmvG/8beflS4jIFc47wPdufhE52skrdN5XvkpEHrN3hJtU0TEnnbsvHMz84t08P2NtssMxJq4aG3JbrapPquqJQF/gdGC4qvZV1Z+q6rwozvsJMERVhxJYNfd255z/UdWjNbAY4lXAuqDz/IPAe78HOtuYKM5vTExdMLQ7ZxzZlb98vJySSpu7YVqvcIbcoqq1qrpFD3wBU8RU9WNV3fuq2K8IvBe8vssJjNhCRLoD+ao603mG8hIwLhaxGBMLIsJ944bgcbl4flGNrU1lWq0mh9zGPQCR94HXVPXleumrgbGqukhERgB/VtUznLyTgVtV9fwG6ryewFUJBQUFhRMmTIjr92BDbpuXlgjJapOpG2p5cYmXa4akc2qvtKjqtiG3savHhtw2LC5DbiPZCLzHe1GIbWxQmTsIjMqSesceBywM2j8WmBK0fzLwfjhx2JBbG0aYyHp8Pr+e9ecPdcjdH+l3u6uiqtuG3MauHhty27B4LI2OiPQVkb1/5WeJSF5Tx6jqGao6JMT2rlPPj4HzgSucAINdhnNrylHMgbewegGbw4ndmERyuYRrhmTgrfPbEuqmVQrnHeE/Bd4AnnKSegHvRHNSERkD3ApcqKqV9fJcwKXAvntKGljGpExEjndGTf0IeDeaGIyJl245Lv73rMP4ZMlWPlz4XbLDMSamwrnS+AVwIrAHQFVXAl2jPO/jQB7wiTO09p9BeacAxaq6pt4xPyOwwu4qYDUwKcoYjImba0/sz9Be7bj7vUWUVniTHY4xMRNOp1Gjqvv+rxcRD4E1qCKmqgNUtbc6w2tV9YagvCJVPT7EMXOc21uHquqNIW5pGZMyPG4XD44fyq7KWu77YEmywzEmZsLpND4Xkd8DWSJyJvBf4P34hmVMy3dk93x+PupQ3vpmE0XLS5IdjjExEU6ncRuwDVgI/A/wIXBnPIMyprX4xWkDGNA1lzveXkR5TV3TBxiT4po1T0NEOgK9VHVB/EKKDVsaPfwybWnseSLqqX/MqlIf939dzcm9PFw7JCPsum2eRuzqsXkaDYv5PA2gCMgHOgIbgLnAI00dlyqbzdOwseeJrifUMX+etFT73jpRJy3cEnbdNk8jdvXYPI2GxWOeRjtV3QNcDDyvqoXAGRF2asa0Sb854zCG9mrHbW8t4Lvd1ckOx5iIhdNpeJy1n74PTIxzPMa0SukeF3/7wdHU1Pr57evz8NsS6qaFCqfTuBeYTOAVr7NF5BACy6MbY5rhkC65/PHCQXy5egf/+Hx1ssMxJiJNdhqq+l9VHaqqP3f216jq+PiHZkzr8/0RvblgWA8e/ng5i7bbaCrT8jT4ute9RCQTuA4YDGTuTVfVa+MYlzGtkojw4PijWLm1jH/ML2PsaZX07pid7LCMCVs4t6f+DXQDzgY+J7D2VFk8gzKmNctO9/DUVYWowvX/nkul1644TMvR5DwNEflWVYeLyAJVHSoiacBkVT0tMSFGxuZphF+mLY09T0Q94R7z9YZy/rlEGNrFzU3DM3C7DnyDsc3TiF09Nk+jYfGYpzHL+ToNGAJ0BtY0dVyqbDZPw8aeJ7qe5swJeGnmOu1760S97c0F6vf7w67H5mk0rx6bp9Gw5s7TaPKZBvC0iHQA7gLeA3KBP0TYqRljglx1fF+27KriyaLVdM3L4DdnHpbskIxpVJOdhqo+43z8HDgkvuEY0/bccvbhbCur4dFPV+JxCb88fWCyQzKmQeGMnsoAxgP9gsur6r2RnlRE7gPGAn6gBLhaVTeLSCcCL3w6FnhBVW8MOqYI6A5UOUlnqaotHWpaPBHhz+OH4vMrD3+yAhG48TTrOExqCuf21LvAbgJrTtXE6LwPqepdACJyE4HbXTcA1QRugw1xtvquUNU5MYrBmJThdgkPXToMgL98vIJdlbWckGOzxk3qCafT6KWqY2J5Ug2sZbVXDs5LnVS1AvhCRAbE8nzGtAR7O478rDSe+WItC7q5OeEkH5lp7mSHZsw+4Qy5fRr4u6oujOmJRe4n8K7v3cBoVd0WlHc1MCLE7alOgA94E/iTNhC8iFwPXA9QUFBQOGHChFDFYibVh3faMMLE1hPt8E5V5aN1dby23MvA9i5+cXQG7TNdYR1rQ26jK2//VvZr9pBbAi9dWgAsAWqB5c7+QmBBQ8cFHT8FWBRiG1uv3O3APfXSrgYer5fW0/maB3wM/KipGNSG3IZVpi0NI0xEPbEa3vngK5/oEXdO0sL7PtGZq7eHdawNuY2uvP1b2Y8Ihtye35zeKkRnFO7y6a8AHwB3N1HfJudrmYi8AowEXoomRmNS2cjuHsaddhw3vDyXH/7rK359xmH8bNShpLnDWcjBmPho7P++rcBFwC3AGGCTqq7fu0VzUhEJHhpyIbCsifIeEensfE4j0KEtiiYGY1qCwwryeO/Gk7hgWA8e+WQFFz/5JSu22io+Jnkau9J4kcBtqenAOcAg4FcxOu+fReRwAkNu1xMYOQWAiKwj8KbAdBEZB5zllJnsdBhuAre+/hWjWIxJabkZHh69bDhnD+7Gne8s4vzHvmBMPzfHneAjK90ekpvEaqzTGKSqRwGIyLPArFidVBtZWl1V+zWQVRir8xvTEp17VHdG9u/IfROX8O68zcx5uIjfn3ck5x3VHRFpugJjYqCx21O1ez+oqi3DaUwK6JybwaOXDef2kZm0z07nxle+ZdwTM5i6vGTvgBFj4qqxK41hIrJ3PoUAWc6+AKqq+XGPzhgT0uEd3bw/7iTenFvMY5+t5JrnZ3N07/aM6lLHKX7F5bIrDxMfTc7TaKlsafTwy7SlseeJqCfRcwLq/MoXm+p4f3UtO6qVgmzhjL5pnNTTQ5ZHbJ5GM8rbv5X9Il4avaVvNk/Dxp4nup5kzQmoqfXpA//5RMc+/oX2vXWiDvnDR3rbmwv06bemHLDsus3TiLxcW/q3QhRLoxtjWoB0j4vje3i47Ycn8u2GUl6auZ63vy2mutbPK6s/5+LhPTl3aPdkh2laOOs0jGmFhvfpwPA+Hbh37GD+9kYRiyoyePiTFTz8yQp65AgX1Szj7MHdOKpnOxt5ZZrFOg1jWrG8zDRO7pXGXaO+x6ZdVXyy+Dtem7GMf36+hiemrqZzbgYnDejEiQM6c9LAznRvl5XskE2Ks07DmDaiZ/ssrj6xP/1q1zPs2BP4dFkJ01du44tV23ln3mYADu2Sw8j+HTmmTweO6duBQzrn2JWIOYB1Gsa0QR1y0rmksBeXFPbC71eWby1jxqrtfLFqOxMXbOHVWRsBaJ+dxjF9OnB07/YM7pHPoB75dMvPtI6kDbNOw5g2zuUSjuyez5Hd8/nJyYfg9yurt5Uzd30p32woZe76Uj5btv8lme2z0xjUPZ9B3fM5ons+h3bJ4ZAuubTLSkvid2ESxeZpxECqzwmwseeJrSfV5wRE8j1V1SnFZX7W7/GzoczPxj1+isv91Pr3l8lPF7rnCN1zXHTLcdEtR+iS5aJTlpDpkYjPHYrN04gdm6dh8zSaXaYtjT1PRD2pPicgVj+b2jqfrtxappMXbdF/FK3Sm1+fpxc/OUOH3TNZ+9468YBt+L0f6wV/n67jH5mk93+wRF/8cq1OWfKdLizepVt3V2mdz9/0CaP4HlK9TcI5d7zYPA1jTEJ43C4GdM1lQNeD/0rdWeFl7fYKiksrKS6tcrZKVu70s+DLdXjr/AeUdwl0ycuga14mBfkZdM3PpKuz3zEn3dnS6JCdTvvs9ER9iyYE6zSMMTG39xd9Yd8OB6QXFRVxyimnsr28hk27qigpq6FkTzUlZTVs3VPN1j01bNpVzbcbdrGjwhuybhHI9kDBnCI65KTTITvQobTLSiMvM438TA95mWnkBX0tqfRTWuElL9ODx15iFZWkdBoich8wlsD7NEqAq1V1s5M3FHiKwDs1/MCxqlotIoXAC0AW8CHwK+cSyhjTgrhcEriSyM9stJy3zs+OihpKK2rZWeFlZ6WX0govOyu8LFq5lqz2+ZRWetm0q4qFm3axu6qW6lp/wxVO+wSArDQ3+VkecjM8ZKd7yEp3k+1su3fU8OmuRc6+h+x09wH5K7fVkbNuJ5keN+keFxkeFxlpLjI8bqrrlFqfH49LWvXosmRdaTykqncBiMhNwB+AG0TEA7wMXKWq80WkE/uXaP8HcD3wFYFOYwwwKeGRG2MSIt3jonu7rJATDovSNjNq1DEHpdf6/JRV11FWXUtZdR17qmrZU13H7HkL6dlvwL68PdW1lNfUUen1Uen1saPcy0ZvHaVlPhbv2kyl13fQLbR95s5sOOgpk3AJZHjcZKS5SHe70Dov7b75PNDBeFyke1xU7Knm5fVzSHMLHreLNJewraSGj0sXkuYKpHlcgscteFyufeU8LiHN7cLjFtJcga97j/c46R6X4HYJbhFG9OuIO8YrHiel01DVPUG7OcDeK4azgAWqOt8ptwNARLoD+ao609l/CRiHdRrGmCBpbte+W2PBMrcvY9RJ/Zs8vqioiFGjRgFQ5/NTVeujyuujwuuj0lvHl1/P4fDBQ/HW+amp81NT59v3ecnylfTq0w+vz8mr9eH1+VlfvIUOnfL2la+p81NRq2zeVUWd30+dT6n1+6mo9LF099b9aT4/dX7F54/8hsqy+8bgdsX27Y5JG3IrIvcDPwJ2A6NVdZuI/JrAG/q6Al2ACar6/0RkBPBnVT3DOfZk4FZVPb+Buq8ncFVCQUFB4YQJE+L6vaT68E4bRpjYelJ9eKctjR55uWT8W/Gr4lfw+aFOwafg86vz1dl30uoU1Nn3KxzR0YWriVtlKTPklsB7vBeF2MbWK3c7cI/z+WZgLdAZyAZmAqcDxwJTgo45GXg/nDhsyK0NI0x0Pak+vNOWRo+8XFv6t0Kih9yqc1UQhleAD4C7gWLgc1XdDiAiHwLHEHjO0SvomF7A5thFa4wxJhxJGXsmIgODdi8EljmfJwNDRSTbeSh+KrBEVbcAZSJyvASGJfwIeDehQRtjjEnOMw0ReRM4nMCQ2vXADaq6ycm7ksAtKwU+VNXfOekj2D/kdhLwSw0jeBHZ5pwjntoReDaTjHrCOaapMo3lN5QXKj1UWmdgexPxxUOqt0k45eLVLslqk1CxJKqeVG8TSK1/K9sBVHXMQaVD3bOyrdnPb55OVj3hHNNUmcbyG8oLld5AWsj7om29TZLZLslqk2S2S6q3STLbpbk/S5saGRvvJ7GecI5pqkxj+Q3lhUqP1c8hFlK9TcIpZ+0Su3qsTRrWrFha7Sq3JjWIyBwNNWzPJI21SWpqKe1iVxom3p5OdgDmINYmqalFtItdaRhjjAmbXWkYY4wJm3UaxhhjwmadhjHGmLDZS5hMQonIKOA+YDGBBSmLkhqQQURcBNokn8BcgReTHFKb5yzKegWB39GDVPWEJIe0j11pmKiJyHMiUiIii+qljxGR5SKySkRuc5IVKAcyCaw1ZuKgmW0yFuhJ4N011iZx0pw2UdXpqnoDMBFIqU7cRk+ZqInIKQQ6gpdUdYiT5gZWAGcS+EU0G7gcWKaqfhEpAB5R1SuSFHar1sw2uRAoVdWnROQNVb0kSWG3as1pE1Vd4uS/DvxED3wHUVLZlYaJmqpOA3bWSx4JrFLVNarqBSYQWBZ/7+vQSoGMBIbZpjSnTQj8sip1yvgSF2Xb0sw2QUT6ALtTqcMAe6Zh4qcnsDFovxg4TkQuBs4G2gOPJyOwNixkmwCPAn937qNPS0ZgbVhDbQJwHfB8wiNqgnUaJl5CvS5MVfUt4K1EB2OAhtukksAvKJN4IdsEQFXvTnAsYbHbUyZeioHeQfv24qzkszZJPS2uTazTMPEyGxgoIv1FJB24DHgvyTG1ddYmqafFtYl1GiZqIvIqgfe5Hy4ixSJynarWATcSeBvjUuB1VV2czDjbEmuT1NNa2sSG3BpjjAmbXWkYY4wJm3UaxhhjwmadhjHGmLBZp2GMMSZs1mkYY4wJm3UaxhhjwmadhmkzRMQnIvOCtn7JjilWRGS4iDzjfL5aRB6vl18kIiMaOX6CiAyMd5ym5bO1p0xbUqWqR4fKEBEhMG/JHyq/Bfg98Kcojv8H8Dvgp7EJx7RWdqVh2iwR6SciS0XkSeAboLeI3CIis0VkgYjcE1T2DudFOVNE5FURudlJ3/cXvIh0FpF1zme3iDwUVNf/OOmjnGPeEJFlIvIfp8NCRI4VkS9FZL6IzBKRPBGZLiJHB8UxQ0SG1vs+8oChqjo/jO/5wqArreUistbJmg6cISL2h6RplP0PYtqSLBGZ53xeC/wGOBy4RlV/LiJnAQMJvONAgPecF+dUEFgTaDiBfzPfAHObONd1BN6FcKyIZAAzRORjJ284MJjAwnQzgBNFZBbwGvADVZ0tIvlAFfAMcDXwaxE5DMhQ1QX1zjUCWFQv7QciclLQ/gAAVX0PZ20j5wU/nzvpfhFZBQwL43szbZh1GqYtOeD2lPNMY72qfuUkneVs3zr7uQQ6kTzgbWcJcUQknAXlzgKGisjet+C1c+ryArNUtdipax7QD9gNbFHV2QB7X7wjIv8F7hKRW4BrgRdCnKs7sK1e2muqemPQ91oUnCkivyPw83giKLkE6IF1GqYR1mmYtq4i6LMAD6jqU8EFROTXOO84CKGO/bd5M+vV9UtVnVyvrlFATVCSj8C/Qwl1DlWtFJFPCLzN7fsErirqq6p37kaJyOnApcAp9bIynbqMaZA90zBmv8nAtSKSCyAiPUWkK4G32V0kIlnO84MLgo5ZBxQ6ny+pV9fPRCTNqeswEclp5NzLgB4icqxTPi/o+cIzwGPAbFWt/7pQCKyOOiCcb1BE+gJPAt9X1fodxGFASq+wapLPrjSMcajqxyJyJDDTeTZdDlypqt+IyGvAPGA9gYfGe/0FeF1ErgI+C0p/hsBtp2+cB93bgHGNnNsrIj8g8NrVLAJ/8Z8BlKvqXBHZQwOv/lTVZSLSTkTyVLWsiW/zaqAT8LbzPW5W1XNFpIDA7aotTRxv2jhbGt2YZhKRPxL4Zf6XBJ2vB1AEHNHQkGAR+Q1QpqrPRHiO3wB7VPXZiAM1bYLdnjImhYnIj4CvgTuamEPyDw58VtJcu4AXozjetBF2pWGMMSZsdqVhjDEmbNZpGGOMCZt1GsYYY8JmnYYxxpiwWadhjDEmbNZpGGOMCdv/B7E5jUT0wM6UAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "matlab.bode(GH, Hz=True)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", " 1.551e-13 s^2\n", "---------------------------------------------------------------\n", "1.126e-27 s^4 + 1.622e-20 s^3 + 1.918e-13 s^2 + 4.191e-07 s + 1\n", "\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOydd3iUVdbAf2cmjSQU6VU6KCKWACpKU1Ts3bXXz7L2VXexrKuuuuva17Xv2thVsWChg2gCqHREIfRepNckpExmzvfHO5EYJpkaZpKc3/O8T+ae99z7Ho5jbt5z7j1XVBXDMAzDCAVXvA0wDMMwag42aRiGYRghY5OGYRiGETI2aRiGYRghY5OGYRiGETI2aRiGYRghkxRvA6qbpk2baocOHWI+bkFBARkZGdXSL5hOZfcDySvKwm3HkkjGNn/F3l/B9Mxf4elF46+KskTxV9OmTZk4ceJEVR16wE1VrdVXVlaWVgfZ2dnV1i+YTmX3A8krysJtx5JIxjZ/VU+fqvTMX+HpReOvirJE8hcwRwP8TrXwlGEYhhEytT48ZdQdtuUV8/OG3azdsY/Ne4vYnl9McamPTZuL+HDdHOqluElLclMvxbm2bixhTfJqMtOSyUx1k5GaxIrdXtpsySMjNYmM1CR8VjHBMH5DjZs0RGQo8E/ADfxHVZ+Os0lGnFBVZq3eyQeLi/nLrGzW7dz3673UJBdNM1NJTXbhKVL27dxHkcdLocdLYYmXfSVeSn3KyOWLDhj3yRlTf9Ou9+0EMlKTyEx1o54iWi6dTnqKm7RkN7t3FjFx58+kJrlJTXaRVu5nWrKb1CQXaclu0pJdNKyXTMuGabRokEay217yjZpJjZo0RMQNvAqcCmwAZovIKFU98P98o9ayr6SUkXM38P70tazYmk+SCwZ0a8w1J7SnV9tGdGmeySHpyYgIADk5OQwaNOA3Y6gqX3+bQ9Zx/Sgo9pJfXEpBSSk/zJ5Hp249KCguJb+4lNylK2jWuh15RaUUFJey9pctAGzPL6G41MuuvT5W5m2lyOOluNRHcakvqP0i0Cwzla4tMunRqgE9Wjegb8cmtGlUL/bOMowYU6MmDaAvsEJVVwGIyAjgPMAmjTpASamPj2at41/fLmd7fglHtmnIsxf3ov7uFQwd0iessUSEFLfQJDOVJpn75QVrkhh0VOtf2znedQwadPj+dk4OgwadUKE96Ne2z6eUeH0Ue3wUlXp/nUyKPF527fOwaXchm/YUsXF3Icu25PH+9LWU+CeaTs0yGNC1Gecc1YpjDz0kTO8YxsFBtAbFbEXkYmCoqv6fv301cJyq3lFB72bgZoAWLVpkjRgxIua25Ofnk5mZGVwxgn7BdCq7H0heURZuO5ZEMnZZn5+3lfLfRSVsK1S6H+Liom4pdG3kQkRqtL+8PmVjvo/FO33kbveyZKeXEh80qyf0burj9C7pNEoNPZQVqj1V6SWyv6oimu9XNHrR+KuiLJH8NXjw4Lmq2vuAG4GWVCXqBVyCk8coa18N/KuqPrbkNjuqdiyJZOyxk77VP346X9sPG6NDns/RnKVb1efzhTVuTfJXXpFHP5uzXq/6zwztMGyMdn1onN778Xxd9MuekPrHewlpTft+xdtfFWWJ5C8qWXJb08JTG4B25dptgV/iZItRzUxZto0/f1/I7uIN3DaoM3cP6UpqkjveZlUrmalJXJTVlouy2jJi7LcsKm3Op3M2MHLeBs48siX3ntqNLs3rx9tMow5T0yaN2UBXEekIbAQuA66Ir0lGrNlb5OGpMYv5eM56WmcIn992Ike3axRvsw46LTNcXDaoJ/ee2o13vlvN29+tZsLCzVxwTFvuO60brS1xbsSBGjVpqGqpiNwBTMRZcvuOqubG2SwjhkxZto0HRv7Mlr1F/H5QZ45J3lQnJ4zyNEpP4d7TunNtvw68MWUlw6evZeyCX/j9wC7cPKAT9VJq99uXkVhENWmISHPgRKA1UAgsxImDBV93GCGqOg4YV13jG/Ehr8jDU2MXM2L2ero0z/z17SInZ3O8TUsYmmSm8vBZPbi2Xwf+Pm4JL05exidz1vPgmYdx1pGtfl1ibBjVSUSThogMBh4AGgM/AluBNOB8oLOIfAY8r6p7Y2WoUXuZ6n+72Ly3iFsHduaeIV1JS7a/niuj7SHpvHrlsVy9agePj17EHR/+yAed1vHkBT3jbZpRB4j0TeNM4CZVXVfxhogkAWfjbMAbGYVtRi0nr8jD38Yt5qNZ6+ncLIORv+/HMbY/IWSO79SEMXeexIez1vHMhCWc8dI0zujg5vgTvTbpGtVGRJOGqv6xinulwJcRW2TUCaYt38awz5y3i1sGduIPQ7rZL7oIcLuEq49vz+lHtODJMYv56qdf+Pmf03jq/J7069I03uYZtZCIcxoiMhDYpao/i8ilwABgJfCaqhbHykCjduG8XSzho1nr6Nwsg89+3892P8eA5vXTePnyY+iWvINPVytX/GcmFxzThofPOpymmanxNs+oRUSa03gV6AWkisgyIBOYAPQD3gGujJmFRq0hd7uXh1+axqY9hdwyoBN/ONXeLmJNz6ZJ/N95/Xk1ewVvTFnJt0u28sAZh/G73u1wuSxRbkRPpG8ag1W1h4ik4eyXaK6qXhF5E/g5duYZtYG8Ig9/H7+ED+cU0alZBp/e2o+s9vZ2UV2kJbu577TunHd0ax7+YiEPfr6AkXM38LcLj6RbC9sYaERHpPWZiwBUtQhYq6pef1sBT4xsM2oBU5Zt4/QXpzJi1jqGdkhi3F39bcI4SHRpXp8RNx/Psxf3YsW2fM785zSenbiEIo833qYZNZhI3zSai8i9gJT7jL/dLCaWGTWaPfs8PDl2EZ/O3UCX5pl89vt+7F31k4WjDjIiwiW923HyYc3527glvJq9ktE/beLSTl4Gxds4o0YS6ZvGv4H6OLmMss9l7f/ExjSjpjJ50RZOfXEKn/+4kdsHd2bMnSdZsjvONMlM5flLj+LDm44jySU8N6eYu0f8yLY8W7NihEekS24fj7UhRs1nV0EJj43O5av5v3BYy/q8fW0fjmzbMN5mGeXo17kp4+7uz7D3vmH8gs1kL9nKA2cczmV9LFFuhEakq6deruq+qt4VmTlGTWXcgk385auF7N7n4Z4hXbltUBdSkuxI00QkLdnNBV1TuOv83jz8xQIe+mIBI+dt4G8XHEn3lpYoN6om0pzGXP/PE4EewMf+9iXl7hl1gK15RTw2KpdxCzbTs00D/nvjcRzeqkG8zTJCoHOzTD666XhGztvIU2MXcdbL07h5QCeOTq45B7MZB59Iw1PvA4jIdTjLbz3+9hvApJhZZyQsPp8yYvZ6nh6/mCKPjz+e3p1bBnQiyW1vFzUJEeHirLb+RPliXstZSbN6Qmq7bQzsZmtajAOJ9v/w1jgJ8DIy/TKjFrNiax6/e2s6D32xgB6tGzDhnv7cPriLTRg1mMYZKTx3iZModwtc+84s7vroR7bmFcXbNCPBiPY8jaeBH0Uk298eCDwW5ZhGglLk8fJazkpez1lBRmoSz1zci0uy2lpJ7lpEv85NeeKkeuT62vBa9kpylu5PlBsGRDlpqOq7IjIeOM4vekBV7QCEWsj0lTt4+IsFrNpewPlHt+bPZ/ewmka1lGSXcM/J3TjnqNa/SZRf0LbajskxahCRrp7qoKprAPyTxFcV7gvQRlU3RG2hEVd27yvhb+MW88mcDbRrXI/hN/RlgMW66wRlifLP523kybGLeHSdh43JS7jr5K52WmAdJtI3jWdFxIUzWcwFtuEcwtQFGAycAjwK2KRRQ1FVRv30C38dvYjdhR5uHdiZu0+xXxZ1DRHhoqy2DD6sOXe/k83rOSsZNf8XHjv3CE7t0SLe5hlxINLVU5eISA+carY3AK2AfcBinKNYn/LXpTJqIOt27OPhLxcwbfl2jm7XiP9deKQto63jNM5I4cYjU7njrN488tVCbho+h1MOa87Q5hayqmtEnNNQ1UXAwzG0BQAReRY4ByjBOZ/jelXdLSIdcCalpX7VGap6a6yfX5cp8nj5akUJ4yZPIdnt4vFzj+Cq49vjtp3Chp/jOjVh7F39eff71bw0eTlTl3nZlLqcWwZ2IjXJ3kLrAom4RvJroKeq9gKWAQ+Wu7dSVY/2XzZhxJCpy7Yx9KWpfLHCw6k9WjD53oFc26+DTRjGASS7Xdw8oDPf3DeQo5u7eeHrZZz+4lSmLNsWb9OMg0DCTRqqOsl/ZCzADKBtPO2p7ewq8nH7B/O45p1ZuET4Y+80XrniWFo2TIu3aUaC06phPW4/Oo3hN/RFRLj2nVm88mMRv+wujLdpRjUizhEYiYmIjAY+VtX/+cNTuThvH3uBP6vqtEr63QzcDNCiRYusESNGxNy2/Px8MjMzq6VfMJ3K7geSV5SVtUt9yuS1pXyxohifCud0TuaMjskU7yuI6N8VCpH4LFH8FWo7llSXv4LpReIvj08Zv9rD6JUluEQ4r0syp7VPpqjC98n8VbUskb5fgwcPnquqvQ+4oaoRXzjnZ1wF/MXfPhToG0K/ycDCANd55XQeBr5g/8SWCjTxf84C1gMNgj0rKytLq4Ps7Oxq6xdMp7L7geQVZdnZ2Tpr9Q49/cUp2n7YGD33ufG6bkdBWPZFSiRjJ4K/wmnHkuryVzC9aPz1ydhv9Mb3Zmv7YWN0yPM5+tpnkyOyLxJqor8qyhLp+wXM0QC/U6PdEf4a4ANOBv4K5AEjgT5VdVLVIVXdF5FrgbOBU/zGo6rFQLH/81wRWQl0A+ZE+W+oM+zIL+Y/C4r5bsJ02jSqx1tXZ5G8dTHtGqfH2zSjltAs3cV/ru3NN4u38OioXP4xu4ilnh956KzDaV7fQp61gWgnjeNU9VgR+RFAVXeJSEo0A4rIUGAYMFBV95WTNwN2qnMWeSegK7AqmmfVFXw+5aPZ63hmwlLyi0r5/aDO3HlyF9JTksjZtiTe5hm1kFMOb8GJXZryp/e+YdyCzXyzeCv3ntaNQ32JGw43QiPaScMjIm7AiVU5v9ijXbj9Ck4o6mt/TaOypbUDgL+KSCngBW5V1Z1RPqvWs2aPlxdf/4Gf1u/m+E6NOad1IVcOPSzeZhl1gLRkNxd2TeEPF/Th0VG5PD56EYfWd9Go8y47J74GE+2k8TJO3qG5iDwFXAz8OZoBVbVLJfKROKEvIwT2FHp4YdJShk8vokmm8s/Ljubco1ozZcqUeJtm1DE6Ns3g/ev7MGHhZh4e+SMXvf4Dl2S1ZdgZh1n9shpItAULPxCRuThlQwQ4X1UXx8QyIyJUlS/nb+SpsUvYWVDMKYcm8fz1A2lYLznephl1GBHhjCNbIVuW8KOnJe98t5qJuZu577TuXHncoVZWvwYRacHCxuWaW4GPyt+zsFF8WJ/n43dvzWDW6p0c1a4R713fh+3Lf7QJw0gY0pKEB4ccziVZ7Xh8dC6Pjsrlo1nr+Ot5PenbsXHwAYy4E81xr4rzdnEosMv/uRGwDugYE+uMkNhT6OHFr5cxfHohDeqV8rcLjuSyPu1wuYSc5fG2zjAOpEvzTIbf0JeJuZt5YsxiLn1zOhcc04YHzziM5g1slVUiE2nBwo7w6/Guo1R1nL99BlDlclojdvh8ysh5G/jHhCXsKChhcNsknr9uEIdkRLWAzTAOCiLC0J6tGNitOa/lrODNKav4etEW7hnSlWv7dSDZQlYJSbT/VfqUTRgAqjoe5/Q+o5pZuHEPF7/xA3/87GcObZzO6DtO4pojUm3CMGoc9VLc3Hdadyb9YQB9OhzCk2MXc8Y/p/HDiu3xNs0IQLSTxnYR+bOIdBCR9iLyMLAjFoYZgdlVUML7ucWc88p3rNu5j+cuOYrPbu1HzzYN422aYURFh6YZvHt9X/5zTW+KS71c8Z+Z3P7hPKtllWBEu+T2cpzDlr7wt6f6ZUaM8fqUEbPX8ezEpewtLOX6fh2559SuNEizJLdRuxjSowUndW3KW1NX8Wr2Cr5dvJU7T+nCjSd1tPLrCUC0S253AnfHyBajEuat28WjX+WyYOMejuvYmLNbF3L1OT3ibZZhVBtpyW7uOqUrFxzThifHLuKZCUv5dM4GHj2nB4O6N4+3eXWaqCYNEcnGvxu8PKp6cjTjGg7b84v5x/glfDp3Ay0apPLy5cdwTq9WtkHPqDO0a5zOm1f3ZsqybTw2Kpfr3p3NaT1a8MjZPaxmWpyINjx1f7nPacBFQGklukaIeH3Ke9+v5vmvl1Hk8XLLwE7ceXJXMlOj/c9lGDWTgd2aMeGe/rzz3Rr+9e1yhrwwhdsGdeGWgZ1IS7aQ1cEk2vDU3Aqi70XE/gyOgpmrdvDoD4VsyF9E/65NefScI+jSvHrq6RtGTSI1yc3vB3Xm/GNa89TYxbw4eRmfzVvPX84+giGHW8jqYBFteKr8Fk4XzjkXLaOyqI6yZW8Rfx+3mC/n/0KTNOGNq7I4/YgW+Is2Gobhp1XDerxyxbFc0Xc7j47K5abhcxjcvRlDm0dbK9UIhWjjHeV3hpcCq4EbozWqLlFS6uO9H1bzz8nL8fiUu07uwhHuXzi9p829hlEV/bo0Zdzd/Xn/hzW8NHk505aVsj5pKbcN7kx6ioVyq4toPXu4qhaVF4iIla0Mke+Wb+fRUQtZua2AUw5rzl/O6UH7Jhnk5GyKt2mGUSNIdrv4v/6dOPeo1tz9Xg6vZK/g83kbeOTsHgzt2dLe1KuBaDf3/RBANj3KMWs9G3cXctsHc7nq7ZmU+pS3r+3N29f1oX2TjHibZhg1kuYN0rilVxqf3HICDeol8/sP5nHNO7NYsTU/3qbVOiKtctsSaAPUE5FjcMJTAA0AWwdXCcWlXkatLGHcNzkA3HdqN24aYKs/DCNW9O3YmDF3nsQHM9fx3KSlDH1pKjee1JE7T7HVh7EiUi+eDlwHtAVeKCfPAx6K0qZaSfaSrTw+Opc1Ozyc0bMlD591OG0PsfnVMGJNktvFtf06cFavVjwzYQlvTl3Fl/M38tCZh3PuUa0tZBUlkVa5fR94X0Qu8p+oZ1TCuh37+OuYXCYv3kqnZhnc3zuNOy7OirdZhlHraZqZyjMXH8XlfQ/lL1/lcveI+Xwwcx2Pn3sEh7dqEG/zaiyRhqeuUtX/AR1E5N6K91X1hQDd6hSFJV5en7KSN6asJMklPHjGYVx/Ykd++G5qvE0zjDrFMYcewpe3n8jHs9fz7MQlnPXyNK4+vj33nto93qbVSCINT5VlbGO+60xEHgNuArb5RQ+VO6/jQZwlvV7gLlWdGOvnR4uqMjF3C0+MWcTG3YWcd3RrHjzjcFo2tINlDCNeuF3CFccdyplHtuSFr5fx3xlrGf3zJs7rAAN8istlIatQiTQ89ab/5+OxNedXXlTV58oLRKQHcBlwBNAamCwi3VTVW002hM3Kbfk8NiqXacu3071FfUbcfDzHd2oSb7MMw/DTKD2Fv57Xk8v6HMqjoxbybu4u5r32PY+f15Oj2zWKt3k1gmh3hDfDeSvoUH4sVb0hOrMCch4wQlWLgdUisgLoSwIs8S0oLuVf367g7e9WkZbk5tFzenD18e1JspPHDCMh6dG6AZ/ccgJ//+gbvlxdxPmvfs+lvdvyp6GH0TTTtppVhageUKQ29M4iPwDTcHaG//oXfzTJcX946jpgLzAHuE9Vd4nIK8AMfy4FEXkbGK+qnwUY42bgZoAWLVpkjRgxIlJzKiU/P5+MjAxmbvby8ZISdhUrJ7VJ4pJuKTRMrfxVNz8/n8zMqqN6wXQqux9IXlEWbjuWRDK2+Sv2/gqmV9f85U7LYNRKD5PWeEhxw4VdUjj50CTc5UJW1eWvirJE8tfgwYPnqmrvA26oasQXMD/CfpOBhQGu84AWgBtn4+FTwDv+Pq8CV5Ub423gomDPysrK0urgv6O+0d+9+YO2HzZGz3p5qs5ZszOkftnZ2VHrVHY/kLyiLNx2LIlkbPNX9fSpSq+u+mv5ljy96j8ztP2wMXraC1P0hxXbQxo7Gn9VlCWSv4A5GuB3arS7XcaIyJla7pzwUFDVIaHoici/gTH+5gagXbnbbYFfwnluLNhb5OGfk5fz7g+FNKhXylMXOPFRtyXSDKNG06V5JsNv6MukRc5Clsv/PYOze7Xi4bMOj7dpCUW0k8bdwEMiUgx4cHaGq6pGvAhaRFqpalnxpQtw3kAARgEfisgLOInwrsCsiC0PE59P+eLHjfx9/BJ2FBQzsG0SL143iEMyUg6WCYZhVDMiwulHtGRgt2a8MWUlr+es5JvFWzmrg4sTTvLacbNEf55G/VgZUo5nRORonOq5a4Bb/M/KFZFPgEU4FXVv14O0cir3lz385atc5q7dxdHtGvHudX3YseJHmzAMo5aSluzmniHduOjYtjw5dhGf5W5hzotTefScIxh8WN0+uyPa1VPHBhDvAdaqakQn+Knq1VXcewonz3FQ2L2vhOcnLeODmWs5JD2FZy7uxcXHtsXlEnJWHCwrDMOIF2XHzb7y2WQ+Xytc/95shhzenEfO7lFnC4xGG556DTgWWOBvHwn8BDQRkVtVdVKU48cFn0/5eM56npmwhD2FHq45oQN/OLUbDeslx9s0wzDiQM+mSdx8/gDe/X41L3+znFNfnMotAzpx26Au8TbtoBPtpLEGuFFVc+HXDXh/BJ4APgdq5KQBMGLWOro2r8/j51mdGsMwICXJxS0DO3P+MW34+7jF/OvbFYycu4ELOioDVetMIcRoJ43DyiYMAFVdJCLHqOqqmuxAl0t49/q+HJKeXGe+CIZhhEaLBmm8dNkxXHFcex4dlcur8/cyP38mj51zBF1bVEeaN7GIdsvyUhF5XUQG+q/XgGX+0/s8MbAvbjTOSLEJwzCMSunbsTGj7ziRq3uksHDjXs745zSeGLOIvUU1+ldfUKKdNK4DVgD3AH8AVvllHmBwlGMbhmEkNEluF6ccmkz2/YO4pHc73vl+NSc/N4XP5m7AF0W1jUQm2iW3hcDz/qsids6iYRh1gsYZKfz9wiO5vG87Hh2Vy/2f/kSXRi6ad9tDzzYN421eTInqTUNEuorIZyKySERWlV2xMs4wDKMm0attI0be2o9nL+7F1n0+znnlOx76YgG7CkribVrMiDYR/i7wKPAiTjjqevafF24YhlHncLmES3q3I2PXCuYWt+C9H9Yw9udN3H96d9rUgpBVtDmNeqr6DU613LWq+hhwcvRmGYZh1GzSk4VHzu7B+Lv706NVAx75ciGP/VDEnDU7421aVEQ7aRSJiAtYLiJ3iMgFQN3eY28YhlGObi3q8+FNx/HqFceS71EufmM69348n617i+JtWkREG566B0gH7sLZ0HcycG20RhmGYdQmRISzerUiaesSFnhb89bUVUxatIW7TulCJ1/NCllFu3pqtv9jPk4+wzAMw6iE1CTh/iHduTirLU+MWcTfxi2hVYaQ2m4b/bs2i7d5IRHRpCEio6q6r6rnRmaOYRhG7adD0wzevq4P3y7ZwgOfzOXqt2cx9IiWDGnqi7dpQYn0TeMEYD3wETATWzFlGIYRNicf1oInT6zHclc7Xvl2Bd8s9rIxeTm3DOxEWnJint0RaSK8JfAQ0BP4J3AqsF1Vp6jqlFgZZxiGUdtJcQu3D+7CN/cN5Jjmbl6cvIwhL0xhUu7msqOtE4qIJg1V9arqBFW9Fjgep5RIjojcGVPrDMMw6gitG9XjtqPT+PCm40hPcXPzf+fy/NxiVm5LrOIaES+5FZFUEbkQ+B9wO/AyTjl0wzAMI0L6dW7K2Lv685eze7Byt5ehL03l6fFLKCiO6Fy7mBNpIvx9nNDUeOBxVV0YpIthGIYRIsluFzec1JEmBWv4bm9j3piyki9+3MBDZx7OuUe1jmsF7kjfNK4GugF3Az+IyF7/lScie2NnnmEYRt2lYarw7CVH8flt/WjRII27R8znd2/NYPGm+P2ajTSn4VLV+v6rQbmrvqpGdcydiHwsIvP91xoRme+XdxCRwnL33ojmOYZhGDWFYw89hC9vO5GnLzyS5VvyOOvlaTz61UL27Dv4Z3dEuyM85qjq78o+i8jzwJ5yt1eq6tEH3yrDMIz44nIJl/U9lKE9W/LC18v474y1jP55E386vTuX9m6Hy3VwQlbR1p6qNsQJ2l2KsxfEMAzDABqlp/DX83oy5s7+dGmWyQOfL+CC175n/vrdB+X5kojrgAFEZADwgqr29rc7ALnAMmAv8GdVnVZJ35uBmwFatGiRNWLEiJjbl5+fT2ZmZrX0C6ZT2f1A8oqycNuxJJKxzV+x91cwPfNXeHrR+KuiLFx/qSozNnn5eGkJu4uV/m2SuKRbCg1Sg791BBt78ODBc8t+/x7w0IN9AZOBhQGu88rpvA7cV66dCjTxf87C2ZHeINizsrKytDrIzs6utn7BdCq7H0heURZuO5ZEMrb5q3r6VKVn/gpPLxp/VZRF6q+8Io/+bdwi7fLQWO356AR957tV6in1Vtkn2NjAHA3wOzUuOQ1VHVLVfRFJAi7EmRzK+hQDxf7Pc0VkJc4KrjnVaKphGEbCk5maxINnHM6lvdvx2KhcHh+9iBGz1jP8xr60aJAW02clXCLczxBgiapuKBOISDNgp6p6RaQT0BWwo2UNwzD8dG6WyfAb+jJp0Ra+mr+RZpmpMX9Gok4al3FgAnwA8FcRKQW8wK2qWrOPwDIMw4gxIsLpR7Tk9CNaVsv4CTlpqOp1AWQjgZEH3xrDMAyjjIRdcmsYhmEkHgm75DZWiMg2YG01DN2Q3248jGW/YDqV3Q8krygL1m4KbA9iX6RE4jPzV/X0qUrP/BWeXjT+qihLFH9tB1DVoQfcCbSkyq6Qlg2/VV39gulUdj+QvKIshHbAZXbx8pn5q3r6VKVn/jp4/qooS3R/qaqFp6JgdDX2C6ZT2f1A8oqyYO3qJJJnmb+qp09Veuav8PSi8VdFWaL7q/aHp4zwEJE5GmgXqBEQ81d4mL/CIxH9ZW8aRkXeircBNQzzV3iYv8Ij4fxlbxqGYRhGyNibhmEYhhEyNmkYhmEYIWOThmEYhhEyNmkYhmEYIZOQtaeMxEBEBgFP4Bx+NUJVc+JqUIIjIkv8xgoAACAASURBVC4cfzXA2ZT1fpxNSmhEpD9wJc7voR6q2i/OJiU8InIo8ArOju1lqvr0wbbB3jTqGCLyjohsFZGFFeRDRWSpiKwQkQf8YgXygTRgQ8Wx6gJh+us8oA3gwfwV1F+qOk1VbwXGAHV2gg3zO9YNGKuqNwA9DrqxYGVE6tqFU2L+WGBhOZkbWAl0AlKAn/xfSJf/fgvgg3jbXgP89QBwi1/ns3jbnuj+Knf/E0I4hbO2XmF+x5oA2cC3wPXxsNfeNOoYqjoVqHgOSV9ghaquUtUSYATO0bs+//1dOMft1jnC8RfO28Uuv4734FmZOITpr7Jwyx5V3XtwLU0cwvTZ9cCjqnoycNbBtdTBJg0DnJDK+nLtDUAbEblQRN4E/osTRzUcAvoL+Bw4XUT+BUyNh2EJSmX+ArgRePegW5T4VOazCcBdIvIGsCYOdlki3ABAAshUVT/H+UVo/JbK/LUP55eg8VsC+gtAVR89yLbUFCr7ji0ELj7YxpTH3jQMcP6KaVeu3Rb4JU621ATMX+Fh/gqfhPWZTRoGwGygq4h0FJEUnDPaR8XZpkTG/BUe5q/wSVif2aRRxxCRj4DpQHcR2SAiN6pqKXAHMBFYDHyiqrnxtDNRMH+Fh/krfGqaz6zKrWEYhhEytTYRLiLnAOekp6ff1LZt25iP7/P5cLnCf1ELpV8wncruB5JXlIXbjiWRjG3+ir2/gumZv8LTi8ZfFWWJ4q+GDRsye/bsiZooZ4QDl+CUpvABvcvJk3F2hi7AeSV7sNy9LL98BfAy/rekYFdWVpZWB9nZ2dXWL5hOZfcDySvKwm3HkkjGNn9VT5+q9Mxf4elF46+KskTyF5WcTx6vnMZC4EIOXMt+CZCqqkfiTBK3iEgH/73XgZuBrv7rwBnQMAzDqFbiEp5S1cUAIgcsRVYgQ0SSgHpACbBXRFrhlBmY7u83HDgfGF9dNm7PL0aA5CQXKW7ncrkCLZ02DMOoO8Q1ES4iOcD9qjrH307G2X18CpAO/EFV3xKR3sDTqjrEr9cfGKaqZ1cy7s04byW0aNEia8SIEWHb9tB3+/gl/7e+cQkkCbhdkCRKstuFWyDZBW6XkOTyfxZI8rcdmZDqhhQ3iNdDZloKKe4ymZDiZv9nF3iLC2lUP520JKFeErgqTK75+flkZmYeYHMgeUVZuO1YEsnYofQJpmP+Ck/P/BWeXjT+qihLJH8NHjx4rqr2riivtjcNEZkMtAxw62FV/aqSbn1xava0Bg4BpvnHqXRHaSBU9S38B7L37t1bBw0aFIblDg8e8gu79pVQUuqjxOvDU6p4vD48Xh/FpT7Wrt9AsxYt8Xh1v47/ctpKcamPPK+P4hIvhSU+ijxeCooFxRPk6QIU/trKTE2iflrZlYynIIn2bRr+KmuQlkyDesls3ryc/j2OpElGCo0zUmiUnsK0qVMo/+/PyckJqx1LIhk7lD7BdCq7H0gern9qor+C6Zm/wtOLxl8VZYnuL6jGSaPsrSBMrgAmqKoH2Coi3wO9gWk4OyLLqPbdkecc1brK+zk52xg06Kiwx83OzqZf/wEUlfgo9HjZV1JKocdLkceZWAo9XubOX0CHLt0oKPGyt9BDXlEpeUX+n8UedpYoCzfuIa/Iw96iUkpKfb+O/+r8Gb9+dglkJEGLeVNonJFCs8xUPHuLWeFeRcuGabRqmMaOQmeiS3bblh3DMIKTaEtu1wEni8j/cMJTxwMvqeomEckTkeOBmcA1wL/iaGfEiAipSW5Sk9w0JDmgTvLWxQzqe2ilY1T8C6HI40wuE6f8QKfDe7GjoISd+cXsLChhwfI1pDXMZEdBCYs37WXDrlImrV38m/HunzqeZpmptDmkHvVKi/nZu5z2TdLp2DSD9k0yaFgvsJ2GYdQ94pLTEJELcH7pNwN2A/NV9XQRycSpeNkDJ0bzrqo+6+/TG3gPJ0E+HrhTqzC+bJ9GmzZtbvrf//4X839DpLHGeMfo8/LycaVlsKtI2VnkY9OeIvZpCjuLlG2FPjbne9ld8ttoYP1kaFPfRbv6LtrWd3FofRdtMl2kuMNbGGA5jfCoiTF681fV8tqQ04j7ASTVfdk+jeyw24Ulpbpk016dsHCTvpGzQod99pOe/+p3evgj47X9sDHaftgY7fjAGD31hRx9YORP+snsdbpya576fL6o/l2R9om3v6qLmrjvwPxVtbw27NNItPCUkQCkJbvp3rI+3VvW/43c51PW7dzHks17WbQpjwUbdjP25018NMsp+39IejJ9OzZmYLfmDOzejDaN6sXDfMMwqhGbNIyQcbmEDk0z6NA0g6E9WwHORLJqez5z1+5izppdfL9iOxNztwDQpXkmg7s34+xerenVtmE8TTcMI0bEK6fxLHAOzua9lThn3e4WkVOBp3HOxC0B/qiq3/r75ACt2L8W9TRV3VrFMyynEUBW3TFnVWVTgfLzNi8LtpeyZKcPr0LzdOGYJj4Gd0inZUboK7Vqu7+qoibG6M1fVcstpxF57anTgCT/538A//B/PgZo7f/cE9hYrk8O5epUhXpZTiM7qna07C4o0Y9nrdMr/z1DO/jzIVf+e4ZOWLhJPaXeoP3rmr+iHTveMXrzV9Vyy2lEiKpOKtecgf/4QlX9sZw8F0gTkVRVLT6Y9hmxo2F6Mpf2acelfdrx5YRv2ZDSjg9mruOW/87l0Mbp3D64Mxce29b2iRhGDSHu52mIyGjgY1X9XwX5xcCtur90SA7QBGfH+EjgSa3E+FiUEQmGhafCp2xsr0/5cauXMas8rNnro2k94cKuKZzQyn1APTLzV80Kt5i/qpZbeKrqENRknGq2Fa/zyuk8DHxBhTLnwBE4uY7O5WRt/D/rA5OAa0Kxw8JT2VG1Y0nFsX0+n367eIue/fI0bT9sjF702ve6YMPusO2pK/6KZR8LT4XXx8JTByE8pUHKiIjItcDZwCl+A8vkbXEmkmtUdWW58Tb6f+aJyIc4daqGV4ftxsFBRBh8WHMGdmvGZ3M38I8JSzjnle+4ZUBn/nBqV1KT3PE20TCMCsQlkCwiQ4FhwLmquq+cvBEwFufwpe/LyZNEpKn/czLOZLPw4FptVBcul3Bpn3Z8e/8gfte7HW9MWcl5r3zP0s158TbNMIwKxGvJ7QogFdjhF81Q1VtF5M/Ag8DycuqnAQU4BzYlA26c0Ne9quqt4hm25DaArCbEnOdvLeWdhSUUeZWruigDOpq/Yt3Hchrh9bGcRjkCxaxq02U5jeyo2rEknLG37CnUC1/7XtsPG6N/H7dYvd7KS5SYv8LvYzmN8PpYTmP/ZescjYSkeYM0PrrpeAa1S+KNKSv5wyfzf1MC3jCM+FBlItyflL4M6I9zMFIhTi5hLDBeVe3/YqPaSElycW2PFPoe0ZlnJixl9z4Pr191LOkpVv3GMOJFpTkNEXkXaAOMAeYAW4E0oBswGMgCHlDVqWE/tPIyIlcCfyyn2gs4VlXni0gW+0ujjwPu1sqMx3Ialclqasw5Z72H93NL6N7YxR+y0kgtV5bd/BV+H8tphNfHchrlCBSz8v8u7lnZPf/9FKBLVTpV9A1YRqSCzpHAqnLtWcAJOOdsjAfOCOVZltPIjqodS6KNOX8xb4N2eGCMXvWfGVrkKQ153Lrqr0j1zF/h6VlOY/9kUuWSVlUtUdUVVelU0XeSqpb6mzP47VGuZVwOfAQgIq2ABqo63f+PGQ6cH8mzjZrL+ce04R8X9mLa8u3c/sE8PF6LjhrGwSZocFhEFgAVw0B7cEJWT6rqjgN7hcUNwMcB5L8DzvN/bgNsKHdvg19m1DEu7dOO4lIvj3yVywMjF/DcJb3ibZJh1CmC7tMQkWdw6j196Bdd5v+5FzhJVc+ppN9koGWAWw+r6ld+nYeB3sCFWs4QETkO+I+qHulv9wH+rvvrUPUH/lTFs632VCWy2hJz/nJFCV+u8HBu52ROa1Vi/gqzj+U0wutjOY1yBIpZ6W9zC99XJgMWBOtfxbjXAtOB9AD3XgQeKtduBSwp174ceDOU51hOIzuqdiyJZczZ5/Ppnz79SdsPG6OPDp8U0Rh1yV/h6Jm/wtOznMaBZPr/8gdARPoCZdNTaeAuVVNZGRH/PRdwCfDr64GqbgLyROR4ccqgXgN8FcmzjdqBiPDkBT0Z1L0Z7+eW8M3iLfE2yTDqBKGEp/oA77B/osgD/g/nvIuzVPWTsB9aSRkR/71BwNOqenyFPr3Zv+R2PHCnVmG8LbkNLKtt4YOiUuWpGQVsKRQe6JNGp0YHFjk0f4WnZ/4KT8/CU5WHkxoCjULVT5TLwlPZUbVjSXWFD76c8K2e9I9v9Ji/TtKVW/NCHqOu+iuYnvkrPD0LT1VARFqIyNvACHU24PUQkRsjnt4MI8Y0TBWG3+BEUK95ZxZb84ribJFh1F5CyWm8B0zEKSMCsAy4p7oMMoxI6Ng0g3ev68OO/BKue2c2eUWeeJtkGLWSUHIas1W1j4j8qKrH+GXzVfXog2JhhFhOI7CsNsacy+v8vK2Uf84rpntjF/dmpZHkEvNXmHrmr/D0LKdxYC4jB+ds7nn+9vHAlGD9goz5BPAzMB/n6NbWfnkTIBvIB14JYMdSf5/5QPNQnmU5jeyo2rGkumLOFXU+m7Ne2w8bo7d9MFc9pV7zV5h65q/w9OpaTiOUcqH3AqOAziLyPdAMuDiEflXxrKo+AiAidwF/AW4FioBHgJ7+qyJXquqcKJ9t1HIuymrLzoISnhq3GLcI57U4+AeNGUZtJeikoarzRGQg0B2nWOBSVY0qYKyqe8s1M/CXKVHVAuA7EekSzfiGcdOATpT6lH9MWMLWVm4GDvSR5LbjYwwjWqoqjX5hVR1V9fOoHizyFM4mvT3AYFXdVu7edUBvVb2jnCwHJ3zlBUbi1L0KaLyVEalcVhtjzlXpjF1VwqfLPGS1cHNLr1RSgpRUrwv+CqZn36/w9CynsT+H8K7/GgvswvlFPRLYCXxeWb9y/SfjHNhU8Tqvgt6DwOMVZNdxYE6jjf9nfZw8yDXBbFDLadSJmHMwnQffnaQdHhijF7/+ve4qKK6yX13wVzA9+36Fp1fXchpVlUa/XlWvxwkd9VDVi1T1IuCIUGYxVR2iqj0DXBXLf3wIXBTCeBv9P/P8ffqGYodhnNYhmX9dfgw/rd/DxW9MZ92OfcE7GYYRkFCCvB3Uqf1Uxhac0/siRkS6lmueCywJop8kIk39n5OBs3HeWgwjJM7u1ZrhN/ZlW14x57zyHdlLt8bbJMOokYSyT+MVoCvOgUiKUxp9hareGfFDRUbiJNZ9wFrg1rI3CRFZAzTAORlwN84pf2uBqUAy4MYJfd2rqt4qnmH7NALIamPMORx/bd3n418/FrMhz8eZhyoXHZ6BSyrPc9RGfwXTs+9XeHqW0wicn7gAp1z5i8AFofRJlMtyGtlRtWNJvHIaFe/vKy7Vuz+ap+2HjdFLXv9B1+0oqFS3NvormJ59v8LTq2s5jUqX3IqI+Duiql8AX1SlYxg1hXopbl783dE0827no2V7OeOf03j83CO48Fg7DNIwglFVTiNbRO4UkUPLC0UkRUROFpH3cQ5SMowah4hwYptkxt/dn8Nb1ee+T3/ixvfnsL3Qzh03jKqoap9GGs753VcCHXHyC/VwJppJwKuqOj+ih4o8gXP+tw/YClynqr/47/UC3sTJa/iAPqpaJCJZ7D9PYxxwd1VvOZbTCCyrjTHnaP3lU2XSmlI+X1ECqlzQNZXT2ifhDlC3qjb4K5iefb/C07OcRuCcRjLOkasxOU8DaFDu813AG/7PSTg1qY7yt5sAbv/nWcAJOLvSxwNnhPIsy2lkR9WOJYmS06hMvn5ngZ773HhtP2yMnvHSVJ29eket9FcwPft+hadX13IaIdVVUFWPqm5S1d2h6IcwXsAyIjgrpX5W1Z/8ejtU1SsirXAmmun+f8xw4PxY2GIYZbQ9JJ17jk3ljauOZWdBCRe/MZ03fipi057CeJtmGAlD0CW31fbgAGVEROQeIAtojlMYcYSqPuM/6vVpVR3i79sfGKaqZ1cytpURqURWG8MH1eGv4lJlzGoP41eV4HIJZ3dKZmiHZEoKC2q8v4Lp2fcrPD0LT8XoIoIyIsD9wGqgKZAOTAdOAfoAk8v16Q+MDsUOC09lR9WOJYkengok+2TsN3rrf+do+2Fj9MSnv9FnPvpafT5fWPZFSk0Mt9TG71cwPQtPBUBE2otI2V/59USkfrA+GlkZkQ04Z3VsV9V9OAnvY/3ytuX6tAV+CcV2w4iGZukuXr8qiw9vOo7M1CRenV/M5f+eweJNe4N3NoxaSChnhN8EfIazogmcX9hfRvPQKsqITAR6iUi6iCQBA4FF6pQxyROR40VEcMJaFScfw6g2+nVuypg7T+KaHiks2ZzHWS9P4/HRuRSW2jYlo24RShmR+TjFAWfq/uNeF6jqkRE/tOoyIlfhhKwUGKeqf/LLe7N/ye144E6twnhbchtYVhtjzgfbX6Rk8PnyErLXl1I/Wbn88DSOb+VGypUjiQU1MUZfG79fwfQsp3FgbmKm/+ePWm5ZbLB+iXJZTiM7qnYsqYk5jara89ft0kFPjdP2w8boZW9O12Wb9wa1NRxqYoy+Nn6/gulZTuNApojIQ0A9ETkV+BQYHe6sZhi1jaPaNeKRE9J48vyeLNq0lzNfnsa/vlmOx2u7yo3aSyiTxgPANmABcAtOcvrP1WmUYdQUXCJcdXx7vr1vIEN7tuL5r5dx4Ws/sHRzXrxNM4xqIax9GiLSGGirqj/H5OEi9wPPAs1Udbtf9iBwI86xrnep6kS/PAdnV3rZTqvTVLXSQxEspxFYVhtjzonkr9mbSxmeW0xhKZzfNZkzOyb/pvR6ONTEGH1t/H4F07OcxoE5jRycOlCNgXXAXOCFYP1CGLcdzmqptUBTv6wH8BOQilPvaiX7y4jk4JwbbjkNizmHpXOw/bUtr0h//z9nb8cV/56uW/cWVWlfZdTEGH1t/H4F07OcxoE0VKfsx4XAu6qaBQwJoV8wXgT+xP4SIuAUMRyhqsWquhpYgR3ratQwmmam8uoVx/LMRb2Ys2YXZ748jR9WbI+3WYYRE0KZNJL8tZ8uBcbE4qEici6wUf01psrRBlhfrr3BLyvjXRGZLyKPSKzXNxpGDBERLu3Tjq/uOJEGaUlc+fZMXstZUfaWbRg1llD2aVwCPAJ8p6q3iUgn4FlVvShIv8lAywC3HgYewslJ7PEf79pbVbeLyKvAdFX9n3+Mt3H2aowUkTaqutG/G30k8D9VHV7Js632VCWy2hhzTnR/FZUq7y4sZuZmLye0cnN9z1RS3MH/5qmJMfra+P0Kpmc5jWq+gCNxztBY479KcXIlLXE29T1YTncicEKAMa4DXgnleZbTyI6qHUvqSk4jED6fT1/5drm2HzZGz33lO92ypzBon5oYo6+N369gepbTqICIpInI7SLymoi8U3aFPp8dMEktUNXmqtpBVTvghKCOVdXNwCjgMhFJFZGOQFdglogkiUhTvz3JwNk4xQ8No0YgItw+uAtvXJXF8i15nPvK9yzZbPWrjJpHKOGpT3FqQ10B/BXnJL/Fqnp3TAwoF57ytx/GOTGwFLhHVceLSAYwFecwKDdOBd17VdVbxbi25DaArDaGD2qav9bt9fLC3GKKvco9x6bRvbE7LLurIt7hltr4/QqmZ+GpA0NBZeVDfvb/TAa+DdYvUS4LT2VH1Y4ldTk8VZENu/bpyc9la9eHx+m4n38J2b5gxDvcUhu/X8H0LDx1IB7/z90i0hNoCHQIY0IzDKMCbRrV47Nb+9GzdQNu+3Ae/52xNt4mGUZIhDJpvCUih+CsoBoFLAKeqVarDKMOcEhGCh/83/GcclhzHvlyIc9PWmpLco2EJ27HvVY3ltMILKuNMeea7i+vT3l/UQlTN5QyoG0S1/ZIwe2SGhmjr43fr2B6ltM4MKeRipMEfwj4S9kVrF8oF87xrsr+MiIdcGpLzfdfb5TTzcIpmrgCeBn/hBfsspxGdlTtWGI5jcrx+Xz63MQl2n7YGL3xvVm6r7i0Rsboa+P3K5heXctpJIUwIX0F7MGpOVUcgn5IiEg74FScPRrlWamqRwfo8jrOhr0ZOJV2h+IcxmQYNR4R4b7TutO8fip/GZXLlf+ZwQ1damcUwKjZhDJptFXVodXw7LLaU0GPbfWXMWmgqtP97eHA+dikYdQyrj6hA00zU7l7xHz+tkM5pm8hbRrVi7dZhvEroezTeAv4l6ouiNlDndpTp6jq3RXKiHQAcoFlwF7gz6o6zX/U69OqOsTfvz8wTFXPrmR8KyNSiaw2xpxro7+W7PTy0txC6iW5uK93Gm3rh7JmJf4x+tr4/QqmZzmN/TmEBcDPOKulPMBSf3sBIRz3irMBb2GA6zxgJk71XHBKiZTlNFKBJro/h7Eepyx7H2ByubH7A6OD2aCW06gTMefa6q/ho77RPk9+rUc+OkG/W74tpD7xjtHXxu9XMD3Laewn4F/xoaL+t4KKiMiROGdl/OQvVNsWmCcifdUpJVLs7z9XRFYC3XBKjbQtN0xb4Jdo7DOMRKddfRef33Y81787m6vfnslDZx7OjSd1xAo8G/GkqnfeLcAFwB9xks4bVXVt2RXpA7WK2lMi0kxE3AD+arpdgVWqugnIE5Hj/SXRryGEXIhh1HTaHpLOF7efyKk9WvDk2MXc+8lPFHkqrZ5jGNVOpTkNEfkYJyw1DTgDWKsxqjdV4Tlr2J/TuAinvlUpznGvj6rqaL9eb+A9oB5OAvxOrcx4bJ9GZbLaGHOuC/7yqTJmlYcvlntoW9/Frb1SaRMgzxHvGH2i+CvWfSynUY5AMSv/7+IF5T4nAfMq003ky3Ia2VG1Y4nlNMIj0NjfLtmiWU9M0m4Pj9N3v1ulPp8vInsspxFeH8tphFZ7qqzmFKpaGuFkZhhGDBncvTnj7x5Av85NeGz0Iq55ZxZrthfE2yyjDlHVpHGUiOz1X3lAr7LPImIHARhGnGhWP5V3ruvDk+f3ZP663Zz20lT+OXk5xaWW6zCqn7jWnhKR+4FngWbq5DRSgDeB3oAPuFtVc/y6OUArnDIj4BwXu7WKsS2nEUBWG2POddlfu4t8fLSkhJmbvTSrJ5zZzsfAjhm4gqywspxGeH0sp1GOQDGrg3EB7XCOc13L/n0atwPv+j83xyld4vK3c3AS5pbTsJhzWDp1wV9Tl23VoS9N1fbDxujQl6bq17mb1ev1VapvOY3w+lhOI7ScRnVTVkak/KtOD+AbAHXeInbjvHUYhlEF/bs2Y+ydJ3FLr1QKikv5v+FzOP2lqXwye72FrYyYEpfwVBVlRG7GKWJ4Oc6byI/Ajao60h+eaoKzFHck8KRWYryVEalcVhvDB+av3/ZJS89g1mYv41d7WJ/no0GK0L9NEv3bJtEywxV07LrmLwtPBeagh6eIrIxIEs4byHyczXvjgPP899r4f9YHJgHXhGKHhaeyo2rHEgtPhUe0/vL5fPrd8m1643uztNODY7X9sDF60Wvf64cz1+pXE78N+7m13V+R6tW18FQoVW4jQiMvI/KHcro/AMv94230/8wTkQ+BvsDw6rLfMGo6IsKJXZpyYpembM0r4ot5G/l4znoe/HwBAnywejpDe7bk1B4taHtIerzNNWoI1TZpVIY61XKbl7UrhKfScUJmBSJyKlCqqotEJAlo5NdJxqmLNflg224YNZXm9dO4ZWBnbh7QicWb8nhj7AyW5Jfw+OhFPD56ER2apHNil6ac1KUp3hI7x8OonLgf9xqgNPpEnOW2G3HyGWtFJAOYCiQDbpwJ415VrTTDZ0tuA8tqY8zZ/BV+nzK9Tfk+Fmz3krvDy9KdXoq8IChtMl10auSmc0MXnRu5aZ0p7CsoqPP+Cuee5TRq6GU5jeyo2rHEchrhcbBj9CWlXp2zZofe8++Jet07M/Woxydq+2FjtP2wMdrjkfF6yt/H6YOf/6zv/7BaZ67aobv3ldRpfwW7ZzkNwzBqNcluF1ntG5PXJYVBg/qiqqzZsY8f1+3ip/W7mbFkPWN++oUPZ+6vKnRIqtB92XQ6Ns2gfZMM8jeX0mLTXjo0yaBeijuO/xqjurBJwzCMgIgIHZtm0LFpBhce25achtsZOHAgm/cWsWRzHks25fHdzyso8iqTcrewo6AEgFfmTwOgcUYK9d2ldF07m1YN69GyYRp7NpWSunIHLRum0SQzhfqpSXY+SA0jXvs0HgNuArb5RQ+p6rhy9w/FOTHwMVV9zi/LYn9p9HE4JUasNHoQeV2IOZu/wu9THTH6fR5l7Y4C8jSNzft87CxStuV72FvqZmeRjwLPAUOSJNAgVaif4lwNUoT6KdAgRchMEdKThIxkIT0J0pOddr0kcLsOnGhqmr8CyRLp+1VZTiOek0Z+2YQQ4P5InGT4zHKTxizgbmAGzqTxsqqOD/as3r1765w5c2Jl+q/k5OQwaNCgaukXTKey+4HkFWXhtmNJJGObv8IbO9Q+VelVl7/2lZQy6uuptOvei817ithZUML2gmJ25pewo6CEHfnF7CgoYXt+MUUeX5X2Z6S4aVAvmQZpydRPSyI9NYmCPTtp37ol9VLcpKe4qZeSRHqKm4xyn+uluElPdn6mJLmYP3cO/U88gRS3i5QkF6lJLlLcLlwVJqXq8ldFWSJ9v0Qk4KSRcOEpETkfWAUUlJO1Ahqo6nR/ezhwPs5hTIZh1ADSU5wd6Sd2aRpUd19JKTsLSsgrKmVvoYe9/p95Rfs/7y3ysLew1P/Tw7YCH1vW7KSwxMu+Ei+FoZ5wOO3bA0TJbiE1yZlYUtwufKXFNJw3xWknuUh2uUhyC26XsGd3ER+sm0OSy2knuYQkt4utW4r5eteCX9tJLmHjhhLmeZb9RnfNGg+rv1+N2yWsWOdh48y1uERw7HQ0cwAACJlJREFUi7B0g4cdczfgdgki4BJx7rmc8KHTBpdr/2e3iP8e9O7QOOBbWTTEc9K4Q0SuAeYA96nqLv/S2mE4pUTuL6fbBudY2DI2+GWGYdRC0lOSSE8J79dTxb+cfT6l0OOfQEq87POUsq/Ey75iL0UeLyVeH/MX5NK5azdKSn0U+6+SUh8lXh/FHh8lXi8lpT7WbdzEIU0yf9Ur9fnweH0UepS9xYpnVyFen49Sn1LqVbw+paDQS+6uzZT6nHapz4en1Id31fIDjV+yaP/nRQt/e2/hT2H54TfDPjEUtyu2CxKqLTwlIpOBlgFuPYwTYtqOU6zwCaCVqt4gIs8Bs1T1k/IhLBHpA/xd/bvMRaQ/8CdVPaeSZ1vtqUpktTFGb/4Kv4/VngqvTyz9lZGRgQJeBa8P9uYXkJ6egQ/Izy8gPT0dH6AK+QX7qJeejk+dtoLzGecIYFV+1S3/uUznsMauSsvk19h9GkAHYKH/8zScWlRrcCrc7gTuwDlHY0m5PpcDb4Yyvu3TyI6qHUtsn0Z41MR9B+avquW2TyNCRKSVqm7yNy/AKWSIqvYvp/MYzpvGK/52nogcj1Ps8BrgXwfVaMMwDCNuq6f+CxyN8wa1Bril3CRSpvMY5VZYiUhv9i+5HQ/cqVUYb0tuA8tqY/jA/BV+HwtPhdfHltzuJ+61p6obEdmGczpgrGkI7KmmfsF0KrsfSF5RFqzdFCffVB1E4jPzV/X0qUrP/BWeXjT+qihLFH9tB1DVoQfcCRSzsiukXMxb1dUvmE5l9wPJK8pCaAeMY8bLZ+av6ulTlZ756+D5q6Is0f2lGt/jXms6o6uxXzCdyu4HkleUBWtXJ5E8y/xVPX2q0jN/hacXjb8qyhLdX7U/PGWEh4jM0UDL7IyAmL/Cw/wVHonoL3vTMCryVrwNqGGYv8LD/BUeCecve9MwDMMwQsbeNAzDMIyQsUnDMAzDCBmbNAzDMIyQSbjS6EbiICKDcApK5gIjVDUnrgYlOCLiwvFXA5z19e/H2aSExl949Eqc30M9VLVfnE1KePwH1L2Cs/lumao+fbBtsDeNOoaIvCMiW0VkYQX5UBFZKiIrROQBv1iBfCCN35amrzOE6a/zcEr2ezB/BfWXqk5T1VuBMUCdnWDD/I51A8aq6g1Aj4NuLNiO8Lp2AQOAY/FXFvbL3MBKoBOQAvzk/0K6/PdbAB/E2/Ya4K8HcOqoAXwWb9sT3V/l7v9/e3cXImUVx3H8+yPKJDcvjCJNMvClCExJuylCSgyKXjHtosIXIooKhfQiCQuCLvIqSxM2KiLULRS8SgszxSwXt1UMDSSVRCEvzHcS9d/FOZvPDrPrM2qzs87vAwPPnGee55znz+z85zxn55w20iJrfd7+Ro8ZMAT4AVgPzOyL9rqn0WQiYiNpyvmi+4A9EfFHRJwBVgBPRkTXmptHgAF1bGbDqCVepN7FkfyaksvGXV1qjFfX7ZajEXGsvi1tHDXGbCawMCIeAh6rb0sTJw2DdEvlz8LzA8AwSc9IWgZ8SbqPaknVeAGrgEckLQY29kXDGlRP8QKYDXxW9xY1vp5i9i3whqRPSDOE150Hwg2g2tJeERGrSB+E1l1P8TpF+hC07qrGCyAiFta5Lf1FT++xncDUejemyD0Ng/QtZnjh+W3AwT5qS3/geNXG8apdw8bMScMA2oFRku6QdB3wHLCmj9vUyByv2jhetWvYmDlpNBlJy4EtwBhJByTNjoizpLXY1wK7gLaI+K0v29koHK/aOF61628x84SFZmZWmnsaZmZWmpOGmZmV5qRhZmalOWmYmVlpThpmZlaak4aZmZXmpGFNQ9I5SZ2Fx4i+btOVImm8pNa8PUPSRxX7N0ia0MvxKySN+r/baf2f556yZnI6IsZV2yFJpN8tna+2vx94C3jvMo5fCswHXroyzbGrlXsa1rQkjZC0S9ISoAMYLmmepHZJOyS9W3jtgrwgzveSlkt6M5f/9w1e0k2S9uXtayR9UDjXy7l8Uj7mG0m7JX2VExaSJkr6SdJ2SVsltUjaJGlcoR2bJY2tuI4WYGxEbC9xzU8Uelq/S9qbd20CJkvyF0nrld8g1kwGSurM23uBucAY0mI2r0qaAowirWUgYI2kB4GTpLl/xpP+ZjqAbRepazZpnYiJkgYAmyWty/vGA3eTJqDbDNwvaSuwEpgeEe2SbgROA63ADGCOpNHAgIjYUVHXBGBnRdl0SQ8Uno8EiIg15DmMJLUBP+by85L2APeUuDZrYk4a1ky63Z7KYxr7I+LnXDQlP37NzweRkkgLsDpPfY6kMhPHTQHGSuqaxnpwPtcZYGtEHMjn6gRGAEeBQxHRDtC1KJGkr4G3Jc0DZgGfV6nrVuBwRdnKiHitcK0bijslzSfF4+NC8V/AUJw0rBdOGtbsTha2BbwfEcuKL5A0h7z+QxVnuXCb9/qKc70eEWsrzjUJ+KdQdI70d6hqdUTEKUnfkVZtm0bqVVQ6XVF3ryQ9DDxLWma06Pp8LrMeeUzD7IK1wCxJgwAkDZN0M2kVvqclDczjB48XjtkH3Ju3p1ac6xVJ1+ZzjZZ0Qy917waGSpqYX99SGF9oBT4E2iOicllQSLOgjixzgZJuB5YA0yKiMkGMBhpiJlVrXO5pmGURsU7SXcCWPDZ9Ang+IjokrQQ6gf2kQeMui4A2SS8A6wvlraTbTh15oPsw8FQvdZ+RNB1YLGkg6Rv/ZOBERGyTdIwelkWNiN2SBktqiYjjF7nMGcAQYHW+xoMR8aikW0i3qw5d5Hhrcp4a3axGkt4hfZgvqlN9Q4ENwJ09/UuwpLnA8YhovcQ65gLHIuLTS26oNQXfnjJrYJJeBH4BFlzkNyRL6T5WUqu/gS8u43hrEu5pmJlZae5pmJlZaU4aZmZWmpOGmZmV5qRhZmalOWmYmVlpThpmZlbav3DEiT9v/RiSAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "G = GH*GL\n", "print(G)\n", "matlab.bode(G, Hz=True)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## ネットワーク図\n", "AnalogDiscoveryを使って、ネットワーク解析をしました。\n", "\n", "\n", "\n", "AnalogDiscoveryとの接続との接続は以下の通りです。\n", "\n", "| AnalogDiscovery | RFフロントエンド |\n", "| :-- | :-- |\n", "| Scope Channel 1 Posigive(Orange) | A段12列(A1入力) |\n", "| Scope Channel 1 Negative(Orange/White) | A段12列(GND) |\n", "| Waveform Generator 1(Yellow) | L段6列(Antena) |\n", "| Ground(Black) | A段12列(GND) |\n", "\n", "ネットワーク解析の結果は、以下の通りです。周波数は100kHzから5MHzまでとしました。\n", "\n", "\n", "\n", "トラ技コラム2図AからAM受信機用アンプとバンドパスフィルタの総合周波数特性を引用します。\n", "\n", "トラ技のコラム2の値より3〜5dBほど低い値ですが、形状はよく捉えていると思います。\n", "\n", "\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.8.8" } }, "nbformat": 4, "nbformat_minor": 2 }