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

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

6.2 模型的检验

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

七、模型推广

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

八、参考文献

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

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

http://baike.http://www.wodefanwen.com//link?url=OShRgznmQPuhTWIv5XGto0b5KpSC2VnF_fQLn97xTW9LVOsTjv3Vz5afEzjdg4_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