淡水养殖池塘水华发生及池水净化处理数学建模论文

淡水养殖池塘水华发生及池水净化处理,网赛论文

淡水养殖池塘水华发生及池水净化处理

摘要

随着社会的的发展，越来越多的人加入到水产养殖行业，而其中池塘养殖产量约占淡水养殖的70%。但是随着淡水生态系统水体污染和富营养化进程的加剧，经常导致有害蓝藻、轮虫等常见的浮游生物高密度发生，很容易诱发大面积水华。水华造成严重的环境污染及水体污染，对养殖业是一个严重的打击。

本文主要采用了MATLAB程序中的相关系数分析，模糊综合评价，单侧检验等方法对淡水养殖池塘水华发生及池水净化处理的相关问题进行了分析，建立相关模型。

针对问题一，首先将题目中要进行分析的量给找出来，同时将他们运用MATLAB进行相关系数分析，在此基础上分析水体、底泥与间隙水中常见主要理化因子之间的关系，并分析原因。分析水体、底泥与间隙水中常见主要理化因子之间的关系，并分析原因。

针对问题二，建立模糊综合评价。首先，对数据指标进行归一化处理，并利用层次分析法和因子分析法确定各指标因素的权重，最后利用确定的权重，建立池塘水体质量的综合评价模型，对池塘水体质量进行分级。

针对问题三，建立单侧检验相关性模型。首先，运用SPSS软件分析理化因子与水华发生的相关性；然后进一步分析，得出结论。

针对问题四，利用MATLAB建立鱼类生长周期体重模型，运用二次函数建立关于体重与生长周期的拟合方程。建立浮游植物密度与时间的关系模型并得到图像。

针对问题五，通过网上查阅资料结合附件资料分析，可以得到有利于池水养殖池塘水体的自净化的方法，并据此提供建立生态养殖模式的方案。

关键词：单侧检验相关系数分析回归分析综合评价

1

淡水养殖池塘水华发生及池水净化处理,网赛论文

一、问题重述

目前在我国水产养殖中，池塘养殖产量约占淡水养殖的70%。近年来，随着淡水生态系统水体污染和富营养化进程的加剧，经常导致有害蓝藻、轮虫等常见的浮游生物高密度发生，很容易诱发大面积水华。水华造成严重的环境污染及水体污染，对养殖业是一个严重的打击。

水华的发生不仅直接影响了养殖对象的正常生长发育，严重时大量排泄废水造成淡水资源污染，还会破坏养殖生态系统的平衡，导致养殖对象的不同程度死亡，造成巨大经济损失。为此我们通过研究淡水养殖池塘相关主要理化因子，主要浮游生物数据及鱼虾生成等数据分析水华发生的原因，控制并预测水华的发生，从而提高养殖产量，减小环境污染等。通过对水华发生的了解，加强大家环保意识。

根据附件1-8完成如下问题：

1）通过附件1中数据分析水体、底泥与间隙水中常见主要理化因子之间的关系，并分析原因。

2）通过附件2中数据对四个池塘水体质量进行评价及分类，分析虾池与鱼池对水体的影响。

3）建立主要理化因子和常见浮游生物致害密度发生关系的模型，给出水华发生时主要理化因子的范围，预测淡水养殖池塘水华发生 (1号池发生轻微水华)。

4）结合附件及以上分析，建立鱼类生长与体重相关模型。在养殖鲢鱼、鳙鱼等的生长过程中可以摄食浮游生物，净化某些藻类，构造一个与1号池相同大小的净化池，通过水循环，并放养鲢鱼或鳙鱼，放养多少才能净化1号池中的藻类，净化效果如何。

5）结合附件及通过查阅资料构建一种生态养殖模式，有利于池水养殖池塘水体的自净化。通过以上养殖从而使淡水养殖减少向江河湖海养殖废水排放。

二、问题分析

对于问题一，对附件1中的六种理化因子与池水、间隙水和底泥的进行分析和整理，并用MATLAB对数据进行相关系数处理，并分析原因。

对于问题二，根据水池取样点报告中，采用均方差法确定各水池的各项指标的权重，并用隶属度函数对水质进行检测，最后用模糊综合评判模型对水质进行评判得出水池的水质类型。

对于问题三，要求建立主要理化因子和常见浮游生物致害密度发生关系的模型，并给出水华发生时主要理化因子的范围，再根据理化因子的范围预测淡水养殖池塘水华发生。考虑对主要理化因子和常见浮游生物密度利用单侧检验建立模型对问题进行相关性和显著性分析。

对于问题四，由附件6里体重和体长的数据，利用MATLAB建立鱼类生长体重模型，运用二次函数建立关于体重体长的拟合方程。根据附件4中浮游植物密度的数据利用MATLAB建立浮游植物密度与时间的关系模型并得到图像。

对于问题五，通过网上查阅资料结合附件资料分析，可以得到有利于池水养殖池塘

2

淡水养殖池塘水华发生及池水净化处理,网赛论文

3

水体的自净化的方法，并据此提供建立生态养殖模式的方案。

为解决上述五个问题，我们将采用Excel 、SPSS 、MATLAB 等软件来协助分析。

三、模型假设

1．假设附件中给出的数据全部真实有效 2．假设不会发生任何天灾影响四个池塘 3. 假设池塘主人不会对池塘疏于管理

四、符号说明

()s k

归一化均方差 P

确定某个因子是否显著

()f DO

隶属度函数

浮游植物密度 s ie X 隶属度 N

浮游植物的个数

k x

每次检测值 m 浮游植物的质量 i x

评价指标 v

池塘1的体积 s i X

第s 个评价点的总评价系数

M

放养鱼的质量 s ie

r 第s 个站点第e 个等级的隶属度

()s B

第s 个站点的评价模型

()s R

第s 个站点的评价矩阵

五、模型的建立及求解

5.1问题一的模型建立及求解

此问题要求对附件1中数据分析水体、底泥与间隙水中常见主要理化因子之间的关系，并分析原因。运用MATLAB 程序对附件1中的数据做相关性分析。（程序代码见附录1-12）

表1：池塘1池水中的主要理化因子的相关系数

总磷 磷酸盐磷

总氮

硝态氮

亚硝态氮

铵态氮

总磷

1

0.95

0.71

-0.0

0.02

0.48

淡水养殖池塘水华发生及池水净化处理,网赛论文

9881 9285 4118 0593 0061

磷酸盐磷

0.95

9881

1 0.62

8607

-0.0

3857

0.05

1013

0.38

0394

总氮0.71

9285

0.62

8607

1 0.07

1894

0.02

2963

0.55

0647

硝态氮

-0.0

4118

-0.0

3857

0.07

1894

1 0.74

2675

0.33

1797

亚硝态氮

0.02

0593

0.05

1013

0.02

2963

0.74

2675

1 0.40

0354

铵态氮

0.48

0061

0.38

0394

0.55

0647

0.33

1797

0.40

0354

1

由以上相关系数矩阵发现所有元素之间有正相关的也有负相关的，以0.5作为衡量标准，可以看出超过0.5的相关系数的相关性较高，而低于0.5的则相关性较低，而相关系数为负的则不相关。

由表中数据可以得到：在池塘1中，总磷与磷酸盐磷、总氮相关系数较大；磷酸盐磷与总磷和总氮的相关系数较大；总氮和总磷、磷酸盐磷以及铵态氮相关系数较高；硝态氮与亚硝态氮的相关系数较大。因此，在池塘1中，总磷、磷酸盐磷、总氮显著相关；亚硝态氮、硝态氮显著相关。

表2：池塘2池水中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.49

6358

0.27

0388

0.07

9883

0.15

4939

-0.1

6755

磷酸盐磷

0.49

6358

1 0.14

1147

0.40

9524

0.46

3739

0.08

2663

总氮0.27

0388

0.14

1147

1 -0.1

1129

-0.2

0602

0.06

7132

硝态氮

0.07

9883

0.40

9524

-0.1

1129

1 0.80

3121

-0.0

2717

亚硝态氮

0.15

4939

0.46

3739

-0.2

0602

0.80

3121

1 -0.0

3923

铵态氮

-0.1

6755

0.08

2663

0.06

7132

-0.0

2717

-0.0

3923

1

由表中我们可以得到：在池塘2中的池水中，硝态氮与亚硝态氮相关系数较大，则硝态氮、亚硝态氮显著相关。

表3：池塘3池水中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.78

8287

0.54

6012

-0.0

8328

0.36

853

0.28

7627

磷酸盐磷

0.78

8287

1 0.57

6525

0.13

0808

0.69

6497

0.49

1296

总氮0.54

6012

0.57

6525

1 -0.1

7743

0.30

2094

0.05

8583

4

淡水养殖池塘水华发生及池水净化处理,网赛论文

硝态氮

-0.0

8328

0.13

0808

-0.1

7743

1 0.52

1686

0.32

549

亚硝态氮

0.36

853

0.69

6497

0.30

2094

0.52

1686

1 0.50

5326

铵态氮

0.28

7627

0.49

1296

0.05

8583

0.32

549

0.50

5326

1

由表中可以得到：在池塘3中的池水中，总磷、磷酸盐磷、总氮相关系数较大；磷酸盐磷、亚硝态氮相关系数较大；硝态氮、亚硝态氮相关系数较大；铵态氮、亚硝态氮相关系数较大。因此，总磷、磷酸盐磷、总氮显著相关；磷酸盐磷、亚硝态氮显著相关；硝态氮、亚硝态氮显著相关；铵态氮、亚硝态氮显著相关。

表4：池塘4池水中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.96

0771

0.76

4333

0.51

4673

0.56

8574

0.85

2377

磷酸盐磷

0.96

0771

1 0.76

0778

0.39

7906

0.48

2119

0.85

3247

总氮0.76

4333

0.76

0778

1 0.13

8018

0.31

6132

0.90

5362

硝态氮

0.51

4673

0.39

7906

0.13

8018

1 0.82

4655

0.28

407

亚硝态氮

0.56

8574

0.48

2119

0.31

6132

0.82

4655

1 0.45

5221

铵态氮

0.85

2377

0.85

3247

0.90

5362

0.28

407

0.45

5221

1

由表可以得到：在池塘4的池水中，总磷、总氮、磷酸盐磷、硝态氮、亚硝态氮、铵态氮相关系数较大；因此，总磷、总氮、磷酸盐磷、硝态氮、亚硝态氮、铵态氮显著相关。

表5：池塘1底泥中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.48

1908

0.59

3845

0.14

6465

-0.2

4423

0.20

8662

磷酸盐磷

0.48

1908

1 0.50

9889

0.23

2456

-0.0

8231

0.13

1286

总氮0.59

3845

0.50

9889

1 -0.2

8262

-0.2

1274

0.45

6345

硝态氮

0.14

6465

0.23

2456

-0.2

8262

1 0.00

3727

-0.2

112

亚硝态氮

-0.2

4423

-0.0

8231

-0.2

1274

0.00

3727

1 -0.1

1637

铵态氮

0.20

8662

0.13

1286

0.45

6345

-0.2

112

-0.1

1637

1

5

淡水养殖池塘水华发生及池水净化处理,网赛论文

由表中可以得到：在池塘1的底泥中，磷酸盐磷、总氮相关系数较大；总氮、总磷相关系数较大。因此，磷酸盐磷、总氮显著相关；总氮、总磷显著相关。

表6：池塘2底泥中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 -0.0

6218

0.35

7413

-0.3

6049

-0.1

1488

-0.0

3401

磷酸盐磷

-0.0

6218

1 0.32

0124

-0.4

6552

-0.0

5499

0.32

1014

总氮0.35

7413

0.32

0124

1 -0.7

6848

-0.5

0011

0.46

6109

硝态氮

-0.3

6049

-0.4

6552

-0.7

6848

1 0.23

9611

-0.5

5745

亚硝态氮

-0.1

1488

-0.0

5499

-0.5

0011

0.23

9611

1 -0.2

3727

铵态氮

-0.0

3401

0.32

1014

0.46

6109

-0.5

5745

-0.2

3727

1

由表中可以得到：在池塘2中，底泥中的主要理化因子大多负相关，即没有相关性，也有少许是正相关，但其相关系数很小。

表7：池塘3底泥中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.03

6951

0.08

8021

-0.1

2545

0.03

9383

0.20

9577

磷酸盐磷

0.03

6951

1 0.80

3188

0.09

345

-0.0

5587

0.32

22

总氮0.08

8021

0.80

3188

1 0.19

2693

-0.2

4988

0.43

1674

硝态氮

-0.1

2545

0.09

345

0.19

2693

1 -0.2

5643

-0.5

5889

亚硝态氮

0.03

9383

-0.0

5587

-0.2

4988

-0.2

5643

1 -0.1

616

铵态氮

0.20

9577

0.32

22

0.43

1674

-0.5

5889

-0.1

616

1

由表中可以得到：在池塘3中，磷酸盐磷、总氮相关系数较大，因此，磷酸盐磷、总氮显著相关。

表8：池塘4底泥中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 -0.0

029

0.16

1555

-0.0

9397

0.30

6827

0.15

4051

磷酸盐磷

-0.0

029

1 -0.0

8431

-0.5

0081

-0.2

9559

0.34

7434

总氮0.16-0.0 1 0.150.100.08

6

淡水养殖池塘水华发生及池水净化处理,网赛论文

1555 8431 8923 0991 367

硝态氮

-0.0

9397

-0.5

0081

0.15

8923

1 0.07

5438

-0.4

9824

亚硝态氮

0.30

6827

-0.2

9559

0.10

0991

0.07

5438

1 0.10

0489

铵态氮

0.15

4051

0.34

7434

0.08

367

-0.4

9824

0.10

0489

1

由表中可以得到：在池塘4中，底泥中的主要理化因子大多负相关，即没有相关性，也有少许是正相关，但其相关系数很小。

表9：池塘1间隙水中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.37

6621

0.53

5111

-0.0

6287

-0.3

899

0.61

7302

磷酸盐磷

0.37

6621

1 0.07

4662

-0.4

1363

-0.2

2696

0.09

0285

总氮0.53

5111

0.07

4662

1 -0.3

1275

-0.6

1391

0.28

4653

硝态氮

-0.0

6287

-0.4

1363

-0.3

1275

1 0.71

2261

-0.0

1146

亚硝态氮

-0.3

899

-0.2

2696

-0.6

1391

0.71

2261

1 0.01

5164

铵态氮

0.61

7302

0.09

0285

0.28

4653

-0.0

1146

0.01

5164

1

由表可以得到：在池塘1中，总磷、总氮、铵态氮相关系数较大；因此，总磷、总氮、铵态氮显著相关。

表10：池塘2间隙水中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.89

5115

0.51

9241

-0.2

6906

-0.3

0245

0.55

173

磷酸盐磷

0.89

5115

1 0.68

9465

-0.3

1863

-0.3

1741

0.77

2288

总氮0.51

9241

0.68

9465

1 -0.3

9573

-0.2

6918

0.68

079

硝态氮

-0.2

6906

-0.3

1863

-0.3

9573

1 0.94

7526

-0.2

7211

亚硝态氮

-0.3

0245

-0.3

1741

-0.2

6918

0.94

7526

1 -0.2

2274

铵态氮

0.55

173

0.77

2288

0.68

079

-0.2

7211

-0.2

2274

1

由表可以得到：总磷、磷酸盐磷、总氮、铵态氮相关系数较大；亚硝态氮、硝态氮相关系数较大。因此，总磷、磷酸盐磷、总氮、铵态氮显著相关；亚硝态氮、硝态氮显著相关。

7

淡水养殖池塘水华发生及池水净化处理,网赛论文

表11：池塘3间隙水中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.69

9411

0.44

4919

-0.3

4321

-0.2

7257

-0.0

6792

磷酸盐磷

0.69

9411

1 0.42

0465

-0.2

9155

-0.0

9431

-0.1

3973

总氮0.44

4919

0.42

0465

1 -0.1

1398

-0.0

6559

0.33

9116

硝态氮

-0.3

4321

-0.2

9155

-0.1

1398

1 0.68

8261

-0.1

626

亚硝态氮

-0.2

7257

-0.0

9431

-0.0

6559

0.68

8261

1 -0.0

9781

铵态氮

-0.0

6792

-0.1

3973

0.33

9116

-0.1

626

-0.0

9781

1

由表可以得到：总磷、磷酸盐磷相关系数较大；硝态氮、亚硝态氮相关系数较大。因此，总磷、磷酸盐磷显著相关；硝态氮、亚硝态氮显著相关。

表12：池塘4间隙水中的主要理化因子的相关系数

总磷磷酸

盐磷总氮硝态

氮

亚硝

态氮

铵态

氮

总磷 1 0.63

594

0.24

7643

-0.3

4721

-0.3

0808

0.28

5353

磷酸盐磷

0.63

594

1 0.22

7267

-0.2

445

-0.4

1077

0.01

3332

总氮0.24

7643

0.22

7267

1 0.36

5797

0.52

3698

0.06

4822

硝态氮

-0.3

4721

-0.2

445

0.36

5797

1 0.84

3259

-0.4

2254

亚硝态氮

-0.3

0808

-0.4

1077

0.52

3698

0.84

3259

1 -0.2

2215

铵态氮

0.28

5353

0.01

3332

0.06

4822

-0.4

2254

-0.2

2215

1

由表可以得到：总磷、磷酸盐磷相关系数较大；总氮、亚硝态氮相关系数较大；硝态氮、亚硝态氮相关系数较大。因此，总磷、磷酸盐磷显著相关；总氮、亚硝态氮显著相关；硝态氮、亚硝态氮显著相关。

5.2问题二模型的建立及求解

5.21模型的准备

问题二，我们对每一个观察站点进行评价，涉及溶解氧、总碱度、盐度、pH值等等因素。根据水池取样点报告中，采用均方差法确定各水池的各项指标的权重。对每一水池各项指标数据按公式求其均方差，将均方差归一化，

公式1

2/1

1

2

))

(

(

1

1

)

(?

?

?

?

?

?

-

-

=∑

=

n

i

ki

k

X

X

n

k

s，3,2,1

=

k，4

求其均方差，将均方差归一化，得权重指标（0.6075, 0.55,0.47,0.51325）。

8

淡水养殖池塘水华发生及池水净化处理,网赛论文

9 隶属度函数的定义

长江水质检测数据带有一定的误差，为了真正反映评价对象属于某一类水的程度，需要确定评价对象的隶属度，6个评价等级，以溶解氧为例，建立隶属度函数

?????≤<-≥=5.775.6,75.075.65.7,1)(1x x x DO f ， ???????<<--≤≤<<-=25.85.7,75.025.85.76,16

5.5,5.05.5)(2x x x x x DO f

???????<<--≤≤<<-=75.66,75.075.665,154,14)(3x x x x x DO f , ???????<<--≤≤<<-=5.55,5.05.55

3,135.2,5.05.2)(4x x x x x DO f

???????<<--≤≤<<-=43,1432,12

1,11)(5x x x x x DO f , ?????<<--≤≤=5.22,5.05.220,1)(6x x x DO f

另三个指标的同样定义隶属度。

5.22 计算单因素评判矩阵

对评价指标i x （1,2,3,4i =），第s 个站点属于第e 个评价等级的评价系数（隶属度）为∑==p

k k e s ie x f x 1)

()(（6,,2,1 =e ），k x 表示每一次检测的值。

对评价指标i x （3,2,1=i ），第s 个站点属于各个评价等级的总评价系数为

公式2

∑==61)()

(e s ie

s i x x

淡水养殖池塘水华发生及池水净化处理,网赛论文

10

对所有检测值就评价指标i x （3,2,1=i ），第s 个站点属于第e 个等级的隶属度为

公式3

)()()

(s i

s ie

s ie

x x r

= 第s 个站点就评价指标i x （3,2,1=i ）对于各类别的水质的单因素评判向量

公式4

),,,()(6)(2)

(1)(s i s i s i s i r r r r =

将第s 个站点的全部指标对于各评价等级的评价向量综合后，得到第s 个站点对于

各评价等级的评价矩阵：

公式5

?

???? ??=)(36)(31

)(16)(11)

(s s s s s r r r R 5.23 模糊综合评价

对第s 个站点做综合评价，模型为

公式6

),,,()

(6)(2)(1)()(s s s s s b b b R A B =?=

根据}{max )(6

1s j j b ≤≤确定第s 个站点近两年的水质情况所属类别。

模糊综合评判结果分析

各因素数据整理放入Excel 文件Data.xls 中，作为程序输入的数据文件，综合评判结果如下：

表13：4个池塘水质模糊综合评判结果

序号

池号 评价结果 （水质类别） 序号

池号

评价结果 （水质类别） 1 1#A 0.7551（1） 5 3#A 0.3652（2） 2 1#B 0.5784（1） 6 3#B 0.7975（1） 3 2#A 0.4708（2） 7 4#A 0.3504（2） 4

2#B

0.4976（2）

8

4#B

0.4413（2）

由表中数据我们可以得到如下结论：1号池塘是1类水质；2号池塘是2类水质；3号池塘是1类水质；4号池塘是2类水质。 5.3问题三模型的建立及求解

5.31模型准备 Ⅰ 单侧检验理论

淡水养殖池塘水华发生及池水净化处理,网赛论文

11

当所设H0为总体参数等于某一定值，而H1为仅从一个方向上偏离此定值者，为单侧检验。单侧检验有如下几个性质：

1．检验两组的差异显著性时，只考虑A ＞B 之意义，不考虑A ＜B 之可能性。 2．检验差值均数的显著性时，只考虑正值的意义，不考虑负值之可能性。 3．检验相关的显著性时，只考虑正相关的意义，不考虑负相关的可能性。

4．为检验多个实验组与一个对照组之差异而进行多重比较时，只考虑在一个方向上的差异，不考虑在另一方向上出现差异之可能性

Ⅱ 显著性检验理论

显著性检验就是事先对总体（随机变量）的参数或总体分布形式做出一个假设，然后利用样本信息来判断这个假设（原假设）是否合理，即判断总体的真实情况与原假设是否显著地有差异。或者说，显著性检验要判断样本与我们对总体所做的假设之间的差异是纯属机会变异，还是由我们所做的假设与总体真实情况之间不一致所引起的。

抽样实验会产生抽样误差，对实验资料进行比较分析时，不能仅凭两个结果（平均数或率）的不同就得出结论，而是要进行统计学分析，鉴别出两者差异是抽样误差引起的，还是由特定的实验处理引起的。

通常我们会引入p 值作为检验样本观察值的原假设可被拒绝的最小显著性差异水平。

公式7 ()0.05P X x ρ=<=为“显著(significant)” 公式8

()0.01P X x ρ=<=为“极显著(extreme significant)”

5.32模型建立与求解

因为问题中给出1号池发生了轻微水华，因此主要对1号水池进行具体分析就可以得出水华发生的具体区间从而解决问题。建立总磷、磷酸盐磷等六种理化因子与浮游生物密度的关系，对1号池塘每周的六种理化因子在A 和B 两个点的值求平均，得到数据表。

表14：理化因子与浮游生物密度总和

周数 浮游生物密度 总磷 磷酸盐磷 总氮

硝态氮 亚硝态氮 氨态氮 1 52.73 10.18905 0.0674 8.8743 1.1641 18.7008 18.7008 2 296.12 8.3721 0.0432 8.3905 1.2571 20.08185 20.08185 3 274.41 8.744 0.0676 7.81895 1.04295 23.6535 23.6535 4 95.94 9.36255 0.06485 8.3797 0.59715 25.23545 25.23545

5 231.88 8.413 0.0769 8.06715 1.4382 22.2213 22.2213

6 172.29 9.90535 0.0606 6.4773 1.755

7 30.3525

30.3525

7 505.65 5.5932 0.0646 2.75755 1.1598 11.26765 11.26765 8 185.53 6.9944 0.06645 4.10125 2.02835 16.26835 16.26835 9 36.24 9.90585

0.077

4.8385 2.31645 21.6749 21.6749 10 262.41 10.88045 0.06955 6.73425 2.5972 17.66835 17.66835

11 466.12 11.86245 0.15585 8.2514 0.80915 27.014 27.014 12 290 8.7985 0.06245 8.11765 0.7943 12.0263 12.0263 13 554.71 12.53835 0.0849 9.5958 1.43 26.76515 26.76515

14 1064.82 13.6464

0.19775

8.3491 1.651 24.0939 24.0939 15

954.71

14.26685 0.2136

7.3688

1.00515

19.7633

19.7633

淡水养殖池塘水华发生及池水净化处理,网赛论文

根据得出的表14，利用SPSS软件对得出的数据表进行单侧检验，分析相关性和显著性。（置信区间为0.95）得出相关性表15和检验统计表16

表：15相关性

浮游

生

物

密

度总

磷

磷

酸盐

磷

总

氮

硝

态氮

亚

硝态

氮

氨

态氮

浮游

生

物

密

度Pearson 相关性 1 .617 .847 .147 -.154 .037 .037 显著性（单侧）.007 .000 .301 .292 .448 .448 N 5 15 15 15 15 15 15

Bootstrapa 偏差0 -.047 -.059 .004 -.014 .002 .002 标准误差0 .284 .199 .225 .234 .202 .202

9 5% 置信区间下

限

1 -.240 .233 -.354 -.593 -.377 -.377

上

限

1 .890 .951 .57

2 .308 .431 .431

由上表可知浮游生物密度与总磷、磷酸盐磷、总氮、亚硝态氮和氨态氮呈现正相关和硝态氮呈现负相关。由表中数据可知磷酸盐磷的相关显著性P<0.01，其余五个理化因子的相关显著性P>0.05,所以浮游生物密度与磷酸盐磷的相关性极为显著与总磷、总氮、亚硝态氮和氨态氮的相关性不显著。

表16：描述性统计量

Statistic Bootstrap a

偏差标准误差95% 置信区间

下限上限

浮游生物密度均值362.9040 .8408 77.4991 213.6780 516.5483 标准差304.44268 -19.41991 68.19416 128.98332 393.91954 N 15 0 0 15 15

总磷均值9.9648333 -.0082094 .6001654 8.8155598 11.1249343 标准差 2.37958599 -.10987416 .40289684 1.44047918 3.03993168 N 15 0 0 15 15

磷酸盐磷均值.091513 -.000010 .012959 .068354 .119619 标准差.0525303 -.0036050 .0131662 .0100794 .0674941 N 15 0 0 15 15

总氮均值7.208147 -.013550 .484701 6.118631 8.060250 标准差 1.9167706 -.0914789 .3932895 .9619177 2.4989239 N 15 0 0 15 15

硝态氮均值 1.403107 -.000320 .145717 1.122757 1.692687 标准差.5755389 -.0313641 .0988979 .3357250 .7258219

12

淡水养殖池塘水华发生及池水净化处理,网赛论文

13

N

15 0 0 15 15 亚硝态氮

均值

21.119153

-.033654 1.376087 18.344239 23.737109 标准差 5.4086829

-.2089932

.8504720

3.4609874

6.8488482

N

15 0 0 15 15 氨态氮 均值

21.119153

-.033654 1.376087 18.344239 23.737109 标准差 5.4086829

-.2089932

.8504720

3.4609874

6.8488482

N

15

15

15

a. Unless otherwise noted, bootstrap results are based on 1000 bootstrap samples

由上表可以得出总磷、总氮等六种理化因子均值的置信区间的上下限。如下表：

表17：理化因子均值的置信区间的上下限

理化因子 置信下限 置信上限 总磷 8.8155598 11.1249343 磷酸盐磷 .068354 .119619 总氮 6.118631 8.060250 硝态氮 1.122757 1.692687 亚硝态氮 18.344239 23.737109 氨态氮

18.344239

23.737109

因为1号水池只发生了轻微水华，因此将1号水池考虑为发生水华的最低临界值，总磷、磷酸盐磷等六种因素在置信区间低于8.8155598、0.068354、6.118631、1.122757、18.344239、18.344239时不发生水华，在置信区间高于11.1249343、0.119619、8.060250、1.692687、23.737109、23.737109时水华严重，位于置信区间内时只发生轻微水华。因此我们推测当理化因子的在水池中的含量大于置信上限时发生严重水华；当理化因子的在水池中的含量小于于置信下限时不发生水华；当理化因子的在水池中的含量位于置信区间内时发生轻微水华。

5.4 问题四模型的建立与求解

对于问题四，要求建立体重与生长关系的模型，分析附件6里体重和体长与生长时间的关系，运用二次函数建立关于体重与生长关系的拟合方程，

2

(1)*.(2)*(3)y p t p t p =++

然后利用MATLAB 得出鱼类生体重与生长的关系的图像如下。（程序见附录13）

淡水养殖池塘水华发生及池水净化处理,网赛论文

14

根据所得图像忽略误差，得到体重与生长关系模型，由图像可知随着时间增加体重增长，并且边际体重增长率递减，并在后期趋于平稳。

题目要求构造一个与1号池相同大小的净化池，通过水循环，并放养鲢鱼或鳙鱼来净化1号池中的藻类。根据附件4中浮游植物密度的数据利用MATLAB 建立浮游植物密度与时间的关系模型，并得出如下图像。程序见附录

14

由附件7“某地区某湖水中各龄鲢鱼、鳙鱼饵料生物的季节性变化”得到鲢鱼与鳙鱼的饵料种类，发现鲢鱼的饵料种类包括需要净化的五种浮游植物，而鳙鱼的饵料种类仅包含硅藻、绿藻两种浮游植物，因此判定投放鲢鱼可达到净化浮游植物的最佳效果。

由浮游植物密度与时间图像，我们看出随时间变化浮游植物有一个最大峰值，查阅生物学中种群的增长与控制相关知识可知，当种群达到最大值时，选取种群最大值一半的时候，作为投放捕食者的最佳时间段，此时的浮游植物增长速率最快。假设放入的鲢鱼能够在最大生长速率时抑制浮游植物生长，即浮游植物生长速率=鱼类生长速率，假如不考虑环境因素和其因素的影响，我们近似地认为在投入鱼类的质量=鱼的饵料质量。

由附件4中的具体数据知在第14周时浮游生物密度达到最大值ρ=1064.82×

淡水养殖池塘水华发生及池水净化处理,网赛论文

10^6个/L，由附件5可知池塘1水体体积v=5.3×10^3×0.6=3.18×10^3L，由此得到浮游生物的个数N= ×v=3.386127×10^12个，查阅资料知微型种类只鉴别到门，按大、中、小级的平均质量计算。极小的(<5 µm)为0.0001 mg/104个；中等的(5 µm~10 µm)为0.002 mg/104个；较大的(10 µm~20 µm)为0.005 mg/104个，根据查阅资料可知，五种藻类的平均质量较大水平，因此浮游生物的质量m=162.794kg，由附件7中浮游生物占饵料的比例得出鱼的饵料质量，即鱼的质量M=237.65kg，因此应该在建立的净化池里放入237.65kg的鲢鱼能够得到净化水质的最佳效果。且在最佳时间段投入能够保持浮游植物生长速率=鱼类生长速率，构成一个水生态循环系统。

5.5问题五的模型建立与求解

通过网上查阅资料结合附件资料分析，可以得到有利于池水养殖池塘水体的自净化的方法为以下几点

(1)利用水生生物吸收利用氮、磷元素进行代谢活动以去除水体中氮、磷营养物质的方法。利用大型水生植物污水处理系统净化富营养化的水体。大型水生植物包括凤眼莲、芦苇、狭叶香蒲、加拿大海罗地、多穗尾藻、丽藻、破铜钱等许多种类，可根据不同的气候条件和污染物的性质进行适宜的选栽。水生植物净化水体的特点是以大型水生植物为主体，植物和根区微生物共生，产生协同效应，净化污水。

(2)投放河蚌、鲢鱼等水生动物净化水体。河蚌具有强大的滤水滤食功能，能有效降低氮、磷等含量。而鲢鱼的饵料成分中浮游植物占68%，可达到清除水华、净化水体的目的。

六、问题模型的评价与检验

6.1 模型的评价

6.11 模型的优点

从题意出发，充分运用MATLAB、SPSS、excel等软件绘制了大量的图表，使得数据分析和处理更加具备科学性和合理性。

在分析问题的过程中建立的合理的模型和方法，主要运用相关系数分析，判断各理化因子之间的差异性；运用隶属度函数、模糊综合和评价，判断出每个水池的水质类型；运用单侧检验，显著性检验，分析了浮游生物与理化因子的相关性；利用拟合检验建立了体重与生长关系模型和物种生态性构建了一个水生态循环系统。运用上述方法，使问题得到解决，文章以图表程序相结合，使问题分析更具体和形象。

6.12模型的缺点

显著性检验在一定程度上会产生抽样误差，由于定量数据较少，定性成分多，在比较分析实验资料时，不能仅凭部分结果的不同就得出结论，鉴别两者差异是抽样误差引起的，还是由特定的实验处理引起的还是需要大量的数据进行验证。

15

淡水养殖池塘水华发生及池水净化处理,网赛论文

在问题处理时，有些问题没有得到更好的解决。在处理数据和求解过程中不可避免各种误差，在一定程度上也影响到模型求解的精确度。

6.2 模型的检验

经初步检验每个问题所建立的模型基本符合建立模型的初衷，有助于合理地分析数据、处理数据，最终得出具有有效的结论。

七、模型推广

依据题目中提供的附件一到附件八，运用 MATLAB 与 SPSS 对提供数据进行数理统计与分析，建立了相关性分析模型，了解了各理化因子之间的相关关系，并用隶属度函数得到了不同池塘水质类型，还运用了单侧检验和拟合，建立了浮游植物关系模型和水生态循环系统。我们建立的模型可以拓展运用到农业水生态循环，鱼类养殖和水质鉴定等方面。

八、参考文献

[1].盛骤.谢式千.潘承毅《概率论与数理统计》单侧检验的相关性分析和置信区间 [教材]-浙江大学

[2].百度百科：富营养化

/link?url=OShRgznmQPuhTWIv5XGto0b5KpSC2VnF_fQLn97xTW9 LVOsTjv3Vz5afEzjdg4_GVm2dkLQVTGP6U0up6amWH_访问时间：2016.7.1

[3].《地表水环境质量标准》.2002.4.28.GB3838-2002

[4].刘峰养殖水体水华发生的原因探讨[M]，东北农业大学，第24-35页,2007年6月20日

[5].应光辉植物种群研究在群落生态学中意义，《现代农业科技》2009年第8期

[6].刘开旻吴小俊基于隶属度函数的模糊支持，江南大学物联网工程学院，2016年4月

[7].王喜刚扶名福基于隶属度函数模糊弹塑性的研究，南昌大学，辽宁科技学院，2015年

[8].杨朝凤单侧检验的相容性及应用，保山高等专科院校，计科系，2007年

[9].张天伟胡凤丽利用MATLAB建立不确定度评定测量分析变量的相关系数，伊犁州计量检定所，新疆，2016年1月

九、附录

附录1

clc;clear;

A1=xlsread('data1.xls',1,'C3:T17');

16

淡水养殖池塘水华发生及池水净化处理,网赛论文

A2=xlsread('data1.xls',1,'C19:T33');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[1:3:n]);

CA11=corrcoef(A11);

附录2

clc;clear;

A1=xlsread('data1.xls',1,'C37:T51');

A2=xlsread('data1.xls',1,'C53:T67');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[1:3:n]);

CA11=corrcoef(A11);

附录3

clc;clear;

A1=xlsread('data1.xls',1,'C71:T85');

A2=xlsread('data1.xls',1,'C87:T101');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[1:3:n]);

CA11=corrcoef(A11);

附录4

clc;clear;

A1=xlsread('data1.xls',1,'C105:T119');

A2=xlsread('data1.xls',1,'C121:T135');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

17

淡水养殖池塘水华发生及池水净化处理,网赛论文

A11=A1(:,[1:3:n]);

CA11=corrcoef(A11);

附录5

clc;clear;

A1=xlsread('data1.xls',1,'C3:T17');

A2=xlsread('data1.xls',1,'C19:T33');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[3:3:n]);

CA11=corrcoef(A11);

附录6

clc;clear;

A1=xlsread('data1.xls',1,'C37:T51');

A2=xlsread('data1.xls',1,'C53:T67');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[3:3:n]);

CA11=corrcoef(A11);

附录7

clc;clear;

A1=xlsread('data1.xls',1,'C71:T85');

A2=xlsread('data1.xls',1,'C87:T101');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[3:3:n]);

CA11=corrcoef(A11);

附录8

clc;clear;

A1=xlsread('data1.xls',1,'C105:T119');

A2=xlsread('data1.xls',1,'C121:T135');

A1=(A1+A2)/2;

[m,n]=size(A1);

18

淡水养殖池塘水华发生及池水净化处理,网赛论文

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[3:3:n]);

CA11=corrcoef(A11);

附录9

clc;clear;

A1=xlsread('data1.xls',1,'C3:T17');

A2=xlsread('data1.xls',1,'C19:T33');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[2:3:n]);

CA11=corrcoef(A11);

附录10

clc;clear;

A1=xlsread('data1.xls',1,'C37:T51');

A2=xlsread('data1.xls',1,'C53:T67');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[2:3:n]);

CA11=corrcoef(A11);

附录11

clc;clear;

A1=xlsread('data1.xls',1,'C71:T85');

A2=xlsread('data1.xls',1,'C87:T101');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[2:3:n]);

CA11=corrcoef(A11);

19

淡水养殖池塘水华发生及池水净化处理,网赛论文

附录12

clc;clear;

A1=xlsread('data1.xls',1,'C105:T119');

A2=xlsread('data1.xls',1,'C121:T135');

A1=(A1+A2)/2;

[m,n]=size(A1);

for i=1:n

a=min(A1(:,i));b=max(A1(:,i));

A1(:,i)=(A1(:,i)-a)/(b-a);

end

A11=A1(:,[2:3:n]);

CA11=corrcoef(A11);

附录13

clc;clear;clear all;

tz=xlsread('data1.xls',6,'B4:F522');

[M N]=size(tz);

t=[1:M];

for i=[2:M-1]

for j=[1:N]

if tz(i,j)==0

tz(i,j)=(tz(i-2,j)+tz(i+2,j))/2;

end

end

end

plot(t,tz(:,2),t,tz(:,N));

xlabel('时间(周)');ylabel('体重(g)');

legend('鲢鱼','鳙鱼',2);

p=polyfit(t,tz(:,2),2);

y=p(1)*t.^2+p(2)*t+p(3);

hold on

plot(t,y,'r');

hold off

附录14

x=xlsread('data1.xls',4,'H5:H64');

t=[1:15]';

plot(t,x(1:15,1));

xlabel('时间(周)');ylabel('浮游植物密度（10^6个/L）'); legend('1号池');

20

本文来源：https://www.bwwdw.com/article/iaj1.html