【MathorCup】第1问论文..docxVIP

第1问1 问题分析氮在水中主要以离子形态存在，主要包括硝态氮、铵态氮和亚硝态氮，也有相当部分的氮以有机质的形式存储在水体有机物中。水中的各种含氮物质通过氨化、硝化、反消化和同化作用等进行相互转换。水体中的磷则以活性磷酸盐为主体，存在部分有机磷和无机磷颗粒。在沉积物中，氮主要以含氮有机物形式存在，并通过微生物作用转变为无机氮向水体扩散。而磷则无机与有机参半，其中能溶解并与水体交换的磷是各种评价的主要研究对象。[1]在了解水体氮磷基本情况的前提下，探究池塘中常见N、P理化因子的关系，应当包含两个维度。一个维度是横向探究各种理化因子自身间的相互影响，相互转换关系。另一个维度是纵向各种理化因子在沉积物、间隙水、水体间的扩散转移。对于各维度的关系，可以在相对忽略另一维度关系影响的情况下，单独讨论。附件一中包含的数据为“连续15周，对1到4号养殖池，分别在A、B两处采集底泥、间隙水、池水三种环境中的样本，并检测包括总氮与总磷在内的6种理化指标的含量值”。因为理化因子间的相互关系具有稳定性，因此抽取1号养殖池数据分析，其结果可以代表整体。对于A、B两个采样点的数据取均值，作为被测理化因子的实际含量。在横向探究理化因子间的自身影响时，采用主成分分析法求三种环境下各理化因子的综合指标，以消除理化因子纵向扩散转移关系对横向关系探究的影响。然后对6项理化因子综合指标做相关性分析，初步判别个理化因子间的相关关系。最后对，具有较强相关性的理化因子做进一步函数拟合，解读其相关关系。纵向探究不同环境中理化因子的扩散转移关系的思路与方法与横向探究一致。2 数据处理与分析2.1 综合指标提取求综合指标时使用主成分分析法，其基本步骤为“数据标准化”，”建立变量相关系数矩阵R“，“求R的特征根??1, ??2, ??3……及其对应的特征向量T1*, T2*, T3*……”，“由特征根及特征向量求主成分表达式与综合评价函数”，“带入标准化的原始数据，求综合评价值”。值得注意地是，从本质上讲，综合评价函数数对原始指标的线性综合。从计算主成分到对之加权，经过两次线性运算后得到综合评价函数。求综合指标的Matlab代码在附录1中给出，这里仅给出计算结果。表1为横向探究中的理化因子综合指标，表2为纵向探究中的池塘环境综合指标。 表1 理化因子综合指标周数总磷磷酸盐磷总氮硝态氮亚硝态氮铵态氮10.1790.16010.1913-0.08890.0724-0.174420.12220.05660.1802-0.16710.0825-0.144530.15330.10320.1324-0.11010.1266-0.190440.15530.08710.18680.00710.1029-0.205550.16170.10480.1828-0.11870.0326-0.184160.1630.09450.1427-0.24470.0622-0.262270.11870.09240.07550.00620.0852-0.096580.16170.13840.1457-0.15420.3066-0.165490.18390.24410.2173-0.34530.8489-0.2543100.16210.14540.1619-0.52430.3524-0.1659110.17490.16120.2281-0.08980.0955-0.1869120.17750.17930.1827-0.0830.0446-0.114130.20790.27230.2446-0.11350.0963-0.2221140.21530.35840.2596-0.01640.0711-0.1896150.2270.44160.2129-0.00770.0752-0.1764表2 池塘环境综合指标周数池水间隙水底泥周数池水间隙水底泥10.2528-0.2109-0.330290.3534-0.0819-0.363820.0714-0.2383-0.311410-0.002-0.2551-0.360930.0974-0.2937-0.3611110.1125-0.4581-0.34140.0968-0.3002-0.328120.2319-0.1433-0.27485-0.0471-0.2742-0.3503130.2933-0.235-0.39966-0.0483-0.3176-0.3335140.4332-0.4387-0.464370.062-0.1647-0.3239150.5473-0.4504-0.392480.139-0.1213-0.36012.2 相关性分析利用malab中的相关性分析函数corr()分别对理化因子综合指标和池塘环境综合指标进行相关性分析，结果分别如表3和表4所示。表3 理化因子综合