531 lines
1.5 MiB
Plaintext
531 lines
1.5 MiB
Plaintext
|
|
{
|
||
|
|
"cells": [
|
||
|
|
{
|
||
|
|
"cell_type": "markdown",
|
||
|
|
"id": "7d7ec9b3-a8f6-4e85-a472-8d5552c53eed",
|
||
|
|
"metadata": {},
|
||
|
|
"source": [
|
||
|
|
"数据处理笔记与记录"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 1,
|
||
|
|
"id": "21bdcaa3-3761-4c92-be8f-e51fc2d99df8",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"from matplotlib import pyplot as plt\n",
|
||
|
|
"import os\n",
|
||
|
|
"import numpy as np\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "markdown",
|
||
|
|
"id": "8834938f-8857-45f2-9ab2-739fdc081fef",
|
||
|
|
"metadata": {},
|
||
|
|
"source": [
|
||
|
|
"# 常用文件变量\n",
|
||
|
|
"1. 参数项\n",
|
||
|
|
"```\n",
|
||
|
|
"-id D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\\TempData\\lee_filter -od D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\\TempData\\out_result -iodf T3 -ang D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\\TempData\\incidence.dat -ofr 0 -ofc 0 -fnr 8064 -fnc 7086 -fr 5.3 -un 0 -th1 100.0 -th2 700.0\n",
|
||
|
|
"```\n",
|
||
|
|
"2. 程序所在地址\n",
|
||
|
|
"```\n",
|
||
|
|
"D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\x64\\Release\n",
|
||
|
|
"```\n",
|
||
|
|
"3. 数据工作空间\n",
|
||
|
|
"```\n",
|
||
|
|
"D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\n",
|
||
|
|
"```\n",
|
||
|
|
"\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 3,
|
||
|
|
"id": "7064c90e-143e-4e8f-88a0-18bdb6a7e27f",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"T3_folder=r\"D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\\TempData\\lee_filter\"\n",
|
||
|
|
"inc_folder=r\"D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\\TempData\\incidence.dat\""
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 4,
|
||
|
|
"id": "824ac42f-8d7a-43a9-b9c5-69f0208faa2f",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
" 驱动器 D 中的卷是 新加卷\n",
|
||
|
|
" 卷的序列号是 D0A2-A693\n",
|
||
|
|
"\n",
|
||
|
|
" D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\\TempData\\lee_filter 的目录\n",
|
||
|
|
"\n",
|
||
|
|
"2024/02/07 03:49 <DIR> .\n",
|
||
|
|
"2024/02/07 03:49 <DIR> ..\n",
|
||
|
|
"2024/01/24 17:36 97 config.txt\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T11.bin\n",
|
||
|
|
"2024/02/07 03:49 76,106,970 T11.bin.enp\n",
|
||
|
|
"2024/02/07 03:49 243 T11.hdr\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T12_imag.bin\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T12_real.bin\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T13_imag.bin\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T13_real.bin\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T22.bin\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T23_imag.bin\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T23_real.bin\n",
|
||
|
|
"2024/01/24 17:39 228,566,016 T33.bin\n",
|
||
|
|
" 12 个文件 2,133,201,454 字节\n",
|
||
|
|
" 2 个目录 471,848,288,256 可用字节\n"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"!dir D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\\TempData\\lee_filter"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "markdown",
|
||
|
|
"id": "ea19f925-cbcf-463b-83d6-3f2627d2c780",
|
||
|
|
"metadata": {},
|
||
|
|
"source": [
|
||
|
|
"# 展示原始数据数据 T3、入射角数据\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 29,
|
||
|
|
"id": "258bc4f2-3bc2-4195-83ac-7619418d21de",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"ename": "ValueError",
|
||
|
|
"evalue": "cannot reshape array of size 46188060 into shape (8064,7086)",
|
||
|
|
"output_type": "error",
|
||
|
|
"traceback": [
|
||
|
|
"\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
|
||
|
|
"\u001b[1;31mValueError\u001b[0m Traceback (most recent call last)",
|
||
|
|
"Cell \u001b[1;32mIn[29], line 15\u001b[0m\n\u001b[0;32m 13\u001b[0m T23\u001b[38;5;241m=\u001b[39mT23\u001b[38;5;241m.\u001b[39mreshape(\u001b[38;5;241m8064\u001b[39m,\u001b[38;5;241m7086\u001b[39m)\n\u001b[0;32m 14\u001b[0m T33\u001b[38;5;241m=\u001b[39mT33\u001b[38;5;241m.\u001b[39mreshape(\u001b[38;5;241m8064\u001b[39m,\u001b[38;5;241m7086\u001b[39m)\n\u001b[1;32m---> 15\u001b[0m incidenceData\u001b[38;5;241m=\u001b[39m\u001b[43mincidenceData\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mreshape\u001b[49m\u001b[43m(\u001b[49m\u001b[38;5;241;43m8064\u001b[39;49m\u001b[43m,\u001b[49m\u001b[38;5;241;43m7086\u001b[39;49m\u001b[43m)\u001b[49m\n",
|
||
|
|
"\u001b[1;31mValueError\u001b[0m: cannot reshape array of size 46188060 into shape (8064,7086)"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"T11=np.fromfile(os.path.join(T3_folder,\"T11.bin\"),dtype=np.float32)\n",
|
||
|
|
"T12=np.fromfile(os.path.join(T3_folder,\"T12_real.bin\"),dtype=np.float32)+np.fromfile(os.path.join(T3_folder,\"T12_imag.bin\"),dtype=np.float32)*1j\n",
|
||
|
|
"T13=np.fromfile(os.path.join(T3_folder,\"T13_real.bin\"),dtype=np.float32)+np.fromfile(os.path.join(T3_folder,\"T13_imag.bin\"),dtype=np.float32)*1j\n",
|
||
|
|
"T22=np.fromfile(os.path.join(T3_folder,\"T22.bin\"),dtype=np.float32) \n",
|
||
|
|
"T23=np.fromfile(os.path.join(T3_folder,\"T23_real.bin\"),dtype=np.float32)+np.fromfile(os.path.join(T3_folder,\"T23_imag.bin\"),dtype=np.float32)*1j\n",
|
||
|
|
"T33=np.fromfile(os.path.join(T3_folder,\"T33.bin\"),dtype=np.float32)\n",
|
||
|
|
"incidenceData=np.fromfile(r\"D:\\codestorage\\testpolsarpro50\\prosarpro_sufarce_oh2004\\testpolsarpro50\\testpolsarpro50\\testpolsarpro50\\TestData\\TempData\\incidence.dat\",dtype=np.float64)\n",
|
||
|
|
"# 数据格式重整化\n",
|
||
|
|
"T11=T11.reshape(8064,7086)\n",
|
||
|
|
"T12=T12.reshape(8064,7086)\n",
|
||
|
|
"T13=T13.reshape(8064,7086)\n",
|
||
|
|
"T22=T22.reshape(8064,7086)\n",
|
||
|
|
"T23=T23.reshape(8064,7086)\n",
|
||
|
|
"T33=T33.reshape(8064,7086)\n",
|
||
|
|
"incidenceData=incidenceData.reshape(8781,10520)"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 30,
|
||
|
|
"id": "f95a7843-3fd8-480a-aa4c-a36950607821",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"text/plain": [
|
||
|
|
"5727.6860119047615"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"execution_count": 30,
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "execute_result"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 27,
|
||
|
|
"id": "8f72330a-2a6f-453b-8aad-18b4a59cd0cf",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T11\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAGVCAYAAAALyQ1dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOy9e3yW1ZXo/33VJJogBDQBSdQEh2Ab0ATagA1YgjMBh8sxcg6XcmkDDpdfI8hlRFoEQU4jDMKIablUSA+GI3AG08OlQhwJSlIxPZBUSJVQSdREISmQIKEmEZ/fH2vtZxEvHRmYovTdn8/7eZPn3c/ea6+991p7r2vA8zyPYAmWYAmWYAmWLyjXXGkAgiVYgiVYguXrW4JMIliCJViCJVi+tASZRLAES7AES7B8aQkyiWAJlmAJlmD50hJkEsESLMESLMHypSXIJIIlWIIlWILlS0uQSQRLsARLsATLl5YgkwiWYAmWYAmWLy1BJhEswRIswRIsX1qCTCJYgiVYgiVYvrQEmUSwBEuwBMs3oLz22msMHTqUzp07EwgE+PWvf/0fvvPqq6/Sq1cvrr/+erp06cLq1asvut8gkwiWYAmWYPkGlMbGRu6++25ycnK+Uv3Kykr+8R//kX79+lFaWspPfvITpk2bxtatWy+q30AwwF+wBEuwBMs3qwQCAfLz83nggQe+tM6cOXPYtm0bb731lv9sypQp/P73v+f111//yn1ddymABkuwBEuw/K2Vjz/+mObm5svSlud5BAKBVs/CwsIICwu75LZff/110tPTWz0bOHAg69ato6WlhZCQkK/UTpBJBEuwBEuwfMXy8ccfE39DDMc5dVnaa9OmDWfPnm31bMGCBTzxxBOX3Pbx48fp2LFjq2cdO3bkk08+4U9/+hO33HLLV2onyCSCJViCJVi+YmlubuY4p3iffNoScUltnaGRW89m8P7779O2bVv/+eW4Rbjy2VuK0y589vlfKkEmESzBEizBcpGlLe1oS5tLbEXIb9u2bVsxictVOnXqxPHjx1s9q62t5brrruOmm276yu0ErZuCJViCJViuwnLPPffw8ssvt3pWUFDAd77zna+sj4AgkwiWYAmWYPlPlNDL9Pnq5ezZs5SVlVFWVgaIiWtZWRnvvfceAHPnzmX8+PF+/SlTpvDuu+8yc+ZM3nrrLdavX8+6deuYPXv2RfUbFDcFS7AES7BcdAnRz6W28dXL//t//4+0tDT//5kzZwLwwx/+kF/96ld8+OGHPsMAiI+P5ze/+Q0zZszg5z//OZ07d2blypUMHz78ovoN+kkES7AES7B8xXLmzBnatWtHA29csk7iDGdpR28aGhr+S3QSl6sEbxLBEizBEiwXXf76N4krVYJMIliCJViC5aJLCBerU/jiNr7+Jai4DpZgCZZgCZYvLcGbRLAES7AEy0WXoLgpWIIlWIIlWL60XLwJ6xe38fUvQXFTsARLsARLsHxpCd4kgiVYgiVYLrpcx6WLi74Z5PebAWWwBEuwBMvXqgTFTd/48otf/IL4+Hiuv/56evXqxb59+y5Lu/9RCkHP83jiiSfo3LkzN9xwA/3796e8vLxVnaamJh5++GFuvvlmIiIiGDZsGNXV1a3qnD59mnHjxtGuXTvatWvHuHHjqK+v/xw82dnZfPe73+XGG28kOjqaBx54gCNHjlxRmFatWsVdd93lBy675557eOmll64YPF9UsrOzCQQCPPLII1cMrieeeIJAINDq06lTpyuOp5qaGsaOHctNN91EeHg4SUlJHDhw4IrBFRcX9zk8BQIBfvzjH19BPIVcps83oHhXYdm0aZMXEhLi/fKXv/T+8Ic/eNOnT/ciIiK8d99995Lb/s1vfuP99Kc/9bZu3eoBXn5+fqvfn3rqKe/GG2/0tm7d6h06dMgbOXKkd8stt3hnzpzx60yZMsWLiYnxXn75Ze/gwYNeWlqad/fdd3uffPKJX2fQoEFe9+7dvd/+9rfeb3/7W6979+7ekCFDPgfPwIEDvdzcXO/w4cNeWVmZN3jwYO+2227zzp49e8Vg2rZtm7dz507vyJEj3pEjR7yf/OQnXkhIiHf48OErAs9nS0lJiRcXF+fddddd3vTp068YnhYsWOAlJiZ6H374of+pra29YvB4nuedOnXKu/32270f/ehH3htvvOFVVlZ6//7v/+798Y9/vGJw1dbWtsLRyy+/7AFeYWHhXx2ehoYGD/AaOOl5tFzSp4GT0lZDwxfOxdelXJVMIiUlxZsyZUqrZ3feeaf32GOPXdZ+PsskPv30U69Tp07eU0895T/7+OOPvXbt2nmrV6/2PM/z6uvrvZCQEG/Tpk1+nZqaGu+aa67xdu3a5Xme5/3hD3/wAG///v1+nddff90DvLfffvsvwlRbW+sB3quvvvq1gcnzPK99+/bec889d8Xh+eijj7yuXbt6L7/8svf973/fZxJXAq4FCxZ4d9999xfCeaXwNGfOHK9v375fir8rPX+e53nTp0/37rjjDu/TTz/9q8PjmEQ9Dd6neJf0qafhG8EkrjpxU3NzMwcOHPhc2r709HR++9vf/pf2XVlZyfHjx1v1HRYWxve//32/7wMHDtDS0tKqTufOnenevbtf5/XXX6ddu3b07t3br9OnTx/atWv3H46hoaEBgA4dOnwtYDp//jybNm2isbGRe+6554rD8+Mf/5jBgwfz93//962eXym4jh49SufOnYmPj2fUqFEcO3bsisKzbds2vvOd7/A//sf/IDo6muTkZH75y19ecTy50tzcTF5eHhMmTCAQCFwxeJov0+ebUK46JvGnP/2J8+fPf2Havs8m4LjcxbX/l/o+fvw4oaGhtG/f/i/WiY6O/lz70dHRf3EMnucxc+ZM+vbtS/fu3a8oTIcOHaJNmzaEhYUxZcoU8vPz+fa3v31FcbRp0yYOHjxIdnb25367EnD17t2bDRs2sHv3bn75y19y/Phxvve973Hy5Mkrhqdjx46xatUqunbtyu7du5kyZQrTpk1jw4YNVwxPF5Zf//rX1NfX86Mf/ehrAc/fQrlqrZu+KG3fxaTs+2v3/dk6X1T/P2onKyuLN998k6KioisOU7du3SgrK6O+vp6tW7fywx/+kFdfffWKwfP+++8zffp0CgoKuP7667+0j78mXPfff7//d48ePbjnnnu44447+F//63/Rp0+fvzo8AJ9++inf+c53+NnPfgZAcnIy5eXlrFq1qlWugiu1xtetW8f9999P586dWz3/a8PTop9LKZf6/l+rXHU3iZtvvplrr732C9P2ffa0cbmLs0z5S3136tSJ5uZmTp8+/RfrnDhx4nPt19XVfekYHn74YbZt20ZhYSGxsbFXHKbQ0FD+7u/+ju985ztkZ2dz991388wzz1wxeA4cOEBtbS29evXiuuuu47rrruPVV19l5cqVXHfddf47V2LuXImIiKBHjx4cPXr0iuHplltu4dvf/narZ9/61rf8PAVXco2/++67/Pu//zsPPfSQ/+xKwRMUN32DS2hoKL169fpc2r6XX36Z733ve/+lfcfHx9OpU6dWfTc3N/Pqq6/6fffq1YuQkJBWdT788EMOHz7s17nnnntoaGigpKTEr/PGG2/Q0NDwuTF4nkdWVhYvvvgie/bsIT4+/orD9EXF8zyampquGDz33Xcfhw4d8jN7lZWV8Z3vfIcxY8ZQVlZGly5drjiempqaeOutt7jllluuGJ5SU1M/Z0JdUVHB7bffDlzZ9ZSbm0t0dDSDBw/2n31d1vdVXf5qKvK/YnEmsOvWrfP+8Ic/eI888ogXERHhVVVVXXLbH330kVdaWuqVlpZ6gLd8+XKvtLTUN6996qmnvHbt2nkvvviid+jQIW/06NFfaI4XGxvr/fu//7t38OBBb8CAAV9ojnfXXXd5r7/+uvf66697PXr0+ELzwKlTp3rt2rXz9u7d28pM8Ny5c36dvzZMc+fO9V577TWvsrLSe/PNN72f/OQn3jXXXOMVFBRcEXi+rFxo3XQl4Jo1a5a3d+9e79ixY97+/fu9IUOGeDfeeKO/Tq8EnkpKSrzrrrv
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T12 amp\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAGVCAYAAAALyQ1dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOy9f3icV3Un/nnzZmbeaMbj8QzSTCwhZKmx3NpKiJPahMRA0t3ELAkPpLuQLri0pE9J2yRA6C7QLj+2PNuEdhva0EJYGugW2A39bpIuga6T0JjgQLB3Y7soDpZTy0aRkpHESKPJjDo/mNzvH59z7r3jJGBjgxMz93nGGr9z3/vj3HvP73NuYIwx6JZu6ZZu6ZZueZ5yxqkeQLd0S7d0S7e8eEuXSHRLt3RLt3TLC5YukeiWbumWbumWFyxdItEt3dIt3dItL1i6RKJbuqVbuqVbXrB0iUS3dEu3dEu3vGDpEolu6ZZu6ZZuecHSJRLd0i3d0i3d8oKlSyS6pVu6pVu65QVLl0h0S7d0S7d0ywuWLpHolm7plm55CZRvfvObuOqqq7B69WoEQYC///u//7HvPPTQQ7jgggsQRRGGh4dx++23H3e/XSLRLd3SLd3yEii1Wg3nnXce/vIv//KY6h8+fBj/5t/8G2zZsgV79+7FH/zBH+DGG2/EXXfddVz9Bt0Ef93SLd3SLS+tEgQB7rnnHrzpTW96wTrvf//78ZWvfAXf+9737LPrrrsO//RP/4RHHnnkmPs680QG2i3d0i3d8vNW6vU6ms3mSWnLGIMgCDqeJRIJJBKJE277kUceweWXX97x7IorrsAdd9yBVquFWCx2TO10iUS3dEu3dMsxlnq9jjVn9aOIhZPSXiqVQrVa7Xj2kY98BB/96EdPuO1isYh8Pt/xLJ/P44c//CF+8IMf4Oyzzz6mdrpEolu6pVu65RhLs9lEEQt4EvcgjeQJtVVBDS+vvhlPPvkk0um0fX4ypAgtR0spal04+vmPKl0i0S3d0i3dcpwljZVII3WCrRD9ptPpDiJxskqhUECxWOx4Njc3hzPPPBO5XO6Y2+l6N3VLt3RLt5yG5aKLLsIDDzzQ8ez+++/HhRdeeMz2CKBLJLqlW7qlW36CEj9Jn2Mv1WoV+/btw759+wDQxXXfvn2YmpoCAHzwgx/Er//6r9v61113Hb7//e/jpptuwve+9z187nOfwx133IHf//3fP65+u+qmbumWbumW4y4x+ZxoG8de/t//+3+49NJL7f9vuukmAMA73vEO/M3f/A2efvppSzAAYM2aNfiHf/gHvPe978Vf/dVfYfXq1bjtttvwq7/6q8fVbzdOolu6pVu65RhLpVLBypUrsYRdJ2yTqKCKldiMpaWln4pN4mSVriTRLd3SLd1y3OVnL0mcqtIlEt3SLd3SLcddYjhem8Lzt/HiL13Ddbd0S7d0S7e8YOlKEt3SLd3SLcdduuqmbumWbumWbnnBcvwurM/fxou/dNVN3dIt3dIt3fKCpStJdEu3dEu3HHc5EyeuLnppoN+Xxii7pVu6pVteVKWrbnrJl0996lNYs2YNoijCBRdcgJ07d56Udn/cFYLGGHz0ox/F6tWrcdZZZ+F1r3sd9u/f31Gn0WjghhtuwMte9jIkk0m88Y1vxPT0dEedxcVFbNu2DStXrsTKlSuxbds2lMvl54zn5ptvxi//8i9jxYoV6Ovrw5ve9CZMTEyc0jF9+tOfxrnnnmsTl1100UX4P//n/5yy8TxfufnmmxEEAd7znvecsnF99KMfRRAEHZ9CoXDK4TQzM4O3v/3tyOVy6OnpwStf+Uo8+uijp2xcQ0NDz4FTEAT4vd/7vVMIp9hJ+rwEijkNy5133mlisZj57Gc/ax5//HHz7ne/2ySTSfP973//hNv+h3/4B/OHf/iH5q677jIAzD333NPx+y233GJWrFhh7rrrLjM+Pm7e+ta3mrPPPttUKhVb57rrrjP9/f3mgQceMHv27DGXXnqpOe+888wPf/hDW2fr1q1mw4YN5tvf/rb59re/bTZs2GCuvPLK54zniiuuMJ///OfNY489Zvbt22fe8IY3mMHBQVOtVk/ZmL7yla+Yr33ta2ZiYsJMTEyYP/iDPzCxWMw89thjp2Q8R5fdu3eboaEhc+6555p3v/vdpwxOH/nIR8z69evN008/bT9zc3OnbDzGGLOwsGBe8YpXmN/4jd8wu3btMocPHzZf//rXzT//8z+fsnHNzc11wOiBBx4wAMyOHTt+5uNZWloyAMwSSsagdUKfJZTY1tLS867Fi6WclkRi06ZN5rrrrut4tm7dOvOBD3zgpPZzNJF49tlnTaFQMLfccot9Vq/XzcqVK83tt99ujDGmXC6bWCxm7rzzTltnZmbGnHHGGWb79u3GGGMef/xxA8B85zvfsXUeeeQRA8AcOHDgR45pbm7OADAPPfTQi2ZMxhizatUq89d//denfDzPPPOMOeecc8wDDzxgXvva11oicSrG9ZGPfMScd955zzvOUwWn97///eaSSy55Qfid6vUzxph3v/vdZmRkxDz77LM/8/EokShjyTwLc0KfMpZeEkTitFM3NZtNPProo8+5tu/yyy/Ht7/97Z9q34cPH0axWOzoO5FI4LWvfa3t+9FHH0Wr1eqos3r1amzYsMHWeeSRR7By5Ups3rzZ1nnVq16FlStX/tg5LC0tAQCy2eyLYkztdht33nknarUaLrroolM+nt/7vd/DG97wBvyrf/WvOp6fqnE98cQTWL16NdasWYNrrrkGk5OTp3Q8X/nKV3DhhRfi3/27f4e+vj6cf/75+OxnP3vK4aSl2Wzii1/8It75znciCIJTNp7mSfq8FMppRyR+8IMfoN1uP++1fUdfwHGyi7b/o/ouFouIx+NYtWrVj6zT19f3nPb7+vp+5ByMMbjppptwySWXYMOGDad0TOPj40ilUkgkErjuuutwzz334Jd+6ZdOKYzuvPNO7NmzBzfffPNzfjsV49q8eTP+9m//Fvfddx8++9nPolgs4tWvfjVKpdIpg9Pk5CQ+/elP45xzzsF9992H6667DjfeeCP+9m//9pTByS9///d/j3K5jN/4jd94UYzn56Gctt5Nz3dt3/Fc2fez7vvoOs9X/8e1c/311+O73/0uHn744VM+ptHRUezbtw/lchl33XUX3vGOd+Chhx46ZeN58skn8e53vxv3338/oih6wT5+luN6/etfb7+PjY3hoosuwsjICP77f//veNWrXvUzHw8APPvss7jwwgvxx3/8xwCA888/H/v378enP/3pjrsKTtUev+OOO/D6178eq1ev7nj+sx5PSz4nUk70/Z9VOe0kiZe97GUIw/B5r+07mts42UU9U35U34VCAc1mE4uLiz+yzuzs7HPan5+ff8E53HDDDfjKV76CHTt2YGBg4JSPKR6P4xd+4Rdw4YUX4uabb8Z5552Hv/iLvzhl43n00UcxNzeHCy64AGeeeSbOPPNMPPTQQ7jttttw5pln2ndOxdppSSaTGBsbwxNPPHHK4HT22Wfjl37plzqe/eIv/qK9p+BU7vHvf//7+PrXv47f+q3fss9O1Xi66qaXcInH47jggguec23fAw88gFe/+tU/1b7XrFmDQqHQ0Xez2cRDDz1k+77gggsQi8U66jz99NN47LHHbJ2LLroIS0tL2L17t62za9cuLC0tPWcOxhhcf/31uPvuu/Hggw9izZo1p3xMz1eMMWg0GqdsPL/yK7+C8fFxe7PXvn37cOGFF+Jtb3sb9u3bh+Hh4VMOp0ajge9973s4++yzTxmcLr744ue4UB88eBCveMUrAJza/fT5z38efX19eMMb3mCfvVj292ldfmYm8p9hURfYO+64wzz++OPmPe95j0kmk+bIkSMn3PYzzzxj9u7da/bu3WsAmFtvvdXs3bvXutfecsstZuXKlebuu+824+Pj5td+7dee1x1vYGDAfP3rXzd79uwxl1122fO645177rnmkUceMY888ogZGxt7XvfA3/md3zErV6403/jGNzrcBJeXl22dn/WYPvjBD5pvfvOb5vD
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T13 amp\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAGVCAYAAAALyQ1dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9fZhkV3Xfi3+KUXUNXa3uox5Vd9M9GloayzNgCWJJQQgsAzIICAJzufcagpEhYK4VDLYsHGPMk5iQBAG5wY5NkBMb7BhIxO8G5AcSWSAbkIWRhK2XWBKMGI/Ubk2PerqY1plWVzPVxXB+f+z13Xvt6pGs0YwZvdR6nuqqPmef/bLO3ut9r12rqqpiAAMYwAAGMIAjwNNOdAcGMIABDGAAj18YMIkBDGAAAxjAw8KASQxgAAMYwAAeFgZMYgADGMAABvCwMGASAxjAAAYwgIeFAZMYwAAGMIABPCwMmMQABjCAAQzgYWHAJAYwgAEMYAAPCwMmMYABDGAAA3hYGDCJAQxgAAMYwMPCgEkMYAADGMATAP7iL/6CV7/61UxPT1Or1fiTP/mTv/eZG264gXPPPZfNmzdzxhln8Hu/93tH3e6ASQxgAAMYwBMAOp0Oz33uc/nYxz72qMrfd999/JN/8k+48MILuf322/mN3/gNfumXfonPfe5zR9VubZDgbwADGMAAnlhQq9W45ppreO1rX/uwZd7znvfwhS98gW9/+9vx2mWXXcb//t//m5tuuulRt3XSsXR0AAMYwACeanDo0CHW19ePS11VVVGr1bJrjUaDRqNxzHXfdNNNXHzxxdm1l7/85XziE5+g1+tRr9cfVT0DJjGAAQxgAI8SDh06xOlPn2GR5eNS38jICKurq9m13/zN3+T973//Mde9uLjI5ORkdm1ycpLvf//7fPe73+UZz3jGo6pnwCQGMIABDOBRwvr6Oosscz/XMErzmOpaocNpq/8H999/P6Ojo/H68dAiBP1airwL/dcfCQZMYgADGMAAjhJGGWOUkWOsJZDf0dHRjEkcL5iammJxcTG7trS0xEknncSWLVsedT2D6KYBDGAAA3gSwgUXXMD111+fXfvyl7/Meeed96j9ETBgEgMYwAAG8Bhg6Dh9Hj2srq5yxx13cMcddwAhxPWOO+5gfn4egPe+97383M/9XCx/2WWX8Xd/93dcccUVfPvb3+aTn/wkn/jEJ/jVX/3Vo2p3YG4awAAGMICjhrp9jrWORw9//dd/zUte8pL4/xVXXAHAm9/8Zv7oj/6IBx54IDIMgNNPP51rr72WX/mVX+E//af/xPT0NL/zO7/D//l//p9H1e5gn8QABjCAATxKWFlZYWxsjIPccsw+iRVWGeN8Dh48+A/ikzheMNAkBjCAAQzgqOGHr0mcKBgwiQEMYAADOGqoc7Q+hSPX8fiHgeN6AAMYwAAG8LAw0CQGMIABDOCoYWBuGsAABjCAATwsHH0I65HrePzDwNw0gAEMYAADeFgYaBIDGMAABnDUcBLHbi56YpDfJ0YvBzCAAQzgcQUDc9MTHj7+8Y9z+umns3nzZs4991xuvPHG41Lv33eEYFVVvP/972d6epqnP/3pvPjFL+buu+/OynS7Xd71rndx6qmn0mw2ec1rXsPevXuzMg8++CCXXnopY2NjjI2Ncemll1KW5Yb+XHnllfzjf/yPOfnkk5mYmOC1r30t99xzzwnt01VXXcVznvOcmLjsggsu4E//9E9PWH+OBFdeeSW1Wo3LL7/8hPXr/e9/P7VaLftMTU2dcDwtLCzwpje9iS1btjA8PMw/+kf/iFtvvfWE9Wt2dnYDnmq1Gr/4i794AvFUP06fJwBUT0K4+uqrq3q9Xv3+7/9+9a1vfav65V/+5arZbFZ/93d/d8x1X3vttdX73ve+6nOf+1wFVNdcc012/0Mf+lB18sknV5/73OeqO++8s3r9619fPeMZz6hWVlZimcsuu6yamZmprr/++uq2226rXvKSl1TPfe5zq+9///uxzCte8YrqrLPOqr7xjW9U3/jGN6qzzjqruuSSSzb05+Uvf3n1h3/4h9Vdd91V3XHHHdWrXvWqatu2bdXq6uoJ69MXvvCF6n/9r/9V3XPPPdU999xT/cZv/EZVr9eru+6664T0px+++c1vVrOzs9VznvOc6pd/+ZdPGJ5+8zd/s/qxH/ux6oEHHoifpaWlE9afqqqq5eXl6pnPfGb1lre8pbrllluq++67r/qzP/uz6m//9m9PWL+WlpYyHF1//fUVUH31q1/9offn4MGDFVAd5EBV0Tumz0EOhLoOHjziu3i8wJOSSTzvec+rLrvssuzazp07q1//9V8/ru30M4kf/OAH1dTUVPWhD30oXjt06FA1NjZW/d7v/V5VVVVVlmVVr9erq6++OpZZWFionva0p1XXXXddVVVV9a1vfasCqptvvjmWuemmmyqg2rVr1yP2aWlpqQKqG2644XHTp6qqqlNOOaX6gz/4gxPen4ceeqg688wzq+uvv7560YteFJnEiejXb/7mb1bPfe5zj9jPE4Wn97znPdVP/MRPPCz+TvT7q6qq+uVf/uVq+/bt1Q9+8IMfen/EJEoOVj+gOqZPycEnBJN40pmb1tfXufXWWzcc23fxxRfzjW984x+07fvuu4/FxcWs7UajwYte9KLY9q233kqv18vKTE9Pc9ZZZ8UyN910E2NjY5x//vmxzPOf/3zGxsb+3jEcPHgQgPHx8cdFnw4fPszVV19Np9PhggsuOOH9+cVf/EVe9apX8dKXvjS7fqL6tXv3bqanpzn99NN5wxvewL333ntC+/OFL3yB8847j//7//6/mZiY4Md//Mf5/d///ROOJ8H6+jqf/vSneetb30qtVjth/Vk/Tp8nAjzpmMR3v/tdDh8+fMRj+/oP4DjeoPofqe3FxUWGhoY45ZRTHrHMxMTEhvonJiYecQxVVXHFFVfwEz/xE5x11lkntE933nknIyMjNBoNLrvsMq655hqe/exnn1AcXX311dx2221ceeWVG+6diH6df/75/PEf/zFf+tKX+P3f/30WFxd5wQtewIEDB04Ynu69916uuuoqzjzzTL70pS9x2WWX8Uu/9Ev88R//8QnDk4c/+ZM/oSxL3vKWtzwu+vNUgCdtdNORju07miP7ftht95c5Uvm/r553vvOd/M3f/A1f//rXT3ifduzYwR133EFZlnzuc5/jzW9+MzfccMMJ68/999/PL//yL/PlL3+ZzZs3P2wbP8x+vfKVr4y/zz77bC644AK2b9/Of/2v/5XnP//5P/T+APzgBz/gvPPO44Mf/CAAP/7jP87dd9/NVVddlZ1VcKLm+Cc+8Qle+cpXMj09nV3/YfenZ59jgWN9/ocFTzpN4tRTT2XTpk1HPLavX9o43qDIlEdqe2pqivX1dR588MFHLLN///4N9bfb7Ycdw7ve9S6+8IUv8NWvfpWtW7ee8D4NDQ3xIz/yI5x33nlceeWVPPe5z+U//sf/eML6c+utt7K0tMS5557LSSedxEknncQNN9zA7/zO73DSSSfFZ07EuxM0m03OPvtsdu/efcLw9IxnPINnP/vZ2bVnPetZ8ZyCEznH/+7v/o4/+7M/4+d//ufjtRPVn4G56QkMQ0NDnHvuuRuO7bv++ut5wQte8A/a9umnn87U1FTW9vr6OjfccENs+9xzz6Ver2dlHnjgAe66665Y5oILLuDgwYN885vfjGVuueUWDh48uGEMVVXxzne+k89//vN85Stf4fTTTz/hfToSVFVFt9s9Yf35qZ/6Ke688854stcdd9zBeeedx8/+7M9yxx13cMYZZ5xwPHW7Xb797W/zjGc844Th6YUvfOGGEOrvfOc7PPOZzwRO7Hz6wz/8QyYmJnjVq14Vrz1e5veTGn5oLvIfIigE9hOf+ET1rW99q7r88surZrNZzc3NHXPdDz30UHX77bdXt99+ewVUH/3oR6vbb789htd+6EMfqsbGxqrPf/7z1Z133ln903/6T48Yjrd169bqz/7sz6rbbrutuuiii44Yjvec5zynuummm6qbbrqpOvvss48
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T12 pha\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAGICAYAAACqSCyrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOydfXyU1ZX4v9N1Qp2sMqiTYqbqUJus7qCGVdJqqIYqUTZITcsOxSASwSAS5UUKRZE0lhpFyosNUiIaS0E0i8alpGqQJlLiS5CSVcYqqTKtDdJMLRNsRs1Y7++Pc57nwqpdLSj82Dmfz3ySmXnm3nPPvfecc889Lx5jjCENaUhDGtKQhv3gC4cbgTSkIQ1pSMORB2nhkIY0pCENafgQpIVDGtKQhjSk4UOQFg5pSEMa0pCGD0FaOKQhDWlIQxo+BGnhkIY0pCENafgQpIVDGtKQhjSk4UOQFg5pSEMa0pCGD8ExhxuBNKQhDWn4/wneffddent7D0lbGRkZfPGLXzwkbR1qSAuHNKQhDWn4hPDuu+8y4Ngge/jLIWmvf//+7Nq164gUEGnhkIY0pCENnxB6e3vZw194gwaOJ/Og2tpHD6fsKaG3tzctHNKQhjSk4WiA48k8aOFwpENaOKQhDWlIw6eGYwDvIWjjyIW0t1Ia0pCGNHxqyDhEr08Gy5cv5+yzz+b444/n+OOP5/zzz+fxxx8/dMP5CEgLhzSkIQ1pOMLhy1/+MnfccQcvvPACL7zwAt/85jf51re+RTQa/cz69KTrOaQhDWlIwyeDffv20bdvX7p5nuP554Nri7/Sl6/R3d3N8ccf/6l/f8IJJ3DXXXcxYcKEg8Lj4+DINnqlIQ1pSMMRCV4O/s5Bfr9v374DPu3Tpw99+vT52F/97W9/4z//8z/p6enh/PPPP0gcPh7SZqU0pCENaTiMcMopp9C3b1/3VV1d/ZHPvfTSS/zzP/8zffr04brrrqOhoYF//dd//czwSp8c0pCGNKThU4OXT3Oh/PFtwBtvvHGAWenjTg3/8i//Qnt7O4lEgkceeYSrr76ap59++jMTEOk7hzSkIQ1p+IRg7xxe4XiOO7i2eJu+nPEP3zlccsklnH766axYseKg8Pg4SJuV0pCGNKTh/0MwxvDee+99Zu2nzUppSEMa0vCp4dPFKXx8G58Mbr75ZoYPH84pp5zC22+/zUMPPURLSwtPPPHEQeLw8ZAWDmlIQxrS8Knh842Q/tOf/sRVV13Fm2++Sd++fTn77LN54oknGDZs2EHi8PFw1JqV7rnnHgYMGMAXv/hFzj33XH79618fknY3b97M5ZdfTnZ2Nh6Ph8cee+yA740x/OAHPyA7O5tjjz2WwsLCDwWqvPfee9xwww2cdNJJZGZmMnLkSP74xz8e8MzevXu56qqrXA+Gq666ikQi8SF8qqurGTx4MMcddxxZWVlcccUVvPrqq4cVp/8tmvPzxuejoLq6Go/Hw7Rp0w4bXj/4wQ/weDwHvPr373/Y6dTZ2cnYsWM58cQT8fl85OXlsW3btsOGVygU+hCdPB4PU6ZMOYx0+nwjpO+77z5isRjvvfceXV1dPPXUU5+pYADAHIXw0EMPGa/Xa+69917z8ssvm6lTp5rMzEzz+9///qDb/uUvf2luueUW88gjjxjANDQ0HPD9HXfcYY477jjzyCOPmJdeesmMHj3anHzyyWbfvn3uM9ddd50JBoNm48aN5je/+Y0ZOnSoOeecc8z777/vPnPZZZeZgQMHmmeeecY888wzZuDAgWbEiBEfwufSSy81dXV1ZseOHaa9vd0UFxebU0891fz1r389bDitX7/eNDY2mldffdW8+uqr5uabbzZer9fs2LHjsODzP6Gtrc2EQiFz9tlnm6lTpx42OlVWVppwOGzefPNN99XV1XXY8DHGmL/85S/mtNNOM+PHjzfPP/+82bVrl3nqqafM7373u8OGV1dX1wE02rhxowFMc3Pz545Pd3e3AUw3e4wheVCvbvZIW93dHzkXhxuOSuGQn59vrrvuugM+O+OMM8z3v//9Q9rP/xQOH3zwgenfv7+544473M/effdd07dvX/PTn/7UGGNMIpEwXq/XPPTQQ+4znZ2d5gtf+IJ54oknjDHGvPzyywYwzz33nPvMs88+awDzyiuv/F2curq6DGCefvrpIwYnY4zp16+fWbly5WHH5+233zY5OTlm48aN5qKLLnKFw+HAq7Ky0pxzzjkfiefhotPs2bPNkCFDPpZ+h3v+jDFm6tSp5vTTTzcffPDB546PFQ5vGUPqoF7dvHVEC4ejzqzU29vLtm3bKCoqOuDzoqIinnnmmc+07127drFnz54D+u7Tpw8XXXSR2/e2bdtIpVIHPJOdnc3AgQPdZ5599ln69u3L1772NfeZr3/96/Tt2/d/HUN3dzcgofVHAk5/+9vfeOihh9xozsONz5QpUyguLuaSSy454PPDhVdHRwfZ2dkMGDCA7373u7z++uuHFZ/169dz3nnn8R//8R9kZWUxaNAg7r333sNOJwd6e3tZvXo111xzDR6P57DhYzjmkLyOZDjqhMOf//xn/va3v/GlL33pgM+/9KUvsWfPns+0b6f9v9f3nj17yMjIoF+/fn/3maysrA+1n5WV9XfHYIxhxowZDBkyhIEDBx5WnD4umvNw0uihhx7iN7/5zUdGoB4OvL72ta+xatUqnnzySe6991727NnDBRdcwFtvvXXY6PT666+zfPlycnJyePLJJ7nuuuu48cYbWbVq1WGj0/7w2GOPkUgkGD9+/BGBz9EMR7boOgjweDwHvDfGfOizI6nv//nMRz3/v7VTUVHBiy++yJYtWw47Th8XzXm48HnjjTeYOnUqTU1Nf7fq1ueJ1/Dhw93/zzrrLM4//3xOP/10fvazn/H1r3/9c8cH4IMPPuC8887j9ttvB2DQoEFEo1GWL1/OuHHjPrbNz2uN33fffQwfPpzs7OwDPv+88enV18HAoalC/dnBUXdyOOmkk/inf/qnD0n7rq6uD2kXhxocT5O/13f//v3p7e1l7969f/eZP/3pTx9qPx6Pf+wYbrjhBtavX09zczNf/vKXDztOGRkZfPWrX+W8886jurqac845h6VLlx42fLZt20ZXVxfnnnsuxxxzDMcccwxPP/00d999N8ccc4z7m8Mxdw5kZmZy1lln0dHRcdjodPLJJ38oHcOZZ57JH/7wB7e9w4EXwO9//3ueeuopJk6c6H52uPBJHaLXkQxHnXDIyMjg3HPPZePGjQd8vnHjRi644ILPtO8BAwbQv3//A/ru7e3l6aefdvs+99xz8Xq9Bzzz5ptvsmPHDveZ888/n+7ubtra2txnnn/+ebq7uz80BmMMFRUVPProo/zqV79iwIABhx2njwKj0ZyHC5+LL76Yl156ifb2dvd13nnnUVpaSnt7O1/5ylcOO53ee+89fvvb33LyyScfNjoVFBR8yBV6586dnHbaacDhXU91dXVkZWVRXFzsfnakrO+jEj63q+/PERxX1vvuu8+8/PLLZtq0aSYzM9PEYrGDbvvtt98227dvN9u3bzeAWbRokdm+fbvrJnvHHXeYvn37mkcffdS89NJLZsyYMR/pVvflL3/ZPPXUU+Y3v/mN+eY3v/mRbnVnn322efbZZ82zzz5rzjrrrI9085s8ebLp27evaWlpOcDdL5lMus983jjNmTPHbN682ezatcu8+OKL5uabbzZf+MIXTFNT02HB5+Ngf2+lw4HXTTfdZFpaWszrr79unnvuOTNixAhz3HHHuev0cNCpra3NHHPMMeZHP/qR6ejoMGvWrDE+n8+sXr36sNHJGGP+9re/mVNPPdXMnj37Q999nvg43kq76DZvYQ7qtYvuI9pb6agUDsYYs2zZMnPaaaeZjIwM82//9m+ua+fBQnNzswE+9Lr66quNMeLqV1lZafr372/69OljLrzwQvPSSy8d0MY777xjKioqzAknnGCOPfZYM2LECPOHP/zhgGfeeustU1paao477jhz3HHHmdLSUrN3794P4fNRuACmrq7Ofebzxumaa65xaR8IBMz
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T13 pha\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAGICAYAAACqSCyrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9f5iX1XXuj98TMoEQBAKdgIiEgBywSKREREupBMWQECLBj5RmgmJpOCiU4NFQLH6JXChmvnhQSg4GMfDFi6J4Jlj8IAaDHcYzQpiDDuIYzIAZERA6QuCRjhBAX98/7rWf0WpaDSQY+97X9b7Emfc8z372s/f6ca97rVUEoMIojMIojMIojHeMT5zpCRRGYRRGYRTGR28UlENhFEZhFEZhvGcUlENhFEZhFEZhvGcUlENhFEZhFEZhvGcUlENhFEZhFEZhvGcUlENhFEZhFEZhvGcUlENhFEZhFEZhvGcUlENhFEZhFEZhvGd88kxPoDAKozAK449pHDt2TMePHz8t1/rUpz6lFi1anJZrne5RUA6FURiFURgfcBw7dkxf+PQ52q9fn5brdezYUfX19R9JBVFQDoVRGIVRGB9wHD9+XPv1a+3Wo2qtz5zStd5Qo87d/00dP368oBwKozAKozA+DqO1PnPKyuGjPgrKoTAKozAK40OPT0oqPg3X+OiOAlupMAqjMArjQ49PnabPBxv33XefvvjFL6p169Zq3bq1Lr30Uj3xxBOn73HeZxSUQ2EURmEUxkd8dO7cWT/4wQ+0ZcsWbdmyRUOGDNFVV12lF1988fd2z6JCP4fCKIzCKIwPNt544w21adNGmTartVqd2rX0b2qjAcqyTK1bt/7Qf9+uXTvNnTtX48ePP6V5/Lbx0Qa9CqMwCqMwPpKjWKcec/Dfv/HGG+/6afPmzdW8efPf+ldvvfWW/vf//t9qbGzUpZdeeopz+O2jACsVRmEURmGcwXHuueeqTZs2+eeuu+563++98MILatWqlZo3b66JEyfq0Ucf1Z/+6Z/+3uZV8BwKozAKozA+9CjWhwko//ZrSLt3734XrPTbvIaePXtq69atOnz4sH7yk5/ouuuuU2Vl5e9NQRRiDoVRGIVRGB9wNMUcXlJrnXVq19IRtVGv3znmcMUVV6h79+5atGjRKc3jt40CrFQYhVEYhfFHOAD95je/+b1dvwArFUZhFEZhfOjx4fIUfvs1Ptj4h3/4B331q1/VueeeqyNHjujhhx/Whg0b9NOf/vQU5/DbR0E5FEZhFEZhfOjxh82Q/td//VeNHTtW+/btU5s2bfTFL35RP/3pTzV06NBTnMNvHx9bWGnhwoX6whe+oBYtWuhLX/qS/s//+T+n5bpPP/20RowYoU6dOqmoqEj//M///K7fA7r99tvVqVMnffrTn9bgwYPfk6jym9/8Rn/3d3+nP/mTP9FnPvMZfeMb39CePXve9Z1Dhw5p7NixOYNh7NixOnz48Hvmc9ddd6l///4666yz9LnPfU4jR47UL3/5yzM6p/8sm/MPPZ/3G3fddZeKioo0derUMzav22+/XUVFRe/6dOzY8Yyv0969e/Xtb39b7du3V8uWLdW3b189++yzZ2xeXbt2fc86FRUVadKkSWdwnf6wGdI//vGP9corr+g3v/mNGhoatH79+t+rYpAk8TEcDz/8MMXFxSxevJhf/OIXfPe73+Uzn/kMu3btOuVrr127lhkzZvCTn/wESTz66KPv+v0PfvADzjrrLH7yk5/wwgsv8Fd/9VecffbZvPHGG/l3Jk6cyDnnnMPPfvYznnvuOb785S9z4YUXcvLkyfw7w4YN44ILLmDjxo1s3LiRCy64gK9//evvmc9XvvIVli5dSm1tLVu3bmX48OF06dKFf/u3fztjc3rsscd4/PHH+eUvf8kvf/lL/uEf/oHi4mJqa2vPyHz+/aiurqZr16588Ytf5Lvf/e4ZW6fvf//79O7dm3379uWfhoaGMzYfgF//+td8/vOfZ9y4cWzevJn6+nrWr1/Pzp07z9i8Ghoa3rVGP/vZz5BERUXFH3w+WZYhiUz7QW+e0ifTfl8ry973XZzp8bFUDhdffDETJ05818969erF9OnTT+t9/r1yePvtt+nYsSM/+MEP8p8dO3aMNm3a8KMf/QiAw4cPU1xczMMPP5x/Z+/evXziE5/gpz/9KQC/+MUvkMTPf/7z/DubNm1CEi+99NJ/OKeGhgYkUVlZ+ZGZE8BnP/tZHnjggTM+nyNHjtCjRw9+9rOfcdlll+XK4UzM6/vf/z4XXnjh+87zTK3T3//93/MXf/EXv3X9zvT7A/jud79L9+7defvtt//g82lSDgdBJ07pk+ngR1o5fOxgpePHj+vZZ5/VlVde+a6fX3nlldq4cePv9d719fXav3//u+7dvHlzXXbZZfm9n332WZ04ceJd3+nUqZMuuOCC/DubNm1SmzZtNGDAgPw7l1xyidq0afOfPkOWZZKcWv9RmNNbb72lhx9+OM/mPNPzmTRpkoYPH64rrrjiXT8/U/PasWOHOnXqpC984QsaM2aMfvWrX53R+Tz22GO66KKLdM011+hzn/uc/uzP/kyLFy8+4+uUxvHjx7V8+XL9zd/8jYqKis7YfNAnT8vnozw+dsrhwIEDeuutt9ShQ4d3/bxDhw7av3//7/Xe6fr/0b3379+vT33qU/rsZz/7H37nc5/73Huu/7nPfe4/fAZA/+N//A/9xV/8hS644IIzOqffls15Jtfo4Ycf1nPPPfe+GahnYl4DBgzQgw8+qHXr1mnx4sXav3+//vzP/1wHDx48Y+v0q1/9Svfdd5969OihdevWaeLEiZoyZYoefPDBM7ZO7xz//M//rMOHD2vcuHEfifl8nMdHW3WdwigqKnrX/wPv+dlH6d7//jvv9/3/7DqTJ0/Wtm3bVFVVdcbn9NuyOc/UfHbv3q3vfve7evLJJ//Drlt/yHl99atfzf/dp08fXXrpperevbuWLVumSy655A8+H0l6++23ddFFF2nOnDmSpD/7sz/Tiy++qPvuu0/XXnvtb73mH2qP//jHP9ZXv/pVderU6V0//0PP53h8TmWcni7Uv7/xsfMc/uRP/kTNmjV7j7ZvaGh4j3VxukdimvxH9+7YsaOOHz+uQ4cO/Yff+dd//df3XP/111//rc/wd3/3d3rsscdUUVGhzp07n/E5fepTn9J5552niy66SHfddZcuvPBCzZ8//4zN59lnn1VDQ4O+9KUv6ZOf/KQ++clPqrKyUv/4j/+oT37yk/nfnIl3l8ZnPvMZ9enTRzt27Dhj63T22We/pxzD+eefr1dffTW/3pmYlyTt2rVL69ev19/+7d/mPztT8zlxmj4f5fGxUw6f+tSn9KUvfUk/+9nP3vXzn/3sZ/rzP//z3+u9v/CFL6hjx47vuvfx48dVWVmZ3/tLX/qSiouL3/Wdffv2qba2Nv/OpZdeqizLVF1dnX9n8+bNyrLsPc8AaPLkyVq1apX+5V/+RV/4whfO+JzebxDZnGdqPpdffrleeOEFbd26Nf9cdNFFKi0t1datW9WtW7czvk6/+c1vtH37dp199tlnbJ0GDhz4Hip0XV2dPv/5z0s6s/tp6dKl+tznPqfhw4fnP/uo7O+P5fiDhb7/gCNRWX/84x/zi1/8gqlTp/KZz3yGV1555ZSvfeTIEWpqaqipqUES8+bNo6amJqfJ/uAHP6BNmzasWrWKF154gb/+679+X1pd586dWb9+Pc899xxDhgx5X1rdF7/4RTZt2sSmTZvo06fP+9L8brjhBtq0acOGDRveRfd788038+/8oed066238vTTT1NfX8+2bdv4h3/4Bz7xiU/w5JNPnpH5/LbxTrbSmZjXzTffzIYNG/jVr37Fz3/+c77+9a9z1lln5fv0TKxTdXU1n/zkJ7nzzjvZsWMH//RP/0TLli1Zvnz5GVsngLfeeosuXbrw93//9+/53R9yPomtVK+Mg+KUPvXKPtJspY+lcgD4X//rf/H5z3+eT33qU/Tr1y+ndp7qqKioQNJ7Ptdddx1gqt/3v/99OnbsSPPmzfnLv/x
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T22\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAGVCAYAAAALyQ1dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9e5zkVXnngb+/FHWhq6ampsruLqbbdi5xhgQGAxgQFRWSRVzBnzG7UReJScxrQxJBJdn1ko2a5LeBJL+YBBMlMWg26gb3tUJWjEExIEFF2GWGOIAMZIax6Zbqbqunuqhq60Jxfn88z3POqR4wjDNxYKjzevX0dNX3e+7nuXyey0mcc45RGZVRGZVRGZUnKccd7Q6MyqiMyqiMyjO3jJjEqIzKqIzKqDxlGTGJURmVURmVUXnKMmISozIqozIqo/KUZcQkRmVURmVURuUpy4hJjMqojMqojMpTlhGTGJVRGZVRGZWnLCMmMSqjMiqjMipPWUZMYlRGZVRGZVSesoyYxKiMyqiMyqg8ZRkxiVEZlVEZlWdB+ad/+icuuugiNm7cSJIk/N3f/d2/+s5tt93GGWecQS6XY8uWLVxzzTWH3O6ISYzKqIzKqDwLSrvd5kUvehF/9md/9rSef/jhh/n3//7fc84557Br1y7e9773cfnll/PZz372kNpNRgn+RmVURmVUnl0lSRJuuOEGXv/61z/lM+9+97v53Oc+x7e+9S3/2aWXXso///M/c8cddzztto4/nI6OyqiMyqg810qn06HX6x2RupxzJEky9Fk2myWbzR523XfccQfnn3/+0GevfvWrufbaa+n3+6TT6adVz4hJjMqojMqoPM3S6XTYfMIUNZaPSH2FQoFWqzX02Qc+8AE++MEPHnbdtVqNycnJoc8mJyd5/PHH+e53v8uJJ574tOoZMYlRGZVRGZWnWXq9HjWWeYQbKJI/rLqatHl+66d55JFHKBaL/vMjoUVYWaulmHVh7effr4yYxKiMyqiMyiGWIuspUjjMWoT8FovFISZxpEq1WqVWqw19tri4yPHHH0+lUnna9Yy8m0ZlVEZlVI7BcvbZZ3PzzTcPffalL32JF7/4xU/bHgEjJjEqozIqo/IDlMwR+nn6pdVqcc8993DPPfcA4uJ6zz33MDs7C8B73/tefu7nfs4/f+mll/Ltb3+bK664gm9961t8/OMf59prr+U3fuM3DqndEdw0KqMyKqNyyCWtP4dbx9Mv/+///T/OPfdc//cVV1wBwFvf+lb++q//mkcffdQzDIDNmzfzhS98gXe96138+Z//ORs3buTqq6/mZ37mZw6p3VGcxKiMyqiMytMszWaT9evXs8Kdh22TaNJiPWexsrLyb2KTOFJlpEmMyqiMyqgccvnhaxJHq4yYxKiMyqiMyiGXNIdqU3jyOp75ZWS4HpVRGZVRGZWnLCNNYlRGZVRG5ZDLCG4alVEZlVEZlacsh+7C+uR1PPPLCG4alVEZlVEZlacsI01iVEZlVEblkMvxHD5c9Owgv8+OXo7KqIzKqDyjyghuetaXj3zkI2zevJlcLscZZ5zB7bfffkTq/deuEHTO8cEPfpCNGzdywgkn8KpXvYr77rtv6Jlut8tll13G8573PPL5PK973euYm5sbeubAgQNccsklrF+/nvXr13PJJZfQaDQO6s+VV17JT/zET7Bu3TomJiZ4/etfz549e45qnz760Y9y6qmn+sRlZ599Nv/wD/9w1PrzZOXKK68kSRLe+c53HrV+ffCDHyRJkqGfarV61Odpfn6et7zlLVQqFcbGxvjxH/9x7r777qPWr02bNh00T0mS8Gu/9mtHcZ7SR+jnWVDcMViuu+46l06n3cc+9jF3//33u3e84x0un8+7b3/724dd9xe+8AX3m7/5m+6zn/2sA9wNN9ww9P1VV13l1q1b5z772c+63bt3uze+8Y3uxBNPdM1m0z9z6aWXuqmpKXfzzTe7nTt3unPPPde96EUvco8//rh/5oILLnCnnHKK+/rXv+6+/vWvu1NOOcVdeOGFB/Xn1a9+tfvEJz7h7r33XnfPPfe41772tW5mZsa1Wq2j1qfPfe5z7u///u/dnj173J49e9z73vc+l06n3b333ntU+rO23HXXXW7Tpk3u1FNPde94xzuO2jx94AMfcCeffLJ79NFH/c/i4uJR649zzi0vL7sXvOAF7ud//ufdnXfe6R5++GH35S9/2f3Lv/zLUevX4uLi0BzdfPPNDnC33nrrD70/KysrDnAr1J2jf1g/K9SlrpWVJ12LZ0o5JpnEmWee6S699NKhz0466ST3nve854i2s5ZJPPHEE65arbqrrrrKf9bpdNz69evdNddc45xzrtFouHQ67a677jr/zPz8vDvuuOPcTTfd5Jxz7v7773eA+8Y3vuGfueOOOxzgHnjgge/bp8XFRQe422677RnTJ+ec27Bhg/urv/qro96fxx57zL3whS90N998s3vlK1/pmcTR6NcHPvAB96IXvehJ+3m05und7363e/nLX/6U83e01885597xjne4rVu3uieeeOKH3h9jEg1W3BO4w/ppsPKsYBLHHNzU6/W4++67D7q27/zzz+frX//6v2nbDz/8MLVabajtbDbLK1/5St/23XffTb/fH3pm48aNnHLKKf6ZO+64g/Xr13PWWWf5Z17ykpewfv36f3UMKysrAJTL5WdEnwaDAddddx3tdpuzzz77qPfn137t13jta1/LT/3UTw19frT69dBDD7Fx40Y2b97Mm970Jvbt23dU+/O5z32OF7/4xfzH//gfmZiY4LTTTuNjH/vYUZ8nK71ej0996lP84i/+IkmSHLX+9I7Qz7OhHHNM4rvf/S6DweBJr+1bewHHkS5W//dru1arkclk2LBhw/d9ZmJi4qD6JyYmvu8YnHNcccUVvPzlL+eUU045qn3avXs3hUKBbDbLpZdeyg033MCP/diPHdU5uu6669i5cydXXnnlQd8djX6dddZZ/M3f/A1f/OIX+djHPkatVuOlL30p9Xr9qM3Tvn37+OhHP8oLX/hCvvjFL3LppZdy+eWX8zd/8zdHbZ7i8nd/93c0Gg1+/ud//hnRn+dCOWa9m57s2r5DubLvh9322mee7Pl/rZ63v/3tfPOb3+SrX/3qUe/T9u3bueeee2g0Gnz2s5/lrW99K7fddttR688jjzzCO97xDr70pS+Ry+Weso0fZr9e85rX+P/v2LGDs88+m61bt/I//sf/4CUveckPvT8ATzzxBC9+8Yv5vd/7PQBOO+007rvvPj760Y8O3VVwtPb4tddey2te8xo2btw49PkPuz99/Tmccrjv/7DKMadJPO95zyOVSj3ptX1rpY0jXcwz5fu1Xa1W6fV6HDhw4Ps+s7CwcFD9S0tLTzmGyy67jM997nPceuutTE9PH/U+ZTIZfuRHfoQXv/jFXHnllbzoRS/iT//0T49af+6++24WFxc544wzOP744zn++OO57bbbuPrqqzn++OP9O0dj7azk83l27NjBQw89dNTm6cQTT+THfuzHhj770R/9UX9PwdHc49/+9rf58pe/zC/90i/5z45Wf0Zw07O4ZDIZzjjjjIOu7bv55pt56Utf+m/a9ubNm6lWq0Nt93o9brvtNt/2GWecQTqdHnrm0Ucf5d577/XPnH322aysrHDXXXf5Z+68805WVlYOGoNzjre//e1cf/313HLLLWzevPmo9+nJinOObrd71Przkz/5k+zevdvf7HXPPffw4he/mIsvvph77rmHLVu2HPV56na7fOtb3+LEE088avP0spe97CAX6gcffJAXvOAFwNHdT5/4xCeYmJjgta99rf/smbK/j+nyQzOR/xCLucBee+217v7773fvfOc7XT6fd/v37z/suh977DG3a9cut2vXLge4D33oQ27Xrl3evfaqq65y69evd9dff73bvXu3e/Ob3/yk7njT09Puy1/+stu5c6c777zzntQd79RTT3V33HGHu+OOO9yOHTue1D3wV37lV9z69evdV77ylSE3wdXVVf/MD7tP733ve90//dM/uYcffth
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T23 amp \n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAGVCAYAAAALyQ1dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9f5xdVXkvjr8P5JyBOUlmk3CSYfLDgZgSNMDVUGOwVlQErEj5eL9XWmyqFa1cFUzBStHPrdzeXqL0Fq2lUH9gbZVe/N6rtNpLI3jVKJefArkNSGgaGEMmmeSQZM8wZ2DmEPbnj73fZ73Xs9eezGRGkzDneb32a++99rOe9ay113p+rbXXLiVJkqANbWhDG9rQhgAcc7gZaEMb2tCGNhy50FYSbWhDG9rQhkJoK4k2tKENbWhDIbSVRBva0IY2tKEQ2kqiDW1oQxvaUAhtJdGGNrShDW0ohLaSaEMb2tCGNhRCW0m0oQ1taEMbCqGtJNrQhja0oQ2F0FYSbWhDG9rQhkJoK4k2tKENbTgK4Mc//jHe+c53oqenB6VSCf/wD/9w0DwbN27EqlWrcNxxx+GUU07BX//1X0+63LaSaEMb2tCGowAajQbOPPNM3HTTTRPCf/rpp/Ebv/EbeOMb34hHH30Un/zkJ3HllVfiW9/61qTKLbU3+GtDG9rQhqMLSqUS7rjjDlx88cWFONdccw2+853v4IknnmilXX755fi///f/4r777ptwWbOmwmgb2tCGNsw0eOGFFzA2NjYttJIkQalU8tI6OjrQ0dExZdr33XcfzjvvPC/t/PPPx6233opms4lyuTwhOm0l0YY2tKENE4QXXngBJx+/CAPYNy30Zs+ejeHhYS/t05/+NK677rop0x4YGMDChQu9tIULF+LFF1/Es88+i5NOOmlCdNpKog1taEMbJghjY2MYwD48gzswF9Up0RpCA0uG/x8888wzmDt3bit9OrwIgvVSOLtg08eDtpJoQxva0IZJwlx0YS5mT5FKKn7nzp3rKYnpgu7ubgwMDHhpe/bswaxZszB//vwJ02mvbmpDG9rQhpchrFmzBnfffbeXdtddd+Gss86a8HwE0FYSbWhDG9pwCFCZpmPiMDw8jE2bNmHTpk0A0iWumzZtwvbt2wEA1157LX73d3+3hX/55Zfj5z//Oa666io88cQT+OpXv4pbb70VH//4xydVbjvc1IY2tKENk4ZydkyVxsThpz/9Kd785je37q+66ioAwHvf+1587Wtfw65du1oKAwBOPvlk3HnnnfiDP/gD/NVf/RV6enrwhS98Af/+3//7SZXb/k6iDW1oQxsmCENDQ+jq6sIgHpjynMQQhtGF1RgcHPyFzElMF7Q9iTa0oQ1tmDT88j2JwwVtJdGGNrShDZOGMiY7pxCmceRDe+K6DW1oQxvaUAhtT6INbWhDGyYN7XBTG9rQhja0oRAmv4Q1TOPIh3a4qQ1taEMb2lAIbU+iDW1oQxsmDbMw9XDR0SF+jw4u29CGNrThiIJ2uOmoh5tvvhknn3wyjjvuOKxatQo/+clPpoXuwX4hmCQJrrvuOvT09OD444/HOeecg8cff9zDGR0dxRVXXIETTzwR1WoVF110EXbs2OHh7N+/H2vXrkVXVxe6urqwdu1axHGc42f9+vX41V/9VcyZMwcLFizAxRdfjCeffPKw8nTLLbfgjDPOaG1ctmbNGvzzP//zYeMnBOvXr0epVMK6desOG1/XXXcdSqWSd3R3dx/2durv78fv/M7vYP78+ejs7MS/+3f/Dg8//PBh46u3tzfXTqVSCR/5yEcOYzuVp+k4CiB5GcLtt9+elMvl5Mtf/nLys5/9LPnYxz6WVKvV5Oc///mUad95553Jpz71qeRb3/pWAiC54447vOef+cxnkjlz5iTf+ta3ks2bNyeXXHJJctJJJyVDQ0MtnMsvvzxZtGhRcvfddyePPPJI8uY3vzk588wzkxdffLGFc8EFFyQrV65M7r333uTee+9NVq5cmVx44YU5fs4///zkb/7mb5LHHnss2bRpU/KOd7wjWbp0aTI8PHzYePrOd76T/K//9b+SJ598MnnyySeTT37yk0m5XE4ee+yxw8KPhQcffDDp7e1NzjjjjORjH/vYYWunT3/608mrX/3qZNeuXa1jz549h42fJEmSffv2Ja94xSuS973vfckDDzyQPP3008n3v//95N/+7d8OG1979uzx2ujuu+9OACQ//OEPf+n8DA4OJgCSQexNEjSndAxib0prcDD4Lo4UeFkqide97nXJ5Zdf7qWtWLEi+aM/+qNpLccqiZdeeinp7u5OPvOZz7TSXnjhhaSrqyv567/+6yRJkiSO46RcLie33357C6e/vz855phjkg0bNiRJkiQ/+9nPEgDJ/fff38K57777EgDJli1bxuVpz549CYBk48aNRwxPSZIkJ5xwQvKVr3zlsPPz3HPPJcuXL0/uvvvu5E1velNLSRwOvj796U8nZ555ZpDPw9VO11xzTfJrv/Zrhe13uN9fkiTJxz72sWTZsmXJSy+99Evnh0oixmDyEpIpHTEGjwol8bILN42NjeHhhx/O/bbvvPPOw7333vsLLfvpp5/GwMCAV3ZHRwfe9KY3tcp++OGH0Ww2PZyenh6sXLmyhXPfffehq6sLq1evbuG8/vWvR1dX10HrMDg4CACYN2/eEcHTgQMHcPvtt6PRaGDNmjWHnZ+PfOQjeMc73oFzzz3XSz9cfG3duhU9PT04+eST8Vu/9Vt46qmnDis/3/nOd3DWWWfhP/yH/4AFCxbgNa95Db785S8f9nYijI2N4Rvf+Abe//73o1QqHTZ+xqbpOBrgZacknn32WRw4cCD42z77A47pBtIfr+yBgQFUKhWccMIJ4+IsWLAgR3/BggXj1iFJElx11VX4tV/7NaxcufKw8rR582bMnj0bHR0duPzyy3HHHXfgVa961WFto9tvvx2PPPII1q9fn3t2OPhavXo1/u7v/g7f+9738OUvfxkDAwM4++yzsXfv3sPWTk899RRuueUWLF++HN/73vdw+eWX48orr8Tf/d3fHbZ2UviHf/gHxHGM973vfUcEPzMBXrarm0K/7ZvML/t+2WVbnBD+weh89KMfxb/8y7/gnnvuOew8nXrqqdi0aRPiOMa3vvUtvPe978XGjRsPGz/PPPMMPvaxj+Guu+7CcccdV1jGL5Ovt7/97a3r008/HWvWrMGyZcvwt3/7t3j961//S+cHAF566SWcddZZuP766wEAr3nNa/D444/jlltu8f5VcLj6+K233oq3v/3t6Onp8dJ/2fw0s2MqMNX8vyx42XkSJ554Io499tjgb/ustTHdwJUp45Xd3d2NsbEx7N+/f1yc3bt35+jX6/XCOlxxxRX4zne+gx/+8IdYvHjxYeepUqngla98Jc466yysX78eZ555Jv7iL/7isPHz8MMPY8+ePVi1ahVmzZqFWbNmYePGjfjCF76AWbNmtfIcjndHqFarOP3007F169bD1k4nnXQSXvWqV3lpp512Wus/BYezj//85z/H97//fXzgAx9opR0uftrhpqMYKpUKVq1alftt3913342zzz77F1r2ySefjO7ubq/ssbExbNy4sVX2qlWrUC6XPZxdu3bhsccea+GsWbMGg4ODePDBB1s4DzzwAAYHB3N1SJIEH/3oR/Htb38bP/jBD3DyyScfdp5CkCQJRkdHDxs/b33rW7F58+bWn702bdqEs846C+95z3uwadMmnHLKKYe9nUZHR/HEE0/gpJNOOmzt9IY3vCG3hPpf//Vf8YpXvALA4e1Pf/M3f4MFCxbgHe94RyvtSOnfL2v4pU2R/xKBS2BvvfXW5Gc/+1mybt26pFqtJn19fVOm/dxzzyWPPvpo8uijjyYAkhtvvDF59NFHW8trP/OZzyRdXV3Jt7/97WTz5s3Jb//2bweX4y1evDj5/ve/nzzyyCPJW97yluByvDPOOCO57777kvvuuy85/fTTg8sD/+N//I9JV1dX8qMf/chbJjgyMtLC+WXzdO211yY//vGPk6effjr5l3/5l+S
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T23 pha \n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAGICAYAAACqSCyrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOy9f5TW5XXufU3IBDJBmEAnIAqZgHPAjgRKxB8lVCpCSCZElFdKM0GxVA4ES/BIKBRfIkuDsvBg5ujBKEYWLorRM4HiqxjJWMCDEOkooGMwgEFEHDpC4KtFENDP+8e17+9ANa0KCcY+d9azgjPPPM/9vX/sH9e+9t5FACqMwiiMwiiMwjhmfOpUT6AwCqMwCqMwPn6joBwKozAKozAK4z2joBwKozAKozAK4z2joBwKozAKozAK4z2joBwKozAKozAK4z2joBwKozAKozAK4z2joBwKozAKozAK4z2joBwKozAKozAK4z3j06d6AoVRGIVRGH9M49ChQzp8+PBJ+azPfOYzatWq1Un5rJM9CsqhMAqjMArjA45Dhw7pS589Q7v125PyeR07dtT27ds/lgqioBwKozAKozA+4Dh8+LB267faqaVqo8+d0Ge9oQPqvPsyHT58uKAcCqMwCqMwPgmjjT53wsrh4z4KyqEwCqMwCuNDj09LKj4Jn/HxHQW2UmEURmEUxocenzlJrw827rrrLn35y19WmzZt1KZNG1144YV67LHHTt7jvM8oKIfCKIzCKIyP+TjzzDN16623qr6+XvX19br44ot16aWX6oUXXvi9fWdRoZ9DYRRGYRTGBxtvvPGG2rZtq0xPq41an9hn6d/UVucryzK1adPmQ/99u3btNGfOHI0ZM+aE5vG7xscb9CqMwiiMwvhYjmKdeMzBf//GG28c99OWLVuqZcuWv/Ov3nnnHf2f//N/dODAAV144YUnOIffPQqwUmEURmEUxikcnTt3Vtu2bfPXLbfc8r7ve/7559W6dWu1bNlS48aN09KlS/Wnf/qnv7d5FTyHwiiMwiiMDz2K9WECyr/7M6SdO3ceByv9Lq+he/fu2rhxo/bv36+f/exnuuqqq7R69erfm4IoxBwKozAKozA+4GiOObyoNjrtxD5Lb6qtenzkmMMll1yibt266e677z6hefyuUYCVCqMwCqMw/ggHoLfffvv39vkFWKkwCqMwCuNDjw+Xp/C7P+ODjX/4h3/Q17/+dXXu3FlvvvmmfvrTn2rVqlX6+c9/foJz+N2joBwKozAKozA+9PjDZkj/67/+q0aNGqXGxka1bdtWX/7yl/Xzn/9cgwYNOsE5/O7xiYWV5s2bpy996Utq1aqVvvKVr+j//t//e1I+98knn9TQoUPVqVMnFRUV6Z/+6Z+O+z2gG2+8UZ06ddJnP/tZDRgw4D2JKm+//bb+7u/+Tn/yJ3+iz33uc/rWt76lV1999bj37Nu3T6NGjcoZDKNGjdL+/fvfM59bbrlFffv21WmnnaYvfOELGjZsmH7961+f0jn9Z9mcf+j5vN+45ZZbVFRUpEmTJp2yed14440qKio67tWxY8dTvk67du3Sd77zHbVv314lJSXq3bu3nnnmmVM2r/Ly8vesU1FRkSZMmHAK1+kPmyH9k5/8RC+//LLefvttNTU1qa6u7veqGCRJfALHT3/6U4qLi5k/fz6/+tWv+N73vsfnPvc5duzYccKfvXz5cqZPn87PfvYzJLF06dLjfn/rrbdy2mmn8bOf/Yznn3+ev/qrv+L000/njTfeyN8zbtw4zjjjDH7xi1/w7LPP8pd/+Zf06tWLo0eP5u8ZMmQI55xzDmvXrmXt2rWcc845fPOb33zPfL72ta+xYMECGhoa2LhxI1VVVXTp0oV/+7d/O2Vzevjhh3n00Uf59a9/za9//Wv+4R/+geLiYhoaGk7JfP79WL9+PeXl5Xz5y1/me9/73ilbpx/84AdUVlbS2NiYv5qamk7ZfAB++9vf8sUvfpHRo0fz9NNPs337durq6ti2bdspm1dTU9Nxa/SLX/wCSaxcufIPPp8sy5BEpt2gt07olWm3PyvL3ncvTvX4RCqH8847j3Hjxh33sx49ejB16tST+j3/Xjm8++67dOzYkVtvvTX/2aFDh2jbti0//vGPAdi/fz/FxcX89Kc/zd+za9cuPvWpT/Hzn/8cgF/96ldI4pe//GX+nnXr1iGJF1988T+cU1NTE5JYvXr1x2ZOAJ///Oe59957T/l83nzzTSoqKvjFL37BRRddlCuHUzGvH/zgB/Tq1et953mq1unv//7v+epXv/o71+9U7x/A9773Pbp168a77777B59Ps3LYCzpyQq9Mez/WyuETBysdPnxYzzzzjAYPHnzczwcPHqy1a9f+Xr97+/bt2r1793Hf3bJlS1100UX5dz/zzDM6cuTIce/p1KmTzjnnnPw969atU9u2bXX++efn77ngggvUtm3b//QZsiyT5NT6j8Oc3nnnHf30pz/NszlP9XwmTJigqqoqXXLJJcf9/FTNa+vWrerUqZO+9KUvaeTIkfrNb35zSufz8MMP69xzz9UVV1yhL3zhC/qzP/szzZ8//5SvUxqHDx/WokWL9Dd/8zcqKio6ZfNBnz4pr4/z+MQphz179uidd95Rhw4djvt5hw4dtHv37t/rd6fP/4++e/fu3frMZz6jz3/+8//he77whS+85/O/8IUv/IfPAOh//I//oa9+9as655xzTumcflc256lco5/+9Kd69tln3zcD9VTM6/zzz9f999+vxx9/XPPnz9fu3bv153/+59q7d+8pW6ff/OY3uuuuu1RRUaHHH39c48aN08SJE3X//fefsnU6dvzTP/2T9u/fr9GjR38s5vNJHh9v1XUCo6io6Lj/Bt7zs4/Td//797zf+/+zz7n22mv13HPPac2aNad8Tr8rm/NUzWfnzp363ve+pxUrVvyHXbf+kPP6+te/nv+7Z8+euvDCC9WtWzctXLhQF1xwwR98PpL07rvv6txzz9WsWbMkSX/2Z3+mF154QXfddZeuvPLK3/mZf6gz/pOf/ERf//rX1alTp+N+/oeez+F4ncg4OV2of3/jE+c5/Mmf/IlatGjxHm3f1NT0HuviZI/ENPmPvrtjx446fPiw9u3b9x++51//9V/f8/mvv/7673yGv/u7v9PDDz+slStX6swzzzzlc/rMZz6js846S+eee65uueUW9erVSzU1NadsPs8884yampr0la98RZ/+9Kf16U9/WqtXr9b/+l//S5/+9KfzvzkVe5fG5z73OfXs2VNbt249Zet0+umnv6ccw9lnn61XXnkl/7xTMS9J2rFjh+rq6vS3f/u3+c9O1XyOnKTXx3l84pTDZz7zGX3lK1/RL37xi+N+/otf/EJ//ud//nv97i996Uvq2LHjcd99+PBhrV69Ov/ur3zlKyouLj7uPY2NjWpoaMjfc+GFFyrLMq1fvz5/z9NPP60sy97zDICuvfZaLVmyRP/8z/+sL33pS6d8Tu83iGzOUzWfgQMH6vnnn9fGjRvz17nnnqvq6mpt3LhRXbt2PeXr9Pbbb2vz5s06/fTTT9k69evX7z1U6C1btuiLX/yipFN7nhYsWKAvfOELqqqqyn/2cTnfn8jxBwt9/wFHorL+5Cc/4Ve/+hWTJk3ic5/7HC+//PIJf/abb77Jhg0b2LBhA5KYO3cuGzZsyGmyt956K23btmXJkiU8//zz/PVf//X70urOPPNM6urqePbZZ7n44ovfl1b35S9/mXXr1rFu3Tp69uz5vjS/8ePH07ZtW1atWnUc3e+tt97K3/OHntO0adN48skn2b59O8899xz/8A//wKc+9SlWrFhxSubzu8axbKVTMa/rr7+eVatW8Zvf/IZf/vKXfPOb3+S0007Lz+mpWKf169fz6U9/mh/+8Ids3bqVf/zHf6SkpIRFixadsnUCeOedd+jSpQt///d//57f/SHnk9hK25WxV5zQa7uyjzVb6ROpHAD+9//+33zxi1/kM5/5DH369MmpnSc6Vq5ciaT3vK666irAVL8f/OAHdOzYkZYtW/IXf/E
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T33\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAGVCAYAAAALyQ1dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOy9e3yW1ZXo/33VJJogBDQBSdQEh2Ab0ATagA1YgjMBh8sxcg6XcmkDDpdfI8hlRFoEQU4jDMKIablUSA+GI3AG08OlQhwJSlIxPZBUSJVQSdREISmQIKEmEZ/fH2vtZxEvHRmYovTdn8/7eZPn3c/ea6+991p7r2vA8zyPYAmWYAmWYAmWLyjXXGkAgiVYgiVYguXrW4JMIliCJViCJVi+tASZRLAES7AES7B8aQkyiWAJlmAJlmD50hJkEsESLMESLMHypSXIJIIlWIIlWILlS0uQSQRLsARLsATLl5YgkwiWYAmWYAmWLy1BJhEswRIswRIsX1qCTCJYgiVYgiVYvrQEmUSwBEuwBMs3oLz22msMHTqUzp07EwgE+PWvf/0fvvPqq6/Sq1cvrr/+erp06cLq1asvut8gkwiWYAmWYPkGlMbGRu6++25ycnK+Uv3Kykr+8R//kX79+lFaWspPfvITpk2bxtatWy+q30AwwF+wBEuwBMs3qwQCAfLz83nggQe+tM6cOXPYtm0bb731lv9sypQp/P73v+f111//yn1ddymABkuwBEuw/K2Vjz/+mObm5svSlud5BAKBVs/CwsIICwu75LZff/110tPTWz0bOHAg69ato6WlhZCQkK/UTpBJBEuwBEuwfMXy8ccfE39DDMc5dVnaa9OmDWfPnm31bMGCBTzxxBOX3Pbx48fp2LFjq2cdO3bkk08+4U9/+hO33HLLV2onyCSCJViCJVi+YmlubuY4p3iffNoScUltnaGRW89m8P7779O2bVv/+eW4Rbjy2VuK0y589vlfKkEmESzBEizBcpGlLe1oS5tLbEXIb9u2bVsxictVOnXqxPHjx1s9q62t5brrruOmm276yu0ErZuCJViCJViuwnLPPffw8ssvt3pWUFDAd77zna+sj4AgkwiWYAmWYPlPlNDL9Pnq5ezZs5SVlVFWVgaIiWtZWRnvvfceAHPnzmX8+PF+/SlTpvDuu+8yc+ZM3nrrLdavX8+6deuYPXv2RfUbFDcFS7AES7BcdAnRz6W28dXL//t//4+0tDT//5kzZwLwwx/+kF/96ld8+OGHPsMAiI+P5ze/+Q0zZszg5z//OZ07d2blypUMHz78ovoN+kkES7AES7B8xXLmzBnatWtHA29csk7iDGdpR28aGhr+S3QSl6sEbxLBEizBEiwXXf76N4krVYJMIliCJViC5aJLCBerU/jiNr7+Jai4DpZgCZZgCZYvLcGbRLAES7AEy0WXoLgpWIIlWIIlWL60XLwJ6xe38fUvQXFTsARLsARLsHxpCd4kgiVYgiVYLrpcx6WLi74Z5PebAWWwBEuwBMvXqgTFTd/48otf/IL4+Hiuv/56evXqxb59+y5Lu/9RCkHP83jiiSfo3LkzN9xwA/3796e8vLxVnaamJh5++GFuvvlmIiIiGDZsGNXV1a3qnD59mnHjxtGuXTvatWvHuHHjqK+v/xw82dnZfPe73+XGG28kOjqaBx54gCNHjlxRmFatWsVdd93lBy675557eOmll64YPF9UsrOzCQQCPPLII1cMrieeeIJAINDq06lTpyuOp5qaGsaOHctNN91EeHg4SUlJHDhw4IrBFRcX9zk8BQIBfvzjH19BPIVcps83oHhXYdm0aZMXEhLi/fKXv/T+8Ic/eNOnT/ciIiK8d99995Lb/s1vfuP99Kc/9bZu3eoBXn5+fqvfn3rqKe/GG2/0tm7d6h06dMgbOXKkd8stt3hnzpzx60yZMsWLiYnxXn75Ze/gwYNeWlqad/fdd3uffPKJX2fQoEFe9+7dvd/+9rfeb3/7W6979+7ekCFDPgfPwIEDvdzcXO/w4cNeWVmZN3jwYO+2227zzp49e8Vg2rZtm7dz507vyJEj3pEjR7yf/OQnXkhIiHf48OErAs9nS0lJiRcXF+fddddd3vTp068YnhYsWOAlJiZ6H374of+pra29YvB4nuedOnXKu/32270f/ehH3htvvOFVVlZ6//7v/+798Y9/vGJw1dbWtsLRyy+/7AFeYWHhXx2ehoYGD/AaOOl5tFzSp4GT0lZDwxfOxdelXJVMIiUlxZsyZUqrZ3feeaf32GOPXdZ+PsskPv30U69Tp07eU0895T/7+OOPvXbt2nmrV6/2PM/z6uvrvZCQEG/Tpk1+nZqaGu+aa67xdu3a5Xme5/3hD3/wAG///v1+nddff90DvLfffvsvwlRbW+sB3quvvvq1gcnzPK99+/bec889d8Xh+eijj7yuXbt6L7/8svf973/fZxJXAq4FCxZ4d9999xfCeaXwNGfOHK9v375fir8rPX+e53nTp0/37rjjDu/TTz/9q8PjmEQ9Dd6neJf0qafhG8EkrjpxU3NzMwcOHPhc2r709HR++9vf/pf2XVlZyfHjx1v1HRYWxve//32/7wMHDtDS0tKqTufOnenevbtf5/XXX6ddu3b07t3br9OnTx/atWv3H46hoaEBgA4dOnwtYDp//jybNm2isbGRe+6554rD8+Mf/5jBgwfz93//962eXym4jh49SufOnYmPj2fUqFEcO3bsisKzbds2vvOd7/A//sf/IDo6muTkZH75y19ecTy50tzcTF5eHhMmTCAQCFwxeJov0+ebUK46JvGnP/2J8+fPf2Havs8m4LjcxbX/l/o+fvw4oaGhtG/f/i/WiY6O/lz70dHRf3EMnucxc+ZM+vbtS/fu3a8oTIcOHaJNmzaEhYUxZcoU8vPz+fa3v31FcbRp0yYOHjxIdnb25367EnD17t2bDRs2sHv3bn75y19y/Phxvve973Hy5Mkrhqdjx46xatUqunbtyu7du5kyZQrTpk1jw4YNVwxPF5Zf//rX1NfX86Mf/ehrAc/fQrlqrZu+KG3fxaTs+2v3/dk6X1T/P2onKyuLN998k6KioisOU7du3SgrK6O+vp6tW7fywx/+kFdfffWKwfP+++8zffp0CgoKuP7667+0j78mXPfff7//d48ePbjnnnu44447+F//63/Rp0+fvzo8AJ9++inf+c53+NnPfgZAcnIy5eXlrFq1qlWugiu1xtetW8f9999P586dWz3/a8PTop9LKZf6/l+rXHU3iZtvvplrr732C9P2ffa0cbmLs0z5S3136tSJ5uZmTp8+/RfrnDhx4nPt19XVfekYHn74YbZt20ZhYSGxsbFXHKbQ0FD+7u/+ju985ztkZ2dz991388wzz1wxeA4cOEBtbS29evXiuuuu47rrruPVV19l5cqVXHfddf47V2LuXImIiKBHjx4cPXr0iuHplltu4dvf/narZ9/61rf8PAVXco2/++67/Pu//zsPPfSQ/+xKwRMUN32DS2hoKL169fpc2r6XX36Z733ve/+lfcfHx9OpU6dWfTc3N/Pqq6/6fffq1YuQkJBWdT788EMOHz7s17nnnntoaGigpKTEr/PGG2/Q0NDwuTF4nkdWVhYvvvgie/bsIT4+/orD9EXF8zyampquGDz33Xcfhw4d8jN7lZWV8Z3vfIcxY8ZQVlZGly5drjiempqaeOutt7jllluuGJ5SU1M/Z0JdUVHB7bffDlzZ9ZSbm0t0dDSDBw/2n31d1vdVXf5qKvK/YnEmsOvWrfP+8Ic/eI888ogXERHhVVVVXXLbH330kVdaWuqVlpZ6gLd8+XKvtLTUN6996qmnvHbt2nkvvviid+jQIW/06NFfaI4XGxvr/fu//7t38OBBb8CAAV9ojnfXXXd5r7/+uvf66697PXr0+ELzwKlTp3rt2rXz9u7d28pM8Ny5c36dvzZMc+fO9V577TWvsrLSe/PNN72f/OQn3jXXXOMVFBRcEXi+rFxo3XQl4Jo1a5a3d+9e79ixY97+/fu9IUOGeDfeeKO/Tq8EnkpKSrzrrrv
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 421.786x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"print(\"T11\")\n",
|
||
|
|
"plt.matshow(T11,cmap=\"gist_rainbow\" ,vmin=0, vmax=1)\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"print(\"T12 amp\")\n",
|
||
|
|
"plt.matshow(np.abs(T12),cmap=\"gist_rainbow\" ,vmin=0, vmax=1)\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"print(\"T13 amp\")\n",
|
||
|
|
"plt.matshow(np.abs(T13),cmap=\"gist_rainbow\" ,vmin=0, vmax=1)\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"print(\"T12 pha\")\n",
|
||
|
|
"plt.matshow(np.angle(T12),cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"print(\"T13 pha\")\n",
|
||
|
|
"plt.matshow(np.angle(T13),cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"print(\"T22\")\n",
|
||
|
|
"plt.matshow(T22,cmap=\"gist_rainbow\" ,vmin=0, vmax=1)\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"print(\"T23 amp \")\n",
|
||
|
|
"plt.matshow(np.abs(T23),cmap=\"gist_rainbow\" ,vmin=0, vmax=1)\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"print(\"T23 pha \")\n",
|
||
|
|
"plt.matshow(np.angle(T23),cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"print(\"T33\")\n",
|
||
|
|
"plt.matshow(T11,cmap=\"gist_rainbow\" ,vmin=0, vmax=1)\n",
|
||
|
|
"plt.colorbar(cmap=\"gist_rainbow\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 24,
|
||
|
|
"id": "99f573d2-5716-43ab-9dfc-0a79c056ef87",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"(array([ 65, 97, 170, 189, 190, 238, 279, 280, 342, 342, 343,\n",
|
||
|
|
" 343, 343, 343, 344, 344, 352, 353, 362, 363, 419, 481,\n",
|
||
|
|
" 482, 484, 484, 488, 777, 805, 805, 806, 813, 835, 836,\n",
|
||
|
|
" 909, 979, 980, 981, 1004, 1009, 1009, 1009, 1009, 1009, 1010,\n",
|
||
|
|
" 1010, 1010, 1010, 1010, 1011, 1014, 1015, 1015, 1015, 1016, 1016,\n",
|
||
|
|
" 1017, 1017, 1017, 1018, 1018, 1019, 1019, 1019, 1020, 1020, 1020,\n",
|
||
|
|
" 1020, 1054, 1168, 1168, 1172, 1220, 1348, 1349, 1550, 1570, 1633,\n",
|
||
|
|
" 1679, 1799, 1800, 1806, 1807, 1808, 1808, 1808, 1808, 1808, 1808,\n",
|
||
|
|
" 1809, 1809, 1809, 1809, 1809, 1809, 1809, 1810, 1810, 1810, 1810,\n",
|
||
|
|
" 1810, 1810, 1810, 1811, 1849, 1854, 1854, 1855, 1855, 1856, 1856,\n",
|
||
|
|
" 1856, 1857, 1857, 1858, 1858, 1859, 1886, 1910, 1910, 1911, 1911,\n",
|
||
|
|
" 1912, 1912, 1916, 1916, 1917, 1917, 1917, 1918, 1921, 1921, 1922,\n",
|
||
|
|
" 1922, 1923, 1923, 1956, 1957, 1957, 2071, 2072, 2072, 2074, 2082,\n",
|
||
|
|
" 2082, 2188, 2197, 2198, 2316, 2317, 2317, 2361, 2362, 2403, 2425,\n",
|
||
|
|
" 2425, 2426, 2426, 2427, 2509, 2510, 2513, 2513, 2513, 2514, 2528,\n",
|
||
|
|
" 2580, 2640, 2669, 2670, 2693, 2694, 2694, 2727, 2728, 2860, 2932,\n",
|
||
|
|
" 2933, 2937, 2940, 2941, 2996, 3006, 3028, 3029, 3029, 3051, 3053,\n",
|
||
|
|
" 3187, 3187, 3188, 3188, 3196, 3460, 3460, 3471, 3472, 3472, 3473,\n",
|
||
|
|
" 3473, 3478, 3479, 3479, 3479, 3479, 3480, 3480, 3491, 3495, 3495,\n",
|
||
|
|
" 3503, 3509, 3510, 3532, 3536, 3536, 3536, 3537, 3537, 3537, 3537,\n",
|
||
|
|
" 3538, 3538, 3538, 3539, 3582, 3621, 3631, 3631, 3682, 3683, 3686,\n",
|
||
|
|
" 3687, 3737, 3738, 3764, 3765, 3765, 3778, 3779, 3800, 3801, 3801,\n",
|
||
|
|
" 3802, 3802, 3836, 3895, 3896, 4000, 4114, 4115, 4115, 4115, 4115,\n",
|
||
|
|
" 4116, 4116, 4116, 4116, 4141, 4349, 4360, 4361, 4518, 4518, 4519,\n",
|
||
|
|
" 4519, 4519, 4520, 4676, 4676, 4691, 4737, 4737, 4738, 4738, 4739,\n",
|
||
|
|
" 4739, 4745, 4745, 4792, 4793, 4794, 4796, 4797, 4798, 4802, 4913,\n",
|
||
|
|
" 4956, 4956, 4980, 4981, 5058, 5058, 5058, 5059, 5059, 5059, 5065,\n",
|
||
|
|
" 5084, 5085, 5085, 5086, 5087, 5123, 5124, 5128, 5165, 5165, 5172,\n",
|
||
|
|
" 5172, 5173, 5173, 5174, 5174, 5191, 5191, 5191, 5192, 5192, 5192,\n",
|
||
|
|
" 5193, 5193, 5248, 5340, 5405, 5405, 5481, 5482, 5637, 5638, 5899,\n",
|
||
|
|
" 5940, 5941, 6180, 6326, 6522, 6551, 6951, 6952, 6952, 6953, 6956,\n",
|
||
|
|
" 6956, 6957, 6957, 6958, 6994, 6995, 7331, 7332, 7385, 7434, 7434,\n",
|
||
|
|
" 7435, 7435], dtype=int64), array([5284, 7047, 5340, 2449, 2449, 1343, 4415, 4415, 1344, 1345, 1344,\n",
|
||
|
|
" 1345, 1978, 1979, 1978, 1979, 1886, 1886, 1321, 1321, 1291, 1343,\n",
|
||
|
|
" 1343, 1094, 1095, 1106, 4732, 877, 878, 877, 883, 869, 869,\n",
|
||
|
|
" 6378, 6132, 6115, 6115, 2450, 2449, 2450, 2451, 2452, 2459, 2451,\n",
|
||
|
|
" 2452, 2458, 2459, 2460, 2459, 2449, 2449, 2458, 2459, 2409, 2410,\n",
|
||
|
|
" 2409, 2410, 2411, 2410, 2411, 2449, 2451, 2452, 2449, 2451, 2452,\n",
|
||
|
|
" 2453, 2437, 4138, 4143, 4139, 1843, 6836, 6836, 4435, 4458, 1210,\n",
|
||
|
|
" 1157, 3457, 3457, 50, 3456, 3435, 3439, 3440, 3447, 3455, 3456,\n",
|
||
|
|
" 3435, 3436, 3437, 3438, 3439, 3440, 3441, 3435, 3436, 3437, 3438,\n",
|
||
|
|
" 3439, 3440, 3441, 3436, 3389, 3351, 3352, 3351, 3352, 3336, 3351,\n",
|
||
|
|
" 3352, 3335, 3336, 3335, 3336, 3336, 3799, 3399, 3400, 3399, 3400,\n",
|
||
|
|
" 3399, 3400, 6614, 6615, 3364, 6614, 6615, 3364, 3390, 3391, 3390,\n",
|
||
|
|
" 3391, 3390, 3391, 3339, 3339, 3363, 5182, 5181, 5182, 5180, 4011,\n",
|
||
|
|
" 4012, 3856, 3855, 3855, 3143, 3143, 3144, 4035, 4035, 2951, 3021,\n",
|
||
|
|
" 3022, 3021, 3022, 3022, 4608, 4608, 4600, 4601, 4613, 4601, 6880,\n",
|
||
|
|
" 6802, 5441, 5482, 5482, 5490, 5489, 5490, 4029, 4029, 5867, 5865,\n",
|
||
|
|
" 5865, 5804, 1212, 1212, 502, 2433, 4994, 4994, 4995, 4765, 4760,\n",
|
||
|
|
" 3833, 3834, 3833, 3834, 1243, 4324, 4325, 4456, 4455, 4456, 4455,\n",
|
||
|
|
" 4456, 672, 672, 4441, 4442, 4443, 4442, 4443, 928, 695, 696,\n",
|
||
|
|
" 7043, 4451, 4451, 4467, 4421, 4422, 4423, 4421, 4422, 4423, 4424,\n",
|
||
|
|
" 4422, 4423, 4424, 4424, 4477, 5111, 361, 362, 1603, 5120, 564,\n",
|
||
|
|
" 564, 735, 735, 340, 340, 341, 4871, 4871, 421, 421, 422,\n",
|
||
|
|
" 421, 422, 343, 4831, 4831, 4079, 5630, 5629, 5630, 5631, 5632,\n",
|
||
|
|
" 3588, 5629, 5630, 5631, 3490, 1183, 2563, 2563, 2120, 2121, 2119,\n",
|
||
|
|
" 2120, 2121, 2120, 748, 750, 6982, 1547, 1548, 1547, 1548, 1547,\n",
|
||
|
|
" 1548, 1534, 1535, 4073, 4073, 4073, 4073, 4073, 4073, 4073, 2950,\n",
|
||
|
|
" 2981, 2982, 2961, 2961, 2982, 2987, 2988, 2982, 2984, 2988, 2975,\n",
|
||
|
|
" 2912, 2911, 2912, 2864, 2864, 3004, 3004, 2986, 2892, 2893, 2933,\n",
|
||
|
|
" 2934, 2933, 2934, 2933, 2934, 2918, 2919, 2920, 2918, 2919, 2920,\n",
|
||
|
|
" 2918, 2919, 6984, 2141, 953, 954, 6850, 6850, 1441, 1441, 726,\n",
|
||
|
|
" 6449, 6449, 3835, 6409, 6479, 5768, 6793, 6792, 6793, 6793, 6792,\n",
|
||
|
|
" 6793, 6792, 6793, 6793, 6810, 6810, 1026, 1026, 5184, 5619, 5620,\n",
|
||
|
|
" 5619, 5620], dtype=int64))\n"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "83c217d4-55b2-4e4a-95cf-814c98db113b",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "1fd36d09-f3c2-493a-8e2f-46c47f9786f3",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "12ef75ae-cdca-4791-b3e8-352a41375b43",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "d7c9af73-b57e-44eb-a362-eceb74beb126",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "433ff571-113b-42c1-b50d-92398ca64d38",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "a0fca093-b2ae-48ee-9d8e-f8d4c037db07",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "c5c4605b-0510-4999-9f33-268f311727f9",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "fdb6324f-8afa-4590-ab86-473f48fdc666",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "4a85daac-8298-48c6-ad2c-e541d9750401",
|
||
|
|
"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.8.18"
|
||
|
|
}
|
||
|
|
},
|
||
|
|
"nbformat": 4,
|
||
|
|
"nbformat_minor": 5
|
||
|
|
}
|