一种养殖环境参数-鱼类生长预测模型建立方法及系统

摘要

本发明涉及一种养殖环境参数‑鱼类生长预测模型建立方法及系统，属于工厂化水产养殖技术领域，首先采集鱼类养殖数据；对鱼类养殖数据进行分类；对分类后的鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数之间的函数关系式；根据函数关系式，建立参数预测模型；对参数预测模型进行参数求解，得到可决系数；根据可决系数对参数预测模型进行修正，得到优化后的参数预测模型；对优化后的参数预测模型进行验证，并将养殖环境参数数据输入至优化后的参数预测模型中，输出鱼类生长效益指标数据，可对鱼类的生长进行精准预测，有利于鱼类养殖企业寻求最佳生产负荷和养殖规模，提高养殖企业的经济效益。

权利要求

1.一种养殖环境参数‑鱼类生长预测模型建立方法，其特征在于，包括：采集鱼类养殖数据；所述鱼类养殖数据包括养殖环境参数数据和鱼类生长指标数据；

具体包括：

采集所述养殖环境参数数据和所述鱼类生长指标数据，并建立数据库；

所述养殖环境参数数据包括水化因子指标数据和生物因子指标数据；所述水化因子指标数据包括养殖温度、氨氮和亚硝态氮；所述生物因子指标数据包括养殖种类、养殖模式、投喂方式、养殖密度和养殖鱼体质量；所述鱼类生长指标数据包括特定生长率；

对所述鱼类养殖数据进行分类；

对分类后的所述鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数之间的函数关系式，具体包括：选择任意种类下的所有所述鱼类养殖数据，将该种类下的所述鱼类生长指标作为因变量，各个所述养殖环境参数作为自变量，分别绘制散点图，并进行数据初步拟合，得到所述鱼类生长指标与各养殖环境参数之间的函数关系式；

根据所述函数关系式，建立参数预测模型，具体包括：

对所述函数关系式添加系数和扰动项，得到所述鱼类生长指标与所有所述养殖环境参数的总函数关系式，初步建立含有多个未知参数的所述参数预测模型；

得到特定生长率和各个养殖环境参数之间的数学函数关系式后，基于互作效应原理对特定生长率与各个养殖之间的函数关系式进行组合，将函数关系式中的参数以不同字母的形式表示，并添加系数和扰动项初步构成包含特定生长率以及所有养殖环境参数的总函数关系式，从而得到参数预测模型，总函数关系式为：其中，SGR表示特定生长率，T表示平均养殖温度，M表示初始鱼体质量，ρ表示初始养殖密度，e为自然常数，a、b、c、d、f为参数；

对所述参数预测模型进行参数求解，得到可决系数；

根据所述可决系数对所述参数预测模型进行修正，得到优化后的参数预测模型；

对所述优化后的参数预测模型进行验证，并将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出鱼类生长效益指标数据。

2.根据权利要求1所述的养殖环境参数‑鱼类生长预测模型建立方法，其特征在于，所述对所述鱼类养殖数据进行分类，具体包括：对所述养殖种类、所述养殖模式和所述投喂方式进行统一，基于此对其它的所述养殖环境参数数据和所述鱼类生长指标数据进行分类。

3.根据权利要求1所述的养殖环境参数‑鱼类生长预测模型建立方法，其特征在于，所述对所述参数预测模型进行参数求解，得到可决系数，具体包括：将分类后的所述鱼类养殖数据分为拟合集和验证集，所述拟合集用于所述参数预测模型的参数求解和模型修正，所述验证集用于模型验证；

采用1stOpt软件使用拟合集数据对所述参数预测模型的参数进行求解，优化算法采用LM算法以及全局优化算法，迭代求解后得到所述可决系数。

4.根据权利要求3所述的养殖环境参数‑鱼类生长预测模型建立方法，其特征在于，所述拟合集和所述验证集的样本比例为3：1。

5.根据权利要求3所述的养殖环境参数‑鱼类生长预测模型建立方法，其特征在于，所述根据所述可决系数对所述参数预测模型进行修正，得到优化后的参数预测模型，具体包括：根据所述可决系数的大小，对所述参数预测模型进行反馈修正，添加相应自变量的一次效应、二次效应，或者所述自变量之间的交互效应；

采用1stOpt软件使用所述拟合集对修正后的所述参数预测模型的参数进行求解，优化算法采用LM算法以及全局优化算法，迭代求解后得到最佳可决系数；

将所述最佳可决系数下的所述参数预测模型作为所述优化后的参数预测模型。

6.根据权利要求3所述的养殖环境参数‑鱼类生长预测模型建立方法，其特征在于，所述对所述优化后的参数预测模型进行验证，并将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出鱼类生长效益指标数据，具体包括：通过Matlab软件将所述验证集中相应的所述养殖环境参数数据的数值代入到所述优化后的参数预测模型中，得到相应的特定生长率；

将所述特定生长率与实际特定生长率进行对比，验证所述优化后的参数预测模型的适用性；

将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出相应的所述鱼类生长效益指标数据，根据所述鱼类生长效益指标数据调节所述养殖环境参数数据，对养殖系统进行合理优化。

7.一种养殖环境参数‑鱼类生长预测模型建立系统，其特征在于，包括：鱼类养殖数据采集模块，用于采集鱼类养殖数据；所述鱼类养殖数据包括养殖环境参数数据和鱼类生长指标数据；具体包括：采集所述养殖环境参数数据和所述鱼类生长指标数据，并建立数据库；

所述养殖环境参数数据包括水化因子指标数据和生物因子指标数据；所述水化因子指标数据包括养殖温度、氨氮和亚硝态氮；所述生物因子指标数据包括养殖种类、养殖模式、投喂方式、养殖密度和养殖鱼体质量；所述鱼类生长指标数据包括特定生长率；

鱼类养殖数据分类模块，用于对所述鱼类养殖数据进行分类；

函数关系式获取模块，用于对分类后的所述鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数之间的函数关系式，具体包括：选择任意种类下的所有所述鱼类养殖数据，将该种类下的所述鱼类生长指标作为因变量，各个所述养殖环境参数作为自变量，分别绘制散点图，并进行数据初步拟合，得到所述鱼类生长指标与各养殖环境参数之间的函数关系式；

参数预测模型构建模块，用于根据所述函数关系式，建立参数预测模型，具体包括：对所述函数关系式添加系数和扰动项，得到所述鱼类生长指标与所有所述养殖环境参数的总函数关系式，初步建立含有多个未知参数的所述参数预测模型；

得到特定生长率和各个养殖环境参数之间的数学函数关系式后，基于互作效应原理对特定生长率与各个养殖之间的函数关系式进行组合，将函数关系式中的参数以不同字母的形式表示，并添加系数和扰动项初步构成包含特定生长率以及所有养殖环境参数的总函数关系式，从而得到参数预测模型，总函数关系式为：其中，SGR表示特定生长率，T表示平均养殖温度，M表示初始鱼体质量，ρ表示初始养殖密度，e为自然常数，a、b、c、d、f为参数；

参数预测模型求解模块，用于对所述参数预测模型进行参数求解，得到可决系数；

参数预测模型修正模块，用于根据所述可决系数对所述参数预测模型进行修正，得到优化后的参数预测模型；

参数预测模型验证和输出模块，用于对所述优化后的参数预测模型进行验证，并将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出鱼类生长效益指标数据。

说明书

一种养殖环境参数‑鱼类生长预测模型建立方法及系统

技术领域

[0001] 本发明涉及工厂化水产养殖技术领域，特别是涉及一种基于数据分析的养殖环境参数‑鱼类生长预测模型建立方法及系统。

背景技术

[0002] 近年来，工厂化循环水鱼类养殖是一种现代工业化生产方式，主要依赖机械设备对水体进行调控，其核心组成是高效生物过滤装置和养殖用水自动监测系统，可以实现养殖环境高度可控。目前生物过滤器主要指标如氨氮负荷率通常设置在不合理的较高水平，这很大程度上增加了能源消耗，而能源消耗问题是制约工厂化循环水养殖盈利的主要原因之一。工厂化循环水养殖的发展趋势之一是进行工程经济性分析，运用工程学、经济学、统计学等方法，寻找最佳的生产负荷和养殖规模，系统地建立养殖环境参数‑鱼类生长指标的响应预测数学模型，探求鱼类如何从生长策略上对养殖环境参数的变化做出响应，模型的建立对于促进养殖系统的优化和养殖效益的提高具有重要意义。

[0003] 目前对鱼类生长的预测大多使用统计学方法进行研究，并集中在回归分析法和数据拟合法两种技术手段上。其中，回归分析法偏重于相关性的探究，它可以依托于现有的数据建立回归模型，从而进行相关性分析，如分析鱼生长效率和不同水环境因子之间的相关性，但其预测能力有很大的局限性，预测精度低。而数据拟合法虽然具有较强的预测能力，但数据拟合法中常用的特殊VonBertalanffy生长方程、一般VonBertalanffy生长方程、Logistic生长方程、Gompertz生长方程、Brody生长方程、指数生长方程等6个动物生长模型并不适用于鱼类生长效益与养殖环境参数关系的研究，无法通过养殖环境参数的变化对鱼类生长进行预测。因此，如何通过养殖环境参数对鱼类的生长进行精准预测是本领域亟需解决的问题。

发明内容

[0004] 本发明的目的是提供一种养殖环境参数‑鱼类生长预测模型建立方法及系统，将鱼类生长对养殖环境参数变化的响应以函数关系式的形式表现出来，通过构建参数预测模型对鱼类生长进行预测和验证，提高了预测精度，为工厂化循环水养殖系统的优化提供理论数据，可用于工厂化循环水养殖系统寻求最佳生产负荷和养殖规模，提高养殖企业的经济效益。

[0005] 为实现上述目的，本发明提供了如下方案：

[0006] 一种养殖环境参数‑鱼类生长预测模型建立方法，包括：

[0007] 采集鱼类养殖数据；所述鱼类养殖数据包括养殖环境参数数据和鱼类生长指标数据；

[0008] 对所述鱼类养殖数据进行分类；

[0009] 对分类后的所述鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数之间的函数关系式；

[0010] 根据所述函数关系式，建立参数预测模型；

[0011] 对所述参数预测模型进行参数求解，得到可决系数；

[0012] 根据所述可决系数对所述参数预测模型进行修正，得到优化后的参数预测模型；

[0013] 对所述优化后的参数预测模型进行验证，并将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出鱼类生长效益指标数据。

[0014] 可选的，所述采集鱼类养殖数据，具体包括：

[0015] 采集所述养殖环境参数数据和所述鱼类生长指标数据，并建立数据库；

[0016] 所述养殖环境参数数据包括水化因子指标数据和生物因子指标数据；所述水化因子指标数据包括养殖温度、氨氮和亚硝态氮；所述生物因子指标数据包括养殖种类、养殖模式、投喂方式、养殖密度和养殖鱼体质量；所述鱼类生长指标数据包括特定生长率。

[0017] 可选的，所述对所述鱼类养殖数据进行分类，具体包括：

[0018] 对所述养殖种类、所述养殖模式和所述投喂方式进行统一，基于此对其它的所述养殖环境参数数据和所述鱼类生长指标数据进行分类。

[0019] 可选的，所述对分类后的所述鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数之间的函数关系式，具体包括：

[0020] 选择任意种类下的所有所述鱼类养殖数据，将该种类下的所述鱼类生长指标作为因变量，各个所述养殖环境参数作为自变量，分别绘制散点图，并进行数据初步拟合，得到所述鱼类生长指标与各养殖环境参数之间的函数关系式。

[0021] 可选的，所述根据所述函数关系式，建立参数预测模型，具体包括：

[0022] 对所述函数关系式添加系数和扰动项，得到所述鱼类生长指标与所有所述养殖环境参数的总函数关系式，初步建立含有多个未知参数的所述参数预测模型。

[0023] 可选的，所述对所述参数预测模型进行参数求解，得到可决系数，具体包括：

[0024] 将分类后的所述鱼类养殖数据分为拟合集和验证集，所述拟合集用于所述参数预测模型的参数求解和模型修正，所述验证集用于模型验证；

[0025] 采用1stOpt软件使用拟合集数据对所述参数预测模型的参数进行求解，优化算法采用LM算法以及全局优化算法，迭代求解后得到所述可决系数。

[0026] 可选的，所述拟合集和所述验证集的样本比例为3：1。

[0027] 可选的，所述根据所述可决系数对所述参数预测模型进行修正，得到优化后的参数预测模型，具体包括：

[0028] 根据所述可决系数的大小，对所述参数预测模型进行反馈修正，添加相应自变量的一次效应、二次效应，或者所述自变量之间的交互效应；

[0029] 采用1stOpt软件使用所述拟合集对修正后的所述参数预测模型的参数进行求解，优化算法采用LM算法以及全局优化算法，迭代求解后得到最佳可决系数；

[0030] 将所述最佳可决系数下的所述参数预测模型作为所述优化后的参数预测模型。

[0031] 可选的，所述对所述优化后的参数预测模型进行验证，并将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出鱼类生长效益指标数据，具体包括：

[0032] 通过Matlab软件将所述验证集中相应的所述养殖环境参数数据的数值代入到所述优化后的参数预测模型中，得到相应的特定生长率；

[0033] 将所述特定生长率与实际特定生长率进行对比，验证所述优化后的参数预测模型的适用性；

[0034] 将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出相应的所述鱼类生长效益指标数据，根据所述鱼类生长效益指标数据调节所述养殖环境参数数据，对养殖系统进行合理优化。

[0035] 一种养殖环境参数‑鱼类生长预测模型建立系统，包括：

[0036] 鱼类养殖数据采集模块，用于采集鱼类养殖数据；所述鱼类养殖数据包括养殖环境参数数据和鱼类生长指标数据；

[0037] 鱼类养殖数据分类模块，用于对所述鱼类养殖数据进行分类；

[0038] 函数关系式获取模块，用于对分类后的所述鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数之间的函数关系式；

[0039] 参数预测模型构建模块，用于根据所述函数关系式，建立参数预测模型；

[0040] 参数预测模型求解模块，用于对所述参数预测模型进行参数求解，得到可决系数；

[0041] 参数预测模型修正模块，用于根据所述可决系数对所述参数预测模型进行修正，得到优化后的参数预测模型；

[0042] 参数预测模型验证和输出模块，用于对所述优化后的参数预测模型进行验证，并将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出鱼类生长效益指标数据。

[0043] 根据本发明提供的具体实施例，本发明公开了以下技术效果：

[0044] 本发明提出了一种养殖环境参数‑鱼类生长预测模型建立方法及系统，采集鱼类养殖数据；鱼类养殖数据包括养殖环境参数指标数据和鱼类生长指标数据；对鱼类养殖数据进行分类；对分类后的鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数之间的函数关系式；根据函数关系式，建立参数预测模型；对参数预测模型进行参数求解，得到可决系数；根据可决系数对参数预测模型进行修正，得到优化后的参数预测模型；利用该优化后的参数预测模型，可对鱼类的生长进行精准预测，解决现有鱼类生长预测方法存在的预测精度低的问题。最终对优化后的参数预测模型进行验证，并将养殖环境参数数据输入至优化后的参数预测模型中，输出鱼类生长效益指标数据。鱼类效益指标数据就是预测出来的鱼类生长性能相关数据，根据鱼类效益指标数据调节养殖环境参数数据，有利于鱼类养殖企业寻求最佳生产负荷和养殖规模，为指导养殖企业对养殖系统现状进行合理优化和改善提供了理论数据，可有效提高养殖企业的经济效益。

附图说明

[0045] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

[0046] 图1为本发明实施例1提供的一种养殖环境参数‑鱼类生长预测模型建立方法的流程图；

[0047] 图2为本发明实施例1提供的草鱼特定生长率对初始养殖密度的响应曲线图；

[0048] 图3为本发明实施例1提供的草鱼特定生长率对初始鱼体质量的响应曲线图；

[0049] 图4为本发明实施例1提供的草鱼特定生长率对养殖温度的响应曲线图；

[0050] 图5为本发明实施例1提供的参数预测模型的适用性验证的结果示意图；

[0051] 图6为本发明实施例2提供的一种养殖环境参数‑鱼类生长预测模型建立系统的结构框图。

具体实施方式

[0052] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0053] 本发明的目的是提供一种养殖环境参数‑鱼类生长预测模型的建立方法及系统，在寻找数据之间相关性的同时，通过严谨的数学计算得到一种鱼类生长对养殖环境参数响应的预测模型，将鱼类生长对养殖环境参数变化的响应以数学关系式的形式表现出来以实现相关性可视化。根据养殖环境参数数据和鱼类生长指标数据构建参数预测模型，用于鱼类养殖过程中的生长预测，对养殖系统进行合理优化，有助于鱼类的生长发育，为养殖企业带来效益。

[0054] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

[0055] 实施例1

[0056] 如图1所示，本实施例示出了一种养殖环境参数‑鱼类生长预测模型建立方法，该方法具体包括以下步骤：

[0057] 步骤S1、采集鱼类养殖数据；所述鱼类养殖数据包括养殖环境参数数据和鱼类生长指标数据；所述养殖环境参数数据包括水化因子指标数据和生物因子指标数据，具体包括：

[0058] 采集所述养殖环境参数数据和所述鱼类生长指标数据，并建立数据库；在建立数据库时，利用各种可以收集鱼类生长以及相应养殖环境参数数据的方法来对数据库进行补充和丰富，数据来源包括且不限于目前已有文献、实地调查和科学实验获取等等。

[0059] 所述水化因子指标数据包括但不限于养殖温度、氨氮和亚硝态氮；所述生物因子指标数据包括但不限于养殖种类、养殖模式、投喂方式、养殖密度和养殖鱼体质量；所述鱼类生长指标数据包括但不限于特定生长率，还可以为增重率、饵料系数、食物转化效率和肥满度等。

[0060] 需要说明的是，本实施例中选择的养殖温度、初始养殖密度和初始鱼体质量这三个养殖环境参数指标数据，以及特定生长率这一个鱼类生长指标数据，是优选出来的对鱼生长影响较大的几种养殖环境参数因子。实际上，养殖环境参数指标还可以包括溶解氧、总氮、氨氮、硝态氮、亚硝酸盐氮、总磷、可溶性活性磷、叶绿素a和化学需氧量等指标。鱼类生长指标包括很多种，如特定生长率、增重率、饵料系数、食物转化效率和肥满度等。本申请选择最具有代表性的特定生长率，特定生长率可直接反映出鱼类生长速度的快慢。因此，本实施例设置的养殖环境参数和鱼类生长指标仅仅是优选的其中几种，并不是固定的、唯一的，不应作为对本发明保护范围的限制，还可以包括其他参数因子，可根据实际情况自行设定，任何通过本发明的方法对养殖环境参数和鱼类生长指标的设定都应在本发明的保护范围之内。

[0061] 下面以草鱼为例进行详细举例说明：

[0062] 本实施例中，通过查阅文献采集草鱼养殖数据共142组并建立数据库，根据养殖数据得到草鱼的特定生长率(SGR)，另外分析养殖环境参数数据，发现各组数据间差异较大的为养殖温度、初始养殖密度、初始鱼体质量三项指标，其他指标如氨氮、亚硝态氮等水平处于安全浓度范围内，未对草鱼生长造成明显影响，故本实施例不加考虑。则有：

[0063] SGR＝(lnW末‑lnW初)/t×100％

[0064] 其中，SGR表示特定生长率，W末表示测量周期末的鱼质量，W初表示测量周期初的鱼质量；t表示一个测量周期的天数。

[0065] 步骤S2、对所述鱼类养殖数据进行分类；具体包括：

[0066] 对所述养殖种类、所述养殖模式和所述投喂方式进行统一，基于此对其它的所述养殖环境参数数据和所述鱼类生长指标数据进行分类。

[0067] 应说明的是，影响鱼类生长效益的因素除了水化因子之外，还有养殖模式、投喂方式以及养殖种类等生物因子因素，为提高建模的成功率，本实施例对数据库中的草鱼养殖数据进行分类，对养殖种类、养殖模式以及投喂方式等进行统一，即控制养殖种类、养殖模式以及投喂方式这三种变量因素保持不变，对包括养殖密度、养殖鱼体质量、养殖温度、氨氮和亚硝态氮等养殖环境参数数据以及包括特定生长率的鱼类生长指标数据进行分类，在此前提下再进行数据建模。例如，按照投喂方式可分有饱食投喂、定量投喂等，进行这种分类的目的在于，因为鱼的不同投喂方式会对鱼类生长产生较大的影响，为消除这种影响而进行分类，比如类别之一为饱食投喂下的草鱼，此类中的具体数据有草鱼生长指标以及对应的养殖环境参数，以此数据分析得到的模型只适用于对饱食投喂下的草鱼的生长情况进行预测，而不适用于其他种类和投喂方式的鱼的生长情况。本实施例中主要考虑养殖密度、养殖鱼体质量、养殖温度这几种对鱼生长影响较大的指标。本发明在对影响鱼类生长情况预测的部分变量因素进行统一后，对水化因子指标数据、生物因子指标数据和鱼类生长指标数据进行分类和分析，而有效地消除了变量因素对预测结果的影响，保证最终预测结果的准确性，提升了预测精度。

[0068] 步骤S3、对分类后的所述鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数指标之间的函数关系式；具体包括：

[0069] 选择任意种类下的所有所述鱼类养殖数据，将该种类下的所述鱼类生长指标作为因变量，各个所述养殖环境参数指标作为自变量，分别绘制散点图，并进行数据初步拟合，得到鱼类生长指标与各养殖环境参数指标之间的函数关系式。

[0070] 本实施例中，对数据库中的数据分类后，以草鱼生长指标如特定生长率为因变量，养殖温度、初始养殖密度等养殖环境参数指标作为自变量，分别绘制出散点图，散点图就是如图2‑图4所示的响应曲线图，曲线两侧的散点分布，并进行数据初步拟合，得到特定生长率与各个养殖环境参数之间的函数关系式及拟合优度。步骤S3的关键点在于进行单养殖环境参数因子分析前，应对其它养殖环境参数进行变量范围控制以获得较高的拟合优度，从而提高模型预测的准确性。其中，变量范围控制即降低其它养殖环境参数对图形的干扰，比如在获得特定生长率和养殖温度之间的散点图时，应保证所选数据中其它因素如初始养殖密度、氨氮浓度等在较小的波动范围区间。

[0071] 如图2所示，对特定生长率和初始养殖密度进行分析，养殖温度控制在25～27℃，草鱼的初始鱼体质量为0～10g，得到草鱼的特定生长率对初始养殖密度的响应：

[0072] SGR＝1.3852ρ‑0.308，R2＝0.5824

[0073] 其中，SGR表示特定生长率，ρ表示草鱼的初始养殖密度，R2表示可决系数，可决系数是度量拟合优度的统计量。

[0074] 如图3所示，对特定生长率和草鱼的初始鱼体质量进行分析，养殖温度控制在25～27℃，初始养殖密度为0～4.2g/L，得到草鱼的特定生长率对初始鱼体质量的响应：

[0075] SGR＝2.9515M‑0.192，R2＝0.5149

[0076] 其中，SGR表示特定生长率，M表示草鱼的初始鱼体质量，R2表示可决系数。

[0077] 如图4所示，对特定生长率和养殖期间内平均养殖温度进行分析，草鱼的初始鱼体质量为0～10g，初始养殖密度为0～4.2g/L，得到草鱼的特定生长率对养殖温度的响应：

[0078] SGR＝0.0017e0.2552T，R2＝0.6122

[0079] 其中，SGR表示特定生长率，T表示草鱼的平均养殖温度，e为自然常数，R2表示可决系数。

[0080] 通过上述步骤得到草鱼的特定生长率与初始养殖密度、初始鱼体质量、养殖温度三项指标之间的相关性以及函数关系式。

[0081] 步骤S4、根据所述函数关系式，建立参数预测模型；具体包括：

[0082] 对所述函数关系式添加系数和扰动项，得到所述鱼类生长指标与所有所述养殖环境参数指标的总函数关系式，初步建立含有多个未知参数的所述参数预测模型。

[0083] 通过上述步骤可以得到特定生长率和各个养殖环境参数之间的数学函数关系。养殖环境参数中水化因子和生物因子通过综合效应影响鱼类的生长，即互作效应。基于互作效应原理对特定生长率与各个养殖之间的函数关系式进行组合，此过程中应将函数关系式中的参数以不同字母的形式表示，并添加系数和扰动项等要素初步构成包含特定生长率以及所有养殖环境参数的函数关系式，从而得到参数预测模型，总函数关系式为：

[0084]

[0085] 其中，SGR表示特定生长率，T表示草鱼的平均养殖温度，M表示草鱼的初始鱼体质量，ρ表示草鱼的初始养殖密度，e为自然常数，a、b、c、d、f为参数。

[0086] 本实施例中，在步骤S3中，可以选取养殖环境参数中关键影响因子进行后续的步骤S4的分析，参照初步拟合结果，并结合相关性分析，选取相关性较高的关键养殖环境参数进行后续建模过程。并且，若养殖环境参数较多，还可以采用主成分分析法或因子分析法对其进行降维，从而更容易确定关键的养殖环境参数，简化计算过程，减少计算量。

[0087] 步骤S5、对所述参数预测模型进行参数求解，得到可决系数；具体包括：

[0088] 步骤S5.1、将分类后的所述鱼类养殖数据分为拟合集和验证集，所述拟合集用于所述参数预测模型的参数求解和模型修正，所述验证集用于模型验证。

[0089] 本实施例中，所述拟合集和所述验证集的样本比例为3：1，实际可根据具体情况自行设定拟合集和验证集的比例以及具体数量。

[0090] 本实施例中，将鱼类养殖数据分为拟合集和验证集，是指将鱼类养殖数据中对应的每一种指标数据都按照比例分为拟合集和验证集，鱼类养殖数据划分完成后，其中拟合集为107组数据，约占总数据量的75％，验证集为35组数据，约占总数据量的25％。同时，为方便对养殖数据进行处理，可对养殖数据进行标准化，以消除各个数据间的量纲的影响，进一步提升模型预测的准确性。

[0091] 步骤S5.2、采用1stOpt软件使用拟合集数据对所述参数预测模型的参数进行求解，优化算法采用LM算法以及全局优化算法，迭代求解后得到所述可决系数。

[0092] 本实施例中，迭代28次达到收敛判断标准后，得到可决系数R2的结果：

[0093]

[0094] 将上述参数代入前述涉及可决系数的公式中，经计算得到可决系数R2＝0.4。

[0095] 本实施例中，由于参数预测模型中的参数较多且为非线性函数，因此应将非线性最小二乘问题转换成线性最小二乘问题，并进行迭代求解。此类函数逼近法常用优化方法为莱文博格‑马夸特算法(Levenberg–Marquardt算法，简称LM算法)，通过1stOpt软件使用拟合集数据对模型参数进行求解并得到模型的可决系数，利用1stOpt软件独特的全局优化算法可以在使用迭代法的同时，使给定合适的初始值成为非必需，方便实用。

[0096] 步骤S6、根据所述可决系数对所述参数预测模型进行修正，得到优化后的参数预测模型；具体包括：

[0097] 步骤S6.1、根据所述可决系数的大小，对所述参数预测模型进行反馈修正，添加相应自变量的一次效应、二次效应，或者所述自变量之间的交互效应；

[0098] 本实施例中的模型拟合效果不佳，则需要对其进行修正，考虑变量可能会存在一次效应，故对其添加一次项：

[0099]

[0100] 其中：SGR表示特定生长率，T表示草鱼的平均养殖温度，M表示草鱼的初始鱼体质量，ρ表示草鱼的初始养殖密度，e为自然常数，a、b、c、d、f、g、h为参数。

[0101] 步骤S6.2、再次采用1stOpt软件使用所述拟合集对修正后的所述参数预测模型的参数进行求解，优化算法采用LM算法以及全局优化算法，迭代求解后得到最佳可决系数；

[0102] 本实施例中，在28次迭代的基础上再次进行迭代，在迭代29次达到收敛判断标准后，计算结束：

[0103]

[0104] 其中：SGR表示特定生长率，T表示草鱼的平均养殖温度，M表示草鱼的初始鱼体质量，ρ表示草鱼的初始养殖密度，e为自然常数。

[0105] 经过步骤S6修正过程后，得到最佳的可决系数R2＝0.79，可决系数R2达到0.79，表明参数预测模型的拟合效果较好。

[0106] 步骤S6.3、将所述最佳可决系数下的所述参数预测模型作为所述优化后的参数预测模型。应说明，拟合集中数据的准确性以及模型修正的合理性直接影响到步骤S6的成功率。

[0107] 步骤S7、对所述优化后的参数预测模型进行验证，并将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出鱼类生长效益指标数据。具体包括：

[0108] 步骤S7.1、通过Matlab软件将所述验证集中相应的所述养殖环境参数数据的数值代入到所述优化后的参数预测模型中，得到相应的特定生长率；

[0109] 步骤S7.2、将所述特定生长率与实际特定生长率进行对比，验证所述优化后的参数预测模型的适用性；

[0110] 步骤S7.3、将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出相应的所述鱼类生长效益指标数据，根据所述鱼类生长效益指标数据调节所述养殖环境参数数据，对养殖系统进行合理优化。

[0111] 需要说明的是，通过拟合集得到的优化后的参数预测模型虽具有较好的拟合优度，但还需要对其适用性进行验证。各个参数的数值在步骤S6中已经确定，只需将验证集中相应的养殖环境参数数值代入到优化后的参数预测模型的公式中，即可得到相应的特定生长率，将其与实际特定生长率进行对比，例如将计算得到的特定生长率和实际的特定生长率通过图像进行对比观察，以此来评价参数预测模型的适用性。

[0112] 如图5所示，x轴表示通过模型拟合得到的SGR值，y轴表示数据库中与之对应的SGR真实值，即实际特定生长率。图5中验证集35组数据均匀地分布在对角线y＝x周围，说明模型拟合能力较好，优化后的参数预测模型具有较好的适用性，且并没有出现过拟合现象。通过本发明的方法，最终得到了精确的养殖环境参数‑草鱼生长预测模型，利用该预测模型可为工厂化循环水养殖寻找最佳的生产负荷和养殖规模。

[0113] 还应该说明，本实施例将鱼类设置为草鱼，以及养殖数据的选取，养殖数据的分类以及拟合集和验证集的分组，初始养殖温度、养殖密度和初始鱼体质量的具体范围，以及其他具体参数的设置等等，仅仅是列举一个实例进行说明，并不是对本发明保护范围的限制，均可根据实际情况自行进行设定。例如，鱼类不仅仅局限于草鱼，还可以是其他任意一种养殖鱼类。同时，养殖环境参数的变化对鱼类生长非常重要，养殖过程中鱼类的生长效益不仅仅被养殖密度、养殖温度等影响，而是受到养殖环境中多种因子综合效应的影响。在收集更多的数据之后，可在本发明参数预测模型的基础上进行补充完善，若变量过多难以操作，还可以通过降维的方法减少自变量，具体可参考主成分分析法或因子分析法，将自变量降维到较少数量后，再进行建模分析，更加简单、方便，计算量更少，预测结果更加准确、可靠。

[0114] 实施例2

[0115] 如图6所示，本实施例提供了一种养殖环境参数‑鱼类生长预测模型建立系统，该系统具体包括：

[0116] 鱼类养殖数据采集模块M1，用于采集鱼类养殖数据；所述鱼类养殖数据包括养殖环境参数数据和鱼类生长指标数据；

[0117] 鱼类养殖数据分类模块M2，用于对所述鱼类养殖数据进行分类；

[0118] 函数关系式获取模块M3，用于对分类后的所述鱼类养殖数据进行相关性分析，得到鱼类生长指标与各养殖环境参数之间的函数关系式；

[0119] 参数预测模型构建模块M4，用于根据所述函数关系式，建立参数预测模型；

[0120] 参数预测模型求解模块M5，用于对所述参数预测模型进行参数求解，得到可决系数；

[0121] 参数预测模型修正模块M6，用于根据所述可决系数对所述参数预测模型进行修正，得到优化后的参数预测模型；

[0122] 参数预测模型验证和输出模块M7，用于对所述优化后的参数预测模型进行验证，并将待预测的所述养殖环境参数数据输入至所述优化后的参数预测模型中，输出鱼类生长效益指标数据。

[0123] 本发明依托于统计学方法，对数据进行收集及深度挖掘，将鱼类养殖过程中关键因素如养殖温度、初始养殖密度和初始鱼体质量等作为响应因子，以鱼类生长效益指标作为响应值进行数据分析，探究养殖环境因子(水质因子和生物因子)与鱼类生长效益之间的相关性，得到鱼类生长效益指标与养殖环境因子两者之间的函数关系式，并对模型进行准确性验证。基于此模型则可通过收集养殖环境因子来对相应鱼类生长效益进行预测。

[0124] 本发明将鱼类生长对养殖环境参数变化的响应以数学关系式的形式表现出来，实现了根据养殖环境参数对鱼类生长效益的预测功能，由于存在验证集，参数预测模型的适用性也得到肯定，其准确性高、应用范围广，可为智慧渔业和数字渔业提供数据支持。参数预测模型的建立，可以为工厂化循环水养殖寻找最佳的生产负荷和养殖规模，降低水处理成本并掌握养殖和苗种生产中的适当放养密度，满足了不同用户和不同养殖对象的生产需求。另外，参数预测模型仅通过采集养殖环境参数即可对鱼类生长效益进行预测，这将极大提高工厂化循环水养殖的自动化水平，可以完善目前养殖过程中缺乏的生产预警机制与控制策略并降低养殖精准生产的数据信息采集难度，有利于我国水产养殖行业的发展。

[0125] 本说明书中各个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。本说明书中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。