LSTM初试遇到障碍，使用较熟悉的TextCNN。

1.基础知识：

Embedding:将词的十进制表示做向量化
起到降维增维的作用
嵌入维度数量（New Embedding维度）的一般经验法则：
embedding_dimensions =  number_of_categories**0.25
也就是说，嵌入矢量维数应该是类别数量的 4 次方根。如词汇量为 81，建议维数为 3。
 

SpatialDropout1D

在模型开始应用，会按一定比例丢弃一个特征图中的多个通道信息，增强特征的独立性，这和dropout是不同的。

正则化：
train_norm = train.apply(lambda x: (x - np.mean(x)) / (np.max(x) - np.min(x)))
 
2.实战
代码1（初版本）
import pickle

from keras.preprocessing.sequence import pad_sequences

from keras_preprocessing.text import Tokenizer

from keras.models import Sequential, Model

from keras.layers import Dense, Embedding, Activation, merge, Input, Lambda, Reshape, LSTM, RNN, CuDNNLSTM, \

    SimpleRNNCell, SpatialDropout1D, Add, Maximum

from keras.layers import Conv1D, Flatten, Dropout, MaxPool1D, GlobalAveragePooling1D, concatenate, AveragePooling1D

from keras import optimizers

from keras import regularizers

from keras.layers import BatchNormalization

from keras.callbacks import TensorBoard, EarlyStopping, ModelCheckpoint

from keras.utils import to_categorical

import time

import numpy as np

from keras import backend as K

from sklearn.model_selection import StratifiedKFold

import pickle

from sklearn.model_selection import train_test_split

from sklearn.feature_extraction.text import TfidfVectorizer

import time

import csv

import xgboost as xgb

import numpy as np

from sklearn.model_selection import StratifiedKFold

my_security_train = './my_security_train.pkl'

my_security_test = './my_security_test.pkl'

my_result = './my_result.pkl'

my_result_csv = './my_result.csv'

inputLen=100

# config = K.tf.ConfigProto()

# # 程序按需申请内存

# config.gpu_options.allow_growth = True

# session = K.tf.Session(config = config)

# 读取文件到变量中

with open(my_security_train, 'rb') as f:

    train_labels = pickle.load(f)

    train_apis = pickle.load(f)

with open(my_security_test, 'rb') as f:

    test_files = pickle.load(f)

    test_apis = pickle.load(f)

# print(time.strftime("%Y-%m-%d-%H-%M-%S", time.localtime()))

# tensorboard = TensorBoard('./Logs/', write_images=1, histogram_freq=1)

# print(train_labels)

# 将标签转换为空格相隔的一维数组

train_labels = np.asarray(train_labels)

# print(train_labels)

tokenizer = Tokenizer(num_words=None,

                      filters='!"#$%&()*+,-./:;<=>?@[\]^_`{|}~\t\n',

                      lower=True,

                      split=" ",

                      char_level=False)

# print(train_apis)

# 通过训练和测试数据集丰富取词器的字典，方便后续操作

tokenizer.fit_on_texts(train_apis)

# print(train_apis)

# print(test_apis)

tokenizer.fit_on_texts(test_apis)

# print(test_apis)

# print(tokenizer.word_index)

# #获取目前提取词的字典信息

# # vocal = tokenizer.word_index

train_apis = tokenizer.texts_to_sequences(train_apis)

# 通过字典信息将字符转换为对应的数字

test_apis = tokenizer.texts_to_sequences(test_apis)

# print(test_apis)

# 序列化原数组为没有逗号的数组，默认在前面填充,默认截断前面的

train_apis = pad_sequences(train_apis, inputLen, padding='post', truncating='post')

# print(test_apis)

test_apis = pad_sequences(test_apis, inputLen, padding='post', truncating='post')

# print(test_apis)

def SequenceModel():

    # Sequential()是序列模型，其实是堆叠模型，可以在它上面堆砌网络形成一个复杂的网络结构

    model = Sequential()

    model.add(Dense(32, activation='relu', input_dim=6000))

    model.add(Dense(8, activation='softmax'))

    return model

def lstm():

    my_inpuy = Input(shape = (6000,), dtype = 'float64')

    #在网络第一层，起降维的作用

    emb = Embedding(len(tokenizer.word_index)+1, 256, input_length=6000)

    emb = emb(my_inpuy)

    net = Conv1D(16, 3, padding='same', kernel_initializer='glorot_uniform')(emb)

    net = BatchNormalization()(net)

    net = Activation('relu')(net)

    net = Conv1D(32, 3, padding='same', kernel_initializer='glorot_uniform')(net)

    net = BatchNormalization()(net)

    net = Activation('relu')(net)

    net = MaxPool1D(pool_size=4)(net)

    net1 = Conv1D(16, 4, padding='same', kernel_initializer='glorot_uniform')(emb)

    net1 = BatchNormalization()(net1)

    net1 = Activation('relu')(net1)

    net1 = Conv1D(32, 4, padding='same', kernel_initializer='glorot_uniform')(net1)

    net1 = BatchNormalization()(net1)

    net1 = Activation('relu')(net1)

    net1 = MaxPool1D(pool_size=4)(net1)

    net2 = Conv1D(16, 5, padding='same', kernel_initializer='glorot_uniform')(emb)

    net2 = BatchNormalization()(net2)

    net2 = Activation('relu')(net2)

    net2 = Conv1D(32, 5, padding='same', kernel_initializer='glorot_uniform')(net2)

    net2 = BatchNormalization()(net2)

    net2 = Activation('relu')(net2)

    net2 = MaxPool1D(pool_size=4)(net2)

    net = concatenate([net, net1, net2], axis=-1)

    net = CuDNNLSTM(256)(net)

    net = Dense(8, activation = 'softmax')(net)

    model = Model(inputs=my_inpuy, outputs=net)

    return model

def textcnn():

    kernel_size = [1, 3, 3, 5, 5]

    acti = 'relu'

    #可看做一个文件的api集为一句话，然后话中的词总量是6000

    my_input = Input(shape=(inputLen,), dtype='int32')

    emb = Embedding(len(tokenizer.word_index) + 1, 5, input_length=inputLen)(my_input)

    emb = SpatialDropout1D(0.2)(emb)

    net = []

    for kernel in kernel_size:

        # 32个卷积核

        con = Conv1D(32, kernel, activation=acti, padding="same")(emb)

        # 滑动窗口大小是2,默认输出最后一维是通道数

        con = MaxPool1D(2)(con)

        net.append(con)

    # print(net)

    # input()

    net = concatenate(net, axis =-1)

    # net = concatenate(net)

    # print(net)

    # input()

    net = Flatten()(net)

    net = Dropout(0.5)(net)

    net = Dense(256, activation='relu')(net)

    net = Dropout(0.5)(net)

    net = Dense(8, activation='softmax')(net)

    model = Model(inputs=my_input, outputs=net)

    return model

# model = SequenceModel()

model = textcnn()

# metrics默认只有loss，加accuracy后在model.evaluate(...)的返回值即有accuracy结果

model.compile(optimizer='adam', loss='sparse_categorical_crossentropy', metrics=['accuracy'])

# print(train_apis.shape)

# print(train_labels.shape)

# 将训练集切分成训练和验证集

skf = StratifiedKFold(n_splits=5)

for i, (train_index, valid_index) in enumerate(skf.split(train_apis, train_labels)):

    model.fit(train_apis[train_index], train_labels[train_index], epochs=10, batch_size=1000,

              validation_data=(train_apis[valid_index], train_labels[valid_index]))

    print(train_index, valid_index)

# loss, acc = model.evaluate(train_apis, train_labels)

# print(loss)

# print(acc)

# print(model.predict(train_apis))

test_apis = model.predict(test_apis)

# print(test_files)

# print(test_apis)

with open(my_result, 'wb') as f:

    pickle.dump(test_files, f)

    pickle.dump(test_apis, f)

# print(len(test_files))

# print(len(test_apis))

result = []

for i in range(len(test_files)):

    # #     print(test_files[i])

    #     #之前test_apis不带逗号的格式是矩阵格式，现在tolist转为带逗号的列表格式

    #     print(test_apis[i])

    #     print(test_apis[i].tolist())

    #     result.append(test_files[i])

    #     result.append(test_apis[i])

    tmp = []

    a = test_apis[i].tolist()

    tmp.append(test_files[i])

    # extend相比于append可以添加多个值

    tmp.extend(a)

    #     print(tmp)

    result.append(tmp)

# print(1)

# print(result)

with open(my_result_csv, 'w') as f:

    #     f.write([1,2,3])

    result_csv = csv.writer(f)

    result_csv.writerow(["file_id", "prob0", "prob1", "prob2", "prob3", "prob4", "prob5", "prob6", "prob7"])

    result_csv.writerows(result)

确定好它的原始文件api序列最大长度：13264587

增加早停机制和模型保存：
EarlyStopping的restore_best_weights参数： 如果要在停止后保存最佳权重,请将此参数设置为True。
 
代码2(加载pkl数据和训练测试部分)：
import pickle

from keras.preprocessing.sequence import pad_sequences

from keras_preprocessing.text import Tokenizer

from keras.models import Sequential, Model

from keras.layers import Dense, Embedding, Activation, merge, Input, Lambda, Reshape, LSTM, RNN, CuDNNLSTM, \

    SimpleRNNCell, SpatialDropout1D, Add, Maximum

from keras.layers import Conv1D, Flatten, Dropout, MaxPool1D, GlobalAveragePooling1D, concatenate, AveragePooling1D

from keras import optimizers

from keras import regularizers

from keras.layers import BatchNormalization

from keras.callbacks import TensorBoard, EarlyStopping, ModelCheckpoint

from keras.utils import to_categorical

import time

import numpy as np

from keras import backend as K

from sklearn.model_selection import StratifiedKFold

import pickle

from sklearn.model_selection import train_test_split

from sklearn.feature_extraction.text import TfidfVectorizer

import time

import csv

import xgboost as xgb

import numpy as np

from sklearn.model_selection import StratifiedKFold

my_security_train = './my_security_train.pkl'

my_security_test = './my_security_test.pkl'

my_result = './my_result1.pkl'

my_result_csv = './my_result1.csv'

inputLen = 5000

# config = K.tf.ConfigProto()

# # 程序按需申请内存

# config.gpu_options.allow_growth = True

# session = K.tf.Session(config = config)

# 读取文件到变量中

with open(my_security_train, 'rb') as f:

    train_labels = pickle.load(f)

    train_apis = pickle.load(f)

with open(my_security_test, 'rb') as f:

    test_files = pickle.load(f)

    test_apis = pickle.load(f)

# print(time.strftime("%Y-%m-%d-%H-%M-%S", time.localtime()))

# tensorboard = TensorBoard('./Logs/', write_images=1, histogram_freq=1)

# print(train_labels)

# 将标签转换为空格相隔的一维数组

train_labels = np.asarray(train_labels)

# print(train_labels)

tokenizer = Tokenizer(num_words=None,

                      filters='!"#$%&()*+,-./:;<=>?@[\]^_`{|}~\t\n',

                      lower=True,

                      split=" ",

                      char_level=False)

# print(train_apis)

# 通过训练和测试数据集丰富取词器的字典，方便后续操作

tokenizer.fit_on_texts(train_apis)

# print(train_apis)

# print(test_apis)

tokenizer.fit_on_texts(test_apis)

# print(test_apis)

# print(tokenizer.word_index)

# #获取目前提取词的字典信息

# # vocal = tokenizer.word_index

train_apis = tokenizer.texts_to_sequences(train_apis)

# 通过字典信息将字符转换为对应的数字

test_apis = tokenizer.texts_to_sequences(test_apis)

# print(test_apis)

# 序列化原数组为没有逗号的数组，默认在前面填充,默认截断前面的

train_apis = pad_sequences(train_apis, inputLen, padding='post', truncating='post')

# print(test_apis)

test_apis = pad_sequences(test_apis, inputLen, padding='post', truncating='post')

# print(test_apis)

def SequenceModel():

    # Sequential()是序列模型，其实是堆叠模型，可以在它上面堆砌网络形成一个复杂的网络结构

    model = Sequential()

    model.add(Dense(32, activation='relu', input_dim=6000))

    model.add(Dense(8, activation='softmax'))

    return model

def lstm():

    my_inpuy = Input(shape=(6000,), dtype='float64')

    # 在网络第一层，起降维的作用

    emb = Embedding(len(tokenizer.word_index) + 1, 5, input_length=6000)

    emb = emb(my_inpuy)

    net = Conv1D(16, 3, padding='same', kernel_initializer='glorot_uniform')(emb)

    net = BatchNormalization()(net)

    net = Activation('relu')(net)

    net = Conv1D(32, 3, padding='same', kernel_initializer='glorot_uniform')(net)

    net = BatchNormalization()(net)

    net = Activation('relu')(net)

    net = MaxPool1D(pool_size=4)(net)

    net1 = Conv1D(16, 4, padding='same', kernel_initializer='glorot_uniform')(emb)

    net1 = BatchNormalization()(net1)

    net1 = Activation('relu')(net1)

    net1 = Conv1D(32, 4, padding='same', kernel_initializer='glorot_uniform')(net1)

    net1 = BatchNormalization()(net1)

    net1 = Activation('relu')(net1)

    net1 = MaxPool1D(pool_size=4)(net1)

    net2 = Conv1D(16, 5, padding='same', kernel_initializer='glorot_uniform')(emb)

    net2 = BatchNormalization()(net2)

    net2 = Activation('relu')(net2)

    net2 = Conv1D(32, 5, padding='same', kernel_initializer='glorot_uniform')(net2)

    net2 = BatchNormalization()(net2)

    net2 = Activation('relu')(net2)

    net2 = MaxPool1D(pool_size=4)(net2)

    net = concatenate([net, net1, net2], axis=-1)

    net = CuDNNLSTM(256)(net)

    net = Dense(8, activation='softmax')(net)

    model = Model(inputs=my_inpuy, outputs=net)

    return model

def textcnn():

    kernel_size = [1, 3, 3, 5, 5]

    acti = 'relu'

    # 可看做一个文件的api集为一句话，然后话中的词总量是6000

    my_input = Input(shape=(inputLen,), dtype='int32')

    emb = Embedding(len(tokenizer.word_index) + 1, 20, input_length=inputLen)(my_input)

    emb = SpatialDropout1D(0.2)(emb)

    net = []

    for kernel in kernel_size:

        # 32个卷积核

        con = Conv1D(32, kernel, activation=acti, padding="same")(emb)

        # 滑动窗口大小是2,默认输出最后一维是通道数

        con = MaxPool1D(2)(con)

        net.append(con)

    # print(net)

    # input()

    net = concatenate(net, axis=-1)

    # net = concatenate(net)

    # print(net)

    # input()

    net = Flatten()(net)

    net = Dropout(0.5)(net)

    net = Dense(256, activation='relu')(net)

    net = Dropout(0.5)(net)

    net = Dense(8, activation='softmax')(net)

    model = Model(inputs=my_input, outputs=net)

    return model

test_result = np.zeros(shape=(len(test_apis),8))

# print(train_apis.shape)

# print(train_labels.shape)

# 5折交叉验证，将训练集切分成训练和验证集

skf = StratifiedKFold(n_splits=5)

for i, (train_index, valid_index) in enumerate(skf.split(train_apis, train_labels)):

    # print(i)

    # model = SequenceModel()

    model = textcnn()

    # metrics默认只有loss，加accuracy后在model.evaluate(...)的返回值即有accuracy结果

    model.compile(optimizer='adam', loss='sparse_categorical_crossentropy', metrics=['accuracy'])

    #模型保存规则

    model_save_path = './my_model/my_model_{}.h5'.format(str(i))

    checkpoint = ModelCheckpoint(model_save_path, save_best_only=True, save_weights_only=True)

    #早停规则

    earlystop = EarlyStopping(monitor='val_loss', min_delta=0, patience=5, verbose=0, mode='min', baseline=None,

                            restore_best_weights=True)

    #训练的过程会保存模型并早停

    model.fit(train_apis[train_index], train_labels[train_index], epochs=100, batch_size=1000,

              validation_data=(train_apis[valid_index], train_labels[valid_index]), callbacks=[checkpoint, earlystop])

    model.load_weights(model_save_path)

    # print(train_index, valid_index)

    test_tmpapis = model.predict(test_apis)

    test_result = test_result + test_tmpapis

# loss, acc = model.evaluate(train_apis, train_labels)

# print(loss)

# print(acc)

# print(model.predict(train_apis))

# print(test_files)

# print(test_apis)

test_result = test_result/5.0

with open(my_result, 'wb') as f:

    pickle.dump(test_files, f)

    pickle.dump(test_result, f)

# print(len(test_files))

# print(len(test_apis))

result = []

for i in range(len(test_files)):

    # #     print(test_files[i])

    #     #之前test_apis不带逗号的格式是矩阵格式，现在tolist转为带逗号的列表格式

    #     print(test_apis[i])

    #     print(test_apis[i].tolist())

    #     result.append(test_files[i])

    #     result.append(test_apis[i])

    tmp = []

    a = test_result[i].tolist()

    tmp.append(test_files[i])

    # extend相比于append可以添加多个值

    tmp.extend(a)

    #     print(tmp)

    result.append(tmp)

# print(1)

# print(result)

with open(my_result_csv, 'w') as f:

    #     f.write([1,2,3])

    result_csv = csv.writer(f)

    result_csv.writerow(["file_id", "prob0", "prob1", "prob2", "prob3", "prob4", "prob5", "prob6", "prob7"])

    result_csv.writerows(result)

可知，增加了早停机制后，约20代程序就被截止，valid不饱和。改进方案呢？

（1）看网络架构

卷积核变小、
层数变多、
去掉spatialdropout、
maxpool变meanpool，掺杂使用
注意初始化使用xavier
（2）数据量增大（小华说的分批次读取数据）
（3）放松早停机制限制条件，使训练较饱和
（4）单模型融合（不同学习率到达的不同局部最优结果汇总做融合）
（5）

尝试参考网上的，前向填充，这个影响大吗？

（6）学习搞xgboost分类器
 
 

阿里云恶意软件检测比赛-第三周-TextCNN的更多相关文章

  1. 确保数据零丢失!阿里云数据库RDS for MySQL 三节点企业版正式商用

    2019年10月23号,阿里云数据库RDS for MySQL 三节点企业版正式商用,RDS for MySQL三节点企业版基于Paxos协议实现数据库复制,每个事务日志确保至少同步两个节点,实现任意 ...

  2. 阿里云 Aliplayer高级功能介绍(三)：多字幕

    基本介绍 国际化场景下面,播放器支持多字幕,可以有效解决视频的传播障碍难题,该功能适用于视频内容在全球范围内推广,阿里云的媒体处理服务提供接口可以生成多字幕,现在先看一下具体的效果: WebVTT格式 ...

  3. 记一次阿里云服务器被用作DDOS攻击肉鸡

    事件描述:阿里云报警 ——检测该异常事件意味着您服务器上开启了"Chargen/DNS/NTP/SNMP/SSDP"这些UDP端口服务,黑客通过向该ECS发送伪造源IP和源端口的恶 ...

  4. 阿里云 ecs win2016 FileZilla Server

     Windows Server 2016 下使用 FileZilla Server 安装搭建 FTP 服务 一.安装 Filezilla Server 下载最新版本的 Filezilla Server ...

  5. 阿里云配置通用服务的坑 ssh: connect to host 47.103.101.102 port 22: Connection refused

    1.~ wjw$ ssh root@47.103.101.102 ssh: connect to host 47.103.101.102 port 22: Connection refused ssh ...

  6. 阿里云入选Gartner 2019 WAF魔力象限，唯一亚太厂商！

    近期,在全球权威咨询机构Gartner发布的2019 Web应用防火墙魔力象限中,阿里云Web应用防火墙成功入围,是亚太地区唯一一家进入该魔力象限的厂商! Web应用防火墙,简称WAF.在保护Web应 ...

  7. SaaS加速器，到底加速了谁？ 剖析阿里云的SaaS战略：企业和ISV不可错过的好文

    过去二十年,中国诞生了大批To C的高市值互联网巨头,2C的领域高速发展,而2B领域一直不温不火.近两年来,在C端流量饱和,B端数字化转型来临的背景下,中国越来越多的科技公司已经慢慢将触角延伸到了B端 ...

  8. 有关阿里云对SaaS行业的思考，看这一篇就够了

    过去二十年,随着改革开放的深化,以及中国的人口红利等因素,中国诞生了大批To C的高市值互联网巨头,2C的领域高速发展,而2B领域一直不温不火.近两年来,在C端流量饱和,B端数字化转型来临的背景下,中 ...

  9. 专访阿里云资深技术专家黄省江：中国SaaS公司的成功之路

    笔者采访中国SaaS厂商10多年,深感面对获客成本巨大.产品技术与功能成熟度不足.项目经营模式难以大规模复制.客户观念有待转变等诸多挑战,很多中国SaaS公司的经营状况都不容乐观. 7月26日,阿里云 ...

随机推荐

  1. Go：排序算法

    一.冒泡排序 package main import "fmt" func BubbleSort(arr []int) { /* 思路:将大的元素一步一步"冒泡" ...

  2. 独立集(bubble) 题解

    问题描述 有一天,一个名叫顺旺基的程序员从石头里诞生了.又有一天,他学会了冒泡排序和独立集.在一个图里,独立集就是一个点集,满足任意两个点之间没有边.于是他就想把这两个东西结合在一起.众所周知,独立集 ...

  3. Java 8新特性（二）：Stream API

    本篇文章继续介绍Java 8的另一个新特性--Stream API.新增的Stream API与InputStream和OutputStream是完全不同的概念,Stream API是对Java中集合 ...

  4. springboot整合druid监控配置

    方式一:直接引入druid 1.maven坐标 <dependency> <groupId>com.alibaba</groupId> <artifactId ...

  5. python元组的概念与基本操作

    元组与列表类似,关于元组同样需要做如下三点: A.概念 1.元组通过英文状态下的圆括号构成“()”.其存放元素与列表一样,可以是不通的数值类型,也可以是不通的数据结构. 2.元组仍然是一种序列,所以几 ...

  6. 20分钟理清Maven构建中的测试相关工具的关系

    如果你用Maven进行系统构建,同时还要同步编写测试用例,获取用例成功与否以及用例覆盖率的相关报告,那么这些工具你肯定接触过不少: JUnit TestNG maven-surefire-plugin ...

  7. muduo源码解析3-currentthread命名空间

    CurrentThread 作用: CurrentThread并不是一个类,而是一个命名空间,在mymuduo内部,目的是提供对于当前线程的管理操作. 内部变量: __thread int t_cac ...

  8. 分布式一致性算法 Paxos、Raft、Zab的区别与联系

    什么是分布式系统? 拿一个最简单的例子,就比如说我们的图书管理系统.之前的系统包含了所有的功能,比如用户注册登录.管理员功能.图书借阅管理等.这叫做集中式系统.也就是一个人干了好几件事. 后来随着功能 ...

  9. js对象数组新增、修改时的验证是否重复的逻辑

    JS代码: // 定义数据集合 const persons = [ { id: 1, name: '张三' }, { id: 2, name: '李四' } ] console.log('') con ...

  10. .net core国际化

    1.背景 公司业务遍及全球各地,对应业务系统国际化就是顺理成章的事情.最近就接手了一批新老系统的国际化任务,这里把一些探索经验.案例记录下来.本身改造和探索过程包括.NET MVC的,以及.NET C ...