PCA-CNN-GRU.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "7c701c67",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 导包\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import glob\n",
    "import os\n",
    "import tensorflow as tf\n",
    "from sklearn.model_selection import train_test_split\n",
    "from sklearn.utils import shuffle\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "%matplotlib inline\n",
    "\n",
    "import os\n",
    "os.environ['CUDA_VISIBLE_DEVICES'] = '/gpu:0'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "6075ea34",
   "metadata": {},
   "outputs": [],
   "source": [
    "train_path = \"./\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "199a7a4c",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['.ipynb_checkpoints', 'data', 'GRU.h5', 'GRU=数据0在第一位.ipynb', 'GRU=数据1在第一位.ipynb', 'GRU=数据2在第一位.ipynb', 'GRU=数据3在第一位.ipynb', 'GRU=数据4在第一位.ipynb', 'GRU=数据5在第一位.ipynb', 'GRU=数据6在第一位.ipynb', 'GRU=数据7在第一位.ipynb', 'GRU=数据8在第一位.ipynb', 'GRU=数据9在第一位.ipynb', 'model.h5', 'Untitled.ipynb']\n"
     ]
    }
   ],
   "source": [
    "print(os.listdir(train_path))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "1c0169b3",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>count</th>\n",
       "      <th>gpsPdop</th>\n",
       "      <th>gpsNumSv</th>\n",
       "      <th>gpsWeek</th>\n",
       "      <th>itow</th>\n",
       "      <th>gpsSec</th>\n",
       "      <th>P_Sol_Status</th>\n",
       "      <th>Pos_Type</th>\n",
       "      <th>AvionicsTemp</th>\n",
       "      <th>MainPowerV</th>\n",
       "      <th>...</th>\n",
       "      <th>rcState</th>\n",
       "      <th>Failure_ad</th>\n",
       "      <th>leftRPM</th>\n",
       "      <th>AirHeight</th>\n",
       "      <th>Failure</th>\n",
       "      <th>WWest</th>\n",
       "      <th>WSouth</th>\n",
       "      <th>Event</th>\n",
       "      <th>Capacity</th>\n",
       "      <th>Jammer</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>140.00000</td>\n",
       "      <td>2.5</td>\n",
       "      <td>11.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.388082e+08</td>\n",
       "      <td>6270.60000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>41.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>9.500000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-3.350000</td>\n",
       "      <td>0.270000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>246.00000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>2.901218e+08</td>\n",
       "      <td>1763.40000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>30.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>180.400000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>4.630000</td>\n",
       "      <td>-3.790000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>192.00000</td>\n",
       "      <td>2.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.385378e+08</td>\n",
       "      <td>6000.20000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>38.0</td>\n",
       "      <td>24.8</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-10.600000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-0.450000</td>\n",
       "      <td>1.830000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>40.00000</td>\n",
       "      <td>3.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.386708e+08</td>\n",
       "      <td>6133.20000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.000000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>0.470000</td>\n",
       "      <td>4.070000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>183.00000</td>\n",
       "      <td>2.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.385288e+08</td>\n",
       "      <td>5991.20000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>38.0</td>\n",
       "      <td>24.8</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-10.700000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-0.370000</td>\n",
       "      <td>2.570000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13823</th>\n",
       "      <td>35.26914</td>\n",
       "      <td>2.5</td>\n",
       "      <td>7.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.390111e+08</td>\n",
       "      <td>6473.46914</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>13.473086</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>0.413852</td>\n",
       "      <td>0.904223</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13824</th>\n",
       "      <td>142.00000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>29.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>1.898014e+08</td>\n",
       "      <td>6300.60000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>29.0</td>\n",
       "      <td>24.3</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>199.400000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-1.600000</td>\n",
       "      <td>-0.010000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13825</th>\n",
       "      <td>93.00000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>28.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>1.892404e+08</td>\n",
       "      <td>5739.60000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>35.0</td>\n",
       "      <td>25.1</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>202.200000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-1.300000</td>\n",
       "      <td>1.920000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13826</th>\n",
       "      <td>189.00000</td>\n",
       "      <td>2.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.385348e+08</td>\n",
       "      <td>5997.20000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>38.0</td>\n",
       "      <td>24.8</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-10.500000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-0.460000</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13827</th>\n",
       "      <td>213.00000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>2.911130e+08</td>\n",
       "      <td>2754.60000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>23.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>179.800000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>4.050000</td>\n",
       "      <td>0.920000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>13828 rows × 65 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           count  gpsPdop  gpsNumSv  gpsWeek          itow      gpsSec  \\\n",
       "0      140.00000      2.5      11.0   2176.0  4.388082e+08  6270.60000   \n",
       "1      246.00000      1.0      31.0   2190.0  2.901218e+08  1763.40000   \n",
       "2      192.00000      2.5      10.0   2176.0  4.385378e+08  6000.20000   \n",
       "3       40.00000      3.0      11.0   2176.0  4.386708e+08  6133.20000   \n",
       "4      183.00000      2.5      10.0   2176.0  4.385288e+08  5991.20000   \n",
       "...          ...      ...       ...      ...           ...         ...   \n",
       "13823   35.26914      2.5       7.0   2176.0  4.390111e+08  6473.46914   \n",
       "13824  142.00000      1.0      29.0   2190.0  1.898014e+08  6300.60000   \n",
       "13825   93.00000      1.0      28.0   2190.0  1.892404e+08  5739.60000   \n",
       "13826  189.00000      2.5      10.0   2176.0  4.385348e+08  5997.20000   \n",
       "13827  213.00000      1.0      31.0   2190.0  2.911130e+08  2754.60000   \n",
       "\n",
       "       P_Sol_Status  Pos_Type  AvionicsTemp  MainPowerV  ...  rcState  \\\n",
       "0               0.0      16.0          41.0        24.7  ...      1.0   \n",
       "1               0.0      50.0          30.0        24.7  ...      1.0   \n",
       "2               0.0      16.0          38.0        24.8  ...      1.0   \n",
       "3               1.0      16.0          39.0        24.7  ...      1.0   \n",
       "4               0.0      16.0          38.0        24.8  ...      1.0   \n",
       "...             ...       ...           ...         ...  ...      ...   \n",
       "13823           0.0      16.0          43.0        24.7  ...      1.0   \n",
       "13824           0.0      50.0          29.0        24.3  ...      1.0   \n",
       "13825           0.0      50.0          35.0        25.1  ...      1.0   \n",
       "13826           0.0      16.0          38.0        24.8  ...      1.0   \n",
       "13827           0.0      50.0          31.0        23.7  ...      1.0   \n",
       "\n",
       "       Failure_ad  leftRPM   AirHeight    Failure     WWest    WSouth  Event  \\\n",
       "0            64.0      0.0    9.500000  3145920.0 -3.350000  0.270000    0.0   \n",
       "1            64.0      0.0  180.400000  3145920.0  4.630000 -3.790000    0.0   \n",
       "2            64.0      0.0  -10.600000  3145920.0 -0.450000  1.830000    0.0   \n",
       "3            64.0      0.0   -2.000000  3145920.0  0.470000  4.070000    0.0   \n",
       "4            64.0      0.0  -10.700000  3145920.0 -0.370000  2.570000    0.0   \n",
       "...           ...      ...         ...        ...       ...       ...    ...   \n",
       "13823        64.0      0.0   13.473086  3145920.0  0.413852  0.904223    0.0   \n",
       "13824        64.0      0.0  199.400000  3145920.0 -1.600000 -0.010000    0.0   \n",
       "13825        64.0      0.0  202.200000  3145920.0 -1.300000  1.920000    0.0   \n",
       "13826        64.0      0.0  -10.500000  3145920.0 -0.460000  2.000000    0.0   \n",
       "13827        64.0      0.0  179.800000  3145920.0  4.050000  0.920000    0.0   \n",
       "\n",
       "       Capacity  Jammer  \n",
       "0           0.0       1  \n",
       "1           0.0       0  \n",
       "2           0.0       1  \n",
       "3           0.0       1  \n",
       "4           0.0       1  \n",
       "...         ...     ...  \n",
       "13823       0.0       1  \n",
       "13824       0.0       0  \n",
       "13825       0.0       0  \n",
       "13826       0.0       1  \n",
       "13827       0.0       0  \n",
       "\n",
       "[13828 rows x 65 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# 读取文件\n",
    "\n",
    "\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/0在第一位.csv')\n",
    "data = pd.read_csv('./data/SMOTE数据拆分/1在第一位.csv')\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/2在第一位.csv')\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/3在第一位.csv')\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/4在第一位.csv')\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/5在第一位.csv')\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/6在第一位.csv')\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/7在第一位.csv')\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/8在第一位.csv')\n",
    "#dataset = pd.read_csv('./data/SMOTE数据拆分/9在第一位.csv')\n",
    "data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "4f23d598",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>count</th>\n",
       "      <th>gpsPdop</th>\n",
       "      <th>gpsNumSv</th>\n",
       "      <th>gpsWeek</th>\n",
       "      <th>itow</th>\n",
       "      <th>gpsSec</th>\n",
       "      <th>P_Sol_Status</th>\n",
       "      <th>Pos_Type</th>\n",
       "      <th>AvionicsTemp</th>\n",
       "      <th>MainPowerV</th>\n",
       "      <th>...</th>\n",
       "      <th>rcState</th>\n",
       "      <th>Failure_ad</th>\n",
       "      <th>leftRPM</th>\n",
       "      <th>AirHeight</th>\n",
       "      <th>Failure</th>\n",
       "      <th>WWest</th>\n",
       "      <th>WSouth</th>\n",
       "      <th>Event</th>\n",
       "      <th>Capacity</th>\n",
       "      <th>Jammer</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>140.0</td>\n",
       "      <td>2.5</td>\n",
       "      <td>11.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>438808200.0</td>\n",
       "      <td>6270.6</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>41.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>9.5</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-3.35</td>\n",
       "      <td>0.27</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>246.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>290121800.0</td>\n",
       "      <td>1763.4</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>30.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>180.4</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>4.63</td>\n",
       "      <td>-3.79</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>192.0</td>\n",
       "      <td>2.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>438537800.0</td>\n",
       "      <td>6000.2</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>38.0</td>\n",
       "      <td>24.8</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-10.6</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-0.45</td>\n",
       "      <td>1.83</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>40.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>438670800.0</td>\n",
       "      <td>6133.2</td>\n",
       "      <td>1.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>0.47</td>\n",
       "      <td>4.07</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>183.0</td>\n",
       "      <td>2.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>438528800.0</td>\n",
       "      <td>5991.2</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>38.0</td>\n",
       "      <td>24.8</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-10.7</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-0.37</td>\n",
       "      <td>2.57</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 65 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "   count  gpsPdop  gpsNumSv  gpsWeek         itow  gpsSec  P_Sol_Status  \\\n",
       "0  140.0      2.5      11.0   2176.0  438808200.0  6270.6           0.0   \n",
       "1  246.0      1.0      31.0   2190.0  290121800.0  1763.4           0.0   \n",
       "2  192.0      2.5      10.0   2176.0  438537800.0  6000.2           0.0   \n",
       "3   40.0      3.0      11.0   2176.0  438670800.0  6133.2           1.0   \n",
       "4  183.0      2.5      10.0   2176.0  438528800.0  5991.2           0.0   \n",
       "\n",
       "   Pos_Type  AvionicsTemp  MainPowerV  ...  rcState  Failure_ad  leftRPM  \\\n",
       "0      16.0          41.0        24.7  ...      1.0        64.0      0.0   \n",
       "1      50.0          30.0        24.7  ...      1.0        64.0      0.0   \n",
       "2      16.0          38.0        24.8  ...      1.0        64.0      0.0   \n",
       "3      16.0          39.0        24.7  ...      1.0        64.0      0.0   \n",
       "4      16.0          38.0        24.8  ...      1.0        64.0      0.0   \n",
       "\n",
       "   AirHeight    Failure  WWest  WSouth  Event  Capacity  Jammer  \n",
       "0        9.5  3145920.0  -3.35    0.27    0.0       0.0       1  \n",
       "1      180.4  3145920.0   4.63   -3.79    0.0       0.0       0  \n",
       "2      -10.6  3145920.0  -0.45    1.83    0.0       0.0       1  \n",
       "3       -2.0  3145920.0   0.47    4.07    0.0       0.0       1  \n",
       "4      -10.7  3145920.0  -0.37    2.57    0.0       0.0       1  \n",
       "\n",
       "[5 rows x 65 columns]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# 查看前5行数据\n",
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "e80698ea",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "RangeIndex(start=0, stop=13828, step=1)"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.index"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "a3f65a12",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<class 'pandas.core.frame.DataFrame'>\n",
      "RangeIndex: 13828 entries, 0 to 13827\n",
      "Data columns (total 65 columns):\n",
      " #   Column          Non-Null Count  Dtype  \n",
      "---  ------          --------------  -----  \n",
      " 0   count           13828 non-null  float64\n",
      " 1   gpsPdop         13828 non-null  float64\n",
      " 2   gpsNumSv        13828 non-null  float64\n",
      " 3   gpsWeek         13828 non-null  float64\n",
      " 4   itow            13828 non-null  float64\n",
      " 5   gpsSec          13828 non-null  float64\n",
      " 6   P_Sol_Status    13828 non-null  float64\n",
      " 7   Pos_Type        13828 non-null  float64\n",
      " 8   AvionicsTemp    13828 non-null  float64\n",
      " 9   MainPowerV      13828 non-null  float64\n",
      " 10  ServPowerV      13828 non-null  float64\n",
      " 11  powerV          13828 non-null  float64\n",
      " 12  APGlobalStatus  13828 non-null  float64\n",
      " 13  GPSStatus       13828 non-null  float64\n",
      " 14  mpStatus        13828 non-null  float64\n",
      " 15  cmd             13828 non-null  float64\n",
      " 16  Lat             13828 non-null  float64\n",
      " 17  Lon             13828 non-null  float64\n",
      " 18  H               13828 non-null  float64\n",
      " 19  tas             13828 non-null  float64\n",
      " 20  Eu[0]           13828 non-null  float64\n",
      " 21  Eu[1]           13828 non-null  float64\n",
      " 22  Eu[2]           13828 non-null  float64\n",
      " 23  acc[0]          13828 non-null  float64\n",
      " 24  acc[1]          13828 non-null  float64\n",
      " 25  acc[2]          13828 non-null  float64\n",
      " 26  pqr[0]          13828 non-null  float64\n",
      " 27  pqr[1]          13828 non-null  float64\n",
      " 28  pqr[2]          13828 non-null  float64\n",
      " 29  InertVNED[0]    13828 non-null  float64\n",
      " 30  InertVNED[1]    13828 non-null  float64\n",
      " 31  InertVNED[2]    13828 non-null  float64\n",
      " 32  StaticP         13828 non-null  float64\n",
      " 33  DynPressure     13828 non-null  float64\n",
      " 34  XMagFiled       13828 non-null  int64  \n",
      " 35  YMagFiled       13828 non-null  int64  \n",
      " 36  ZMagFiled       13828 non-null  int64  \n",
      " 37  Compass         13828 non-null  float64\n",
      " 38  aileronLeftCmd  13828 non-null  float64\n",
      " 39  elevatorCmd     13828 non-null  float64\n",
      " 40  rudderCmd       13828 non-null  float64\n",
      " 41  throttleCmd     13828 non-null  float64\n",
      " 42  latCmd          13828 non-null  float64\n",
      " 43  lonCmd          13828 non-null  float64\n",
      " 44  pedCmd          13828 non-null  float64\n",
      " 45  throtCmd        13828 non-null  float64\n",
      " 46  Aileron         13828 non-null  float64\n",
      " 47  Elevator        13828 non-null  float64\n",
      " 48  Throttle        13828 non-null  float64\n",
      " 49  Rudder          13828 non-null  float64\n",
      " 50  Action          13828 non-null  float64\n",
      " 51  photoNum        13828 non-null  float64\n",
      " 52  XaccBias        13828 non-null  float64\n",
      " 53  YaccBias        13828 non-null  float64\n",
      " 54  ZaccBias        13828 non-null  float64\n",
      " 55  rcState         13828 non-null  float64\n",
      " 56  Failure_ad      13828 non-null  float64\n",
      " 57  leftRPM         13828 non-null  float64\n",
      " 58  AirHeight       13828 non-null  float64\n",
      " 59  Failure         13828 non-null  float64\n",
      " 60  WWest           13828 non-null  float64\n",
      " 61  WSouth          13828 non-null  float64\n",
      " 62  Event           13828 non-null  float64\n",
      " 63  Capacity        13828 non-null  float64\n",
      " 64  Jammer          13828 non-null  int64  \n",
      "dtypes: float64(61), int64(4)\n",
      "memory usage: 6.9 MB\n",
      "None\n"
     ]
    }
   ],
   "source": [
    "print(data.info())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "15ca9bd3",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 去除缺失数据\n",
    "data.dropna(axis=0, how='any', inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "32b92e4d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>count</th>\n",
       "      <th>gpsPdop</th>\n",
       "      <th>gpsNumSv</th>\n",
       "      <th>gpsWeek</th>\n",
       "      <th>itow</th>\n",
       "      <th>gpsSec</th>\n",
       "      <th>P_Sol_Status</th>\n",
       "      <th>Pos_Type</th>\n",
       "      <th>AvionicsTemp</th>\n",
       "      <th>MainPowerV</th>\n",
       "      <th>...</th>\n",
       "      <th>rcState</th>\n",
       "      <th>Failure_ad</th>\n",
       "      <th>leftRPM</th>\n",
       "      <th>AirHeight</th>\n",
       "      <th>Failure</th>\n",
       "      <th>WWest</th>\n",
       "      <th>WSouth</th>\n",
       "      <th>Event</th>\n",
       "      <th>Capacity</th>\n",
       "      <th>Jammer</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>140.00000</td>\n",
       "      <td>2.5</td>\n",
       "      <td>11.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.388082e+08</td>\n",
       "      <td>6270.60000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>41.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>9.500000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-3.350000</td>\n",
       "      <td>0.270000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>246.00000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>2.901218e+08</td>\n",
       "      <td>1763.40000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>30.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>180.400000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>4.630000</td>\n",
       "      <td>-3.790000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>192.00000</td>\n",
       "      <td>2.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.385378e+08</td>\n",
       "      <td>6000.20000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>38.0</td>\n",
       "      <td>24.8</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-10.600000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-0.450000</td>\n",
       "      <td>1.830000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>40.00000</td>\n",
       "      <td>3.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.386708e+08</td>\n",
       "      <td>6133.20000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.000000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>0.470000</td>\n",
       "      <td>4.070000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>183.00000</td>\n",
       "      <td>2.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.385288e+08</td>\n",
       "      <td>5991.20000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>38.0</td>\n",
       "      <td>24.8</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-10.700000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-0.370000</td>\n",
       "      <td>2.570000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13823</th>\n",
       "      <td>35.26914</td>\n",
       "      <td>2.5</td>\n",
       "      <td>7.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.390111e+08</td>\n",
       "      <td>6473.46914</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>24.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>13.473086</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>0.413852</td>\n",
       "      <td>0.904223</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13824</th>\n",
       "      <td>142.00000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>29.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>1.898014e+08</td>\n",
       "      <td>6300.60000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>29.0</td>\n",
       "      <td>24.3</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>199.400000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-1.600000</td>\n",
       "      <td>-0.010000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13825</th>\n",
       "      <td>93.00000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>28.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>1.892404e+08</td>\n",
       "      <td>5739.60000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>35.0</td>\n",
       "      <td>25.1</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>202.200000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-1.300000</td>\n",
       "      <td>1.920000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13826</th>\n",
       "      <td>189.00000</td>\n",
       "      <td>2.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2176.0</td>\n",
       "      <td>4.385348e+08</td>\n",
       "      <td>5997.20000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>38.0</td>\n",
       "      <td>24.8</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-10.500000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>-0.460000</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13827</th>\n",
       "      <td>213.00000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>2190.0</td>\n",
       "      <td>2.911130e+08</td>\n",
       "      <td>2754.60000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>23.7</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>179.800000</td>\n",
       "      <td>3145920.0</td>\n",
       "      <td>4.050000</td>\n",
       "      <td>0.920000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>13828 rows × 65 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           count  gpsPdop  gpsNumSv  gpsWeek          itow      gpsSec  \\\n",
       "0      140.00000      2.5      11.0   2176.0  4.388082e+08  6270.60000   \n",
       "1      246.00000      1.0      31.0   2190.0  2.901218e+08  1763.40000   \n",
       "2      192.00000      2.5      10.0   2176.0  4.385378e+08  6000.20000   \n",
       "3       40.00000      3.0      11.0   2176.0  4.386708e+08  6133.20000   \n",
       "4      183.00000      2.5      10.0   2176.0  4.385288e+08  5991.20000   \n",
       "...          ...      ...       ...      ...           ...         ...   \n",
       "13823   35.26914      2.5       7.0   2176.0  4.390111e+08  6473.46914   \n",
       "13824  142.00000      1.0      29.0   2190.0  1.898014e+08  6300.60000   \n",
       "13825   93.00000      1.0      28.0   2190.0  1.892404e+08  5739.60000   \n",
       "13826  189.00000      2.5      10.0   2176.0  4.385348e+08  5997.20000   \n",
       "13827  213.00000      1.0      31.0   2190.0  2.911130e+08  2754.60000   \n",
       "\n",
       "       P_Sol_Status  Pos_Type  AvionicsTemp  MainPowerV  ...  rcState  \\\n",
       "0               0.0      16.0          41.0        24.7  ...      1.0   \n",
       "1               0.0      50.0          30.0        24.7  ...      1.0   \n",
       "2               0.0      16.0          38.0        24.8  ...      1.0   \n",
       "3               1.0      16.0          39.0        24.7  ...      1.0   \n",
       "4               0.0      16.0          38.0        24.8  ...      1.0   \n",
       "...             ...       ...           ...         ...  ...      ...   \n",
       "13823           0.0      16.0          43.0        24.7  ...      1.0   \n",
       "13824           0.0      50.0          29.0        24.3  ...      1.0   \n",
       "13825           0.0      50.0          35.0        25.1  ...      1.0   \n",
       "13826           0.0      16.0          38.0        24.8  ...      1.0   \n",
       "13827           0.0      50.0          31.0        23.7  ...      1.0   \n",
       "\n",
       "       Failure_ad  leftRPM   AirHeight    Failure     WWest    WSouth  Event  \\\n",
       "0            64.0      0.0    9.500000  3145920.0 -3.350000  0.270000    0.0   \n",
       "1            64.0      0.0  180.400000  3145920.0  4.630000 -3.790000    0.0   \n",
       "2            64.0      0.0  -10.600000  3145920.0 -0.450000  1.830000    0.0   \n",
       "3            64.0      0.0   -2.000000  3145920.0  0.470000  4.070000    0.0   \n",
       "4            64.0      0.0  -10.700000  3145920.0 -0.370000  2.570000    0.0   \n",
       "...           ...      ...         ...        ...       ...       ...    ...   \n",
       "13823        64.0      0.0   13.473086  3145920.0  0.413852  0.904223    0.0   \n",
       "13824        64.0      0.0  199.400000  3145920.0 -1.600000 -0.010000    0.0   \n",
       "13825        64.0      0.0  202.200000  3145920.0 -1.300000  1.920000    0.0   \n",
       "13826        64.0      0.0  -10.500000  3145920.0 -0.460000  2.000000    0.0   \n",
       "13827        64.0      0.0  179.800000  3145920.0  4.050000  0.920000    0.0   \n",
       "\n",
       "       Capacity  Jammer  \n",
       "0           0.0       1  \n",
       "1           0.0       0  \n",
       "2           0.0       1  \n",
       "3           0.0       1  \n",
       "4           0.0       1  \n",
       "...         ...     ...  \n",
       "13823       0.0       1  \n",
       "13824       0.0       0  \n",
       "13825       0.0       0  \n",
       "13826       0.0       1  \n",
       "13827       0.0       0  \n",
       "\n",
       "[13828 rows x 65 columns]"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "4a6fb86e",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\SUN\\anaconda3\\envs\\LSTM\\lib\\site-packages\\seaborn\\_decorators.py:43: FutureWarning: Pass the following variable as a keyword arg: x. From version 0.12, the only valid positional argument will be `data`, and passing other arguments without an explicit keyword will result in an error or misinterpretation.\n",
      "  FutureWarning\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "Text(0.5, 1.0, 'Number of  messages')"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEWCAYAAACXGLsWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAZoUlEQVR4nO3dfbRddX3n8ffH8OhDhAyBQkIJ1agFfGrSiA+rjsVKrFUYFzixVTItmsowVtca24F2pj41DtNqW1FhSlUI9QFTHAd0BoWVCo7KgBdlDAFZpIIkJpIIUgJqJOE7f+zfXRySS/YN5Nx7w32/1jrr7PPd+7fP75zcnM/Zv73P3qkqJEnalSdNdgckSVOfYSFJ6mVYSJJ6GRaSpF6GhSSpl2EhSeplWOgJLclFSf5ikp47SS5M8pMk109GH6Q9xbDQhEpyR5K7kjxloPaWJFdPYreG5WXAbwFzq2rRZHdGejwMC02GfYB3THYndleSGbvZ5Cjgjqp6YBj9kSaSYaHJ8FfAu5IctOOMJPOSVJJ9BmpXJ3lLm/53Sb6R5G+S3Jvk+0le0urrkmxKsnSH1R6S5KokW5Jck+SogXU/p827J8mtSd4wMO+iJOcn+d9JHgBeMUZ/j0hyeWu/NslbW/104OPAi5Pcn+S9fW9Kkvck+cckn2p9XZ3kWUnObq9rXZJXDSz/9CSfSLIxyQ+T/MVooCV5Znut/5Lkx0k+1+pp792mNu+7SY5r816T5DtJ7mvP9Z4d+ndakh8kuTvJf2lbia9s856U5Kwk/9zmr0wyq807oL2mu9u/2beSHNb3fmhqMSw0GUaAq4F3Pcb2LwK+C/wr4DPAJcCvA88E3gR8NMlTB5b/PeD9wCHAjcCnAdpQ2FVtHYcCbwTOS3LsQNvfBZYDTwO+PkZfPgusB44ATgE+kOSEqvoE8Dbg2qp6alW9e5yv7bXAPwAHA98BvkL3/3QO8D7g7waWXQFsa6/7hcCrgLe0ee8HrmzrmQt8pNVfBfwG8CzgIODfAne3eQ8Ap7X6a4AzkpwMkOQY4Dy69/Jw4OmtT6P+CDgZeHl7L34CfKzNW9qWP5Lu3+xtwM/G+X5oijAsNFn+HHh7ktmPoe3tVXVhVW0HPkf3IfS+qtpaVVcCv6D7AB31v6rqa1W1Ffgzum/7RwK/QzdMdGFVbauqbwOfp/vQH3VZVX2jqh6qqp8PdqKt42XAf6qqn1fVjXRbE29+DK9p1P+pqq9U1TbgH4HZwDlV9SBdKM5LclD7Zv5q4J1V9UBVbQL+BljS1vMg3TDYEa1vXx+oPw14DpCquqWqNgJU1dVVtbq91u/SBeHLW7tTgC9W1der6hd0/36DJ5b7Q+DPqmp9e5/fA5zSthAfpAuJZ1bV9qq6oaruexzvkSaBYaFJUVU3AV8CznoMze8amP5ZW9+OtcEti3UDz3s/cA/dt9+jgBe1oZF7k9xL9835l8ZqO4YjgHuqastA7Qc88hv37trxdfy4heLoY+he21HAvsDGgb7/Hd0WEsCfAAGuT7ImyR8AVNU/AR+l+9Z/V5ILkswESPKiJF9NsjnJv9BtARwy8FoH38ef8vAWCa0/Xxjoyy3AduAwui2lrwCXJNmQ5C+T7PvY3yJNBsNCk+ndwFt55Ifr6M7gJw/UBj+8H4sjRyfa8NQsYAPdh981VXXQwO2pVXXGQNtdnZZ5AzArydMGar8M/PBx9nc81gFbgUMG+j6zqo4FqKofVdVbq+oIum/95yV5Zpt3blUtAI6lG47647bOzwCXA0dW1dOB/04XOAAb6YazAEhyIN3WwmB/Xr3De3lAVf2wqh6sqvdW1THAS+i26E4bxpui4TEsNGmqai3dMNIfDdQ2033YvinJjPaN+BmP86l+O8nLkuxHN5Z/XVWto9uyeVaSNyfZt91+PcmvjrP/64BvAv+17cR9HnA6bZ/IMLWhoyuBDyWZ2XYwPyPJywGSnJpk9MP9J3Sht729vhe1b/YPAD+n2wKAbnjqnqr6eZJFdPtrRl0KvDbdwQT7Ae/l4SCBLliWpx08kGR2kpPa9CuSPLftfL+PblhqO9qrGBaabO8DnrJD7a1033bvpvv2+83H+RyfoduKuQdYQDfURBs+ehXdOP8G4EfAfwP23411vxGY19p/AXh3VV31OPs7XqcB+wE30wXCpXQ7n6Hb4X9dkvvpthbeUVW3AzOBv2/L/4DuPf5ga/Pvgfcl2UK3T2Ll6BNV1Rrg7XT7TTYCW4BNdFs3AB9uz3Nla/9/6Q5EgG7L8FK6oLgFuAb41J56EzQx4sWPJO2uNpx3LzC/hZCe4NyykDQuSV6b5MntkOMPAquBOya3V5oohoWk8TqJbrhtAzAfWFIOTUwbDkNJknq5ZSFJ6rVP/yJ7p0MOOaTmzZs32d2QpL3KDTfc8OOq2unMCk/YsJg3bx4jIyOT3Q1J2qsk+cFYdYehJEm9DAtJUi/DQpLUa2hhkeTZSW4cuN2X5J1JZqW72Mxt7f7ggTZnp7uAzK1JThyoL0h3IZi1Sc5NkrGfVZI0DEMLi6q6tapeUFUvoDsfz0/pzp1zFrCqquYDq9rj0YurLKE7F9BiurNkjl7G8nxgGd0Pgea3+ZKkCTJRw1AnAP9cVT+g+xXoilZfQXd1LVr9knYBm9uBtcCiJIcDM6vq2vZr0YsH2kiSJsBEhcUSuqtuARw2cGWujTx8sZY5PPJCM+tbbU6b3rG+kyTLkowkGdm8efMe7L4kTW9DD4t27vvX0V0icpeLjlGrXdR3LlZdUFULq2rh7NmP5WqdkqSxTMSWxauBbw9c9vKuNrREu9/U6usZuKIZ3VW5NrT63DHqkqQJMhG/4H4jDw9BQXeBlKXAOe3+soH6Z5L8Nd31fucD11fV9iRbkhwPXEd3wZePDLvTC/744mE/hfZCN/zV1Lga6J3ve+5kd0FT0C//+eqhrXuoYZHkycBv0V0DeNQ5wMokpwN3AqdCdyWuJCvprvq1DThz4EL1ZwAXAQcCV7SbJGmCDDUsquqnPPKi7lTV3XRHR421/HJg+Rj1EeC4YfRRktTPX3BLknoZFpKkXoaFJKmXYSFJ6mVYSJJ6GRaSpF6GhSSpl2EhSeplWEiSehkWkqRehoUkqZdhIUnqZVhIknoZFpKkXoaFJKmXYSFJ6mVYSJJ6GRaSpF6GhSSpl2EhSeplWEiSeg01LJIclOTSJN9LckuSFyeZleSqJLe1+4MHlj87ydoktyY5caC+IMnqNu/cJBlmvyVJjzTsLYsPA1+uqucAzwduAc4CVlXVfGBVe0ySY4AlwLHAYuC8JDPaes4HlgHz223xkPstSRowtLBIMhP4DeATAFX1i6q6FzgJWNEWWwGc3KZPAi6pqq1VdTuwFliU5HBgZlVdW1UFXDzQRpI0AYa5ZfErwGbgwiTfSfLxJE8BDquqjQDt/tC2/Bxg3UD79a02p03vWN9JkmVJRpKMbN68ec++GkmaxoYZFvsAvwacX1UvBB6gDTk9irH2Q9Qu6jsXqy6oqoVVtXD27Nm7219J0qMYZlisB9ZX1XXt8aV04XFXG1qi3W8aWP7IgfZzgQ2tPneMuiRpggwtLKrqR8C6JM9upROAm4HLgaWtthS4rE1fDixJsn+So+l2ZF/fhqq2JDm+HQV12kAbSdIE2GfI63878Okk+wHfB36fLqBWJjkduBM4FaCq1iRZSRco24Azq2p7W88ZwEXAgcAV7SZJmiBDDYuquhFYOMasEx5l+eXA8jHqI8Bxe7RzkqRx8xfckqRehoUkqZdhIUnqZVhIknoZFpKkXoaFJKmXYSFJ6mVYSJJ6GRaSpF6GhSSpl2EhSeplWEiSehkWkqRehoUkqZdhIUnqZVhIknoZFpKkXoaFJKmXYSFJ6mVYSJJ6GRaSpF5DDYskdyRZneTGJCOtNivJVUlua/cHDyx/dpK1SW5NcuJAfUFbz9ok5ybJMPstSXqkidiyeEVVvaCqFrbHZwGrqmo+sKo9JskxwBLgWGAxcF6SGa3N+cAyYH67LZ6AfkuSmskYhjoJWNGmVwAnD9QvqaqtVXU7sBZYlORwYGZVXVtVBVw80EaSNAGGHRYFXJnkhiTLWu2wqtoI0O4PbfU5wLqBtutbbU6b3rG+kyTLkowkGdm8efMefBmSNL3tM+T1v7SqNiQ5FLgqyfd2sexY+yFqF/Wdi1UXABcALFy4cMxlJEm7b6hbFlW1od1vAr4ALALuakNLtPtNbfH1wJEDzecCG1p97hh1SdIEGVpYJHlKkqeNTgOvAm4CLgeWtsWWApe16cuBJUn2T3I03Y7s69tQ1ZYkx7ejoE4baCNJmgDDHIY6DPhCO8p1H+AzVfXlJN8CViY5HbgTOBWgqtYkWQncDGwDzqyq7W1dZwAXAQcCV7SbJGmCDC0squr7wPPHqN8NnPAobZYDy8eojwDH7ek+SpLGx19wS5J6GRaSpF6GhSSpl2EhSeplWEiSehkWkqRehoUkqZdhIUnqZVhIknoZFpKkXoaFJKmXYSFJ6mVYSJJ6GRaSpF6GhSSpl2EhSeplWEiSehkWkqRehoUkqZdhIUnqZVhIknoNPSySzEjynSRfao9nJbkqyW3t/uCBZc9OsjbJrUlOHKgvSLK6zTs3SYbdb0nSwyZiy+IdwC0Dj88CVlXVfGBVe0ySY4AlwLHAYuC8JDNam/OBZcD8dls8Af2WJDVDDYskc4HXAB8fKJ8ErGjTK4CTB+qXVNXWqrodWAssSnI4MLOqrq2qAi4eaCNJmgDjCoskq8ZTG8PfAn8CPDRQO6yqNgK0+0NbfQ6wbmC59a02p03vWB+rn8uSjCQZ2bx58zi6J0kaj12GRZIDkswCDklycNvfMCvJPOCInra/A2yqqhvG2Zex9kPULuo7F6suqKqFVbVw9uzZ43xaSVKffXrm/yHwTrpguIGHP7jvAz7W0/alwOuS/DZwADAzyaeAu5IcXlUb2xDTprb8euDIgfZzgQ2tPneMuiRpguxyy6KqPlxVRwPvqqpfqaqj2+35VfXRnrZnV9XcqppHt+P6n6rqTcDlwNK22FLgsjZ9ObAkyf5JjqbbkX19G6rakuT4dhTUaQNtJEkToG/LAoCq+kiSlwDzBttU1cWP4TnPAVYmOR24Ezi1rWtNkpXAzcA24Myq2t7anAFcBBwIXNFukqQJMq6wSPIPwDOAG4HRD/DRI5N6VdXVwNVt+m7ghEdZbjmwfIz6CHDceJ5LkrTnjSssgIXAMe3QVUnSNDPe31ncBPzSMDsiSZq6xrtlcQhwc5Lrga2jxap63VB6JUmaUsYbFu8ZZickSVPbeI+GumbYHZEkTV3jPRpqCw//ano/YF/ggaqaOayOSZKmjvFuWTxt8HGSk4FFw+iQJGnqeUxnna2q/wn85p7tiiRpqhrvMNTrBx4+ie53F/7mQpKmifEeDfXageltwB1015+QJE0D491n8fvD7ogkaeoa78WP5ib5QpJNSe5K8vl2FTxJ0jQw3h3cF9KdQvwIuqvUfbHVJEnTwHjDYnZVXVhV29rtIsBL0UnSNDHesPhxkjclmdFubwLuHmbHJElTx3jD4g+ANwA/AjYCpwDu9JakaWK8h86+H1haVT8BSDIL+CBdiEiSnuDGu2XxvNGgAKiqe4AXDqdLkqSpZrxh8aQkB48+aFsW490qkSTt5cb7gf8h4JtJLqU7zccbGONa2ZKkJ6bx/oL74iQjdCcPDPD6qrp5qD2TJE0Z4z7rbFXdXFUfraqPjCcokhyQ5Pok/y/JmiTvbfVZSa5Kclu7HxzeOjvJ2iS3JjlxoL4gyeo279wk2d0XKkl67B7TKcrHaSvwm1X1fOAFwOIkxwNnAauqaj6wqj0myTHAEuBYYDFwXpIZbV3nA8uA+e22eIj9liTtYGhhUZ3728N9263ozla7otVXACe36ZOAS6pqa1XdDqwFFiU5HJhZVddWVQEXD7SRJE2AYW5Z0H7tfSOwCbiqqq4DDquqjQDt/tC2+Bxg3UDz9a02p03vWB/r+ZYlGUkysnnz5j36WiRpOhtqWFTV9qp6ATCXbivhuF0sPtZ+iNpFfaznu6CqFlbVwtmzPXWVJO0pQw2LUVV1L3A13b6Gu9rQEu1+U1tsPXDkQLO5wIZWnztGXZI0QYYWFklmJzmoTR8IvBL4Ht2pzpe2xZYCl7Xpy4ElSfZPcjTdjuzr21DVliTHt6OgThtoI0maAMP8FfbhwIp2RNOTgJVV9aUk1wIrk5wO3AmcClBVa5KsBG6mu3TrmVW1va3rDOAi4EDginaTJE2QoYVFVX2XMc4fVVV3Ayc8SpvljPHL8KoaAXa1v0OSNEQTss9CkrR3MywkSb0MC0lSL8NCktTLsJAk9TIsJEm9DAtJUi/DQpLUy7CQJPUyLCRJvQwLSVIvw0KS1MuwkCT1MiwkSb0MC0lSL8NCktTLsJAk9TIsJEm9DAtJUi/DQpLUy7CQJPUaWlgkOTLJV5PckmRNkne0+qwkVyW5rd0fPNDm7CRrk9ya5MSB+oIkq9u8c5NkWP2WJO1smFsW24D/WFW/ChwPnJnkGOAsYFVVzQdWtce0eUuAY4HFwHlJZrR1nQ8sA+a32+Ih9luStIOhhUVVbayqb7fpLcAtwBzgJGBFW2wFcHKbPgm4pKq2VtXtwFpgUZLDgZlVdW1VFXDxQBtJ0gSYkH0WSeYBLwSuAw6rqo3QBQpwaFtsDrBuoNn6VpvTpnesj/U8y5KMJBnZvHnzHn0NkjSdDT0skjwV+Dzwzqq6b1eLjlGrXdR3LlZdUFULq2rh7Nmzd7+zkqQxDTUskuxLFxSfrqr/0cp3taEl2v2mVl8PHDnQfC6wodXnjlGXJE2QYR4NFeATwC1V9dcDsy4HlrbppcBlA/UlSfZPcjTdjuzr21DVliTHt3WeNtBGkjQB9hniul8KvBlYneTGVvtT4BxgZZLTgTuBUwGqak2SlcDNdEdSnVlV21u7M4CLgAOBK9pNkjRBhhYWVfV1xt7fAHDCo7RZDiwfoz4CHLfneidJ2h3+gluS1MuwkCT1MiwkSb0MC0lSL8NCktTLsJAk9TIsJEm9DAtJUi/DQpLUy7CQJPUyLCRJvQwLSVIvw0KS1MuwkCT1MiwkSb0MC0lSL8NCktTLsJAk9TIsJEm9DAtJUi/DQpLUa2hhkeSTSTYluWmgNivJVUlua/cHD8w7O8naJLcmOXGgviDJ6jbv3CQZVp8lSWMb5pbFRcDiHWpnAauqaj6wqj0myTHAEuDY1ua8JDNam/OBZcD8dttxnZKkIRtaWFTV14B7diifBKxo0yuAkwfql1TV1qq6HVgLLEpyODCzqq6tqgIuHmgjSZogE73P4rCq2gjQ7g9t9TnAuoHl1rfanDa9Y31MSZYlGUkysnnz5j3acUmazqbKDu6x9kPULupjqqoLqmphVS2cPXv2HuucJE13Ex0Wd7WhJdr9plZfDxw5sNxcYEOrzx2jLkmaQBMdFpcDS9v0UuCygfqSJPsnOZpuR/b1bahqS5Lj21FQpw20kSRNkH2GteIknwX+NXBIkvXAu4FzgJVJTgfuBE4FqKo1SVYCNwPbgDOrantb1Rl0R1YdCFzRbpKkCTS0sKiqNz7KrBMeZfnlwPIx6iPAcXuwa5Kk3TRVdnBLkqYww0KS1MuwkCT1MiwkSb0MC0lSL8NCktTLsJAk9TIsJEm9DAtJUi/DQpLUy7CQJPUyLCRJvQwLSVIvw0KS1MuwkCT1MiwkSb0MC0lSL8NCktTLsJAk9TIsJEm9DAtJUq+9JiySLE5ya5K1Sc6a7P5I0nSyV4RFkhnAx4BXA8cAb0xyzOT2SpKmj70iLIBFwNqq+n5V/QK4BDhpkvskSdPGPpPdgXGaA6wbeLweeNGOCyVZBixrD+9PcusE9G06OAT48WR3YirIB5dOdhe0M/8+R707e2ItR41V3FvCYqx3oHYqVF0AXDD87kwvSUaqauFk90Mai3+fE2NvGYZaDxw58HgusGGS+iJJ087eEhbfAuYnOTrJfsAS4PJJ7pMkTRt7xTBUVW1L8h+ArwAzgE9W1ZpJ7tZ04tCepjL/PidAqnYa+pck6RH2lmEoSdIkMiwkSb0MC+2Sp1nRVJXkk0k2JblpsvsyHRgWelSeZkVT3EXA4snuxHRhWGhXPM2Kpqyq+hpwz2T3Y7owLLQrY51mZc4k9UXSJDIstCvjOs2KpCc+w0K74mlWJAGGhXbN06xIAgwL7UJVbQNGT7NyC7DS06xoqkjyWeBa4NlJ1ic5fbL79ETm6T4kSb3cspAk9TIsJEm9DAtJUi/DQpLUy7CQJPUyLKRxSnL/ZPdBmiyGhbSXSrJXXBZZTwyGhbQbkjw1yaok306yOslJrT4vyfeSfDzJTUk+neSVSb6R5LYki9py70myIsmVSe5I8vokf9nW9eUk+7blFiS5JskNSb6S5PBWvzrJB5JcA7xj0t4ITTuGhbR7fg78m6r6NeAVwIeSjJ5w8ZnAh4HnAc8Bfhd4GfAu4E8H1vEM4DV0p3v/FPDVqnou8DPgNS0wPgKcUlULgE8CywfaH1RVL6+qDw3pNUo7cTNW2j0BPpDkN4CH6E7Zflibd3tVrQZIsgZYVVWVZDUwb2AdV1TVg60+A/hyq48u92zgOOCqlkMzgI0D7T83hNcl7ZJhIe2e3wNmAwvaB/4dwAFt3taB5R4aePwQj/y/thWgqh5K8mA9fM6d0eUCrKmqFz9KHx543K9C2k0OQ0m75+nAphYUrwCOGsJz3ArMTvJigCT7Jjl2CM8jjZtbFtI4tCOPtgKfBr6YZAS4Efjenn6uqvpFklOAc5M8ne7/6d8CnvFXk8azzkrjkOT5wN9X1aLJ7os0GRyGknokeRvwWeA/T3ZfpMniloUkqZdbFpKkXoaFJKmXYSFJ6mVYSJJ6GRaSpF7/HwA7pDTmJ/0HAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# 样本分布\n",
    "# 以图方式表示\n",
    "sns.countplot(data[\"Jammer\"])\n",
    "plt.xlabel(\"Jammer\")\n",
    "plt.title(\"Number of  messages\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "1d11fa0d",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 重新排序\n",
    "# df = data.sample(frac=1).reset_index(drop=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "b2d1c15b",
   "metadata": {},
   "outputs": [],
   "source": [
    "df = data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "5de9e1fd",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 空值检查\n",
    "# df[df.isnull().values==True]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "e80e8402",
   "metadata": {},
   "outputs": [],
   "source": [
    "x = df.iloc[:,:-1]\n",
    "y = df.iloc[:,-1]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "eb1804e0",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 划分数据集\n",
    "from sklearn.model_selection import train_test_split\n",
    "x_train,x_test,y_train,y_test = train_test_split(x,y,test_size=0.3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "1f5865fa",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 数据归一化处理\n",
    "from sklearn.preprocessing import StandardScaler\n",
    "\n",
    "scaler = StandardScaler()\n",
    "scaler.fit(x_train)\n",
    "\n",
    "X_train = scaler.transform(x_train)\n",
    "X_test = scaler.transform(x_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "7e3c1495",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(9679, 64)\n",
      "(9679,)\n",
      "(4149, 64)\n",
      "(4149,)\n"
     ]
    }
   ],
   "source": [
    "print(X_train.shape)\n",
    "print(y_train.shape)\n",
    "print(X_test.shape)\n",
    "print(y_test.shape)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "d114a85e",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 对数据处理\n",
    "from keras.utils import np_utils\n",
    "\n",
    "X_train = X_train.reshape((-1,1,64))\n",
    "Y_train = np_utils.to_categorical(y_train)\n",
    "X_test = X_test.reshape((-1,1,64))\n",
    "Y_test = np_utils.to_categorical(y_test)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "5a6a5548",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(9679, 1, 64)\n",
      "(9679, 2)\n",
      "(4149, 1, 64)\n",
      "(4149, 2)\n"
     ]
    }
   ],
   "source": [
    "print(X_train.shape)\n",
    "print(Y_train.shape)\n",
    "print(X_test.shape)\n",
    "print(Y_test.shape)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "0ffef768",
   "metadata": {},
   "outputs": [],
   "source": [
    "#导入保存模型\n",
    "from keras import models \n",
    "from keras.models import load_model\n",
    "\n",
    "\n",
    "model = load_model('model.h5')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "e15e6533",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\SUN\\anaconda3\\envs\\LSTM\\lib\\site-packages\\tensorflow\\python\\keras\\engine\\sequential.py:455: UserWarning: `model.predict_classes()` is deprecated and will be removed after 2021-01-01. Please use instead:* `np.argmax(model.predict(x), axis=-1)`,   if your model does multi-class classification   (e.g. if it uses a `softmax` last-layer activation).* `(model.predict(x) > 0.5).astype(\"int32\")`,   if your model does binary classification   (e.g. if it uses a `sigmoid` last-layer activation).\n",
      "  warnings.warn('`model.predict_classes()` is deprecated and '\n"
     ]
    }
   ],
   "source": [
    "# 预测\n",
    "# predict_y = (model.predict(X_test) > 0.5)*1 \n",
    "# predict_y = np.around(predict_y) # np.around\n",
    "predict_y = model.predict_classes(X_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "2095b024",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([0, 0, 1, ..., 0, 0, 0], dtype=int64)"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predict_y"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "4d3653c0",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "5089     0\n",
       "1696     0\n",
       "3071     1\n",
       "10219    0\n",
       "6156     0\n",
       "        ..\n",
       "11102    0\n",
       "11052    0\n",
       "3311     0\n",
       "286      1\n",
       "13463    0\n",
       "Name: Jammer, Length: 4149, dtype: int64"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "y_test"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "fffb1a15",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "              precision    recall  f1-score   support\n",
      "\n",
      "           0     0.9955    0.9965    0.9960      2017\n",
      "           1     0.9967    0.9958    0.9962      2132\n",
      "\n",
      "    accuracy                         0.9961      4149\n",
      "   macro avg     0.9961    0.9962    0.9961      4149\n",
      "weighted avg     0.9961    0.9961    0.9961      4149\n",
      "\n"
     ]
    }
   ],
   "source": [
    "from sklearn.metrics import accuracy_score,f1_score,confusion_matrix,classification_report\n",
    "\n",
    "print(classification_report(y_test,predict_y,digits=4))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "c719fb3b",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[2010    9]\n",
      " [   7 2123]]\n"
     ]
    }
   ],
   "source": [
    "from sklearn.metrics import confusion_matrix\n",
    "\n",
    "a = confusion_matrix(predict_y, y_test)\n",
    "print(a)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "03b6bdab",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "PCA测试集准确率: 0.9961\n",
      "\n",
      "\n",
      "PCA测试集精确率: 0.9958\n",
      "\n",
      "\n",
      "PCA测试集召回率: 0.9967\n",
      "\n",
      "\n",
      "PCA测试集F1值: 0.9962\n"
     ]
    }
   ],
   "source": [
    "#模型预测\n",
    "from sklearn.metrics import accuracy_score\n",
    "from sklearn.metrics import precision_score\n",
    "from sklearn.metrics import recall_score\n",
    "from sklearn.metrics import f1_score\n",
    "\n",
    "\n",
    "accuracy_pca = accuracy_score(predict_y, y_test)\n",
    "print(\"PCA测试集准确率: {0:.4f}\".format(accuracy_pca))\n",
    "\n",
    "precision_pca =precision_score(predict_y, y_test)\n",
    "print(\"\\n\\nPCA测试集精确率: {0:.4f}\".format(precision_pca))\n",
    "\n",
    "recall_pca =recall_score(predict_y, y_test)\n",
    "print(\"\\n\\nPCA测试集召回率: {0:.4f}\".format(recall_pca))\n",
    "\n",
    "f1_pca =f1_score(predict_y, y_test)\n",
    "print(\"\\n\\nPCA测试集F1值: {0:.4f}\".format(f1_pca))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "452ccc83",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAT8AAAEHCAYAAADYj0FrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAcaUlEQVR4nO3df5RcZZ3n8fenujsdQhJISIj5BQkhiRJWgmQRRRwdVIK6BN0Bk93hh3BEWXDgrDs74MyqM3Oy6xkHXVmVERRBxwEDCAQGYRARBicIgYmQEAIJBNIkJCYoiQT653f/uLdDJVR339up6qqu+rw493TVc+99nqfgnC/Pj3ufRxGBmVmjKVS7AmZm1eDgZ2YNycHPzBqSg5+ZNSQHPzNrSM3VrkAxtYwKjTy42tWwHI6dM6XaVbAcXnhhI9u3b9f+5NE09vCIrtczXRuv//aeiFi4P+VVSm0Fv5EH03rsZ6pdDcvhV7/462pXwXI48d0L9juP6Hqd1rlnZrr2jVXfnrDfBVZITQU/MxsOBBr+I2YOfmaWj4BCU7Vrsd8c/MwsP+3XsGFNcPAzs5zc7TWzRuWWn5k1HOGWn5k1IrnlZ2YNqg5me4d/29XMhlg64ZHlGCgnabqk+yWtlbRG0iVp+nhJ90p6Nv07ruieyyWtl7RO0ilF6cdJejI9d6XUf/PUwc/M8hFJtzfLMbAu4AsR8Q7gBOAiSUcBlwH3RcRs4L70O+m5xcA8YCHwHUm9zdCrgAuA2enR72t1Dn5mll+ZWn4RsSUiHk8/7wLWAlOBRcD16WXXA6ennxcBN0ZEe0Q8D6wHjpc0GRgbESsiWZ7+h0X3lOQxPzPLKddzfhMkrSz6fnVEXF0yV2kGcCzwa2BSRGyBJEBKOjS9bCrwcNFtbWlaZ/p53/Q+OfiZWX6FzLO92yNiwNUUJI0GbgEujYid/QzXlToR/aT3ycHPzPIp87u9klpIAt+PI+KnafJWSZPTVt9kYFua3gZML7p9GrA5TZ9WIr1PHvMzs5zKOtsr4PvA2oj4etGp5cA56edzgNuL0hdLapU0k2Ri45G0i7xL0glpnmcX3VOSW35mll/5HnI+ETgLeFLSqjTti8BXgWWSzgdeBM4AiIg1kpYBT5HMFF8UEd3pfRcC1wEHAD9Ljz45+JlZfmV6vS0iHqL0eB3AyX3csxRYWiJ9JXB01rId/Mwsn+zP8NU0Bz8zy68OXm9z8DOznLyen5k1Knd7zazheD0/M2tM7vaaWaNyt9fMGpJne82s4cjdXjNrVO72mlkjGmCF+GHBwc/McklWsXfwM7NGI/peimAYcfAzs5xEoeAJDzNrQO72mllDcvAzs8ZTJ2N+w7/jbmZDSggp2zFgXtK1krZJWl2U9hNJq9JjY+/y9pJmSHq96Nw/FN1znKQnJa2XdKUyFO6Wn5nlVsYJj+uAb5FsMg5ARHyq97OkK4BXi67fEBHzS+RzFXAByZ6+dwELGWAPD7f8zCy3crX8IuJB4JU+yhBwJnDDAHWZDIyNiBURESSB9PSBynbwM7N8lOOACZJWFh0X5CjpJGBrRDxblDZT0r9LekDSSWnaVJJ9e3u1pWn9crfXzHLLMdu7PSIWDLKYJezd6tsCHBYROyQdB9wmaR6lp19ioMwd/Mwsl94Jj4qWITUDnwSO602LiHagPf38mKQNwBySlt60otunAZsHKsPdXjPLrVxjfv34EPB0ROzpzkqaKKkp/XwEMBt4LiK2ALsknZCOE54N3D5QAQ5+ZpaPQAVlOgbMSroBWAHMldQm6fz01GLeOtHxfuAJSb8BbgY+FxG9kyUXAt8D1gMbGGCmF9ztNbNBKFe3NyKW9JF+bom0W4Bb+rh+JXB0nrId/MwsN7/eZmYNZygmPIaCg5+Z5Tf8Y5+D31CJ9lfpXHcb0fEaSDS97V00T3030fk6nU/fTLzxKhp5EC1v/xPUcgDRuZvOtTfRs2szTZPm03LkqXvy6tm1mc5nlkNPJ4Xxs2k+4pS6+D9xrXlm3TrO+i973rTi+eef4399+W/4/CWXVq9StUDu9loeKtB8xEcojJ5MdLXTseoaCgcfQfe2VRQOnknz9PfRtekhutp+RcvMD0GhmebDP0jP7m3Ea7/dK6vO9XfRMvtjaMw0Otf8Ez2/W0/T+NlV+mH1a87cufz6sVUAdHd3M+vwqZx2+ieqW6kaUQ+LmQ7/XzBMaMQYCqMnJ5+bW9EBE4iOnfTseIamSccA0DTpGHp2rEuuaRpB4aDDoLD3/5+iYxd0t1MYOx1JNB365j1WOff/4j5mHjGLww8/vNpVqQ3ZX2+rWW75VUHPG7+n57WXaRkzjej4AxoxBkgCZHS+1u+90b4LtY7d812tY5KAaBV1009u5MxPlXwqoyHVQ7e3oi0/SQslrUvX2LqskmUNF9HdQefam2g54hTU3DqYHMpeJ+tfR0cH/3zncj75J2dUuyo1IevbHbUeICsW/NLXUL4NnAocBSyRdFSlyhsOoqebzqeW0TTxaJomvAMAjRi9p+UWHbtQy4H95qHWsUT7zjfzbN+1p+VolXHP3T9j/rHvYtKkSdWuSs1w8Ovf8cD6iHguIjqAG4FFFSyvpkUEnc/egUZNpHnae/akF8bPoXvrbwDo3vobCofM6TcfjRgDTa307GwjIuje9hsKh8ytaN0b3bKf3OAu7z7K9XpbNVVyzG8qsKnoexvw7n0vStf3Stb4aj2ogtWprti5iZ5tT6BRh9L++HcBaJ7xxzRPP5HOtTfT/vIq1DqWlne82bV645FvQnc79HTTveNpRhz9pxQOnEjLkR+l85nboaeLwrgjKYw7slo/q+7t3r2bX/z8Xr71ne9Wuyo1pdZbdVlUMvhlWmMrIq4GrgYojJlStwNahYMOY+RJXyp5bsQ7zy6ZPvL4S0rnNWYKrcddWLa6Wd9GjRrFS1t3VLsatcXP+Q2oDZhe9D3TGltmVtsE1EHsq+iY36PAbEkzJY0gWaJmeQXLM7MhUR+zvRVr+UVEl6SLgXuAJuDaiFhTqfLMbOjUeFzLpKIPOUfEXSTbyJlZvRAUanwmNwu/3mZmuYgk+GU5Bsyr9KblX5H0UtHm5B8tOnd5+tLEOkmnFKXn3rTcwc/McpOyHRlcR7LB+L6+ERHz0+OupEwdRTJ3MC+95zu9e3rw5qbls9OjVJ57cfAzs9yGYtPyEhYBN0ZEe0Q8T7Jfx/HetNzMhkbGVp/2b9PyiyU9kXaLx6VppV6cmIo3LTezoZA851fRTcuvAv6W5KWIvwWuAM6j7xcnvGm5mQ2FbJMZgxURW/eUJF0D3Jl+7evFCW9abmZDo5IPOadjeL0+AfTOBC8HFktqlTSTZGLjkcFuWu6Wn5nlk30md+Cskk3LP0AyNtgGfBn4gKT5JF3XjcBnASJijaRlwFNAF3BRRHSnWV1IMnN8AMmG5d603MzKK+eYX7/62LT8+/1cvxRYWiLdm5abWeX59TYza0j18Hqbg5+Z5eP1/MysEdXLen4OfmaWU+2v1ZeFg5+Z5VYHsc/Bz8zyc8vPzBqO6mQxUwc/M8vNLT8za0h1EPsc/MwsP7f8zKzxlHFhg2py8DOzXFTh9fyGioOfmeVWqIOmn4OfmeVWB7HPwc/M8pEXNjCzRlUHQ37ew8PM8ivXHh7p1pTbJK0uSvuapKfTrStvlXRwmj5D0uuSVqXHPxTdc5ykJyWtl3SlMhTu4GdmuYhkwiPLkcF1wMJ90u4Fjo6IdwLPAJcXndsQEfPT43NF6VcBF5BsajS7RJ5v4eBnZrkVlO0YSEQ8CLyyT9q/RERX+vVh9t6W8i3S3d7GRsSKiAjgh8DpA/6GgatnZlYkY5c37XlOkLSy6LggZ2nnsfdObDMl/bukBySdlKZNJdm7t1dbmtYvT3iYWW45Jnu3R8SCwZWhvyTZovLHadIW4LCI2CHpOOA2SfNIeuL7ioHyd/Azs1x6x/wqWoZ0DvBx4OS0K0tEtAPt6efHJG0A5pC09Iq7xtOAzQOV4W6vmeUmZTsGl7cWAn8BnBYRu4vSJ0pqSj8fQTKx8VxEbAF2STohneU9G7h9oHLc8jOzXMq5mKmkG4APkIwNtgFfJpndbQXuTccNH05ndt8P/I2kLqAb+FxE9E6WXEgyc3wAyRhh8ThhSQ5+ZpZbubq9EbGkRPL3+7j2FuCWPs6tBI7OU7aDn5nlVgcveAw85ifpR1nSzKxxlOsNj2rK0vKbV/wlHXA8rjLVMbNal8z2VrsW+6/P4CfpcuCLwAGSdvYmAx3A1UNQNzOrRaqPxUz77PZGxP+JiDHA1yJibHqMiYhDIuLyvu4zs/pXD93eLM/53SnpQABJfyrp65IOr3C9zKxG9XZ7y/FubzVlCX5XAbslHQP8T+AFkheHzaxBNUrLryt9vWQR8M2I+CYwprLVMrNapoxHLcsy27srnfw4Czgpne1tqWy1zKxWSfWxgVGWlt+nSF4mPi8iXiZZKuZrFa2VmdW0QkGZjlo2YPBLA94tJO/aAWwHbq1kpcystlVyYYOhkuUNj88ANwPfTZOmArdVsE5mVsNEtiXsa71rnKXbexFwIrATICKeBQ6tZKXMrIZlbPXVeOzLNOHRHhEdvdPWkprJsErqYBw7Zwq/+sVfVyJrq5Bx//HialfBcmhf92JZ8qn1x1iyyNLye0BS72tuHwZuAu6obLXMrJYVMh61LEv9LgN+CzwJfBa4KyL+sqK1MrOaJaCpoExHLcvS7f18+mDzNb0Jki5J08ysAdV4XMskS8vvnBJp55a5HmY2TCSTGeV5vU3StZK2SVpdlDZe0r2Snk3/jis6d7mk9ZLWSTqlKP04SU+m565UhsL7DH6Slki6g2SfzOVFx/3AjgF/lZnVrTIubHAdsHCftMuA+yJiNnBf+h1JRwGLSdYYXQh8p3dDI5I1CC4g2dRodok836K/bu+/keyTOQG4oih9F/DEQBmbWf0q12RvRDwoacY+yYtINjUCuB74JclubouAG9MtLJ+XtB44XtJGYGxErEjqph8CpzPAJkZ9Br+IeIFkBZf39JeBpBUR0e81ZlY/BDRnj34TJK0s+n51RAy0GPKkdDtKImKLpN7niqcCDxdd15amdaaf903vVzk2MBpZhjzMbBjJ0fLbHhELylVsibToJ71f5Qh+FXng2cxqkyr/6tpWSZPTVt9kYFua3gZML7puGrA5TZ9WIr1ftf4copnVoAq/3racN58yOQe4vSh9saRWSTNJJjYeSbvIuySdkM7ynl10T5+yLGxwcfFUc6lLBsrDzOpLuWZ7Jd0ArADmSmqTdD7wVeDDkp4FPpx+JyLWAMuAp4C7gYsiojvN6kLge8B6YAMDTHZAtm7v24BHJT0OXAvck67s3OusDHmYWZ1I9vAoT5snIpb0cerkPq5fCiwtkb4SODpP2VnW8/srkubl90kebn5W0v+WNCs9v7qf282s3giaCtmOWpapemlL7+X06ALGATdL+rsK1s3MapQy/lPLBuz2SvozkkHH7SR96j+PiE5JBeBZkh3dzKxB9G5dOdxlGfObAHwyfeh5j4jokfTxylTLzGpZQwS/iPhSP+fWlrc6ZjYc1MNipuV4yNnMGkgjdXvNzN4kan6h0iwc/MwsF7f8zKxh1cGQn4OfmeUlCjX+DF8WDn5mlotwy8/MGpGguQ4G/Rz8zCwXt/zMrGFVeDHTIeHgZ2a51UHsc/Azs3xEfSwBXw+/wcyGUnk3LZ8raVXRsVPSpZK+IumlovSPFt1TcuPyvNzyM7NcBDSVbyXndcB8gHQD8peAW4FPA9+IiL/fq+y9Ny6fAvxc0pyi5ewzc8vPzHJTxiOnk4EN+y6ft489G5dHxPMke3Ycn78oBz8zG4QK7d62GLih6PvFkp6QdG3RJmpTgU1F12TaoLwUBz8zyynbeF865jdB0sqi44KSOUojgNOAm9Kkq4BZJF3iLcAVewp/q0HtHe4xPzPLJeds7/aIWJDhulOBxyNiK0DvXwBJ1wB3pl/72rg8N7f8zCy3gpTpyGEJRV1eSZOLzn0C6N0lsuTG5YP5DW75mVk+Ku8y9pJGkWxO/tmi5L+TNJ+kS7ux91xErJHUu3F5F3tvXJ6Lg5+Z5VLuh5wjYjdwyD5pZ/VzfcmNy/Ny8DOz3LyBkZk1pOEf+hz8zGwQ6qDh5+BnZvmU8/W2anLwM7OchOqg4+vgZ2a51UHDz8HPzPJJHnUZ/tHPwc/M8hncogU1x8HPzHJz8LNcvnXlN/nBtdcQEXz6vM/w+UsurXaVrB/RsYvOF+8jOneDoOmQeTRPPIbu36+n6+VHiDd+x4g5Z1AYdSgA3bs20bV5BUQ3qInmKe+lacw0ADo23EF0vgb0UDhwCs3T3o80PF+t92yv5bJm9Wp+cO01/Ou/PcKIESM47WMLOfWjH+PI2bOrXTXriwo0TzmRwqiJRHcHHc8sozBmOho5npYZp9K56Zd7X940khFHfAy1HEjP6zvoeO4OmuadC0DLjFNQ0wgigs6Nd9Pz+w00jRu+/+3rYbZ3eP6vZxh6+um1HH/8CYwaNYrm5mZOev8fcfvtt1a7WtYPtRxIYdTE5HPTCNQ6juh8jcLI8RRGjnvL9YVRE1HLgcn1I8dDTxfR073n/kQPRM+Q1L+SKrSY6ZBy8Bsi8+YdzUMPPciOHTvYvXs3d//sLto2bRr4RqsJPe076Xl9O4VRk7Jd/+oGCgdMRIWmPWkdG5bTvvoHUGihcPCsSlV1SCjjP7Wsot1eSdcCHwe2RcTRlSyr1r39He/gC//jL/j4wg9z4OjRvPOdx9Dc7FGH4SC6O+jceDctU99X1ILrW8/rO+javIKWWaftlT5i1mlETxedL9xLzx9eomnM9D5yqG0CCrUd1zKpdMvvOmBhhcsYNs4973xWPPo4P7//QcaNH8+RRw7fMZ9GEdFN58a7aRo3h6YMrbXo+AOdG39Gy2EfotB60FvOq9BM00Ez6Xn1+UpUd2hkXMg052KmQ66iwS8iHgReqWQZw8m2bdsAePHFF7n9tp9y5uIlVa6R9Sci6HzxftQ6juZD5w98fVc7Hc/dSfPk91AY/eZCxNHdkc70QkQP3TtfQK0HV6jWQ6NCu7cNKfe7htCSM/8zr7yyg5bmFv7vld9m3Li3Dppb7YjXttDzu3Vo5CG0P30jAM1TToCeHjpfehC6XqfjuTspHDCBEbNOo3v7k0THq3S9/ChdLz8KJF1dgI7n7koegaGHwuhpNE0YvqNASbe31kPbwKoe/NLdnC4AmH7YYVWuTWXd98t/rXYVLIfC6CmMnH9RyXNNBx/xlrTmty2g+W2l9+ppnXtGWetWbeUMfZI2AruAbqArIhZIGg/8BJhBsoz9mRHxu/T6y4Hz0+v/LCLuGUy5VZ/tjYirI2JBRCyYOGFitatjZlmUv9/7wYiYX7TT22XAfRExG7gv/Y6ko0j2951HMp/wHUlNpTIcSNWDn5kNP0PwqMsi4Pr08/XA6UXpN0ZEe0Q8D6wHjh9MARUNfpJuAFYAcyW1STq/kuWZ2dAoKNtBtk3LA/gXSY8VnZ8UEVsA0r+HpulTgeIHZNvStNwqOuYXEZ7ONKtH2Rt1WTYtPzEiNks6FLhX0tM5S47MtSnibq+Z5ZIM55Wv2xsRm9O/24BbSbqxW3s3Lk//bksvbwOKnw6fBmwezO9w8DOzfDK+15vlaRhJB0oa0/sZ+AiwGlgOnJNedg5we/p5ObBYUqukmcBs4JHB/IyqP+piZsNPGR91mQTcmu4D3Az8U0TcLelRYFk6T/AicAZARKyRtAx4CugCLoqI7sEU7OBnZvmVKfpFxHPAMSXSdwAn93HPUmDp/pbt4GdmOdX+e7tZOPiZWS7D4b3dLBz8zCy/Ooh+Dn5mllutL1SahYOfmeVWB0N+Dn5mltMw2J8jCwc/M8vN3V4zazjCLT8za1B1EPsc/MxsEOog+jn4mVluHvMzs4ZUD/v2OviZWX4OfmbWaHoXMx3uHPzMLB8/5GxmjaoOYp+Dn5kNQh1EP+/hYWY5JYuZZjkGzEmaLul+SWslrZF0SZr+FUkvSVqVHh8tuudySeslrZN0ymB/hVt+ZpZLmRcz7QK+EBGPpxsZPSbp3vTcNyLi7/cqWzoKWAzMA6YAP5c0ZzD7eLjlZ2b5KeMxgIjYEhGPp593AWvpfxPyRcCNEdEeEc8D60m2uszNwc/Mcsuxb+8ESSuLjgv6zFOaARwL/DpNuljSE5KulTQuTZsKbCq6rY3+g2WfHPzMLLcc+/Zuj4gFRcfVpfPTaOAW4NKI2AlcBcwC5gNbgCt6Ly1xewzmN3jMz8zyUXlfb5PUQhL4fhwRPwWIiK1F568B7ky/tgHTi26fBmweTLlu+ZnZIJRn0E/JbuXfB9ZGxNeL0icXXfYJYHX6eTmwWFKrpJnAbOCRwfwCt/zMLJcyL2Z6InAW8KSkVWnaF4ElkuaTdGk3Ap8FiIg1kpYBT5HMFF80mJlecPAzs0EoV+yLiIf6yO6ufu5ZCizd37Id/MwsN7/ba2YNyau6mFlDcsvPzBpO0TN8w5qDn5nl5m6vmTWm4R/7HPzMLL86iH0OfmaWn8f8zKzhiGwLldY6v9trZg3JLT8zy60OGn4OfmaWnx91MbPG44eczawRlXlJq6px8DOz3NztNbOG5JafmTWkOoh9Dn5mNgh1EP0c/Mwst3oY81PEoLa8rAhJvwVeqHY9KmACsL3albBc6vW/2eERMXF/MpB0N8m/nyy2R8TC/SmvUmoq+NUrSSsjYkG162HZ+b9Z/fO7vWbWkBz8zKwhOfgNjaurXQHLzf/N6pzH/MysIbnlZ2YNycHPzBqSg1+FSVooaZ2k9ZIuq3Z9rH+SrpW0TdLqatfFKsvBr4IkNQHfBk4FjgKWSDqqurWyAVwH1ORDuVZeDn6VdTywPiKei4gO4EZgUZXrZP2IiAeBV6pdD6s8B7/KmgpsKvrelqaZWZU5+FVWqbe//WyRWQ1w8KusNmB60fdpwOYq1cXMijj4VdajwGxJMyWNABYDy6tcJzPDwa+iIqILuBi4B1gLLIuINdWtlfVH0g3ACmCupDZJ51e7TlYZfr3NzBqSW35m1pAc/MysITn4mVlDcvAzs4bk4GdmDcnBz8wakoOf5SbpA5LuTD+f1t9SXZIOlvTfhq52Ztk4+Nke6RJcuUTE8oj4aj+XHAw4+FnNcfBrEJJmSHpa0vWSnpB0s6RRkjZK+pKkh4AzJH1E0gpJj0u6SdLo9P6F6f0PAZ8syvdcSd9KP0+SdKuk36THe4GvArMkrZL0tT7q9iNJi4q+/1jSaZX892Hm4NdY5gJXR8Q7gZ282SJ7IyLeB/wc+CvgQxHxLmAl8N8ljQSuAf4TcBLwtj7yvxJ4ICKOAd4FrAEuAzZExPyI+PM+7vse8GkASQcB7wXu2q9fajYAB7/GsikifpV+/kfgfennn6R/TyBZcfpXklYB5wCHA28Hno+IZyN5H/If+8j/j4GrACKiOyJezVKpiHgAOFLSocAS4Jb0vWizimmudgVsSO37Infv99fSvwLujYglxRdJml/i3nL7EfBfSVa+Oa/CZZm55ddgDpP0nvTzEuChfc4/DJwo6UiAdExwDvA0MFPSrKJ7S7kPuDC9t0nSWGAXMCZD3a4DLgXwyjc2FBz8Gsta4BxJTwDjSbuovSLit8C5wA3pNQ8Db4+IN4ALgH9OJzxe6CP/S4APSnoSeAyYFxE7SLrRq/ua8EjL3prW7wf78wPNsvKSVg1C0gzgzog4utp1KUXSKOBJ4F1ZxwrN9odbflZ1kj5E0rX+fw58NlTc8rMhI+k/kExsFGuPiHdXoz7W2Bz8zKwhudtrZg3Jwc/MGpKDn5k1JAc/M2tI/x9pe7MTYeQO8QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "classes = list(set(y_test))\n",
    "classes.sort()\n",
    "confusion = confusion_matrix(predict_y, y_test)\n",
    "plt.imshow(confusion, cmap=plt.cm.Blues)\n",
    "indices = range(len(confusion))\n",
    "plt.xticks(indices, classes)\n",
    "plt.yticks(indices, classes)\n",
    "plt.colorbar()\n",
    "plt.xlabel('predict_y')\n",
    "plt.ylabel('y_test')\n",
    "for first_index in range(len(confusion)):\n",
    "    for second_index in range(len(confusion[first_index])):\n",
    "        plt.text(first_index, second_index, confusion[first_index][second_index])\n",
    " \n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "4f8f636d",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0 1]\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEWCAYAAAB42tAoAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAW5UlEQVR4nO3df5BlZX3n8fdnekBBflkyGB3GwLIYHF2G4AgaYwRiFIiGuGsSUEOJloQoifvLgk1FjUU2Go1uygKCE0VMmUhtImvQBVlXFH8gymBgYGTRWTAw4sZBfumI6Azf/eOeYS7N7dO3+/a5fafn/arq6nvOec55vueZO/fbz3nOeW6qCkmSZrJssQOQJE02E4UkqZWJQpLUykQhSWplopAktTJRSJJamSi0W0rymiT/a4hyFyV52zhiGock30nykub1nyT52GLHpMlnotBEaj7QHkryoyT/kuQjSfZZqONX1d9W1UuHKHdWVZ23UPX2S1JJtjbn+N0k708y1UVd0ihMFJpkr6iqfYCjgecBf9y/McnyRYlqYa1pzvHFwO8Ar1/keKTHMVFo4lXVd4Ergec0f4W/Ocm3gW8DJHl5khuT3J/k2iRH7tg3yaoklyXZkuQHSc5v1r8uyZeb10ny35J8P8kDSTYkeU6z7ZIkf9p3vDcm2ZTk3iSXJ3l637ZKclaSbye5L8kFSTLkOW4CvgIc1Xe8+ZzXYUmubtbdk+Rvkxww50aX+pgoNPGSrAJOBv6pWfWbwLHA6iRHAxcDvwc8BfggcHmSJzSXcT4N/DNwCLASuHRAFS8FfgV4JnAAvb/sfzAgjhOAdwG/DTytOe70472cXu9nTVPuZUOe4xHAi4BNzfJ8zytNjE8HngWsAv5kmBikmZgoNMk+meR+4MvANcCfNevfVVX3VtVDwBuBD1bV16pqe1V9FHgYeD5wDL0PzLdW1daq+klVfXlAPT8D9gWOAFJVt1bV9waUew1wcVV9o6oeBv4L8IIkh/SVeXdV3V9VdwKfp6+HMINvJNkK3Ap8AbiwWT+v86qqTVX12ap6uKq2AO+nd1lLmjcThSbZb1bVAVX181X1piYxANzVV+bngf/UXJ65v0ksq+h9kK4C/rmqtrVVUlVXA+cDFwD/kmRdkv0GFH06vb/id+z3I3o9j5V9Zf5f3+sfA/sAJNnYDFr/KMmL+soc3ZT5HXq9pCeNcl5JDkpyaTM4/iDwMeDAtvOXZmOi0K6of8rju4D/2iSUHT97V9XHm23PGGbQu6o+UFXPBZ5N7xLUWwcUu5veBzgASZ5E77LQd4c4/rOrap/m50vTtlVV/Xfgq8DbRzyvd9FrnyOraj/gtfQuR0nzZqLQru6vgbOSHNsMSj8pya8n2Rf4OvA94N3N+icmeeH0AyR5XrP/HsBW4CfA9gF1/R1wRpKjkjyB3qWwr1XVdxboXN4NnJnk50Y4r32BHwH3J1nJ4IQnzYmJQru0qlpP73r++cB99AaDX9ds2w68AvjXwJ3AZnqXeKbbj94H8330Li39APiLAXV9Dngb8Al6H9SHAacu4LncTG8s5q0jnNc76V3OegD4n8BlCxWfdl/xi4skSW3sUUiSWnWWKJJc3DzAdMsM25PkA83DSxua+8YlSROmyx7FJcCJLdtPAg5vfs4E/qrDWCRJ89RZoqiqLwL3thQ5Bfib5tbA64ADkjytq3gkSfOzmJOqreSxD05tbtY97onYJGfS63WwbK/9nrt8/4PGEqA0aYZ6ICI7fg3/+ERmXJhj3Y+WfTSIudU9w8Z51T1EJXM6bkvhDHg1n+PMfNyZV85W97LAz+2/Fxtu/MY9VbVi+Np3WsxEMeisBt6CVVXrgHUAa37x6Lry8195TOEdd27138C143VR05an7fOY/WYqO+A4Mxx/2LprWtnWuvtiHXR8auZtbXUzU9liYDvNFO8wdU+/u27Gsi1198c7qJ1mPO+Z2m/Ae2B63YPiHbbux9TRss/0bYPinfnf87F1P/a9MPz7+XHxzqFu+s9h1vffY2N97LbZ/z1a3wtD1D3z/7MR3wut77/Z4pq57qHeCy3xPrztEe64ZyvnnfaLnHLUykdnFZirxUwUm+lNRbDDwfSefG21x9Qynn7AXp0FJUlLxR33bOX4v/gCjzxSsxdusZi3x14OnN7c/fR84IEZJmKTJM3DVHOda/uIiaKzHkWSjwPHAQcm2Qy8A9gDoKouAq6gN3X0JnqTp53RVSyStDta1nQFJjZRVNVps2wv4M1d1S9Ju7vlTabYPuIMHD6ZLUlL1I4exbZdeIxCktShHT2KXXkwW5LUoYUazDZRSNISNTVlopAktXi0R+FgtiRpkKll9igkSS1MFJKkVk2e8PZYSdJgSZhaFm+PlSTNbCqxRyFJmtnUsvCIdz1JkmYytSwOZkuSZmaikCS1MlFIklpNLXMwW5LUYireHitJamGPQpLUyttjJUmt7FFIklo5hYckqVVvCo9HRjqGiUKSlrDecxSjHcNEIUlLmIPZkqRWDmZLklo5mC1JatXrUTiYLUmaQW8Kj9GOYaKQpCXMHoUkqdXUsrB9tCEKE4UkLWW95yjsUUiSZuADd5KkVn4fhSSp1dSUg9mSpBZTCSN2KEwUkrSULZ/022OTnJjktiSbkpw7YPv+ST6V5KYkG5Oc0WU8krS7WbZsgh+4SzIFXACcBKwGTkuyelqxNwPfrKo1wHHA+5Ls2VVMkrS7mfTvozgG2FRVt1fVT4FLgVOmlSlg3yQB9gHuBbZ1GJMk7Vampib79tiVwF19y5ubdf3OB54F3A3cDLylqh53SknOTLI+yfotW7Z0Fa8kLTlTmewH7jJg3fSx95cBNwJPB44Czk+y3+N2qlpXVWurau2KFSsWOk5JWrJ6D9xN7nMUm4FVfcsH0+s59DsDuKx6NgF3AEd0GJMk7VYmPVFcDxye5NBmgPpU4PJpZe4EfhUgyVOBXwBu7zAmSdqtLF8Wto/4VajLFyiWx6mqbUnOBq4CpoCLq2pjkrOa7RcB5wGXJLmZ3qWqc6rqnq5ikqTdzULcHttZogCoqiuAK6atu6jv9d3AS7uMQZJ2ZxP/wJ0kaXEtcwoPSVKbqWWDbkCdGxOFJC1hJgpJUisThSSp1XIThSSpzbKYKCRJLZZPmSgkSS3sUUiSWjlGIUlqtcxEIUlqY49CktTK5ygkSa0czJYktfLSkySplYPZkqRW9igkSa3sUUiSWtmjkCS1mvKuJ0lSG5+jkCS1MlFIklo5mC1JauVgtiSplVN4SJJa+Q13kqRW3h4rSWrlXU+SpFYmCklSKxOFJKmViUKS1MrBbElSK3sUkqRWJgpJUquJTxRJTkxyW5JNSc6docxxSW5MsjHJNV3GI0m7m4VIFMsXII6BkkwBFwC/BmwGrk9yeVV9s6/MAcCFwIlVdWeSg7qKR5J2R5PeozgG2FRVt1fVT4FLgVOmlXk1cFlV3QlQVd/vMB5J2u1M+l1PK4G7+pY3N+v6PRN4cpIvJLkhyemDDpTkzCTrk6zfsmVLR+FK0tIz6T2KQdHVtOXlwHOBXwdeBrwtyTMft1PVuqpaW1VrV6xYsfCRStISlYRRc0VnYxT0ehCr+pYPBu4eUOaeqtoKbE3yRWAN8K0O45Kk3cqovYqhehRJXpjks0m+leT2JHckuX2W3a4HDk9yaJI9gVOBy6eV+UfgRUmWJ9kbOBa4da4nIUma2aiJYtgexYeB/wDcAGwfZoeq2pbkbOAqYAq4uKo2Jjmr2X5RVd2a5DPABuAR4ENVdctcT0KSNLNRB7SHTRQPVNWVcz14VV0BXDFt3UXTlt8LvHeux5YkDWdcPYrPJ3kvcBnw8I6VVfWNkWqXJHVuXIni2Ob32r51BZwwUu2SpM5NLRvtBtehEkVVHT9SLZKkRTM14oMQw971tH+S9+946C3J+5LsP1rVkqRxWD5ij2LYvS8Gfgj8dvPzIPCRkWqWJI3FiHli6DGKw6rq3/UtvzPJjaNVLUkah3H1KB5K8ss7FpK8EHhopJolSWMxrik8fh/4aDMuEeBe4HWjVS1JGoex3B5bVTcCa5Ls1yw/OFKtkqSx6fT22CSvraqPJfmP09YDUFXvH6l2SVLnRr09drYexZOa3/uOVo0kabF02qOoqg82v985Ui2SpEUzNeJg9rAP3L0nyX5J9kjyuST3JHntaFVLksZhXLfHvrQZwH45vS8beibw1pFqliSNxSEH7j3S/sMmij2a3ycDH6+qe0eqVZI0Nu951ZqR9h/2OYpPJfk/9B6ye1OSFcBPRqpZkrRLGKpHUVXnAi8A1lbVz4CtwCldBiZJmgyzPUdxQlVdneTf9q3rL3JZV4FJkibDbJeeXgxcDbxiwLbCRCFJS95sz1G8o/l9xnjCkSRNmmGfo/izJAf0LT85yZ92FpUkaWIMe3vsSVV1/46FqrqP3q2ykqQlbthEMZXkCTsWkuwFPKGlvCRpiRj2OYqPAZ9L8hF6g9ivBz7aWVSSpIkx7PdRvCfJBuAl9L646LyquqrTyCRJE2HYHgXArcC2qvrfSfZOsm9V/bCrwCRJk2HYu57eCPwD8MFm1Urgkx3FJEmaIMMOZr8ZeCHwIEBVfRs4qKugJEmTY9hE8XBV/XTHQpLl9Aa1JUlL3LCJ4pokfwTsleTXgL8HPtVdWJKkSTFsojgH2ALcDPwecAXwx10FJUmaHLPe9ZRkGbChqp4D/HX3IUmSJsmsPYqqegS4KckzxhCPJGnCDPscxdOAjUm+Tu9LiwCoqt/oJCpJ0sQYNlG8s9MoJEkTq/XSU5InJvn3wG8BRwBfqaprdvzMdvAkJya5LcmmJOe2lHteku1JXjXXE5AkdWu2MYqPAmvp3e10EvC+YQ+cZAq4oNlvNXBaktUzlPtzwLmjJGkCzXbpaXVV/RuAJB8Gvj6HYx8DbKqq25v9LwVOAb45rdwfAJ8AnjeHY0uSxmS2HsXPdryoqm1zPPZK4K6+5c3NukclWQm8Erio7UBJzkyyPsn6LVu2zDEMSdIoZksUa5I82Pz8EDhyx+skD86ybwasmz7tx18C51TV9rYDVdW6qlpbVWtXrFgxS7WSpIXUeumpqqZGOPZmYFXf8sHA3dPKrAUuTQJwIHBykm1V9ckR6pUkLaC5fB/FXF0PHJ7kUOC7wKnAq/sLVNWhO14nuQT4tElCkiZLZ4miqrYlOZve3UxTwMVVtTHJWc321nEJSdJk6LJHQVVdQW8Cwf51AxNEVb2uy1gkSfMz7OyxkqTdlIlCktTKRCFJamWikCS1MlFIklqZKCRJrUwUkqRWJgpJUisThSSplYlCktTKRCFJamWikCS1MlFIklqZKCRJrUwUkqRWJgpJUisThSSplYlCktTKRCFJamWikCS1MlFIklqZKCRJrUwUkqRWJgpJUisThSSplYlCktTKRCFJamWikCS1MlFIklqZKCRJrUwUkqRWJgpJUisThSSplYlCktSq00SR5MQktyXZlOTcAdtfk2RD83NtkjVdxiNJmrvOEkWSKeAC4CRgNXBaktXTit0BvLiqjgTOA9Z1FY8kaX667FEcA2yqqtur6qfApcAp/QWq6tqquq9ZvA44uMN4JEnz0GWiWAnc1be8uVk3kzcAVw7akOTMJOuTrN+yZcsChihJmk2XiSID1tXAgsnx9BLFOYO2V9W6qlpbVWtXrFixgCFKkmazvMNjbwZW9S0fDNw9vVCSI4EPASdV1Q86jEeSNA9d9iiuBw5PcmiSPYFTgcv7CyR5BnAZ8LtV9a0OY5EkzVNnPYqq2pbkbOAqYAq4uKo2Jjmr2X4R8HbgKcCFSQC2VdXarmKSJM1dqgYOG0ystWvX1vr16xc7DEnapSS5Yb5/iPtktiSplYlCktTKRCFJamWikCS1MlFIklqZKCRJrUwUkqRWJgpJUisThSSplYlCktTKRCFJamWikCS1MlFIklqZKCRJrUwUkqRWJgpJUisThSSplYlCktTKRCFJamWikCS1MlFIklqZKCRJrUwUkqRWJgpJUisThSSplYlCktTKRCFJamWikCS1MlFIklqZKCRJrUwUkqRWJgpJUisThSSplYlCktTKRCFJatVpokhyYpLbkmxKcu6A7UnygWb7hiRHdxmPJGnuOksUSaaAC4CTgNXAaUlWTyt2EnB483Mm8FddxSNJmp8uexTHAJuq6vaq+ilwKXDKtDKnAH9TPdcBByR5WocxSZLmaHmHx14J3NW3vBk4dogyK4Hv9RdKcia9HgfAw0luWdhQd1kHAvcsdhATwrbYybbYybbY6Rfmu2OXiSID1tU8ylBV64B1AEnWV9Xa0cPb9dkWO9kWO9kWO9kWOyVZP999u7z0tBlY1bd8MHD3PMpIkhZRl4nieuDwJIcm2RM4Fbh8WpnLgdObu5+eDzxQVd+bfiBJ0uLp7NJTVW1LcjZwFTAFXFxVG5Oc1Wy/CLgCOBnYBPwYOGOIQ6/rKORdkW2xk22xk22xk22x07zbIlWPGxKQJOlRPpktSWplopAktZrYROH0HzsN0RavadpgQ5Jrk6xZjDjHYba26Cv3vCTbk7xqnPGN0zBtkeS4JDcm2ZjkmnHHOC5D/B/ZP8mnktzUtMUw46G7nCQXJ/n+TM+azftzs6om7ofe4Pf/Bf4VsCdwE7B6WpmTgSvpPYvxfOBrix33IrbFLwFPbl6ftDu3RV+5q+ndLPGqxY57Ed8XBwDfBJ7RLB+02HEvYlv8EfDnzesVwL3Anosdewdt8SvA0cAtM2yf1+fmpPYonP5jp1nboqqurar7msXr6D2PshQN874A+APgE8D3xxncmA3TFq8GLquqOwGqaqm2xzBtUcC+SQLsQy9RbBtvmN2rqi/SO7eZzOtzc1ITxUxTe8y1zFIw1/N8A72/GJaiWdsiyUrglcBFY4xrMQzzvngm8OQkX0hyQ5LTxxbdeA3TFucDz6L3QO/NwFuq6pHxhDdR5vW52eUUHqNYsOk/loChzzPJ8fQSxS93GtHiGaYt/hI4p6q29/54XLKGaYvlwHOBXwX2Ar6a5Lqq+lbXwY3ZMG3xMuBG4ATgMOCzSb5UVQ92HNukmdfn5qQmCqf/2Gmo80xyJPAh4KSq+sGYYhu3YdpiLXBpkyQOBE5Osq2qPjmWCMdn2P8j91TVVmBrki8Ca4ClliiGaYszgHdX70L9piR3AEcAXx9PiBNjXp+bk3rpyek/dpq1LZI8A7gM+N0l+Ndiv1nboqoOrapDquoQ4B+ANy3BJAHD/R/5R+BFSZYn2Zve7M23jjnOcRimLe6k17MiyVPpzaR6+1ijnAzz+tycyB5FdTf9xy5nyLZ4O/AU4MLmL+lttQRnzByyLXYLw7RFVd2a5DPABuAR4ENVteSm6B/yfXEecEmSm+ldfjmnqpbc9ONJPg4cBxyYZDPwDmAPGO1z0yk8JEmtJvXSkyRpQpgoJEmtTBSSpFYmCklSKxOFJKmViUIaoJl59sYktzSzjh6wwMf/TpIDm9c/WshjSwvNRCEN9lBVHVVVz6E3ydqbFzsgabGYKKTZfZVm4rQkhyX5TDPJ3peSHNGsf2qS/9F838FNSX6pWf/JpuzGJGcu4jlI8zaRT2ZLkyLJFL2pHz7crFoHnFVV305yLHAhvYnmPgBcU1WvbPbZpyn/+qq6N8lewPVJPrGE5+LSEmWikAbbK8mNwCHADfRmG92H3pdE/X3fzLRPaH6fAJwOUFXbgQea9X+Y5JXN61XA4YCJQrsUE4U02ENVdVSS/YFP0xujuAS4v6qOGuYASY4DXgK8oKp+nOQLwBO7CFbqkmMUUouqegD4Q+A/Aw8BdyT5LXj0+4d3fD/554Dfb9ZPJdkP2B+4r0kSR9D76klpl2OikGZRVf9E73uYTwVeA7whyU3ARnZ+5eZbgOOb2UlvAJ4NfAZYnmQDvdlLrxt37NJCcPZYSVIrexSSpFYmCklSKxOFJKmViUKS1MpEIUlqZaKQJLUyUUiSWv1/38Azi+I1IIsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "from sklearn.metrics import precision_recall_curve\n",
    "precision, recall, t =precision_recall_curve(y_test, predict_y) #y_score是预测的概率值，y_test是真实值标签值\n",
    "print(t) #t是阈值\n",
    "plt.xlabel('Recall')\n",
    "plt.ylabel('Precision')\n",
    "plt.ylim([0.0, 1.0])\n",
    "plt.xlim([0.0, 1.0])\n",
    "plt.plot(recall, precision)\n",
    "plt.title(\"Precision-Recall\")\n",
    "plt.show()"
   ]
  }
 ],
 "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.6.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}