数学建模社区-数学中国

标题: 东三省二等奖论文 [打印本页]

作者: 刘李 时间: 2013-1-23 20:18
标题: 东三省二等奖论文
深圳市人口与医疗需求预测
摘要
本文是一个如何准确预测深圳市人口数量和年龄结构，并结合深圳市疾病别的患病情况，来预测深圳市未来十年医疗床位需求的问题。以概率论和数理统计为基础，我们建立了基于灰色预测和神经网络的组合模型并结合人口动力发展方程来做出比较精确的预测。
针对问题(1)，即深圳市未来十年的总人口预测，我们从《深圳市统计年鉴》上搜集的部分数据，利用灰色预测建模所需信息少、方案简单的特点和神经网络具有较强的非线性映射能力的特征，提出一种基于灰色预测和神经网络的人口预测方案。首先对人口规模的灰色预测GM（1,1）等维新息模型预测，然后再用优化的BP神经网络对GM（1,1）模型的预测值进行组合预测以作为其最终的预测值，以获得更高的预测精度。
针对问题(2)，即深圳市全市的医疗床位需求预测，我们根据历年的出生率、死亡率、总和生育率等指标，提出了人口动力发展方程，来预测深圳市未来的人口年龄结构。
针对问题(3)，即不同疾病在不同机构的床位需求，我们选取了高血压和分娩这两种具有代表性的疾病，按卫生部推荐的卫生服务需求量法进行卫生资源配置标准的测算，我们倾向于从提高和促进人群健康水平的高度和卫生资源公平与效率的角度，确定区域卫生资源的数量、质量、结构、层次、布局以及发展速度与规模，并结合了深圳市的人均可支配收入、产业结构的调整所带来的影响以及居民对不同医疗机构的选取指标来作出预测，而不是仅仅站在卫生机构生存和发展的立场上加以考虑。
关键词：GM(1,1)灰色预测  人工神经网络  遗传模拟退火算法(GASA)  卫生服务需求法

1．问题的重述
1.1  问题背景
深圳是我国经济发展最快的城市之一，30多年来，卫生事业取得了长足发展，形成了市、区及社区医疗服务系统，较好地解决了现有人口的就医问题。
根据深圳市第六次人口普查的数据，“深圳常住人口十年来年均增长3．98％，低于上一个十年（1990－2000）6．34％的增长，但仍快于全国（0．57％）、全省（1．90％）的平均水平。”
过去十年，深圳人口仍增长较快，主要原因有两个：一是人口流入量较大。十年来，深圳经济保持较高增速，提供了较多的就业岗位，吸引了大量市外人员到深圳工作和生活。深圳常住人口中，约有798万人是非户籍人口，占常住总人口的77％；二是户籍人口增长较快。十年间，深圳采取多项措施，降低了入户门槛，加快了户籍人口增长步伐。普查登记的户籍人口达251．03万人（未剔除流出市外半年以上人口），与第五次全国人口普查的121．48万人相比，10年共增加了130万人，增长1．1倍，年平均增加13万人，年平均增长7．53％。老年人口占总人口比例呈正增长趋势。因为深圳市是我国改革开放以后首批开始发展的新型城市，建市至今将近30年，人口逐渐凸显老龄化趋势。人口的增长，尤其是人口老龄化的加剧，必然导致现有的医疗资源的不足。国内外有关研究表明，人均医疗资源的需求与不同年龄的人群有关，一般情况下，60岁以上年龄组的医疗需求是60岁以下人群的3～5倍。由于老年人口是发病率和患病率最高的人群，所以老年人口的增加，意味着医疗卫生需求大幅增加。未来的医疗需求与人口结构、数量和经济发展等因素相关，合理预测能使医疗设施建设正确匹配未来人口健康保障需求，是保证深圳社会经济可持续发展的重要条件。人口预测是否科学准确，直接关系到其总体规划方案是否合理和实用。
对深圳市未来人口的准确预测，能够为深圳市医疗资源的合理、有效配置提供科学依据，这对加速推进深圳市的经济发展和全市消费有效需求的利用有着极为重要的现实意义。因此，根据已有数据，运用数学建模的方法，对深圳市人口数量和结构做出分析和预测，进而分析深圳市的医疗资源（主要是床位）的需求是一个重要问题。
1.2  问题提出
过去十年，深圳市的人口发展出现了一些新的特点，例如：老龄化进程加速，出生人口性别比持续升高，具体表现为男多女少，家庭规模偏小以及乡村人口城镇化等因素，这些都影响着深圳市人口的增长。
关于人口问题已有多方面的研究，并积累了大量数据资料。现在得到了《深圳市人口统计年鉴》上收集到的部分数据，其中包含1979至2010年的全市常住人口的数量、死亡率、生育率、男女出生比例、育龄妇女生育率的千分比( ‰) ；以及从国家卫生部公开信息网站、深圳市卫生和人口计划生育委员会网站与深圳市统计局官网上查询的其他相关数据。
根据已知数据或搜索相关文献和补充新的数据，解决以下的问题：从深圳市的实际情况和人口增长的上述特点出发，建立深圳人口增长的数学模型，并由此对深圳市未来十年人口增长的短期做出预测，并以此为依据预测深圳市不同医疗机构的资源（主要是床位数）的合理配置。
预测深圳市某一种疾病在不同类型医疗机构的床位需求就要预测该种疾病总的床位需求，而决定床位需求的就是该种疾病的患病人口，也就是深圳市该种疾病的患病率。而患病率又和人口结构有着密切的关系，因为每一种疾病都会有一个高发人群年龄段，例如，心脏病、高血压等疾病多发生在老年人中，而小儿麻痹发生在少年儿童中，因此改年龄段人口的比例严重影响着该种疾病患病率。因此，通过分析深圳市历年该种疾病患病率与人口结构的关系，回归出患病率函数，预测未来几年深圳市某种疾病患病率，并且通过预测人口求出患病人数，进而求出该种疾病的总床位需求。
同时根据高血压在老年人中的发病率较高这一现象，我们在已知高血压这一疾病的两周患病率和慢性病患病率的数据的前提下，对该疾病的床位需求进行预测是很有必要的。
2. 模型的假设
( 1) 考虑所研究的人口为常住人口；
( 2) 深圳市各区人口体质相同，即同一年度各区患病率相同，且各区相对封闭，本区人口不会外出跨区就医。
( 3) 人口数随时间变化是连续的，而且充分光滑。
( 4) 农村人口一旦迁入城镇或者城镇化, 其人口行为和特征即与城镇人口相同，即忽略城镇人口与迁入城镇人口或城镇化人口的差别.
( 5) 人口数据质量高，无误报和漏报等。
( 6) 在分析老年人口问题时，是以65岁作为老年人口的起点年龄；
( 7) 病床使用率不变，即未来的病床使用率沿用现在的标准病床使用率的数据，并且对现在病床使用率进行评价。
3．相关名词解释
人口数  指一定时点，一定地区范围内有生命的个人总和。
出生率  指在一定时期内( 通常为一年) 一定地区的出生人数与同期内平均人数( 或期中人数) 之比，用千分率表示。其计算公式为：

死亡率  指在一定时期内( 通常为一年) 一定地区的死亡人数与同期内平均人数( 或期中人数) 之比，用千分率表示。其计算公式为：

婴儿死亡率  指年内一定地区未满1岁婴儿死亡人数与同年出生的活产数之比，一般用‰表示。
人口自然增长率  是指在一定时期内( 通常为一年) 人口自然增加数( 出生人数减死亡人数) 与该时期内平均人数( 或期中人数) 之比，用千分率表示。计算公式为：

预测时期  短期( <10年)，中期( 10-25年)，长期( >25年) ；
出生人口性别比  是活产男婴数与活产女婴数的比值，通常用女婴数量为100时所对应的男婴数来表示。正常情况下，出生性别比是由生物学规律决定的，保持在103～107之间。
人口抚养比  指人口总体中处于供养年龄( 一般指15岁以下和64岁以上) 的人口与处于“经济活动”年龄( 15-64岁) 人口的比率。用百分比表示，计算公式为：

总和生育率：一定时期( 如某一年) 各年龄组妇女生育率的合计数，说明每名妇女按照某一年的各年龄组生育率度过育龄期，平均可能生育的子女数，是衡量生育水平最常用的指标之一。
总诊疗人次数: 所有诊疗工作的总人次数。诊疗人次数按挂号数统计，包括：①病人来院就诊的门诊、急诊人次；②出诊人次数；③单项健康检查及健康咨询指导人次；④未挂号就诊、本单位职工就诊及外出诊疗不收取挂号费的，按实际诊疗人次统计。患者一次就诊多次挂号，按实际诊疗次数进行统计，不包括根据医嘱进行的各项检查、治疗、处置工作量。
床位数：年底固定实有床位(非编制床位)，包括正规床、简易床、监护床、正在消毒和修理床位、因扩建或大修而停用的床位，不包括产科新生儿床、接产室待产床、库存床、观察床、临时加床和病人家属陪侍床。
出院人数：所有住院后出院的人数，包括治愈、好转、未愈、死亡及其他人数。其他人数指正常分娩、未产出院、住院经检查无病出院、未治出院及健康人进行人工流产或绝育手术后正常出院者。
实际占用总床日数：医院各科每日夜晚12点实际占用病床数(即每日夜晚12点住院人数)总和。包括实际占用的临时加床在内。病人入院后于当晚12点前死亡或因故出院的病人, 作为实际占用床位1天进行统计,同时亦应统计“出院者占用总床日数”1天,入院及出院人数各1人。
出院者占用总床日数：出院者(包括正常分娩、未产出院、住院经检查无病出院、未治出院及健康人进行人工流产或绝育手术后正常出院者)住院日数的总和。

潜在需求增长率：由于居民收入增长而导致的对医疗需求的增长，用百分比表示，计算公式为：

医院病床配置推荐标准：

病床使用率：

病床周转次数：

病床工作日：

出院者平均住院日：

居民年实际住院率：

居民两周就诊率：调查前两周内居民因病或身体不适到医疗机构就诊的人次数与调查人口数之比。
居民两周未就诊率：调查前两周内居民患病而未就诊的人次数与两周患病人次数之比。
居民两周患病率：

年平均床开放日数：

年净流入住院人数：

期望寿命  又称平均期望寿命。指0岁时的预期寿命。一般用“岁”表示。即在某一死亡水平下，已经活到岁年龄的人们平均还有可能继续存活的年岁数。
高、低限计算依据：《标准》中配置的低限为实际需求的计算结果，高限是按实际需求加上可能转化的潜在需求计算结果。
城镇居民可支配收入：将家庭总收入扣除交纳的个人所得税和个人交纳的各项社会保障支出后的收入，等于工资性收入、经营净收入、财产性收入、转移性收入之和，再减去交纳的个人所得税和个人交纳的社会保障支出。目前城镇居民可支配收入中只包括城镇居民以现金形式获得的收入。
人均可支配收入：居民可支配收入的算术平均数，即样本家庭可支配收入总额除以其家庭人口总数计算得出。
4. 模型的建立与求解
4.1  预用模型介绍
传统的人口预测方法主要有逻辑方法、常微分方程方法和动态预测方法等。这些方法在人口预测领域起到了一定的作用，但采用这些方法时都要对数据进行模型假设。由于真实模型往往是非线性的，如果在一些简单的模型假设下就进行数据模拟，常常不能达到较好的模拟效果。灰色系统理论中的灰色预测GM（1,1）模型因其所需信息少、运算方便、建模精度较高而被广泛应用于各种预测领域。近年来该模型已应用于人口规模预测，但其精度不高。人工神经网络具有高度的非线性映射能力，它能以任意精度逼近任意非线性函数，因此，它比较适合于一些复杂问题的建模。
4.1.1  GM(1,1)预测模型的基本原理
GM(1,1)灰色系统
所谓灰色系统是指既含有已知信息，又含有未知信息的系统，是由邓聚龙教授在1986年提出的。灰色理论自诞生以来，发展很快，由于它所需因素少，模型简单，特别是对于因素空间难以穷尽，运行机制尚不明确，又缺乏建立确定关系的信息系统，灰色系统理论及方法为解决此类问题提供了新的思路和有益的尝试。
灰色预测方法是根据过去及现在已知的或非确知的信息，建立一个从过去引申到将来的GM模型，从而确定系统在未来发展变化的趋势，为规划决策提供依据。在灰色预测模型中，对时间序列进行数量大小的预测，随机性被弱化了，确定性增强了。此时在生成层次上求解得到生成函数，据此建立被求序列的数列预测，其预测模型为一阶微分方程，即只有一个变量的灰色模型，记为GM(1,1)模型。
(2) GM(1,1)预测模型的基本原理
GM(1,1)模型是灰色预测的核心，它是一个单个变量预测的一阶微分方程模型，其离散时间响应函数近似呈指数规律。建立GM(1,1)模型的方法是：
设为原始非负时间序列，为累加生成序列，即
            (4-1)
GM(1,1)模型的白化微分方程为：
                                             (4-2)
式(2) 中，为待辨识参数, 亦称发展系数；为待辨识内生变量，亦称灰作用量。设待辨识向量，按最小二乘法求得式中

于是可得到灰色预测的离散时间响应函数为：
               (4-3)
为所得的累加的预测值，将预测值还原即为：
         (4-4)
4.1.2  BP神经网络模型的基本原理
4.1.2.1 神经网络的定义简介
神经网络是由多个神经元组成的广泛互连的神经网络，能够模拟生物神经系统真实世界及物体之间所做出的交互反应。人工神经网络处理信息是通过信息样本对神经网络的训练，使其具有人的大脑的记忆，辨识能力，完成名种信息处理功能，它不需要任何先验公式，就能从已有数据中自动地归纳规则，获得这些数据的内在规律，具有良好的自学习，自适应，联想记忆，并行处理和非线性形转换的能力，特别适合于因果关系复杂的非确定性推理、判断、识别和分类等问题. 对于任意一组随机的，正态的数据，都可以利用人工神经网络算法进行统计分析，做出拟合和预测。
基于误差反向传播(Back propagation)算法的多层前馈网络(Multiple-layer feedforward network, 简记为BP网络)，是目前应用最成功和广泛的人工神经网络。
4.1.2.2  BP模型的基本原理
学习过程中由信号的正向传播与误差的逆向传播两个过程组成。正向传播时，式作用于输入层，经隐含层处理后，传入误差的逆向传播阶段，将输出误差按某种子形式，通过隐含层向输入层逐层返回，并“分摊”给各层的所有单元，从而获得各层单元的参考误差或称误差信号，以作为修改各单元权值的依据。权值不断修改的过程，也就是网络学习过程。此过程一直进行到网络输出的误差准逐渐减少到可接受的程度或达到设定的学习次数为止。BP网络模型包括其输入输出模型，作用函数模型, 误差计算模型和自学习模型。
BP网络由输入层，输出层以及一个或多个隐含层节点互连而成的一种多层网，这种结构使多层前馈网络可在输入和输出间建立合适的线性或非线性关系，又不致使网络输出限制在-1和1之间。见图(4-1)

                        O          O       O             O

( 大于等于一层) W(1)…W(L)

4.1.2.3  BP神经网络的训练
BP算法通过“训练”这一事件来得到这种输入，输出间合适的线性或非线性关系。“训练”的过程可以分为向前传输和向后传输两个阶段：
1．向前传输阶段：
①从样本集中取一个样本，将输入网络；
②计算出误差测度和实际输出；
③对权重值各做一次调整，重复这个循环，直到。
2．向后传播阶段——误差传播阶段：
①计算实际输出与理想输出的差；
②用输出层的误差调整输出层权矩阵；
③ ；
④用此误差估计输出层的直接前导层的误差，再用输出层前导层误差估计更前一层的误差。如此获得所有其他各层的误差估计；
⑤并用这些估计实现对权矩阵的修改，形成将输出端表现出的误差沿着与输出信号相反的方向逐级向输出端传递的过程。
网络关于整个样本集的误差测度：

4.1.2.4  基于GASA混合学习策略的BP神经网络优化
鉴于BP算法易陷入局部极小，遗传算法（Genetic Algorithm，GA）虽有较强的全局搜索性能，但在实际应用中容易产生早熟收敛的问题，而改善其性能需增大种群数目或采用高级的进化策略，因此本文结合GA并行搜索结构和模拟退火（Simulated Annealing，SA）的概率突跳特性，提出一种神经网络权值学习的GASA混合学习策略以避免网络陷于局部极小，并将这种混合学习策略用于NARMA（p,q）递归网络和用作组合预测的多层前馈网络的学习。GASA混合策略在神经网络中学习步骤如下(流程图见图(4-2))：

图(4-2)  GASA混合策略在神经网络中学习步骤
(1)确定神经网络的结构，并用浮点数编码表示该网络。
鉴于浮点数编码更接近问题的本身，且在 v解的质量和算法效率方面均优于二进制编码，所以用基因串表示第层和层的权值向量，为第层的阈值向量，用浮点数编码来表示层BP网络。基于训练误差的单目标学习常会出现过拟合现象，影响网络的推广性，因此，采用训练样本集的训练误差和测试样本集的测试误差的综合误差指标作为网络性能的评价目标，适应度函数定义为：
     (4-5)
式中为学习误差，计算公式为

式中: 和分别为训练样本集和测试样本集中的样本数，才和，分别为训练样本或测试样本集第个样本输入对应网络输出层第个结点的期望输出和实际输出。
选择操作是按照优胜劣汰、适者生存的原则选择交配个体以产生子个体的过程。本文采用截断选择法，在该方法中，个体按适应度由大到小排序，只有截断阈值之上的优秀个体能够被选择作为父个体，在该阈值之下的个体不能产生子个体。
浮点数编码遗传算法(float-encoding genetic algorithm, FGA)的常用交叉算子有多种，如简单交叉、算术交叉、启发交叉、混合交叉和平坦交叉等。本文采用混合交叉算子。设，在交叉算子作用下产生两个新个体和，混合交叉算子是对的每个分量对，设，在区间中按均匀分布选择两个随机数作为，这里为事先给定的常数，一般取。
FGA常用变异算子有均匀变异、边界变异和高斯变异等。本文采用均匀变异的。
(2)确定种群大小、进化最大代数、选择操作时的淘汰率、变异概率、初温和退温速率，初始化种群，令。
(3)对中的各个个体进行SA搜索：
(3.1)利用SA状态产生函数产生新的个体。
(3.2)计算新、旧个体的综合误差之差。
(3.3)计算接受概率。
(3.4) 若，用新个体取代旧个体；否则，旧个体不变。
(4) 对中的个体按适应度由大到小进行排序，然后以淘汰率淘汰掉适应度值低的个体，被淘汰掉的个体数为，再将未被淘汰掉的优良个体保留到下一代种群中。
(5) 次从的优良个体中随机选择两个个体进行交叉操作，每次产生一个新个体填补到中以保持种群规模不变。
(6) 以变异概率对步骤(5)产生的新个体进行变异操作,至此产生种群。利用退温函数进行退温，交换训练样本集和测试样本集。
(7)若最优个体的综合误差满足要求精度或已达到进化最大代数，则终止算法并输出结果；否则，令，转到步骤(3)
4.2  针对问题（1）
4.2.1 预测深圳是未来十年全市总人口所采用模型----GM(1,1)模型和BP神经网络模型的组合预测
尽管目前时序预测的方法很多，但它们的应用条件、建模机理各异，都存在一定的局限，因此，在实际应用中传统的单一预测方法难以获得较为满意的预测结果。如果将各种预测方法进行组合，充分利用各自包含的一些有用的信息，扬长避短，将能产生较好的预测效果。自J.M.Bates和C.W.Granger(1969)首次提出组合预测以来，组合预测的方法和应用被广泛研究，但这些方法的大部分都侧重于预测误差方法的最小。文献 [1]用各种预测方法如平均值法、递归最小二乘法、模糊逻辑系统、神经网络法等对煤气日用量的无隐含层前馈网络和单隐层前馈网络的预测值进行了组合预测的研究，研究结果表明神经网络的组合预测效果优于其它方法。鉴于此，本文将人口规模的GM（1,1）等维新息模型的预测值和BP神经网络进行组合预测并求出其最终的预测值。用于人口预测的灰色预测和神经网络的组合模型的基本方法与步骤为：
(1)以人口的历史统计数据建立GM(1,1)等维数为的等维新息模型，其灰色预测值为。
(2)将灰色预测值作为用于组合预测的多层前馈网络的输入，其对应的实际值作为该网络的输出构成匹配对( , )然后将这些匹配对分为两部分，A部分用于网络的训练，B部分用于测试该网络的性能。
(3)如同递归网络，用GASA混合学习策略对A部分样本进行组合预测网络FNN的学习，然后作一步预测，从而得到人口规模的最终预测值。
4.2.2  未来十年深圳市总人口预测结果分析

表4-1  2011-2020年深圳市总人口预测

单位人口预测值（万人）
2011 1083.018
2012 1131.64
2013 1182.52
2014 1235.3
2015 1289.743
2016 1345.691
2017 1403.042
2018 1461.729
2019 1521.705
2020 1582.941

图(4-3)深圳市未来十年总人口预测图
按此预测，从以上图表中可以知道深圳市人口增长的中短期趋势为：深圳市人口自然增长率继续保持平稳下降，人口总量在中短期内继续增长，增速较为平稳。人口总量在未来10年还将净增500万人左右，总人口将于2020年达到1582.941万人。
单独的GM(1,1)模型和神经网络模型都能非常精确的预测出在短期内深圳市未来人口总数，在中期预测中，也能较好的反映出人口变化趋势。在长期预测时，GM(1,1)的预测数据偏差非常大，神经网络虽然变化趋势较符合深圳市未来人口变化情况，但偏差还是较大。而组合模型在这三个时期内，都能够很好的预测出深圳市未来人口总量的及其变化情况，显然组合模型的预测效果最佳。关于GM(1,1)长期预测的局限性将在模型评价中讨论。
4.3  针对问题（2）
4.3.1 预测深圳市未来十年人口年龄结构所采用模型----人口发展动力方程（部分指标值的预测见附表1）
4.3.1.1  人口的年龄结构
将人口按年龄大小逐年分成个年龄组，模型要讨论在不同时间人口的年龄分布，对时间也加以离散化，其单位与年龄组的间隔相同。时间离散化为 .设在时间段第年龄组的人口总数为，定义向量，根据分年龄的人口结构递推公式进行预测，模型的数学表达如下：

其中：

式中：为第年常住人口0岁出生婴儿数，为第年的第岁的常住人口数，为第年到第年迁入的总人口，为第年到第年迁入的岁人口，  为第年容纳的最大迁入人口，为第t年到第t+1年的迁出的总人口，为第年到第年的迁出的岁的人口，为年的迁入率（在一定时期内可以认为不变），为年的迁出率（在一定时期内可以认为不变），为第年常住人口出生婴儿当年死亡率，为常住人口中妇女总和生育率，即社会人中平均意义下一个妇女在整个育龄时期的生育总数（，即为生育年龄的上下限），为常住人口中第年年龄为的女性性别比，为第年第岁的常住人口死亡率，表示生育模式，反映某一地区某一个育龄妇女生育状态分布；表示迁入模式，反映某一地区迁入人口随迁入年龄的状态分布，表示迁出模式，反映某一地区迁出人口随迁出年龄的状态分布。
假定一：
三个模式在一定时期内稳定，生育模式采取对数正态分布，迁入模式和迁出模式均采取分布，
即：

其中；取总和生育率，为最低生育年龄，的取值范围为(1.7,2.5)，的取值范围为(0,1)，在本文中取2.4，取0.5，拟合效果比较理想。
其中期望生育年龄为，经计算等于26.49，比较符合实际。

其中为迁入峰值年龄，为迁入人群年龄下限，为可控变量，根据拟合可以得出其值约为0.67；

其中为迁出峰值年龄，为迁出人群年龄下限，为可控变量，根据拟合可以得出其值约为1.4；
假定二：
迁入人口的阻尼程度与总人口的阻尼程度成正比；迁出人口的逆阻尼程度与总人口的阻尼程度成正比。
即：

人口的阻尼程度利用经验公式：
;
4.3.1.2  人口的死亡率
人口的死亡率的计算公式为：

其中死亡率表示第岁的人口所占总人口的比例，表示第岁人口的分年龄层死亡率。
其中分年龄的死亡模式计算公式：

因而，预测未来十年的人口分年龄死亡模式为：

在这里，
(1)我们假定环境因子不变，设定为1；
(2)收入因子的计算公式为：，其中人均可支配收入的增长速率以最近五年的平均速率来计算。用反衰减因子法预测未来的人均可支配收入。即：

则，未来的人均可支配收入的计算公式为：
(3)居民健康水平：
4.3.1.3  人口的性别比
性别比的预测公式：

其中第三产业的增长速率数据可以从深圳市国民经济和社会发展统计公报获知。
4.3.1.4  各项预测指标调整的方法
1．先由已知的各指标值的数据，算出其均值和标准差，分别记为记为，我们定义一个相似随机因子，计算公式为。
2．对已预测出的数据进行随机操作：记以预测数的数据所组成的数列为，
①求出已预测出的数据的均值和标准差，分别记为和；
②调整后的数列为，则，调整后的预测值所组成的数列的计算公式为：

4.3.2 未来十年深圳市人口年龄结构预测结果分析
4.3.2.1  深圳市年龄结构变化趋势

图(4-4)  深圳市未来十年年龄别人口预测图

图(4-4-1)  深圳市未来十年中青年人口发展趋势

图(4-4-2)  深圳市未来十年老年人口发展趋势
按此预测，有以上图表可知，深圳市未来十年年龄别人口的发展趋势是：中青年的流动性大，人口结构相对不稳定，易受外部影响（政府的调整、经济的吸引力等）；而老年人口在缓慢增加，导致老龄化在平稳中逐年加剧。
但是人口老龄化将导致抚养比不断提高，尤其是乡老年人口变化显著，老龄化速度仅仅是深圳市在10年就增加了0.7个百分点。
老年人口比重的不断提高和老年人口绝对数的大量增加，在如何搞好老年人口的衣食住行和老有所养、老有所乐、老有所为等方面是社会不容忽视的问题。
4.3.2.2  深圳市人口分年龄死亡率变化
人口死亡水平是造成人口数量和结构变化的重要因素之一，随着深圳市经济的迅速发展和医疗卫生条件的进一步改善，深圳市人口的健康水平有了明显改善, 人口的死亡率总体上较低。

图(4-5)  深圳市未来十年人口年龄别死亡率变化预测图

图(4-5-1)  深圳市未来十年中青年人口死亡率变化预测图

图(4-5-2)  深圳市未来十年老年人口死亡率变化预测图
按此预测，由以上图表可知，深圳市未来十年青少年的死亡率在逐年下降，这与居民生活水平的提高、生活质量的改善有关，从而导致患病率降低；小中青年和青壮年的死亡率逐年递增，与外来务工人员的不断流入，和从事的工作（体力、高风险等）性质有关。
在刚进入老年时期时，死亡率不增反降，与居民的生活水平的提高和刚退休生活压力小有关，而后上升幅度较大，是由于居民的年龄增大和身体素质的下降。
4.3.2.3  深圳市性别比与第三产业关系

图(4-6)  深圳市人口性别比与第三产业关系图
图中居上的一条线代表人口性别比，居下的一条线表示第三产业的增长速度，根据深圳市国民经济和社会发展统计公报上的数据，在近十年中深圳市产业结构发生着巨大的改变，第一产业的增长速率在逐年减少，而第二、第三产业的增长速率却在逐年增加，因而大量的外来务工人员进入深圳市，导致了深圳市性别比、人口结构、分年龄的死亡模式发生了变化。
由图可以知道，第三产业的发展在一定程度上影响性别比（即模型中性别比与第三产业的关系式）。
4.3.3  预测深圳市全市医疗床位需求----卫生服务需求量法
在如何合理配置医疗卫生资源的问题上，有两点大家是有共识的：第一要以区域内居民的卫生服务需要和需求为基础；第二，规划是面向未来的，因此必须以某种方式较准确地反映未来的需要和需求。预测未来医疗服务的需要和需求的变化，不能忽略了一些重要的社会经济指标，多个社会经济因素对卫生服务需要和需求产生影响，这些因素包括人口(数量、年龄结构、性别等)、收入、价格、疾病谱及医疗保障制度等等。因此，在制订医疗资源配置规划时，必须对这些因素的未来状态做一个恰当的估计，并实实在在地纳入影响卫生服务需要和需求的计量中去。
同时在分析本市医疗资源历史发展及现有状况的基础上，考虑病人的流向、潜在需求、流动人口及人口增长等因素的影响，来进一步确定医疗资源的数量、结构、布局以及发展速度与规模。
4.3.3.1  卫生服务需求量法基本公式
医院床位需求量：
   医院床位需要量：

4.3.4  未来十年深圳市全市医疗床位需求预测结果分析
表4-2  2011-2020年深圳市全市医疗床位需求预测
单位全市床位需求（张）
2011 33949
2012 39063
2013 45795
2014 55200
2015 69759
2016 92470
2017 128464
2018 191571
2019 310318
2020 498966

图(4-7)深圳市未来十年全市医疗床位需求预测图
4.4  针对问题（3）
4.4.1  深圳市不同医疗机构的需求权重
通过综合患者偏好、医院规模（主要以床位数据为准）、医师比例等因素，求出了综合医院、专科医院以及街道医院三种医院床位需求权重。
表4-3  不同医疗机构床位需求权重
综合医院专科医院街道医院
消费者偏好 58% 33% 9%
医院规模 55% 40% 5%
医师比例 60% 33% 7%
权值 0.58 0.35 0.07
4.4.2  疾病别的床位需求预测
4.4.2.1  因病住院人数的计算
设，利用衰减因子法：，则对于深圳市未来的，有；其中，住院人数计算公式为：


图(4-8)住院人数与就诊人数的比值随年份的变化趋势图
4.4.2.2  因患高血压所需的床位数
针对深圳市老龄化的日趋严重，以及高血压在老年人口中的高发病率这一现状，随着经济的增长和社会福利制度的完善，老年人患病后的就诊需求会更加迫切，因此，我们预测仅仅因为高血压这一疾病深圳市全市就要增加的床位数总数，其计算公式：

4.4.2.3  因妊娠、分娩和产褥期小计所需的床位数
根据深圳市国民经济和社会发展统计公报的数据，第三产业的增长速率在逐年增加，因而人口性别比也在变化，同时根据深圳市卫生统计年鉴(数据来源于深圳市卫生与人口计划生育委员会网站)，在深圳市医院住院病人疾病顺位里，从2003年—2010年妊娠、分娩和产褥期小计的比重连续8年排在第一位，所以对妊娠、分娩和产褥期小计的床位需求预测是很有必要的，其计算公式为：

4.4.3 高血压在不同的医疗机构的床位需求的预测结果
表4-4  2011-2020年高血压在不同医疗机构的床位预测
单位高血压所需床位数（张）综合医院专科医院街道医院
2011 99.986 57.99 34.99 6.99
2012 100.054 58.03 35.02 7.00
2013 100.133 58.08 35.05 7.01
2014 100.212 58.12 35.07 7.015
2015 100.290 58.16 35.10 7.02
2016 100.369 58.21 35.13 7.025
2017 100.457 58.26 35.16 7.03
2018 100.546 58.31 35.19 7.038
2019 100.634 58.36 35.22 7.04
2020 100.722 58.41 35.25 7.05
4.4.4 妊娠、分娩和产褥期小计在不同的医疗机构的床位需求的预测结果
表4-5  2011-2020年分娩在不同医疗机构的床位预测
单位分娩所需床位数（张）综合医院专科医院街道医院
2011 7276 4220.1 2546.6 509.3
2012 7357 4266.9 2574.8 515.0
2013 7297 4232.1 2553.8 510.8
2014 7266 4214.1 2543.0 508.6
2015 7127 4133.7 2494.5 498.9
2016 6966 4040.4 2438.2 487.6
2017 7181 4165.1 2513.4 502.7
2018 7181 4165.1 2513.4 502.7
2019 7441 4315.8 2604.4 520.9
2020 10681 6194.9 3738.3 747.7
5. 模型检验
5.1  人口预测
现有深圳市1979—2010年的人口总数,、出生率和死亡率数据。用C语言编写的灰色GM(1,1)模型和改进的BP神经网络模型的程序，首先根据深圳市1979-2009年的数据，利用此组合模型进行预测，将所预测得到的2010年深圳市人口总数与已有数据进行比较，并验证组合模型的短期预测精度。
通过计算，可知GM(1,1)预测的数据与实际数值的绝对误差大于组合模型得到的仿真数据与实际数值的绝对误差，相对误差大于组合模型得到的仿真数据与实际数值的相对误差，均方差大于组合模型得到的仿真数据与实际数值的均方差，由此可见此组合模型在人口总数预测方面仿真数据的精确性较单独的GM(1,1)模型要高得多。
5.2  床位数量预测
5.2.1  卫生服务需求量法的精度和优势
预测未来医疗服务的需要和需求的变化，不能忽略了一些重要的社会经济指标。前面的分析己经证明，多个社会经济因素对卫生服务需要和需求产生影响，这些因素包括人口(数量、年龄结构、性别等)、收入、价格、疾病谱及医疗保障制度等等。因此，在制订医疗资源配置规划时，我们对这些因素的未来状态做一个恰当的估计，并实实在在地纳入影响卫生服务需要和需求的计量中去。同时在分析本市医疗资源历史发展及现有状况的基础上，考虑病人的流向、潜在需求、流动人口及人口增长等因素的影响，来进一步确定医疗资源的数量、结构、布局以及发展速度与规模，这样结合了实际情况的预测值比仅仅依靠历年床位数据所得出的预测值更可行，更贴合实际情况。
5.2.2  评价医院床位供需平衡状况
床位供需平衡评价采用供需比方法，即将当地实际医院床位供给数与医院床位需要（需求）数进行对比，求出比值再进行评价。
(1) 供需比值在范围内，为供需平衡；
(2) 供需比值在范围外，但在范围外，为基本平衡；
(3) 供需比值超出范围为供需不平衡，比值为短缺或不足，比值为过剩。
5.2.3  对医疗机构的划分
在这个问题上，我们依据经济学和消费心理学中的相关知识，综合患者偏好、医院规模、医师比例这三个指标求出不同医疗机构的需求权重，而不是根据医院科室的设置来划分，这样的处理更贴合居民看病消费的实际需求。
6. 模型的评价和推广
6.1  灰色GM(1,1) 模型评价及推广
6.1.1  灰色GM(1,1)模型优点
灰色GM(1,1)预测模型在计算过程中主要以矩阵为主，它与C语言的结合解决了它在计算中的问题。由C语言编制的灰色预测程序简单实用，容易操作，预测精度较高。
6.1.2  灰色GM(1,1)模型的缺点
该模型是指运用曲线拟合和灰色系统理论对我国人口发展进行预测的方法，因此它对历史数据有很强的依赖性，而且GM (1,1)的模型没有考虑各个因素之间的联系。因此，误差偏大，尤其是对中长期预测，例如对中国人口总数变化情况做长期预测时，误差偏大，脱离实际。下面我们来讨论GM(1,1)模型的适用范围。
GM(1,1)模型的白化微分方程：

其中为发展系数，
有命题当GM(1,1)发展系数时,GM(1,1)模型无意义。
通过分析，可得下述结论：
当时，GM(1,1)可用于中长期预测；
当时，GM(1,1)可用于短期预测，中长期预测慎用；
当时，GM(1,1)作短期预测应十分谨慎；
当时，应采用残差修正GM(1,1)；
当时，不宜采用GM(1,1)；
在对深圳市人口总数变化情况的预测过程中， =－0.040597，因此可做短期和中长期的预测。
如果要考虑到多因素的联系和影响，此时我们不妨建立GM( 1, n) 模型。 GM( 1, N) 模型能模拟系统发展的动态过程，不但吸收了传统的灰色模型的建立，而且建立了多中改进的灰色模型，提高了预测精度。具体参考[5]
6.2  BP神经网络模型的优化、评价及推广
6.2.1  学习因子h 的优化
采用变步长法根据输出误差大小自动调整学习因子，来减少迭代次数和加快收敛速度。
，为调整步长，在0～1之间取值
6.2.2  隐含层节点数的优化
  隐含层节点数的多少对网络性能的影响较大，当隐含层节点数太多时，会导致网络学习时间过长，甚至不能收敛；而当隐含层节点数过小时，网络的容错能力差。利用逐步回归分析法并进行参数的显著性检验来动态删除一些线形相关的隐节点，节点删除标准：当由该节点出发指向下一层节点的所有权值和阈值均落于死区（通常取±0.1、±0.05等区间）之中，则该节点可删除。最佳隐含层节点数可参考下面公式计算：

表示输入节点数；表示输出节点数；表示介于1～10的常数。
6.2.3  输入和输出神经元的确定
利用多元回归分析法对神经网络的输入参数进行处理，删除相关性强的输入参数，来减少输入节点数。
6.2.4  运用BP神经网络进行人口预测的可行性
一个好的人口预测模型首先应符合人口基本理论和数学建模要求，这是选择模型的关键。其次要保证模型数据可得，一致和可比性，在数据预测检验阶段应充分拟合原始数据，特别是有波动的数据，因为波动性数据往往蕴藏了系统重要信息。具有P个输入，Q个输出的BP网络可看作P维欧氏空间到Q维欧氏空间的一个非线性映射，含一个隐含层的BP网络可以逼近任何连续函数。BP网络的非线性映射与人口非线性特征是相近的，对于未知的人口动力学系统可通过神经网络来学习系统的表征量并进行预测，已有定理从数学上保证了网络用于时间序列预测的可行性。网络的输入向量和结构参数可由相空间重构技术予以确定，同时BP网络的动态学习能力对人口波动性数据具有较强的识别和拟合能力，克服了以往人口预测模型剔除奇异值或划分阶段处理的缺陷。因此从模型映射性质, 输入变量的选择和波动性数据拟合等方面表明BP神经网络模型用于人口预测是可行的。用神经网络对人口历史数据进行分析拟合，是人口预测的有效方法。与传统的人口预测方法相比，将基于BP神经网络的时间序列预测方法用于人口预测，避免了繁琐的常规建模过程。神经网络模型良好的适应和自学习能力，使预测系统计算简单、灵活，运用计算机强大的组合能力，可以更好地实现人口分类预测和管理，大大提高人口预测效率和预测精度。
6.2.5  BP神经网络的不足：
BP神经网络需要大量的样本数据用来训练和测试，当样本数量不够时，预测的误偏差很大。
针对这些问题，我们提出GM(1,1)和BP神经网络的组合模型，这样利用灰色系统理论具有所需要的样本数据少，原理简单，运算方便，短期预测精度高等优点，也发挥了神经网络并行计算，容错能力强，自适应能力强等优点。
6.3  供需比评价法对医院床位的预测的优缺点
6.3.1  供需评价法对医院床位的预测的优势
卫生服务的公益福利性质也决定了不可能使市场机制在卫生资源的配置中起基础性作用。因为在医疗卫生领域市场机制的局限性突出表现在这个机制不能解决医疗卫生服务分配不公的问题，不能解决人人享有基本医疗保健问题，这个局限性是市场固有的，是无法克服的。所以在医疗卫生领域，资源的配置只有在政府的宏观调控下，辅以市场调节，二者有机结合的资源配置方式。宏观调控的基本内容就是根据医疗需求，对医疗资源的总量、功能、结构和布局进行调控，合理配置医疗资源，保持总供给与总需求的基本平衡。
影响床位数的因素多而复杂，如人口数、人口年龄结构、经济水平、人均收入水平、居民健康状况、居民卫生服务需求、卫生服务利用、卫生服务能力等。其中，居民健康状况、卫生服务需求和利用是综合性指标，代表着医疗卫生服务的工作量，从而决定着床位的客观需要量。本文利用卫生服务调查中居民健康状况、居民卫生服务需要(需求)与利用以及卫生统计年鉴中卫生服务利用效率等数据，对深圳市全市床位的需要(需求)量进行测算，然后与实际床位供给数进行比较，求出供需比值，进行供需评价。此法具有以下的优点：
①按照深圳市居民的实际需要(需求)量来测算床位需要(需求)量，科学性和实用性较好；
②方法简单可行，可操作性强；
③结果是定量的而非定性的，在床位资源供需平衡状态评价方面，提供了一个客观的量化标准，为总量控制提供依据。
6.3.2  供需评价法对医院床位的预测的劣势
本法仅限于总量评价，而无法对各级各类医院床位的供需平衡状况进行客观的评价。
6.3.3  对医疗资源预测的评价和改进
该模型考虑的是常住人口，而根据深圳市的人口现状，有大量的外来务工人员在深圳市务工，因而可以对该模型做进一步的改进，即把深圳市人口的中非户籍人口远远高于户籍人口这一现状考虑进去，并考虑流动人口的经济水平和实际就诊需求，对医疗需求的预测会更完善。
7. 我们的一些建议
1.今后20年前后我国人口趋于年青化的发展态势，之后人口老年化的趋势逐渐显著，且成年群体的负担会越来越重，深圳市也不例外。鉴于此, 要加强发展基层的养老保险等社会保障事业，现在的年轻群体可以投入这样的组织来保障老年的生活，同时也减轻年轻一代的负担。
2.继续实施计划生育政策，促使人口健康、稳定发展，并利用教育，宣传等形式来。引导男女性别比的协调发展，如果按照上述的态势发展，势必会造成一系列严重的社会问题.
3.发挥深圳市当地的地区优势，合理调整甚至是的产业结构，创造就业机会，发展自己的经济，活跃市场，给上述的群体创造更多的就业和再就业机会。
结语：
此组合模型将灰色预测与神经网络相结合，充分利用了灰色预测建模所需信息少、方法简单的特点和神经网络具有较强的非线性映射能力的特性，采用多层前馈网络预测模型和灰色预测GM(1,1)等维新息模型的预测值进行组合预测，同时针对BP算法易陷入局部极小的缺陷和GA易产生早熟收敛的问题，提出了一种GASA混合策略的神经网络权值学习方法，并将这种混合策略应用于BP神经网络中进行学习，从而提高了预测的精度。此组合模型简单、实用、特别是在原始数据非常有限的情况下，可以得到相当精度的拟合数据，并能得到可靠的预测数据。
由于神经网络的功能之强大, 型式之多样，若能将其它网络形式同灰色预测相结合，则有可能进一步提高预测精度，并能在人口预测方面发挥更大的作用。
在参考宋健人口发展方程的基础上，我们提出了更切合深圳市人口特色的人口动力发展方程，在预测出深圳市未来十年的总人口的基础上，预测出了未来十年的人口结构，并基于此提出了生育模式和死亡模式。
在预测床位的数量这一问题上，我们结合了深圳市的产业结构发展特点、外来务工人员较多、非户籍人口远远大于户籍人口的特点、深圳市人均可支配收入、居民健康状况、卫生服务需求等种种因素，以医疗服务需求法来预测深圳市未来十年全市的床位需要量、需求量以及疾病别的床位需求量，得出的数据比单纯使用历年数据所预测出的结果要精确得多。
在医疗资源稀缺的情况下，对于医疗资源配置的方式主要包括一下三种：以需求为基础的模式，以供给为基础的模式，以供给和需求相结合为基础的模式。但是随着人类社会的发展和进步、疾病谱的改变、医学模式的转变和人均期望寿命的延长等等，如何更好地体现医疗服务的公平性和可及性以及更好地配置医疗资源成为困扰宏观调控的又一课题。
8. 参考文献
[1]Khotanzad A，Elaragal H and Lu Tsun_Liang. Combination of Artificial Neural-Network Forecasters for Prediction of Natural Gas Consumption[J]. IEEE Transactions on Neural Networks，2000,11(2):464-473.
[2]赖红松，祝国瑞，董品杰.基于灰色预测和神经网络的人口预测[J].经济地理，2004，24（2）：197-201；
[3]汪璇，谢德体，吕家恪，武伟.集成GASA混合学习策略的BP神经网络优化研究[J].(2007)，03(12)；
[4]邓聚龙.灰预测与灰决策[M].武汉：华中科技大学出版社，2002；
[5]深圳市卫生和人口计划生育委员会[EB\OL]. http://www.szhpfpc.gov.cn/view?f ... 032=menunews&nt ；
[6]深圳市统计局[EB\OL].http://www.sztj.com/main/xxgk/tjsj/tjgb/gmjjhshfzgb/
[7]政务公开站点——卫生部政府信息公开专题[EB\OL].http://www.moh.gov.cn/sofpro/cms/previewjspfile/zwg
[8]林安弟.吐鲁番地区医院床位医生配比标准研究[D].；
[9]刘丽娜，我国医院卫生院床位配置情况及预测研究[D].山东大学.2007。

附录:
附表1  2011-2020年深圳市各指标值预测
单位老龄化（老年人口比例）（%）人均可支配收入增长速率（%）第一产业增速（%）第三产业增速（%）
2011 0.0178 0.1270 -21.6090 11.08
2012 0.0183 0.1150 -16.5183 10.62
2013 0.019 0.0912 -19.2546 10.96
2014 0.0198 0.1012 -15.9456 11.14
2015 0.0206 0.1037 -18.3001 11.99
2016 0.0214 0.1046 -22.9454 13.08
2017 0.0222 0.0829 -20.4666 11.65
2018 0.0231 0.0983 -20.0182 11.65
2019 0.0240 0.0887 -11.2366 10.16
2020 0.0249 0.0943 -3.0914 9.37

(1) GM(1,1) 灰色系统C语言源程序清单：
#include<stdio.h>
#include<math.h>
#define N 100
void main()
{
int t0,tn,t,n,i,m;
float e[N],q[N],z,z11,z12,z21,z22,c1,c2,a0,p,u,yy,y,s1,s2,g0,g;
float y1[N-1],h1[N],h0[N],h4[N],r[N-1],z1[N-1],a[N],a1[N];
float z2,v,s11,s12,s21,s22,ee,l1,l2,c3,o;

z=0;
z11=0;
z12=0;
z21=0;
z22=0;
a0=0;
p=0;
u=0;
printf("===========================GM（1，1）===========================\n");
printf("请输入起始年份：");
scanf("%d",&t0);
printf("请输入终止年份：");
scanf("%d",&tn);
n=tn-t0+1;
for(i=0;i<n;i++)
{
printf("请输入%d个数据：",i+1);
scanf("%f",&a);
, ^: v/ M4 M, ?1 w) u }8 W5 G! q: M6 L
m=1;
$ z6 z0 k% q' L! x0 |9 F while(m>0)
4 ~4 {* r5 v# c4 O; q2 ` {
/ o7 j5 j- z, D# |1 U' J8 o. F$ I for(i=1;i<n;i++)
2 X0 u! u, f( `& ~ {
. j, c  y4 ~- e; x r[i-1]=a[i-1]/a;( I: I% W. `: h, o- W0 Q# i8 H3 I
}/ W$ {7 m) w( p' D* a; s
for(i=0;i<n-1;i++)
4 s" ?# T; p$ j: }! K- J/ b, L {1 g1 V! [/ u/ }( ~8 D9 ^
// printf("r[%d]:%f\n",i,r);
& c/ l% a0 t# L' X+ _ }
6 s* Q& @( M& F% q: F" J8 o+ W ) I4 n/ A# j/ W( J6 ~
m=0;6 K& _+ p2 ~. `8 c, a5 O8 J+ z
v=-2/(float)(n+1);/ R0 S$ z0 `, F" ~9 m. @; T
l1=pow(2.7183,v);
- j! @8 |. d2 r l2=pow(2.7183,-v);
, u; J8 @. W; N  V" V0 u; s; Z% E- L  }4 H2 [0 M1 K
for(i=0;i<n-1;i++)/ Y5 E) b  g6 f+ N% x+ \
{ " {: D  e4 v! ~. w7 t6 q, e
if(r>l2||r<l1)( W' x4 `# t5 j3 k5 d& Q3 q6 w
{; @  r% c" b# P- f" p
printf("r[%d]:%f\n",i,r);
; L% Y/ V+ u$ R9 v printf("11");
& m2 g4 \, b+ z6 F% n m=1;; A$ A, ]! X' p, p+ M3 h
printf("需要调整，请输入调整值：");
. r6 D0 A" T2 f& }  \ scanf("%f",&c3);; a4 B8 r/ U) y" h
for(i=0;i<n;i++)' T! l8 u. Q+ D$ o
{2 z2 w! n1 h* X% H* v1 B3 p. q9 _+ F& F" |
a=a+c3;4 n, }/ A* b7 L; u* }$ b8 Z. d: |
}2 v; w! y4 ^6 j/ M
break;6 k+ z5 o8 I7 {& E0 P+ G& g
) E# k: H! I2 \( H
}0 ]7 E6 y( s+ ^! |$ w0 r: q
}
5 N9 c' n) ^8 D+ o }/ }. E1 e3 t- Y1 s0 f1 S
a1[0]=a[0];  z6 s2 U0 r/ _* P- q
for(i=1;i<n;i++)+ G5 l$ J, _; ^; Q, {' _' k
{
1 D' R) V% W+ V5 n! G0 x2 D a1=a1[i-1]+a;* F3 t& G$ E* d3 L5 I" W& T$ y
// printf("a1 %f \n",a1); //累和
& N* c  w& G- T8 }" Y6 e, h" z/ U }
, v7 T7 S- C$ X% e* j" r$ ] % o% Y9 y2 u7 f
for(i=1;i<n;i++)
, `8 X. ?% W5 i {
0 v8 n# V0 |8 L0 @1 A y1[i-1]=a;
. x- b# R5 s+ |1 o }5 U' ?: ~( L; q" A
for(i=0;i<n-1;i++)
6 l5 I. q# G2 J! | {' ~3 a/ w7 k+ G/ D8 b/ ^
// printf("y %f \n",y1);//缺原始数据第一个元素: Z0 _# d; P( }, O. e7 {. K% x( K
}
- S8 o- D4 y# n& w1 x- j for(i=0;i<n-1;i++)6 A$ P* M( l5 W7 {! c7 F1 K7 {
{5 f/ `2 J3 V, h% P
z1=(a1+a1[i+1])*-1/2;* R) J3 s6 \2 z: X7 a3 x1 g
}
7 E3 n8 [+ y. Y3 U9 o3 T: L for(i=0;i<n-1;i++); l, i/ M0 Y, o% H6 y; Q: N7 `
{
' J9 t! E+ Z* m, u  b5 P) @// printf("z:%f \n",z1);// 求1/2（相邻两数和）
( C. m$ W* h% V$ N' C }8 @4 v& ~& y- m8 o: ^0 A2 b
for(i=0;i<n-1;i++)
5 Q( U6 [1 V1 \9 n {
/ @6 }6 I7 L! D' ?" a& f" D6 F z11+=z1*z1;+ P' c4 w" z+ F) `0 y) w
z12+=z1;' ^, u2 N; T% P/ v) ]9 h1 `
}% P- C- G( Y; M" W
z22=n-1;' K' b6 a( J0 D' E
z21=z12;& {3 _# k/ W9 `9 G* {
// printf("z矩阵：\n");
, N8 v: w, J. l6 I2 C// printf("%f %f\n",z11,z12);7 Y( Z2 U- N- c' |* _  P
// printf("%f %f",z21,z22);
8 I" Z: h' l5 a8 H6 s+ ]6 F) s z2=z11*z22-z12*z21;
1 _- M. t9 f8 q) \8 `+ ^// printf("行列式：%f",z2);6 Y) R( c  t. A, O) r3 I; `! a- ?
s11=z22/z2;
5 e( k7 }" h1 G( i- X7 [ s12=-z12/z2;5 m. Z5 Q1 y3 L; |. Q: y
s22=z11/z2;8 O8 \% g5 l) R$ y
s21=-z21/z2;
0 Y  v2 F7 D1 ?2 D) P4 n// printf("\n");1 M& h$ x' S; T2 N% z, ?: o
// printf("%f  %f\n",s11,s12);# g: y- @- C) M6 Z$ d& S  H
// printf("%f  %f",s21,s22);4 z4 k( g+ L, r; q2 Z
c1=0;( c" M/ E4 ~7 R1 J
c2=0;' l: F2 W0 m! }: J
for(i=0;i<n-1;i++)9 _! R: n% h: g  o& w. x6 k, h
{
& x) }  ~6 x/ O4 c5 _; ^) \ c2=c2+y1;9 ?" z* _' E' c( H+ j* X3 Q. q- n
c1=c1+z1*y1;
' u" b6 g+ O: ?4 Q }
  P: x7 L) L1 H. Q+ c// printf("c1:%f c2:%f",c1,c2);
2 m9 y3 w$ }" J  F, v6 m8 T* N+ |) V& F
a0=s11*c1+s12*c2;" R8 |$ O5 T, c  _" ^; b
u=s21*c1+s22*c2;( ]3 L* q1 |3 O5 W3 k) ^
// printf("a=%f u=%f\n",a0,u);
, R- M8 C1 Z+ j: q2 S; I$ m for(i=0;i<n;i++)
5 b& y0 @4 y! h9 B/ |, k {
, ~; Q! x+ G  x) I% b" d3 @0 m6 L h1=(a[1]-u/a0)*exp(-a0*i)+u/a0;
5 l3 O8 ^+ M' X+ o }
4 J5 o- n6 j5 j+ b& k' e' A+ ?! t h1[0]=a[0];
: N5 X) x& Z1 |3 Z h0[0]=a[0];
" |& n) Z6 |4 e6 k for(i=1;i<n;i++)
* E: _' V' o3 Z8 e {6 y$ B2 l( ^0 K! i$ u
h0=h1-h1[i-1];
' p3 r' x0 [8 B }0 Y3 Q9 S# v/ ^8 ^, s+ w
for(i=0;i<n;i++)' G- u, o7 {0 @
{
) _% T7 l" M( e" x: Y' t! K// printf("预测：%f\n",h0); //预测值！* L5 c- ~: f  r6 H% u1 N
}' H! \1 T* L1 |1 _. s7 d
for(i=0;i<n;i++)
- U' @+ x+ m2 w  k {
" r2 N, [& C+ o( L$ ]3 L; n" @* i e=a-h0;
8 q, G/ f+ Z5 u  j& q! d }: `$ T8 t* O. [2 @; m3 H( d
for(i=0;i<n;i++)
3 I( B" t; X% F6 s3 \; v {+ T+ q( w3 e; r- {
q=e/a;& r& W$ ~4 Y8 ^% e( }" i8 _5 }
}
, @. E* X7 p8 `: k) V# W& w5 c for(i=0;i<n;i++)9 G: l, S% u3 ~
{
+ r  R, s( G) j* e// printf("e(%d)=%f\t\tq(%d)=%f\n",i,e,i,q); //--------e为残差！ q相对误差！、
' Q) u$ W% b! j1 [* y }
& g' C% l9 s2 S! k8 s( I6 X, \ yy=0;
1 ^% P1 X! ?( x" Q4 P+ @- ?+ { for(i=0;i<n;i++)4 g! F) C' [4 s1 L
{' H7 Y' f1 ?" t5 M9 \  K2 }
yy+=a;( T/ X8 ~' s! T: q/ U3 j9 f6 ~
}
! _6 b5 q# q: @6 G3 S; Z yy=yy/n;+ V1 c" @( J7 L4 F
y=0;
" c6 p9 x& p+ N, q4 e8 E for(i=0;i<n;i++)
/ R% c2 B5 C  g {
7 i3 i/ F4 n4 e9 G/ C6 u* D3 U, m7 Q* { y+=e;
! O. ^4 S5 d2 `. r6 T }, I% E* v4 s4 D) x  T4 e
y=y/n;
: U6 L& s7 P9 E( D8 @0 F$ B// printf("平均数：y：%f e:%f\n",yy,y);//原始数据的平均数; m- e: |2 X$ e. j7 v9 |( _" y+ V
ee=0;
6 y# c$ w; {- w- B6 ]& y for(i=0;i<n;i++)3 d/ ~- B8 A" `  p
{
# S# l: L7 g& n2 w ee+=(a-yy)*(a-yy);0 u7 J+ o8 |& u" B1 F1 U7 B
}
4 l# b' {  }9 E: ^2 S s1=ee/(float)n;/ Q* C, {# {( b2 Q7 S0 Y
// printf("%f ,%d\n",s1,n);8 f0 O+ c' e% T1 ~8 ?# O
s1=sqrt(s1);
* S+ }6 M( Y/ a0 A2 L: G  S$ U/ L. j// printf("s1:%f\n",s1);//---------------------s10 m3 _/ G* ]: f7 ~3 y( o
ee=0;+ F+ F4 g8 u, |1 `6 n% \/ ^
for(i=0;i<n;i++)
9 s& E0 U. B1 x3 |/ S {* `8 K0 Z1 E- |
ee+=(e-y)*(e-y);
& J# O) m* {1 L5 A, U% n }
0 y! H6 ^4 l! ], W, s s2=ee/(float)n;
* w  u, o" Z- l' c% C3 r# B// printf("%f ,%d\n",s2,n);
. o2 ]5 }- C. X) V2 l  D s2=sqrt(s2);6 ?. |; v# u5 p  E% J& \8 L
// printf("s2:%f\n",s2);//---------------------s2
2 S  W' k2 c7 g1 Z- P ee=s2/s1;
7 u! @* z6 Q* `' i8 t, n4 ~// printf("均方差比值：%f",ee);//-------------ee均方差比值
+ P! k4 W6 Z% r6 m  k o=0;
; `4 q6 F$ c' r! I2 M- j: M for(i=0;i<n;i++)% I, T: j" n  e* {" R
{! p) l1 q# K7 [
if((y-0.6745*s1)<e<(y+0.6745*s1))/ f5 Z' @4 o5 S2 l# l: s
{4 x- f* W& e0 |7 k" z# ~+ R
o++;
3 p: p( B* X% R) D }
7 G  H+ I+ t& n% M. n# |9 l3 F: o, f: S8 Z- A4 m% i. p3 Y3 o
}" Q7 u7 w* y( u) h
o=o/(float)n;7 J2 x( X% k! R) k8 a
// printf("小误差概率为：%lf\n",o);//----------o 小误差概率- h& y% [+ O6 F
printf("请输入预测终止年份：");
# B) p* |+ ]/ W$ q$ S& B scanf("%d",&t);
& b* V, ?9 h- j. }; a4 A9 s t=t-tn;! H% g$ [& }! T* V6 V5 w5 A
g0=(a[1]-u/a0)*exp(-a0*n)+u/a0;: ~6 r- d3 A2 j# z! I5 H, a4 K
h4[0]=g0-a1[n-1];3 w5 N7 J/ b7 G
for(i=0;i<t-1;i++)" n6 b; m& C8 a9 b' _
{
: z/ V5 n& O! E9 Y1 E/ R g=(a[1]-u/a0)*exp(-a0*(n+i+1))+u/a0;* \' G% r$ L# F$ H& M9 Y! I
h4[i+1]=g-g0;9 y* g+ K3 A% A8 V5 b$ i
g0=g;
, ]  i8 B* \; p( c% ^ }//  预测1 a7 P8 C! d$ R- \+ U  L
printf("===============================结果==============================\n");
/ u; r' \. F( d/ G! D* {& R, Y printf("序号：");
7 ^8 A& T. y% g% R2 _ for(i=0;i<n+t;i++)
: \. s! u% c2 B$ I  V: @8 w {
& F  G5 B9 u# X: p$ e6 m2 N printf("%8d",i+1);: f7 U5 m& J2 [
}# e# s5 \" |7 z7 B) m  U' b% J
printf("\n");
/ y* ^- Q3 E: Q3 m4 I8 E printf("年份：");* b" Q- }: M+ P3 B5 y
for(i=0;i<n+t;i++)
+ B  r7 i$ F4 V4 `+ Q5 t {
7 Z8 U0 [( \4 U printf("%8d",t0+i);
- }: P1 ]/ b) P" {* x! r7 e/ N }
7 X6 N1 o" A, W& K' C; E. L3 h! c1 j printf("\n");
& n' a% M  S( a0 J printf("原始：");' ?0 c4 ^! _' J6 r( ], y5 J  M8 L
for(i=0;i<n;i++)
1 C; m4 h( N0 ^+ y" D; H {
. r# d; R7 n1 V printf("%8.4f",a);# D" w3 E& c4 \. n) q1 N: ^
}/ D% ?8 ?( c/ p6 p+ u/ h/ A
printf("\n");
7 f5 {% A7 E6 E& j2 K: X printf("预测：");
* e$ T% I8 B) m" i for(i=0;i<n;i++)
7 `4 ~& X+ U! [# ^1 l( [ {& b3 Z: C0 }- k  l
printf("%8.4f",h0);; i4 w: x. n. g8 }7 P& _
}
1 _: B9 z) @2 E/ o$ j& k for(i=0;i<t;i++)& B6 o) h2 |# h" [( z
{
8 m4 W! ^" t3 c- F% C$ D printf("%8.4f",h4);
' F( P& Z6 T1 @. E }
0 j4 D+ T6 s. V7 c7 c9 V printf("\n");6 U5 K5 N9 n2 y+ n( h
printf("残差：");7 i: O: d* w# _1 H# Z# l  s
for(i=0;i<n;i++)
9 U2 _: X* ^% p* I7 _ {
5 p6 w0 i+ p$ A printf("%8.4f",e);$ \# z6 u; g/ u6 Q
}3 E3 y0 ^) i% v* [( h4 |  K' Z9 y
printf("\n");" R* I  M% y! \# Y) e
printf("均方差比值：");1 a  X% ]2 g" b7 r6 V
printf("\t\t%f",ee);: R! |7 `$ w" P
printf("\n");
2 o" ?( f( G- a2 Y printf("小概率：");
) Q7 F3 g2 }, A1 _0 o7 ]7 H printf("\t\t%f",o);/ Q: O8 u1 U! G* a6 N
printf("\n");
1 I' t9 E# c2 V) o$ U5 ? printf("==============================精度评价===========================");. f$ W2 ~: w7 l& x* I

  |0 i# E  i/ N2 M% D printf("\n");
  J8 R! `4 Y7 A; D7 y0 B) q9 i* L, A if(ee<0.35&&o>0.95)
* \2 ?2 o9 h  G4 e% t printf("\t模型预测精度好！");
: m/ Q, y: i5 c: d else if(ee<0.5&&o>0.8)+ {* t  r: m) G- ]0 x& K8 Y8 t7 m3 p
printf("\t模型预测精度一般！");0 ^; d! t+ {! _7 O' e
else if(ee<0.65&&o>0.7)
0 O+ U7 v1 f8 n+ N printf("\t模型预测精度勉强合格！");
% K2 A' J# A" ?/ y" y* W else printf("\t模型预测精度不好");0 O* I  r# ^7 {
printf("\n");4 z$ G0 q4 `  D  `4 H
printf("==============================预测结束===========================");
6 |0 I3 j1 W$ _# I- ?+ N printf("\n");5 S& q/ A% ^) ?/ j: b. ?/ G9 y
printf("\n");0 i5 O7 d- v) S) u) L& W
printf("\n");0 _1 A; f7 x, Z: p. h; T% k
printf("\n");; u& P/ k3 l6 X, }

3 l' V0 i. V0 q2 b
, m  [$ P; m4 P! O+ o7 {; c scanf("%f",&o);
+ p/ d: `; k/ S$ k" w- T+ f0 s1 z! s$ ?) e- T( f
}
* \  P& A$ w( v(2)基于GASA混合学习策略的BP神经网络的C语言源程序清单：8 Y3 g( T- N4 \( W% m0 U
#include<stdio.h>
  |8 r- H0 O6 [% b6 E#include<stdlib.h>. h+ o9 G) K; Q! V) T3 N( K* @
#include<math.h>- q5 K" T5 I, J7 ]4 d3 L
#define M 10
3 q7 U. ^1 _. L8 N& K6 j2 R#define N 16
+ D# U- s3 b# Q#define m 10
" C4 p. {7 G# S0 W5 H/ d# ^#define jingdu 2000
. U* e3 X5 v2 X#define jingdu1 1000
6 M1 P7 ]( I7 L* f9 l#define jingdu2 0.1+ q8 [3 x+ A* Q3 s, e- H
  V: s+ K. q1 {7 i+ ~" O! T5 ]
const double e=2.71828;
0 ~- N8 v3 g" c7 k$ `
1 m/ W: I+ |1 qfloat in[N];/ F% k' F0 f* a% [
float out[N];8 H$ E) ^, @2 {5 S( C7 M0 X3 W  [
float get[m];
7 j6 M2 Z( t8 ~& t3 n. E; G1 M! x- g+ S

/ b8 i( ?7 U  O$ a+ l$ r* e" Mfloat one[8][M][N];
; A9 @& y; P( ~float two[8][N][m];8 v& k) I  h: B  |3 c+ T- E
typedef struct0 y2 |6 f. k% t
{1 I6 c1 U1 ]3 |' G" f) _; @  b3 H) Q
float yi[M][N];
) h+ a; R) M) x5 l float er[N][m];
  Z* N  @* D5 l7 K6 e2 E. } float fi;
) e5 [+ z/ e6 N8 R float Fi;
" y  b! h0 C( Q3 ?; u}dui;2 U; J* v8 C. c/ b3 J
dui du[64];5 Z% i; H3 l  N; f2 J; o
dui du1[64];
2 d9 |$ \! M" J" F3 e( wdui du2;: k+ |, |4 |' l  e5 J7 J
dui rao;* r5 h9 D0 A$ m2 |! U9 N" I8 @
dui gd;* x! ]+ a) `- y: C

# v* m' ~% G& i, [& D; }float T=100;; h' _: Q$ J; ?( V4 f* j$ {
; s2 ]( k* V3 `
typedef struct( {7 U) u/ f0 T" D4 |
{
$ c7 Z- b7 }& O" X! P! I4 R- l float year;8 f; ?3 V& }1 f; u7 S. S( \( v6 z
float num;
  B- {" ]6 H9 H) @' Y}Data;) N" L* B. H/ k2 s" ?, V5 ~
Data data[10];
: I# ^3 s1 f! n; T+ \* u2 g
4 M4 T  S8 G2 W( ~" @7 A( o4 l+ e$ A( `6 q1 s5 Q3 t4 _- ?
void chushihua()  L) U$ k$ O( W3 c" Q# }" Q
{) p" _' f+ p1 i3 x4 X
for(int i=0;i<8;i++)
; [/ T+ U) j" D) {% U, L1 @ {! y1 c4 ?' E. I, o; \
for(int j=0;j<M;j++)9 ^# P1 K% I* d) V0 v: Y7 u
{
4 z' ]$ e% s) m2 m for(int k=0;k<N;k++)- e( V" K' L! y0 b9 M
{
0 \$ n0 ~3 a2 [6 K( o int h=rand()%700;9 b# K1 R* @2 x# H) X
int w=rand()%10000;
8 e$ Y' O* L% X1 Y float q;
" P' z  X. ]# K q=(float)h/10000;
# U' C9 p. j$ V% ~/ N q=q+(float)w/10000/10000;. m2 o! T3 a# ~% B. y
one[j][k]=0.0375-q;//+ K" V: _! Q) e" V% L6 c& X9 q7 X
h=0;
- D/ U/ H5 m& M; m }
- A/ l' x4 G# e3 m5 ?) O }( d4 c7 j! L; l* m/ `" {" b
}
8 J: N$ {( w, O3 L4 S8 i
3 Q$ g* b: h. m3 W for(i=0;i<8;i++)
1 d6 g/ f9 Y0 \2 \ {
9 d) ^" A5 u$ I9 D3 z, e for(int j=0;j<N;j++)1 W6 B; W8 s; ]4 O: N: m/ M5 J
{4 `9 @4 s* d/ n! ?
for(int k=0;k<m;k++)# ?9 Y& Q* `) ?: \7 O
{
9 z5 x& R2 G4 _# F int h=rand()%24;  S8 `- q) a$ N' u
int w=rand()%1000;9 D! k! T$ R' ?: ]: H
float q;1 w7 x+ W8 t* o8 ?( [% G; d+ l
q=(float)h/1000+(float)w/1000/1000;8 R: V3 e5 U  s" D$ W. \. A
two[j][k]=q;! a0 n1 |( G! Q  T* |% s
h=0; //+ y' e9 f, v4 I  R0 R% b* l  S( |$ u
}
3 T& y- N; \% S. } }% r- a6 a5 _, Y; ?
}
6 f" ^+ H5 P$ c- K% @) ~. L for(i=0;i<64;i++)
' K% x# E. w$ K1 t3 W/ m { int h=i/8;
: r# [4 a' E% x. g for(int j=0;j<M;j++)8 ^" f# c* I) G4 [" ^
{" o2 ]; E1 ?0 c: |5 G/ j5 b0 e
for(int k=0;k<N;k++)5 w& q8 N- B6 a9 Z" E% H# d
{4 _1 M) H7 _$ s
du.yi[j][k]=one[h][j][k];4 Y* u- }! Z# A
}
6 Q* d% b3 L1 Y0 S0 N* j, `( s }3 Z" {; q5 x, T6 M1 }* w
int z=i%8;: d4 M) n+ |( e9 g0 Z
for(j=0;j<N;j++)
( A) _3 i/ R1 C: {1 C* X& u {
" d5 e6 ^$ C: P  v+ F for(int k=0;k<m;k++)4 K" t' O: I( ~: z2 k3 X  ~# H
{& U; ~4 I3 |. z( e' m
du.er[j][k]=two[z][j][k];
) P( x. f$ B! a% `; G# l% | }
9 i; m4 [' N5 l3 s/ [ }
3 O) @2 a9 D# J8 ~0 x; o* k8 K5 ~. p4 k5 u2 n! ]$ Z& n' z
}- B; Y- k5 k0 z3 ]" M4 S

3 S' F6 J5 U9 \0 q9 y+ U& l: C/* for(i=0;i<10;i++)+ _0 Z& R7 c- {' F4 z  E" \8 ^; Q/ [
{
% c6 S! t! v# n data.year=i;' E  C4 S! K3 d
printf("请输入第%d个数据：",i+1);
3 z- f; A+ E) O+ b& e: E scanf("%f",&data.num);
* _2 c, M+ v) m& Z! E8 ? }*/
/ T6 A/ M( F, n data[0].year=1;
& g3 R2 m0 l: V data[0].num=0.1106;
3 A& t" C4 Z! X5 L2 p data[1].year=2;8 Y1 g( D! n+ H% y* e/ n4 W4 \9 L
data[1].num=0.1879;
8 g9 }' g, s& d* I" X. Z# |, S' q data[2].year=3;+ E( I, Y+ `( y2 U* Q
data[2].num=0.3153;' C; U# ~/ s3 S6 x4 f
data[3].year=4;  P" @4 ]0 C' i" e9 N7 o7 _
data[3].num=0.4591;6 b5 J; u: @' J4 p5 ?# g
data[4].year=5;
& a7 Q6 X( R$ b' B0 E data[4].num=0.5274;
% F( [. F0 u4 v1 Z, J data[5].year=6;
2 e9 }5 O2 D/ N, ] data[5].num=0.5907;  r- a- s) f$ \. V! f. P- J
data[6].year=7;
9 D* O2 K7 V1 Z- ^$ r. _ data[6].num=0.6748;" U% b4 h- Q( n5 [' q" n
data[7].year=8;6 ~9 I- m% p8 J% ^) h8 f( T* O$ d
data[7].num=0.7733;4 m& E9 B* `3 {0 G7 r8 ]
data[8].year=9;
0 F0 P& E/ ?  v1 e data[8].num=0.8713;  n" _, U! ~$ O! o9 U, N
data[9].year=10;
1 J% X1 l4 Y# A: f2 g4 e data[9].num=1.0000;. J* f: t" S2 F7 b: }1 ^" V9 Z
2 ?& g! x+ a, v( C
}8 _5 E! I. u& C+ [6 I

. `! {! Z8 A1 h1 G% z$ _- w- T
0 j) q$ w; R4 x% F, o, U; H: u
* [% \; i2 ~5 S
- d2 K0 g1 b$ V' ?

: ~8 ]+ a2 [1 V$ c4 q% J6 n, T" g/ b- E
void yuce(dui h)+ H  i6 `: d9 v. d" S; W. ?3 k" ~
{! u: }# C2 G, k0 N, p
float z,g;
* b9 L; V3 Q" w; |4 P% f2 ~0 U% f for(int i=0;i<N;i++)
5 u; S0 h5 o- @" |# M0 y( z0 @0 G9 V { g=0;4 h* E5 V5 n% N+ S% P
for(int j=0;j<M;j++)
8 A( d# s- b/ E3 U {
$ {. o% C5 _& E- A5 H g=g+(j+1)*h.yi[j];! j( S9 [, Y# k$ }' ]6 W: }
}: [  ~4 L% h" I0 D: ?- A2 c$ p
6 H7 W. _) f$ `; `7 d( @" Z
in=g;' u  y* \. `' _% z: k* N+ m
}
$ l0 T# G. z9 ^: d; J- v2 q% I6 M+ f" F2 o+ [# V5 q
for(i=0;i<N;i++)
) a& z7 [) V0 |: v# m+ I/ v/ W {
: D) ?! E0 Y$ f$ I# b z=in;
) p% m) r7 o, J z=pow(e,z);1 S2 c6 K/ _0 }
z=1/z;
& E5 k7 ~$ O+ ~$ X z=1+z;$ S7 e" J1 \1 W
z=1/z;
/ S% y4 P- z9 r* P8 T out=z;3 f# N/ E0 H3 v& m; n
}
. W1 Z( w/ D- }/ E* u* F
3 u7 W* a% \8 y, J" f for(i=0;i<m;i++)
- U" S( i- G  i) n { g=0;+ _7 n6 B& N1 R1 r/ ]
for(int j=0;j<N;j++)
" H9 L" e, b) V5 \6 e8 I7 q {
' e, D  ?5 l7 h: b  j g=g+out[j]*h.er[j];6 q7 n) c9 F# K$ @/ a7 e

8 H$ B' i' h( v5 e1 Z- \ }0 d$ l2 U) m# [/ q
get=g*10;
( h1 B+ Z# Y+ P* T. F' K0 n }
* I/ H! K( l) Z( w3 j5 f9 a+ |+ E; l for(i=0;i<m;i++)
: x( d: m+ ]( ~+ t4 f  k {% Q2 T8 h% y. F# h5 j$ ^3 X
printf("get:%f\n",get);3 o9 y1 s$ y/ u; T. A
}; h" O; t! F+ `2 G6 \3 ~
scanf("%f",&z);# ^; z; E4 V7 |% K* K5 L
}& R7 s; D4 f* Y; {8 S- l1 O
float shiyingdu(dui h)- x! l4 q7 [) j: ^0 J
{
! y" v9 j) X: [ float Fi=0;) O1 Z0 C: E" [
float E=0;: Y* w! w/ b9 n8 T
float z=0;6 q, {& Q4 U8 Q0 ]5 i
float g=0;
$ U/ L- O% D4 @# k2 j) O' d float in[N];
0 |* g! r  n: |0 m" m( P, t8 R float get[m];: }2 o& y* n8 r4 j4 S# c* o
float out[N];$ z. \6 I- y! A5 G: v7 {6 U) l2 G
for(int i=0;i<N;i++)+ e- I8 v) g3 S. m
{
# ]. `3 L0 T; B0 R, p for(int j=0;j<M;j++)
) @1 G3 U1 f2 N( J) R {
8 i& ^/ o! e+ h) N g=g+data[j].year*h.yi[j];' Z; V! W0 X1 z! _% V0 s
}- g* Z/ V+ Q- {, @
in=g;" W) X- v) S: ?, e# }6 e
5 B7 N6 _' h9 P% u- e- c
}" s" H* D1 u4 O7 |8 y

/ H' N6 ~3 f; D% a
: o% }) {$ Z! O! ?3 h for(i=0;i<N;i++)  V6 R( @, J; o3 M8 p
{9 E. A2 r! ]/ ?9 e/ g# o' [
z=in;
( `! k* r% c! t! p! {. `2 e z=pow(e,z);* \+ R8 c; e0 U* U3 H; j
z=1/z;
5 E- ?# Z* G+ x- F* J* }4 d! } z=1+z;4 I( w# r6 `4 e8 u) C1 `7 R
z=1/z;# M0 V- {$ o( G$ K: ]
out=z;; w, I9 {& \. @# {
}4 Y4 v* Q5 _3 a
7 i; Z, X. X0 I$ |( p
for(i=0;i<m;i++)- ~9 ?4 u- J- j. R$ Z# c8 c( k9 Z
{ g=0;
) X& \. _: P) \; p for(int j=0;j<N;j++)& ?2 P6 L+ s; V6 \) L9 t
{
) ?! `/ B) _$ K8 s: C g=g+out[j]*h.er[j];' H5 p- K* Y; G  B% \
}* w0 W* G2 b: P$ F- v( F9 b
get=g;5 }0 G: f  m) H& h
}
" g, ^* l: {9 \( w: U8 |4 }4 `0 A2 p0 ~( e- B3 T5 n; U! C
+ l& k) y# q9 B
E=0;& S# ?+ C4 q0 g$ p5 t# Z* |. A, ?
for(i=0;i<m;i++)& p: \, \  z( b7 ~/ S& j
{
6 p8 Z% _7 t# b5 _- B" l E=E+(get-data.num)*(get-data.num)/2;9 L2 v; D: k, S2 L( H6 m9 }
}
/ L! [5 `' ?- W5 L # T% y' r% M: `) _( c, b5 }
if(E<0.01)
6 k5 N2 F0 e& Q% N) s8 @ {2 U" |1 m. D- G) K, e, ]
yuce(h);* B8 }  L% E% A; K" h; g, C
exit(0);
3 c- L0 d( G7 ^5 ?' T% U2 a }5 C- t# N% Q" G
return(E);! Z5 M; L0 W  {) C% h9 }
}0 Y' A6 B9 ^$ T# X" e, O

5 z& g% p# \& ~" {6 }5 E# F0 h% ?) i$ F' e9 d( y
void gailv()
) o$ |6 d/ ?) I{9 ^3 g, m/ r% E! o% r
4 X# y1 Q# h3 w; t
float l=0;/ q+ `* ^4 v, f/ O' i
float max=0;- p  J2 v" Z; K7 z' J' R& b
float min=1000;
" y' ?" D$ P: ^0 V; c3 ^8 f; E2 T, a, h
for(int i=0;i<64;i++)4 v2 Z6 l5 k( i9 X. W) b
{
" @! \1 v# q3 Z% T: t5 v8 q8 F if(du.Fi>max)+ p1 }' u8 @1 U% s' _0 y# q, d" W
max=du.Fi;, X' T0 B4 `+ x6 Y& H
}7 F. S, O  X4 P; t% \: i# r
for(i=0;i<64;i++)
3 e# y/ T( U) K& \ {% J! s4 i# D( N6 v7 t/ ^& b
if(du.Fi<min) 5 }) \$ U$ ]1 U- P0 @4 {
min=du.Fi;+ I4 {# @5 b# j. ^( r5 c! K, t2 y
}
' x+ K0 M% q2 p4 ^: t3 U for(i=0;i<64;i++)
: B! [; v+ ^3 p) N/ N {0 P% x: W/ H; I( v/ M* j5 }" F  N
du.Fi=(du.Fi-min+0.000001)/(max-min);
1 h$ f7 Y. q- ]6 {2 d }
+ D9 e% G! a0 s' V1 y# [6 a' p( o; X
for( i=0;i<64;i++)
  W( D" T" k8 T' L {* W/ J9 h+ ?# U2 _6 ]! x+ x
7 o1 C, [) D! z$ U$ k# x- V0 G/ @
l=l+du.Fi;
! K9 A! L8 I! W7 h }5 i+ i. V& `7 Z! r
for(i=0;i<64;i++)
( w% v1 n5 s6 H( K: g- X( v {
/ @/ I; ]( {9 @/ v, s! H' { du.fi=du.Fi/l;2 T) R: \) ]# ]9 ?7 G0 u
}
0 d) B$ D4 A1 o+ n* v0 v6 q3 ^) I# A( N% k2 }
}% V; @* H0 X/ b  O0 M/ E8 q

$ R( g, g! v1 f. ~, t0 I9 G+ }& ~( s$ c6 a) z) M+ L& @# t( w1 U
. B4 R- K! j- v7 X9 F
% w+ Y& Z/ p8 Z, g/ B9 J

# G; {6 `) E, U  I, w
; i) W2 f! Y0 q2 M' I$ E; B6 z$ x" t# G  \7 E
. F1 {6 y; X% g9 ]
0 y- t" F8 x% Q& r3 T

( b. N; l, S7 Y' Z1 Q0 _$ Hvoid lunpan()0 ?7 x+ C+ R& V% n% |* ]
{int t=0;! m5 ]- m7 w8 e
for(int k=0;k<64;k++)" d/ _+ T2 B+ V
{* ~5 i# @) @8 j- z9 c
float p=0;
4 w8 d" P" ~: m6 Z1 ^5 C & b9 K3 G) T3 E: h+ s$ J' `
int h=rand()%1000;
8 ^0 U1 F$ `  M' @ int w=rand()%1000;# @( ]# V: C3 k* c3 V% ^
float q;! I" `5 Z) B% `# h3 U4 I4 X9 M
q=(float)h/1000;: m* R" U$ R$ ~8 J& ?! O
q=q+(float)w/1000/1000;
* E/ x, Q. {# Y1 s
+ o: |5 ^) f: y4 @" L for(int i=0;i<64;i++)
. q0 ?) F' K  ]# T5 u" N {0 D$ ?6 E& x. t0 J  \( P6 {
p=p+du.fi;" s; F* Q; @) I: b) R
if(q<p)! E0 W; V7 v- ~& N
{
* i5 O+ t/ j% }0 d1 b  ^ for(int w=0;w<M;w++)% F3 W1 F- Y* F8 u
{# B& y+ j. L# @) Q( ~5 K
for(int r=0;r<N;r++): r2 y6 L8 _& T' @
{
% \+ y7 O/ A2 s% k. f du1[t].yi[w][r]=du[i-1].yi[w][r];* K2 q6 O0 B! m: y9 w  B8 @
}, R7 D* `1 s8 L6 p. i' P8 J
}$ s: w9 s# C1 J2 Q
for(w=0;w<N;w++)
) ^3 U0 e/ K: F8 i3 Q: ~+ a# X+ U" m( O {
0 [2 K5 Q9 U1 S for(int r=0;r<m;r++)
  L2 w- t2 I5 T8 n, j {3 Q4 b6 ^' W9 _* }7 S
du1[t].er[w][r]=du[i-1].er[w][r];
! W+ x8 `7 R. M% Z0 I }% T1 D+ t, G. ~( R4 `
}; z1 z0 o$ ^6 }  [. E$ q/ d1 ?  w
t++;: p* M& t5 A* h5 n7 M2 r
break;
0 y  w. P! ]  h% U* W2 e+ L }
2 z" x: V7 V4 m% o }- L( s% G. l4 z" E4 M0 B- ]2 ]4 h" p

  v1 e! y8 \8 n, V+ ^4 R$ V% Y  w6 k; R% `; H0 w! x
} . \1 J, e9 v+ r# J0 |4 Y
}
4 |- T' K3 L7 M: o8 n( |) L! W) [+ b/ [, J( ^  z5 }. h. f

" K# N/ r; M7 Y( f) V- s5 M; F& Y. K; U- K6 E: L7 ~! U: M: ?
8 i3 n; f' x+ m  V9 q

# v: e) Q5 M% k0 D3 i( r% @# s: h8 e( L8 m' w  |5 d, v0 X
6 h: C5 K* D2 \( t# i" o; A" |2 e- U
void jiaocha()8 W6 G/ K, Z; w- v+ T+ S
{) I: |1 H) Y% v( t
for(int b=0;b<64;b++)5 V( s& X! k) O( \9 E) x2 m, h7 F
{- B6 p; m+ _+ ~  h; e
int p=0;
7 \! T" B$ `  n+ c3 ~, i int q=0;
! @7 u) g- [1 L# E' O& O3 T3 g0 C int o=0;# b' d8 T+ W. _$ a6 R5 l# @5 w
; Y: g* Y/ l/ b  |8 Z. V( Y
int s1,s4;
. d+ i& ^; j/ V float s2[M][N],s3[N][m],s5[M][N],s6[M][N];
6 U4 I, V' o: u2 K$ H p=rand()%64;
0 ^1 q% |3 C8 P6 j8 I8 v q=rand()%64;: ^) p1 O6 a- h$ b9 t/ ^
o=rand()%10;
/ e9 s2 b2 I3 n o=(float)o/10;9 X: i& R9 a- U
if(o<0.2)
/ V1 o2 ~' ~8 h/ z. H9 }& X, k" M/ Q' T {2 V' B. ]" K5 A. t4 a- d) \. K
s1=rand()%10;: O3 `: `7 C1 b" a
s1=(float)s1/10;* Z% M. |- d  x9 X* S
s4=rand()%10;
6 W8 t  g. ]) ^! e6 i9 G) j s4=(float)s4/10;0 I: O* M+ X* s! N( S# n

. b: X, F! T+ G: [ for(int w=0;w<M;w++)$ n) G) C+ f! E7 K" ^# r1 I
{- F# \$ u( l* ~/ Z
for(int j=0;j<N;j++)5 y% i/ `. X; ^9 G- O% o
{
* F" R8 v) @( G. q0 J& Y3 u) t s2[w][j]=du1[p].yi[w][j]*s1+du1[q].yi[w][j]*(1-s1);3 v4 W% x/ u9 V0 D$ _8 n5 E, k
}& {' n: k8 n0 n( z1 F* B$ |1 D
}
# g/ e% h: \; c& `$ { for(int i=0;i<N;i++)
6 p6 {* i) t! D* E {7 D, ?: e( C2 a8 M+ w
for(int j=0;j<m;j++)" s. T; s- p& B5 v8 u
{
$ u2 Z" D! ^, h, U/ T0 ` s3[j]=du1[p].er[j]*s4+du1[q].er[j]*(1-s4);
# J) x" }: u% j" {/ [ }9 N% T1 b& L6 ]! x( _/ r
}; O# z  e' J  }

: {' c1 s. c/ ]- b for(i=0;i<M;i++)
/ s/ W" s/ j8 c {
, |# b4 C+ y1 ^. H* A for(int j=0;j<N;j++)0 n$ l2 x* v+ V  G$ A$ d9 s
{  L2 H" ?- F! J# S& o8 ~' s+ D/ R- H- c
s5[j]=du1[p].yi[j]*(1-s1)+du1[q].yi[j]*s1;- p8 B, }3 o6 C  ]/ F% l
}6 f1 x2 i: b5 ?! e3 b
}  d/ x" a9 e4 Q
for(i=0;i<N;i++)( I7 z( ^3 L: R5 r
{
/ j; N( o% E+ s" m for(int j=0;j<m;j++)
$ P6 j4 u( P  j& j) B {
! w9 t7 Z1 p, x6 j s6[j]=du1[p].er[j]*(1-s4)+du1[q].er[j]*s4;4 z2 M. m; }# M$ x' S1 @
}) `& T( M3 I! J
}
0 X5 j& |9 L" X" Q9 c6 t+ V/ a; E7 y2 g# u! T; T+ v

( z- r' Q0 M7 _) A6 Q' E
3 G4 [( V4 M2 }1 H/ o for(i=0;i<M;i++)7 g0 d/ E: r; \- r7 `4 y
{
+ g) R+ f, O6 }  N  _: C6 w* S for(int j=0;j<N;j++)
1 I5 l# }" C0 V$ g0 J' q5 _ {
5 }& @" u/ X0 l  ]& o9 i du1[p].yi[j]=s2[j];4 F& V5 j. F& c0 p( U- [
}3 F. C* Z" M& g% e
}
1 G' Y" Q( ~5 b, s9 @! i for(i=0;i<N;i++)4 M7 y8 @. J- I: Y
{* `4 M/ z2 X: Z4 M, X
for(int j=0;j<m;j++)
6 F% |4 q+ d0 K7 i8 b {
' B) J- \% x$ g4 z; v, [ du1[p].er[j]=s3[j];
6 u/ ]# E; x& d: u, t }' P; ]0 {6 X% e' u
}) s. W+ ^5 T  I* \& l  ^

: h, l+ W% C$ ^, s for(i=0;i<M;i++)
# |/ N" g5 Y3 g* L$ e+ Q/ e9 r {
, K, s" z3 `& r4 D6 v. b! Z for(int j=0;j<N;j++)( s+ ~; w0 ], [  m/ L
{
7 E2 m% ?+ x0 j1 @ du1[q].yi[j]=s5[j];( s8 t4 g8 ^, j, |$ B) ~& Q
}
2 B2 l) u4 {" O2 a, a" B! n* r }
! E, r8 C4 b- e4 Z3 I for(i=0;i<N;i++)
4 b, r4 f& s, g- X {
% x; \0 P0 g+ F3 F# X- U* g3 p for(int j=0;j<m;j++)0 ]5 a" ?* H1 \) x; d2 z! B2 p
{9 ?$ U  [+ `0 {+ A
du1[q].er[j]=s6[j];
( C- H' E: J7 E  t2 z }/ ]" X6 J9 X+ S5 B2 E9 T, Y; j' s/ q
}) F4 H/ g2 w7 o* f) ~$ Y) O
}0 f8 v- e- X3 n  y" I3 z
}" f* g/ g5 f6 r1 x7 Q
}" ~& s! r; s. A. `8 |4 {2 c

- N: U" S# E/ q# @% u1 N, V% r0 w- t0 I: D! V8 A+ U( j+ f! V) d3 M6 ?
void bianyi()4 x; F1 W. A6 O$ C: B3 P
{
3 Q/ ]; G( C/ Q2 w7 j int p,o;
. {4 u8 A4 I% r9 F; C9 V3 B float r1[M][N],r2[N][m];8 t9 \. v- z. w6 g7 w1 G
p=rand()%64;
! c: F4 A" g6 ?6 L, S  _ o=rand()%10;; H2 R7 r7 g2 p0 y  g  w
o=(float)o/10;, E) D# p* }# Y; H# R+ l: k
if(o<0.01)
/ x" i8 ^+ j/ L3 F, g/ n {
9 ~& f& p: I6 k1 a for(int i=0;i<M;i++)
# ^) s6 b* Z* G9 ? {. |" W' H, U) |! t2 h$ ]' P
for(int j=0;j<N;j++)
3 Q5 x' V5 s/ R4 ^" \) O$ h {% M9 c, E; |& Y+ ?- f& O
r1[j]=rand()%10/1000;
& [/ g4 U% z1 e# h1 C, H4 `, C9 g6 K# p }! K9 b& K5 r2 m  a' c6 d
}
- n6 o5 m& i% ] for(i=0;i<N;i++)
# Y) ~; \1 ^" L" N- |" z {
) i1 c; w* w1 U; e# n3 k for(int j=0;j<m;j++): }# e6 G, j+ q7 p/ T1 g' ?% T
{
$ g6 ~5 o" R# G r2[j]=rand()%10/100;
" a/ f% n; m/ a. M }
% {6 p0 K! T7 l) s# }$ w# D/ c }
+ M' N2 F- l% J5 ?) K
; e! D4 R6 _$ j2 c# u" c for(i=0;i<M;i++)
% U' [& c& s/ @, I' s( ] {
) p: f7 e& d5 p8 B for(int j=0;j<N;j++)8 Z* C& Z2 j! |2 i+ h
{
9 S, N, m+ M: o. p- K1 M du1[p].yi[j]=du1[p].yi[j]+r1[j];
8 K  U: Z6 h& m5 K5 ? }5 Q5 D; |$ E$ c
}
$ ~7 h" R+ N8 n; C# D/ r1 D7 A for(i=0;i<N;i++)4 g$ X, V1 x% E
{, _0 c9 N+ J) k+ N( K4 Y6 S
for(int j=0;j<m;j++): n+ u2 N7 t* `" U
{
. z6 r4 ^- R% a; e9 ^" c9 T du1[p].er[j]=du1[p].er[j]+r2[j];8 b& v7 @9 J6 R1 n# G& O( ^
}
/ O2 y% \$ e/ D5 F0 a9 e' f! t }" p: `# G0 J8 z! h! B

8 s8 H# ]5 Q' `7 ~) ]; a }0 _. r, K! l7 F2 e" g3 O' W) D0 w
}
7 [5 W8 U( h) B9 K- {
5 @: h: z; y" A' N6 T0 U5 O. q  f  t6 O5 e6 ~
3 K( C1 A/ r; ?, N; J
void change()
" [4 S" ?- j1 V% W5 b{. r* J9 p4 ~+ Z" S
for(int i=0;i<64;i++)
9 J' ~# s2 i+ V3 q {
+ k% I* |- C" h) L% V for(int j=0;j<M;j++)
' o5 x: Z" S. g1 d/ i0 O8 i {
: A- j, q+ ]2 ]$ P& x for(int k=0;k<N;k++)
& E: ~6 M7 ]% x" i {
" \: B' A, f( ~4 }8 T du.yi[j][k]=du1.yi[j][k];$ c; m& s' B9 T3 K$ m
}5 \+ P) J& w6 i  U3 A' m* s$ l- s
}
' t% d. E4 v! Q2 L6 e6 a  ]5 l( x for(j=0;j<N;j++); a3 ]8 S: B4 L! `7 m+ Q- J: |. q& G
{
. z% X$ j' k1 P  m( j for(int k=0;k<m;k++)
' S! J& [. C5 O1 D* v# i3 p {2 J- u* a8 C" i$ V# x8 k6 C% {2 D
du.er[j][k]=du1.er[j][k];) c( M  r& i/ V" @4 R6 }
}
7 b4 B  ~7 O3 t& {/ J9 t6 D! ? }
* J/ U! `+ @* l; H: p! _ }
# M" b& Q5 I8 p, H1 _  U}% v) p* v' k) Q: F/ z+ I0 V/ s

! A8 S$ _* Z: D5 ?. R; v
) V4 @+ f4 C7 L+ Q: I% F- D) Z) G3 G% P2 u% L% g* S

+ a) B/ C! s0 p! W( _void raodong()
  `, i2 v3 q2 Z& @- I, |$ M: C9 D{
  H0 ?' H! v4 f/ J1 K for(int i=0;i<M;i++)
8 k6 z1 |' G+ O6 J {2 ~' C. W- r# @' `1 ]/ X% E7 U
for(int j=0;j<N;j++)
. Q' {8 }* O. S0 C3 W( z {
% a& G- P1 i( \: ` rao.yi[j]=du2.yi[j]+(0.1-rand()%2000/10000)/10;
8 f7 [6 s6 ]% C2 R$ \2 Z  V& L }9 m9 h. f# T  _/ M$ P3 G; p
}# U2 }- W  h9 i$ `; g* R7 V
for(i=0;i<N;i++)
- P9 X7 c3 [3 E3 s" ?& r {
" L4 p9 K( o. P9 U; E$ D2 O for(int j=0;j<m;j++)
  B$ P1 |! }1 J9 @6 R {
0 T4 e! V0 }' r9 t6 ^ rao.er[j]=du2.er[j]+(20-rand()%20)/10/100;
0 x, O4 Z9 ]% y& g! Q1 d- x }
9 d6 E* |1 l+ Z5 u }4 N) y6 I2 h# V# l
}
: g; I9 J8 z0 w8 C
! A1 Q/ ?1 u1 h0 Y% k3 j$ l  e) S0 z* H1 b. _9 X0 J/ }$ ^
. g% u6 R+ }1 x( a" z% a8 ~

+ a( o( B! M' z9 M( v3 s* G+ `1 y8 `8 ~% F6 H: x
, D; W  }+ _) s* L) P
int qushe(dui d,dui r,float T)( R9 v; o( v( C
{
1 u- W: \- i. d# K/ Q1 u8 p4 W: N float w=0;
; K8 i7 m$ G- P1 U! h; S$ Z1 s float p=0;
$ O# ]5 k. R: `& h' K0 B8 q float h=0;" N9 A+ ^. t* A4 _2 D
w=d.Fi-r.Fi;
2 y7 b& g( [+ ^- g: f0 X8 a w=w/T;- [5 o9 d! T4 D, f  w
p=pow(e,w);1 G2 ^) L  K2 H
h=rand()%1000/1000;
, R8 D8 C5 c% }5 R: m9 }* h# e if(p>h)
1 ?: ^9 @9 |% \2 K. H {7 t! C* {0 U; F+ _4 }3 T; r1 K) e
return(1);6 D7 V* J  X, _& w5 Q& F
printf("\n概率：%f,%f\n",p,h);' X  O, I1 [9 M+ h: B
}, @: H7 W9 R" N% s3 C
else
9 e' l/ b! V6 Q! L, ]* l7 H {+ f" k( E9 ?: t8 Y2 p
return(0);
4 ]0 i5 B, J' U0 l }6 S, V& b+ R# q. R; ?& D- b
}
4 a- K4 ~* K9 R. K/ n$ U" C$ a2 U2 r$ Q$ R7 x! X( K4 \2 S5 A

* l( v  c0 F! W4 c$ ]) q
2 E) w$ C2 U6 k) e) v# e- @' K+ G" R0 t) ^- ]- h

+ j4 E7 T5 b5 c1 Svoid yichuan()8 m% \, U) c* F: X1 B
{
" s) O& I3 V2 W1 V+ o for(int i=0;i<64;i++)
! t; G, [0 J0 D* D0 g& v* r {; b  \' U/ {! M: j2 O; G% e0 }
for(int j=0;j<M;j++)
0 \7 Y: d. S- d+ K {- e! v) Q& l$ I0 K$ f$ F# \
for(int k=0;k<N;k++)4 D1 W8 X: z; D& |( R
{
  z1 e& C8 M. z3 f! A% m; I  b gd.yi[j][k]=du.yi[j][k];2 k5 v- D  y2 t9 M( f1 }2 [
}
/ n4 e5 \& t' f }/ p- {+ m  j* a" Y( [# O

6 E* r% e! c; f/ @. h" _ for(j=0;j<N;j++)) r- o- [2 T6 D
{
8 a" I) ^) F) X for(int k=0;k<m;k++)
" `  K2 D0 k4 c' V" t5 H( W8 Q {
" v4 J7 ]: d- o% }5 t$ {  M gd.er[j][k]=du.er[j][k];3 U! n& ]! k* l  M' |
}
: u* x3 t* Z+ n6 ~( P }
5 w+ N) Y* V' M- a2 a& U0 N: X: g9 w5 u9 O

* p% o! K; O2 m6 o5 B du.Fi=5-shiyingdu(gd);* @' k- k& b7 O6 W' |
; @5 a! N* t- V0 [
}
8 F+ p/ d; Q! h- ^8 t
$ h% l- K1 c2 v gailv();
1 E8 g* r6 k- |5 Q8 t lunpan();
1 \1 }3 l" D$ U" Z jiaocha();
* K; u* {- `+ B& `* M bianyi();5 W9 a$ X  V, g& h0 t% E4 M8 r' `" t
change();
1 A2 K" t: n, _* Z( l% @
2 \/ s& P4 v  m/ Q+ ?; s+ Y# d for( i=0;i<64;i++)4 N2 J) E# |8 ^& F4 ]/ @3 l3 \. {0 c
{1 f5 `2 b9 S5 y3 I" V* i
for(int j=0;j<M;j++)) A. R) i/ H' b% \/ R
{5 [4 F, x! Q, D
for(int k=0;k<N;k++)' H, S. {. H% }; }) j% C
{
' l$ ]7 E$ [4 V  d: e; j gd.yi[j][k]=du.yi[j][k];
" ^, l8 H* C/ Z) {  \ }8 f' D3 {% y6 q. C, I
}+ C& f) I2 i5 q8 Q) ]$ D
- C! E  `7 o' J) J7 Q7 P
for(j=0;j<N;j++)6 i  S3 Y: Z  L$ u* d
{
! f7 Y- H; v0 C' h# r8 h# j for(int k=0;k<m;k++)
- J9 r* H4 Z2 {1 h {( u! k5 P6 E( U5 A, _& O! y$ v9 ?
gd.er[j][k]=du.er[j][k];) w% R# g# c& u8 K  w0 @" z
}
8 [( y% \" z. Y2 A  h }
& V6 r( \* i3 a( d+ G9 l
4 R( v0 ~4 B: E6 ?" E' F1 m, d4 T6 I% o& U
if(shiyingdu(gd)<0.01)
$ I  A2 |# G+ e# n& n' L( Y {$ S$ A, @! {3 d) q; b! R3 J  v
for(int j=0;j<N;j++)
0 G5 m4 R- j  j( ~+ u {. D" k( n$ w0 G* O; w2 y2 Y5 l$ Y
for(int k=0;k<M;k++)" f/ i& V( x3 y6 {
{1 \# u+ M; \  X& e. M
printf("%f ",gd.yi[j][k]);9 P( l9 D+ R$ b0 V
}
- ?! g3 q* B! b9 v }
& c9 u0 J% L: V6 R% B7 g printf("\n");% _: H4 r. J9 ?; ]$ j6 e$ O
for(j=0;j<M;j++)
! d/ I' e( v3 n+ L+ {& J {2 S5 A3 d2 G# m
for(int k=0;k<m;k++)! L# P  t& B+ p" O' F  m( S0 N+ T
{
/ @- t/ T- a' m printf("%f ",gd.er[j][k]);, L( `) B* O' R0 W
}) }! X4 [8 S2 ^& f" p! r6 R
}
' X" {6 ?2 K! `  h( G4 V' Z printf("\n");
+ d# [3 x! L1 f+ q; ^ exit(0);- f$ a$ t; x6 a
}
% U0 ]% h" w$ D$ `% H  N- f1 x& G4 {: @1 ^# C9 u3 j+ d5 z
}
. d9 {4 j  d7 _/ g' J+ i7 h$ A) N. U. |; T8 w; A; X

' C) S: J  F4 R  o" q; N}
% V- z1 G0 l4 z3 c  k; K0 ~" j1 W% C# j/ J: J
void tuihuo()2 j( u, }2 j4 w/ w+ i' A
{1 h5 z  E! t! ]! \
for(int i=0;i<10;i++)" b) E: d1 f) K
{0 c% d/ J8 t2 k/ l$ B! [% O5 N, L
int sj;
6 ^0 g( N' W! u8 C' i' X sj=rand()%63;
5 ]& z* u& x$ w  ^% d) u0 d# D7 t2 f( k. ]) g# C4 q
for(int j=0;j<M;j++)* l* y, D2 f6 _1 }) G. }
{
- b! m! P! F+ L$ e for(int k=0;k<N;k++)7 D! r( z, f% q- [8 l
{) l/ w" M% m4 u0 u6 s6 R/ s
du2.yi[j][k]=du[sj].yi[j][k];
& T% h! l% S# f- L# I' L; [& l$ h }3 H+ u* V- j( R( r: c
}
% [- Y( i8 ]" q7 F7 h& a0 ~) `" V8 m1 y" O% }4 c
for(j=0;j<N;j++)8 [. _! W! Z6 H# @7 ~
{+ l8 `' ^: U7 D$ W: B, ^; E
for(int k=0;k<m;k++)
2 Y& @: h& W" j  N {" |, t7 ~, s" ]5 ~. y: K
du2.er[j][k]=du[sj].er[j][k];
0 A2 W. s" b( J. Z3 w }
: k, ?% D! i9 o  z6 v& _) {9 s6 ? }
7 V7 S' j, Q6 k$ ~* C; j/ [ float Ei,Ej;
& ?$ a' p% F1 T4 Q/ c+ M: x Ei=shiyingdu(du2);
( p; j- x5 n6 P: Z raodong();
  d0 t4 ^3 _+ X3 C  x9 u9 V- t3 ] Ej=shiyingdu(du2);
2 O% d- i7 l( y2 M7 J while(Ei>Ej)
% ~8 J. B/ |  s  z' v {& o, V1 ]! t) g, L
raodong();' U/ b( y# ?5 s- ^4 x5 `$ J
Ej=shiyingdu(du2);9 l. C3 a9 V/ I& s# B3 u; _4 y9 U
}
+ G# I; B. S$ }1 v if(Ej<0.01)) H& Y: L. C* g8 ^
{5 |9 v* m  D1 b7 ?: u) k8 k" x, Y: L
yuce(du[sj]);
6 i0 r0 S1 D1 ^" x- | exit(0);& Q" u; O; Y0 \) Z- v4 ?) Y# P
}
5 A, w5 x: b) j& u; Y/ T if(T<0.01)
# o! t$ N' ~0 s. v$ _; [ {
, g+ f/ P" o3 t- X3 C yuce(du[sj]);; V6 V% s! [, c" Z1 k: z
exit(0);
. J- j" H3 [8 [% B# g3 ~5 ` }# i) \7 j+ N. q! [

& t' b% E# O9 N9 W7 ~" P; r! ~7 b% p8 a3 T$ \
if(qushe(du2,rao,T)==1), I0 Y! _4 a; P* {$ ?+ V( @$ H
{
  s3 Z8 s- T$ ~# m; R$ X3 R for(int j=0;j<M;j++)
* J) G. {2 ?3 y: U {
* X% j" T2 U, ^# D- ~' z+ V for(int k=0;k<N;k++)7 i( K# C" i- w0 V/ i: X- M
{  M% {- j; v5 S' A" i9 h" [
du[sj].yi[j][k]=du2.yi[j][k];9 e. k+ ^5 `, D: i
}
/ y+ P' m+ V5 t$ a  X }
% v& w+ L! d( W( F3 G$ P% F. \$ p0 {6 {7 D
for(j=0;j<N;j++)9 `2 Y$ v6 c0 @) i0 e. c8 W
{3 d4 E& {* U( e1 I
for(int k=0;k<m;k++), f4 l  O+ E" i4 y4 O
{) \. c4 x2 u& Q/ d$ M' m( P
du[sj].er[j][k]=du2.er[j][k];4 l0 Y" u* T, u5 w. ]: @& ]; _5 P
}% G9 @0 v9 |5 Z* c) c! v
}
: `+ n; T& _) T- [/ N2 } }
! ]$ x* C( r  ~1 J
; @- w: u( k0 S) q9 `; i  ?' d
. Z& \3 e% s) Z( J8 B9 h }
$ c7 n( a; V! Y5 U/ D
) D$ X( P1 ]; a  ~3 d7 G' n( u7 u; @7 h! _0 a9 q% ~7 R
}
- ]5 @2 v. a5 u
0 l' d0 z0 p' [  Q& Y6 K# P$ o
; W3 [: Y2 a" y% _& s9 B( x, a* X/ _

3 H: }* F' T5 h
- [! p3 l# [  c# I: n* \& L% Z% s) U5 t. p, ~  n; K" C' Z
void main()3 F2 u$ R7 i5 k1 t6 \9 A* \& A
{ - [9 ^3 ]& m& j3 g' u
chushihua();
) _1 f8 U. J8 T0 i: \) Z int n=1;# L& C: t4 O9 D' P! D4 E! {  \
int jishu=0;
' }# ?8 R9 H% d8 A) n5 l, \ for(int i=0;i<2000;i++)( k% a/ g3 M  b! i  w& G$ F; l
{
# @2 S* e" E6 r8 @2 a: j: W5 G
! B( \3 W! ?0 [
& c) O/ P0 m) B7 g yichuan();
8 b% S4 x2 Y1 T( l tuihuo();5 V$ K" J% a. i: q3 M3 D6 o$ x
8 f# r# |% ~7 v* f
float p=0;
  N5 g5 v5 I( v) ]& L% Y p=log(n)/log(e);& ?5 G  [. R* z/ }+ V* e, E
n++;/ \( q( a" b3 _9 d
p=1+p;( x7 ?; F4 e1 i/ C. ?: I5 n- ?
T=100/p;3 @4 \. i: d' K
printf(" 次数:%d,温度：%f\n",jishu++,T);
+ k2 H8 R% ~7 H- d }
7 D5 B( I2 ]8 l6 q% F yuce(du[0]);9 v( @5 u+ j4 g* W& _" M
printf("\n\n");9 T8 z5 y; P+ {+ E7 X& G
printf("\n");. E* b) Q# N- H
exit(0);) h/ Z! d: _- Y
  `3 A) s# J  r' T) I# _8 M
}5 A5 y' u7 d8 V0 _7 J, I
(3)人口动力发展方程C语言源程序清单：
0 w5 ~& r# E* @; |- W* W' s#include <stdio.h>+ q* Q! Y6 X5 G( _  f8 P( Y% [" t
#include <stdlib.h>( ^$ B+ J9 M( [: x
#include <string.h>
3 j% V9 a# \9 B. g9 p$ E. {9 N#include <ctype.h>
: k# p* f4 O( L4 ?+ T/ w
8 @0 X6 J+ @+ ^/ G7 w#define N 100//年龄分段3 X. [% U4 ?) Q- ]7 B7 w' w. |% F
#define M 10//未来几年3 N/ J/ q+ M1 H( N/ y* X7 q
void main()
/ |# j3 J$ ]6 ^{ # R5 I  B- o4 J$ X4 }
FILE *fp = fopen("实验结果.txt","w");3 R; _8 \% [. [9 N

' l" y( C. B. I) f" g: \6 R float num[11][N],out[N],in[N],Nin[N],Nout[N],ch[M][N],gd[N+1],m,t[M-1];
. t) y$ p, X: d float c1[M][35],c2[35],c3,c4[N],h[9];
9 N, r2 U( {, n; t float he;/ B* `  I3 m+ Y6 Z3 S$ q
int t0,tn,jh;( j/ D  Y2 {; @* O9 }
int i;
/ G% T9 v+ m1 y: O1 B. {
9 }8 ]: J2 Q5 G: s+ K7 k t0=1;
& l' A. i0 Y' X3 Q$ C; v tn=10;" ?* L) z( g1 e
( D& [% w$ S, K* O
for(i=0;i<N;i++)& j) P9 Y5 K" F
{7 e. F! D% ?, l
printf("请输入第%d年%d岁的人口:",t0,i);
( U# m: l$ \2 A% F$ r+ n scanf("%f",&num[0]);% P) U$ W4 ^6 h& d
//printf("%f",num[0]);! [1 J' l' t2 {+ @" I
}+ G& l% `  Y0 h% w

% T2 V% Y2 r+ M6 ^0 |/ t' o0 K0 A& P$ `1 P5 s% r0 k9 _* _

( f' k( d% `5 J3 n* K! H! [ for(i=0;i<N;i++)4 {! c! o9 f. E+ |
{
, O( F( A- U4 \8 ~9 P! C  M out=0;& G: a( v9 P7 h0 r
in=0;
+ y& }5 F$ @5 ]5 r+ q' A }
1 l2 Q, R& ^1 g/ k8 J7 t% y0 s* n3 J0 a) Y. U8 q* L( n0 L! [2 u
for(i=0;i<30;i++)* F) w! O9 D( S6 o
{
, S" X/ \0 |* v$ k3 ?" C9 X' Q( U printf("请输入%d岁迁入比例：",i+20);( u9 J6 P# R* O+ A5 y
scanf("%f",&out[i+20]);
1 s+ A( v" E4 y/ N7 v9 z( J9 @# z6 T // printf("%f",out);9 u6 L" }, }& @9 P% y: F" q( f
}8 Z1 x$ r* t# J: V  `

& W% j0 t; C9 h1 U for(i=0;i<30;i++)
' V. L. i# g; I6 h! z {+ K/ z- z# @' y5 n6 b! d1 [2 C
printf("请输入%d岁迁出比例：",i+20);
! p9 }( J& e3 k7 B0 H2 z! d scanf("%f",&in[i+20]);
- {5 y* T" H6 ^; s+ A, ` // printf("%f",in);
' J0 ]0 j: J7 X/ U6 L# g& I }) [& w" m' b. ~& Q( F0 x

  _( Q2 f: G- A" g7 M/ k: d for(i=0;i<M;i++)& [% S5 O  N0 P2 Q  k/ L9 \
{
. q& r/ C: o8 s/ w9 y' W, ? printf("请输入%d年迁入人数：",t0+i);
6 I9 c. N# U1 Y7 A4 O+ {: r) G scanf("%f",&Nin);: \# `0 Q" J6 ?, M# s+ V1 a* j
// printf("%f",Nin);
  q, q6 _: E" X$ W% Q" d" n  H$ y }
( i# P5 ?: E/ E7 S for(i=0;i<M;i++)0 L$ B+ ?& b" p4 O# |& h
{$ n5 Y7 l+ I; f4 T. P5 y
printf("请输入%d年迁出人数：",t0+i);" x# u+ F! @# g- ]
scanf("%f",&Nout);4 M- R1 Q# H( ~& X& q" K
// printf("%f",Nout);
" s* ]& X* e' B1 H" O" G }
, ]" B% c  o/ T& ]0 F5 Z$ G+ u+ \' R
for(i=0;i<M;i++)
% C( M2 z, B0 Z! }+ Y# l  n! ~ {8 _# p# j; l: o
for(int j=0;j<N;j++)" D( f; G5 g7 E, o; {5 p
{
" i& y/ M5 ]9 u* d+ T$ u: f; { printf("请输入%d年%d岁死亡率：",t0+i,j);
) ~# ~# O" }# q% c scanf("%f",&ch[j]);
- p0 D% J) P) B  h. S7 \ // printf("%f",t0+i,j+1,ch[j]);
, F# e6 ]- }) Q5 }$ f }6 L6 X' x" w7 t0 `, S
}
# w4 {9 \5 C' ?+ b5 H( m7 o/ `+ |5 f , l8 }  e) c$ h
3 B8 O* ?1 O# }$ \
for(i=0;i<10;i++), N: z& H. n! v
{ for(int j=0;j<35;j++)
+ ^: J3 r( d' l% ^, j! r {& U$ _& X! N( h8 I* g5 W0 p
printf("请输入第%d年%d岁妇女比：",i+1,j+15);
3 B8 Q& M( h. V% \' V1 t( \# X8 o scanf("%f",&c1[j]);; x7 o" x' |3 ?
}5 B% b4 s1 {* }) n6 ~4 T, \

; w: N& C4 n2 @8 g% v; _4 K; Z }
' ]* X& B3 v9 x2 w% U! K# ?2 B" i* Y2 R
) p' Q. X: Q( L  H. J2 z( A
for(i=0;i<35;i++)
5 F, f, C& Z4 _, h& ~: [6 e2 T {8 h$ i8 S1 H9 {; U
printf("请输入生育模式率");/ ^9 b1 t  a  b
scanf("%f",&c2);: X  ~/ e. p7 W$ e$ A; @& R
}( u7 T' d3 T/ P; Q" _9 \$ z+ M: L

' o. Y: w2 f9 H0 e, a( d6 \
& t# W8 |6 L7 q6 \3 H7 i
8 Q1 _& u1 n/ w' K* D2 I" n3 G7 Y printf("请输入总和生育率：");/ o9 M  \& s( ^4 X$ w$ B
scanf("%f",&c3);: z8 `- _2 p7 \' }. o

, \1 B8 W, J% @) V+ P  z4 }4 B. c" R
- r! O! G+ G+ K% t
for(i=0;i<10;i++)
* A/ i, g$ K2 S- u& g/ [3 I {* s5 D( C" X6 H0 H" |2 M
printf("请输入新生婴儿存活率：");
" X6 W- E( J; @$ [7 F scanf("%f",&c4);0 ^2 h( _6 L" ^6 n  s9 o
}
3 V. G" h4 h( s# V+ N
. X( n- W# m$ I  D: ?8 H( I he=0;  b. I) z0 w8 c7 \
for(i=0;i<10;i++)3 z" V% k& s+ \* Y4 N+ k
{$ B% o7 b! r/ O# @; a1 h
for(int j=15;j<51;j++)
6 B8 k7 C! y3 j% s {
# I# V3 V* o# D4 B0 } he=he+num[j]*c1[j-15]*c2[j-15];
& U0 L! r1 K2 I) d, P9 a6 R }
0 P3 r% F- u2 ~/ z% F& Z. D
% l! ^- q) G# \7 O, e5 r' w h=he*c3*c4;
, D1 m! W) z) I1 `$ P% o: z  ?
' [5 @2 M  |& P* K }
8 v1 W9 Y6 M: h' k8 W/ W$ Z/ ?# R
. V! ]3 D# Q0 s- U for(i=0;i<M;i++)
) N5 |: F- k9 m- w {# Z/ F% C5 Z$ @  g$ m
for(int j=0;j<N;j++)
# p: N1 m4 j5 Y' M" p+ M {
8 z  [' r) Y+ \/ c( I m=in[j]*Nin-out[j]*Nout;! d* y0 H' r9 k$ v
gd[j+1]=num[j]*(1-ch[j])+m;
5 e! S3 u0 t% u
1 l8 U; _7 M* @6 O }
6 k: `/ T2 ]$ ~- r8 g# x gd[0]=h;6 ^% g+ c8 J) V# ~' E) `! J
0 w; m1 K) w. l. R' m

7 d* ?  B- b" K/ \8 B for(j=0;j<N;j++)3 Y  u8 R- z+ j: t- U7 O7 F: v
{4 J/ [% x3 p5 l! ]$ ?+ u
num[i+1][j]=gd[j];
( ^2 R5 a: v+ b+ c }
. p2 b) ~+ O: b/ b1 ~; K }  ^2 ~" [( N! U) a4 R  w: i- A& P
  _. y- T0 h4 S" b$ u
printf("\n");. N, `9 H9 m9 q8 |
for(i=0;i<M+1;i++)3 y0 z) z  p# I: t% M1 t
{
" c" O& a) d" D% W for(int j=0;j<N;j++)
% F+ ?8 R) D1 E7 [( S9 V {
# L) Y: Z9 V) U9 G4 w0 F0 I4 E  x printf("%f ",num[j]);0 r& n. f) _& k! W4 M
fprintf(fp,"%f  ",num[j]);! c/ z8 I, I3 s! o4 `
}
) d7 s. j; `; C* K0 q& J printf("\n");0 r7 C0 }: t0 ]
fprintf(fp,"\n");5 {, f5 |4 N# ^) r' I
fprintf(fp,"\n");. ]) F. z- p' m0 m( p: i
fprintf(fp,"\n");
# p% `5 }- ^3 S( d }  l2 W  L/ A, U: |3 a% C

' Y% r' N- J. b1 R! ` scanf("%f",&num[0][0]);
4 F, B! Q9 f# {! }
: w4 h; `( L, Q6 S}
! C$ H8 y2 u( h  W9 b- K
& `* |$ Q8 H% U4 @8 x; a: v# }
作者: 刘李 时间: 2013-1-24 13:32
哇~居然通过啦！好开心呢
作者: woxingzhe 时间: 2013-4-16 15:30
楼主好人啊
作者: 淡语夕林 时间: 2013-9-13 15:56
gooooooooooooooooooooooooooooooooooooooooooooooooooooooood

欢迎光临数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5