《系统动力学实验报告》由会员分享，可在线阅读，更多相关《系统动力学实验报告（25页珍藏版）》请在金锄头文库上搜索。
1、系统动力学实验报告姓名：徐键班级:管科131班学号：5504113023学院：管理学院、背景：高塘乡德邦牧业有限公司是一家大型种猪养殖场，在高速发展的同 时存在两个急需解决的问题：1、养殖场猪粪尿污染环境；2、高塘乡已建的300 余户用沼气池大部分因缺乏原料致使沼气池闲置，农户买化肥、农药种植粮食、 蔬菜，农作物受到污染。二、基于顶点赋权分析确定规划实现的管理对策：略三、基于逐树入仿真技术建立仿真入树模型建立流位流率系：(年出栏L1(t)头)年出栏变化量R1(t)头/年),规模养殖利润L2(t)万元)，规模 养殖利润年变化量R2(t)万元/年),(日均存栏L3(t (头)，日均存栏年变化 量R
2、3 (t)(头/年)，(年猪尿量L4(t)(t猪尿年变化量R4(t)(年)，(场猪 尿年产沼气量L5(t)(m3)场猪尿产沼气年变化量R5(t)(m3年),(年猪粪量 L6(t)(t)猪粪年变化量(t) (年),(户猪粪年产沼气量L7(t)(m3)户猪粪产沼 气年变化量L7(t)(m3年)据实际意义，将流位流率系分为两部分 第一部分一一生产销售利润流位流率系(年出栏L1(t)(头)，年出栏变化量R1(t)(头/年),(规模养殖利润L2(t)(万元)，规模 养殖利润年变化量R2(t)(万元/年)，(日均存栏L3(t (头)，日均存栏年变化 量 R3(t)(头/年)第二部分生物质资源开发流位流率系
3、(年猪尿量L4(t)(t)潴尿年变化量R4(t)(t/年)(场猪尿年产沼气量L5(t)(mA3), 场猪尿产沼气年变化量R5(t)(mA3/年)，(年猪粪量L6(t)(t)，猪粪年变化量(t)(t/ 年)，(户猪粪年产沼气量L7(t)(mA3),户猪粪产沼气年变化量L7(t)(mA3/年) 第一部分逐枝建树逐树仿真建立生产销售利润子模型(一) 年出栏年变化量R1(t)(头/年)仿真流率基本入树T1(t)1逐枝建立的R1(t)(头/年)前期流率基本入树T1(t)见图3.1O v年岀栏Ll (D头)岀栏皐变兀头拜)豔离蠶悴疇划岀量表函数B1A11 (t)量表函数B3A11 (t)日均存栏保 障因子
4、岀栏年基本变化 量All (t)(头 年)+2016 年至 2020 年规划基本变化量表函数B2A11 (t)图3.1年出栏变化量R1(t)(头/年)前期流率基本入树T1(t)2建立年出栏变化量R1(t)(头/年)流率基本入树T1(t)各变量方程：略3取DT=1.仿真区间20052020年，对T1(t)向后看设调控参数进行仿真检然后，对T1(t)参数调控进行仿真，获得仿真结果T1(t)曲线图3.2和图3.3年岀栏LI (t)(头)年岀栏Ll (t)(头):Current 图3.2 20052020年年出栏L1(t)(头)曲线岀栏年基本变化量All (t)(头/年)岀栏年基本变化量All (t)
5、(头，年)11 : CimeDt 图3.3 20052020年出栏年变化量R1(t)(头/年)曲线4年出栏仿真定量结果可靠性定性评价分析(1 )符合历史年实际变化规律；(2 )符合规划值规律(2015年年规划是年出栏15000头，2020年年规划是年 出栏20000头)；(3 )符合种猪价五年一个周期的低-涨-涨-跌-跌规律(二) 规模养殖利润年变化量流率R2(t)仿真基本入树T2(t)1规模养殖利润年变化量流率R2(t)仿真基本入树T2(t)见图3.47ZT规模养殖利扌|蓉化量R2 (t)(万元.年)L历史年刑润年基本变彳绰A21 (t)(万兀年)4 价格制约下平均每头利润画数A22 (t)
6、(元.头)历史年利润塞未变量芸函数月A21 (t)弋岀栏年变化量R1(头/年)2010-2015:可I利润表函数2016-2020年每头利润表函数E2A22(t)图3.4前期规模养殖利润年变化量流率R2(t)基本入树T2(t)2.基本入树T2(t)中变量方程:略3.T1(t)与T2(t)组合仿真检验T1(t)与T2(t)两树构成的模型仿真结果曲线分别为图3.5和图3.6规模养殖年利润L2 (t)(万元)Time (Year)规模养殖年利润L2(万元)+规模养瞽魏规模养殖年利润L2(t)访元)曲线规模养殖利润变化量R2 (t)(万元/年)规瞬36规模养殖年利润变化量R2(t)(万元/年)曲线4年
7、出栏仿真定量结果可靠性定性评价分析(1) 符合历史年实际变化规律;(2) 符合种猪价五年一个周期的低-涨-涨-跌-跌规律(三) 日均存栏生产年变化量流率R3(t)仿真基本入树T3(t)1日均存栏生产年变化量流率R3(t)仿真基本入树T3(t)+猪尿量猪尿沼气量日均存栏L3(t)X)利润促进 因子IS规模养殖 年利环壮2(0)( 万元)丄量嬴进因子M34t)者尿召气促进3(0)未图3.7前期日均存栏生产年变化量流率R3(t)仿真基本入树T3(t)2. 建立变量仿真方程：略3. T1(t)，T2(t)，T3(t)组合仿真检验进行T1(t)，T2(t)，T3(t)组合仿真，得到以下2组曲线图3.8和
8、图3.9日均存栏13 (t)(头)日均存栏L3 (t)头):Current 图3.8日均存栏L3(t)(头)曲线日均存栏年变化量Rd （t）（头/年）Jr-/、/./2501.0007500200520112013201520172019Time (Year)11S均存栏年变化量朗（t）（头.年）11 : Current图3.9日均存栏年变化量R3（t）（头/年）曲线从三个变化率同步符合市场变化规律可靠性中可以得到两个曲线，日均存栏（生产）.出栏（即销售）.利润年变化量曲线比较和日均存栏（生产）.年出栏（即销售）年利润量曲线比较日均存栏年变化量R3 （头，年）Current 岀栏年变化量R （
9、t）（头.年）Cunent 规模养殖利润变化量肥（t）（万元$年）”：Crnrent 图3.10日均存栏（生产）.出栏（即销售）.利润年变化量曲线比较40： 0006:0001.000Time (Y 已 ai)年岀栏LI. (t) (头):Current 日均存栏L3 (t) 头 j : Current 规模养殖年利润L2 (t) 万元;| : Current 图3.11日均存栏（生产）.年出栏（即销售）.年利润量曲线比较4三棵树模型可靠性评价分析（1 ）符合历史年实际变化规律；（2 ）符合规划值规律（2015年年规划是年出栏15000头，2020年年规划是年出栏20000头）；（3 ）符合种
10、猪价五年一个周期的低-涨-涨-跌-跌规律（4）仿真定量结果实现了符合实际两同步，三同步反馈规律（5）使用了还原论与整体论的有效结合建模方法建立的反馈子模型在Vensim下，三流率基本入树已构成生产及销售两个正反馈环，已构成产销利润正反馈环，见图3.12历史2005年至2009出栏年基本变化-量表年岀栏岀栏年变化量R1決，年）量Al 1(f)头年)201。至2015规划出栏年基本变化一 一+量表函数B34H6年至202屛规划基本变化量表函数B2A11历史年利蜩日均存栏L3(t)X):心年）规模养殖年利润L2（t）X万元）歸）刑润促进进M34（t） 录留气促进 4IZI33 猪尿沼气SL51猪尿量
11、L41（t）:bSi(o#A22（t）（兀头），: i V一2010至；0巧每头2016至202闹头利乳薙敗数BhiJ利润表函数B2A22 A22图3.10三流率基本入树已构成产销及产销利润两个正反馈环反馈模型 第二部分生物质资源开发四棵流率基本入树逐步建树逐树仿真检验(四) 猪尿年变化量R4(t)流率仿真基本入树T4(t)(1) 前期猪尿年变化量R4(t)流率仿真基本入树T4(t)见图3.13图3.11猪尿年变化量R4(t)流率仿真基本入树T4(t)(2 )年猪尿年变化量R4(t)流率仿真基本入树T4(t)变量仿真方程：略T1(t)，T2(t)，T3(t)组合仿真检验(3) T1(t)，T2
12、(t)，T3(t),T4(t)组合仿真检验1) 猪尿沼液开发促进因子M32(t)的子因果链图见图3.14表函数BK11猪尿沼液钾开发 促进因子M323+年可开发猪尿招液中的钾Kl(t)(吨)v年猪尿量L4 (t) 历史年开发楮尿 中钾K11 (吨)猪尿沼液开发 促进因子M32猪尿沼液 ：磷开发促进1因子M322:h +猪尿沼液氮开发促进因子M3210)h +年可开发猪i 尿沼液中的i 磷Pl(吨”+历史年开发猪尿 中磷P11(吨) 猪尿含钾比K10/ f ?后起时间T2!猪尿含:磷比P10表函数BP 11年可开发猪尿沼液中的氮Nl(t)(吨)历史年开发猪尿 中氮K11(吨)n4f表函数BN11
13、后起时间T1v年猪尿量 L4 (t)(t) 猪尿含氮比N1Q图3.12猪尿沼液开发促进因子M32(t)子因果链图2) 建立猪尿沼液开发促进因子M32(t)仿真方程：略进行T1(t)，T2(t)，T3(t),T4(t)组合仿真，得到图3.15的仿真反馈曲线结果年猪尿量L4Time (Year)年者尿量L斗(t)(t) : Current 图3.13 20052020年猪尿年量L4(t)3) 仿真反馈变化结果显示：仿真结果与基本发展规律一致(五) 场猪尿产沼气年变化量R5(t)仿真流率基本入树T5(t)(1)逐枝建立的前期场猪尿产沼气年变化量R5(t)( mA3/年)流率基本入树T5(t)见图3.16图3.14场猪尿产沼气年变化量R5(t)( mA3/年)流率基本入树T5(t)(2) 变量发展方程：略(3) T1(t)，T2(t)，T3(t),T4(t)，T5(t)组合仿真检验打者尿招气促进因子M33 (t)1 )场猪尿沼气促进因子M33(t)的子因果链图见图3.17期猪尿年产招气L5 (t)(m3) I图3.17前期反馈参数的反馈函
