MATLAB之数学建模：深圳市生活垃圾处理社会总成本分析
MATLAB之数学建模：深圳市生活垃圾处理社会总成本分析
注：MATLAB版本--2016a,作图分析部分见《MATLAB之折线图、柱状图、饼图以及常用绘图技巧》
一.现状模式下的模型
%第一题：建立总成本分析模型/年：按现状分析
% 总成本=直接成本 +经济技术成本 + 社会成本
function dataPro = Total_Cost_Analysis(year)
%垃圾每年预测表：2017-2030
table = [ 6.4450e+06 6.8317e+06 7.2416e+06 7.6761e+06 7.9832e+06 8.3025e+06 8.6346e+06 8.9800e+06 9.3392e+06 9.6193e+06 9.9079e+06 1.0205e+07 1.0511e+07 1.0827e+07]; %垃圾总量每年数值（2017-2030）
rubbish_quantity = table(year-2016);
%将时间分期处理：2017-2020，2021-2025，2026-2030
switch year
case { 2017,2018,2019,2020} %近期
rubbish_num_burn=215*10^4;
rubbish_num_landfill = rubbish_quantity-rubbish_num_burn;
class_cost = 0;
handle_cost = rubbish_num_landfill*60+ rubbish_num_burn*100;
transport_cost = 0.5*rubbish_num_landfill*60+0.5*rubbish_num_landfill*70+...
0.5*rubbish_num_burn*60+0.5*rubbish_num_burn*70 ;
social_cost =132*rubbish_quantity;
technology_cost=1300*10^4; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 100*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity/102.49*10^8 ;
case {2021,2022,2023,2024,2025} %中期
rubbish_num_burn =215*10^4;
rubbish_num_landfill = rubbish_quantity-rubbish_num_burn;
class_cost = 0;
handle_cost = rubbish_num_landfill*60+ rubbish_num_burn*150;
transport_cost = rubbish_quantity*100;
social_cost =2*132*rubbish_quantity;
technology_cost=1300*10^4; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 50*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity*(1/104.15+1/100.1+1/317.46)*10^8 +rubbish_quantity/1.54*10^4; %不同时间 定值
case {2026,2027,2028,2029,2030} %远期
rubbish_num_burn =215*10^4;
rubbish_num_landfill = rubbish_quantity-rubbish_num_burn;
class_cost = 0;
handle_cost = rubbish_num_landfill*60+ rubbish_num_burn*180;
transport_cost = rubbish_quantity*100;
social_cost =2*132*rubbish_quantity;
technology_cost=1300*10^4; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 0*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity*(1/104.15+1/50.05+1/222.22)*10^8 +rubbish_quantity/1.54*10^4 ; %不同时间 定值
otherwise
msgbox('亲，请重新输入年份：');
end
%设施投资:
equipment_cost = 1.56*10^8;
%输出，分析:dataPro为数据集合
profit =profit *0.15;
direct_cost = class_cost + transport_cost + equipment_cost + handle_cost;
total_cost = direct_cost+technology_cost +social_cost+subsidy-profit ;
%dataPro（11）: 分类，收运，设施，处理，技术，社会，补贴，收益，直接，总，均
dataPro = [ class_cost,transport_cost,equipment_cost,handle_cost, ...
technology_cost,social_cost,subsidy,profit,direct_cost,total_cost,total_cost/rubbish_quantity];
end
二. 模式一
%模式一：总成本=直接成本 +经济技术成本 + 社会成本
function dataPro = Total_Cost_Analysis_model1(year)
%垃圾每年预测表：2017-2030
table = [ 6.4450e+06 6.8317e+06 7.2416e+06 7.6761e+06 7.9832e+06 8.3025e+06 8.6346e+06 8.9800e+06 9.3392e+06 9.6193e+06 9.9079e+06 1.0205e+07 1.0511e+07 1.0827e+07]; %垃圾总量每年数值（2017-2030）
rubbish_quantity = table(year-2016);
%将时间分期处理：2014-2020，2021-2025，2026-2030
switch year
case {2016,2017,2018,2019,2020}
rubbish_num_burn=215*10^4; %近期
rubbish_num_landfill = rubbish_quantity-rubbish_num_burn;
transport_cost = 0.5*rubbish_quantity*60+0.5*rubbish_quantity*70;
handle_cost = rubbish_num_landfill*60+rubbish_num_burn*100;
social_cost = 132*rubbish_quantity;
technology_cost=1300*10^4; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 100*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity/102.49*10^8 ; %不同时间 定值
case {2021,2022,2023,2024,2025}
transport_cost = rubbish_quantity*100;
handle_cost = rubbish_quantity*150;
social_cost =8*132*rubbish_quantity;
technology_cost=1300*10^4; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 50*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))* rubbish_quantity*(1/52.1+1/100.1 )*10^8 +rubbish_quantity/0.77*10^4;
case {2026,2027,2028,2029,2030}
transport_cost = rubbish_quantity*100;
handle_cost = rubbish_quantity*180;
social_cost =8*132*rubbish_quantity;
technology_cost=1300*10^4; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 0*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity*(1/52.07+1/50.05 )*10^8 +rubbish_quantity/0.77*10^4;
otherwise
msgbox('亲，请重新输入年份：');
end
%分类费用
class_cost = 0;
%设施投资:
equipment_cost = 0 ;
%输出，分析
profit =profit *0.15;
direct_cost = class_cost + transport_cost + equipment_cost + handle_cost;
total_cost = direct_cost+technology_cost +social_cost+subsidy-profit ;
%dataPro（11）: 分类，收运，设施，处理，技术，社会，补贴，收益，直接，总，均
dataPro = [ class_cost,transport_cost,equipment_cost,handle_cost, ...
technology_cost,social_cost,subsidy,profit,direct_cost,total_cost,total_cost/rubbish_quantity];
end
三. 模式二
%模式二：源头分类收集+湿垃圾生物处理+干垃圾焚烧+中心城区干垃圾转运
function dataPro = Total_Cost_Analysis_model2(year)
%垃圾每年预测表：2017-2030
table = [ 6.4450e+06 6.8317e+06 7.2416e+06 7.6761e+06 7.9832e+06 8.3025e+06 8.6346e+06 8.9800e+06 9.3392e+06 9.6193e+06 9.9079e+06 1.0205e+07 1.0511e+07 1.0827e+07]; %垃圾总量每年数值（2017-2030）
rubbish_quantity = table(year-2016);
switch year
case {2016,2017,2018,2019,2020}
rubbish_num_burn=215*10^4; %近期
rubbish_num_landfill = rubbish_quantity-rubbish_num_burn;
class_cost = 0;
transport_cost = 0.5*rubbish_quantity*60+0.5*rubbish_quantity*70;
handle_cost = rubbish_num_landfill*60+rubbish_num_burn*100;
social_cost =132*rubbish_quantity;
technology_cost=1300*10^4+10^8; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 100*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity/102.49*10^8 ; %不同时间 定值
case {2021,2022,2023,2024,2025}
class_cost = 10.6*10^8;
transport_cost = 0.4*rubbish_quantity*60+0.6*rubbish_quantity*100;
handle_cost = rubbish_quantity*150;
social_cost =1.2*132*rubbish_quantity; % (year-2017)
technology_cost=1300*10^4+0.7*10^8; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 50*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity*(1/72.31+1/396.83+1/100.1)*10^8 +rubbish_quantity/1.28*10^4 ; %不同时间 定值
case {2026,2027,2028,2029,2030}
class_cost = 10.6*10^8;
transport_cost = 0.5*rubbish_quantity*60+0.5*rubbish_quantity*100;
handle_cost = 0.5*rubbish_quantity*180+ 0.5*rubbish_quantity*200;
social_cost =1.2*132*rubbish_quantity; % (year-2017)
technology_cost=1300*10^4+0.4*10^8; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 0*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))* rubbish_quantity*(1/77.98+1/277.78+1/50.05)*10^8 +rubbish_quantity/1.28*10^4;
otherwise
msgbox('亲，请重新输入年份：');
end
%设施投资:
equipment_cost = 0;
%输出，分析
profit =profit *0.15;
direct_cost = class_cost + transport_cost + equipment_cost + handle_cost;
total_cost = direct_cost+technology_cost +social_cost+subsidy-profit ;
%dataPro（11）: 分类，收运，设施，处理，技术，社会，补贴，收益，直接，总，均
dataPro = [ class_cost,transport_cost,equipment_cost,handle_cost, ...
technology_cost,social_cost,subsidy,profit,direct_cost,total_cost,total_cost/rubbish_quantity];
end
四. 模式三
%模式三：混合收集+末端分类+湿垃圾生物处理+干垃圾焚烧+中心城区干垃圾转运
function dataPro = Total_Cost_Analysis_model3(year)
%垃圾每年预测表：2017-2030
table = [ 6.4450e+06 6.8317e+06 7.2416e+06 7.6761e+06 7.9832e+06 8.3025e+06 8.6346e+06 8.9800e+06 9.3392e+06 9.6193e+06 9.9079e+06 1.0205e+07 1.0511e+07 1.0827e+07]; %垃圾总量每年数值（2017-2030）
rubbish_quantity = table(year-2016);
switch year
case {2016,2017,2018,2019,2020}
rubbish_num_burn=215*10^4; %近期
rubbish_num_landfill = rubbish_quantity-rubbish_num_burn;
transport_cost = 0.5*rubbish_quantity *60+0.5*rubbish_quantity *70 ;
handle_cost = rubbish_num_landfill*60+rubbish_num_burn*100;
social_cost =132*rubbish_quantity;
technology_cost=1300*10^4+10^8; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 100*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity/102.49*10^8 ; %不同时间 定值
case {2021,2022,2023,2024,2025}
handle_cost = 0.5*rubbish_quantity*150+ 0.5*rubbish_quantity*200;
transport_cost = 0.4*rubbish_quantity*60+0.6*rubbish_quantity*100;
social_cost =132*rubbish_quantity;
technology_cost=1300*10^4+0.7*10^8; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 50*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity*(1/86.87+1/317.46+1/100.1)*10^8 +rubbish_quantity/1.54*10^4; %不同时间 定值
case {2026,2027,2028,2029,2030}
handle_cost = 0.5*rubbish_quantity*180+ 0.5*rubbish_quantity*200;
transport_cost = 0.4*rubbish_quantity*60+0.6*rubbish_quantity*100;
social_cost =132*rubbish_quantity;
technology_cost=1300*10^4+0.4*10^8; % 湿处理分期：10^8,0.7*10^8，0.4*10^8
subsidy = 0*rubbish_quantity; %前期 100，中期50，后期取消，成本计算取负
profit = (10^(-4))*rubbish_quantity*(1/86.80+1/222.22+1/50.05)*10^8 +rubbish_quantity/1.54*10^4; %不同时间 定值
otherwise
msgbox('亲，请重新输入年份：');
end
%设施投资:
equipment_cost = 0;
%分类
class_cost = 0;
%输出，分析
profit =profit *0.15;
direct_cost = class_cost + transport_cost + equipment_cost + handle_cost;
total_cost = direct_cost+technology_cost +social_cost+subsidy-profit ;
%dataPro（11）: 分类，收运，设施，处理，技术，社会，补贴，收益，直接，总，均
dataPro = [ class_cost,transport_cost,equipment_cost,handle_cost, ...
technology_cost,social_cost,subsidy,profit,direct_cost,total_cost,total_cost/rubbish_quantity];
end
五. 垃圾总量预测
%垃圾总量预测
rubbish_table = zeros(1,14);
rubbish_table(1,1) = 541.14*10^4; %2014年垃圾产量541.14万吨
for year = 2015:2020
rubbish_table(1,year-2013) = rubbish_table (1,year-2014)*(1+0.06);
end
for year = 2021:2025
rubbish_table(1,year-2013) = rubbish_table (1,year-2014)*(1+0.04);
end
for year = 2026:2030
rubbish_table(1,year-2013) = rubbish_table (1,year-2014)*(1+0.03);
end
六.各模式数据汇总
% 总成本=直接成本 +经济技术成本 + 社会成本
%数据收集data_model（从2017-2030年：现状，模式一，模式二，模式三）
clear;close;clc;
data_model0 = zeros(14, 11);
data_model1 = zeros(14, 11);
data_model2 = zeros(14, 11);
data_model3 = zeros(14, 11);
for year = 2017 : 2030
dataPro0 = Total_Cost_Analysis(year ); %现状
dataPro1 = Total_Cost_Analysis_model1(year );
dataPro2 = Total_Cost_Analysis_model2(year );
dataPro3 = Total_Cost_Analysis_model3(year );
for i = 1:11
data_model0(year-2016,i) = dataPro0(i);
data_model1(year-2016,i) = dataPro1(i);
data_model2(year-2016,i) = dataPro2(i);
data_model3(year-2016,i) = dataPro3(i);
end
end
七.最优模式评选
%优选模式计算：分类0.3，设施0.5，收运1.5，处理1，技术1，社会1.7，收益-2
%原理较复杂，优选模式以远期成本最优，并且设定不同成本的比重，所得结果为每吨垃圾的成本
%
table = [9.6193e+06 9.9079e+06 1.0205e+07 1.0511e+07 1.0827e+07]; %垃圾总量每年数值（远期2025-2030）
rubbish_quantity = sum(table,2);
x=0.3*sum(data_model0(10:14,1))+0.5*sum(data_model0(10:14,3))+1.5*sum(data_model0(10:14,2))+...
1*sum(data_model0(10:14,4))+1*sum(data_model0(10:14,5))+1.7*sum(data_model0(10:14,6))+...
2*sum(data_model0(10:14,8));
table0 =x/rubbish_quantity %现状模式 /rubbish_quantity
x=0.3*sum(data_model1(10:14,1))+0.5*sum(data_model1(10:14,3))+1.5*sum(data_model1(10:14,2))+...
1*sum(data_model1(10:14,4))+1*sum(data_model1(10:14,5))+1.7*sum(data_model1(10:14,6))+...
2*sum(data_model1(10:14,8));
table1 = x/rubbish_quantity %模式一
x=0.3*sum(data_model2(10:14,1))+0.5*sum(data_model2(10:14,3))+1.5*sum(data_model2(10:14,2))+...
1*sum(data_model2(10:14,4))+1*sum(data_model2(10:14,5))+1.7*sum(data_model2(10:14,6))+...
2*sum(data_model2(10:14,8));
table2 = x/rubbish_quantity %模式二
x=0.3*sum(data_model3(10:14,1))+0.5*sum(data_model3(10:14,3))+1.5*sum(data_model3(10:14,2))+...
1*sum(data_model3(10:14,4))+1*sum(data_model3(10:14,5))+1.7*sum(data_model3(10:14,6))+...
2*sum(data_model3(10:14,8));
table3 = x/rubbish_quantity %模式三
MATLAB之数学建模：深圳市生活垃圾处理社会总成本分析的更多相关文章
- Matlab与数学建模
一.学习目标. (1)了解Matlab与数学建模竞赛的关系. (2)掌握Matlab数学建模的第一个小实例—评估股票价值与风险. (3)掌握Matlab数学建模的回归算法. 二.实例演练. 1.谈谈你 ...
- 卓金武《MATLAB在数学建模中的应用》 第2版
内容介绍 本书的作者都具有实际的数学建模参赛经历和竞赛指导经验.书中内容完全是根据数学建模竞赛的需要而编排的,涵盖了绝大部分数学建模问题的matlab求解方法.本书内容分上下两篇.上篇介绍数学建模中常 ...
- “GANs”与“ODEs”：数学建模的终结？
在本文中,我想将经典数学建模和机器学习之间建立联系,它们以完全不同的方式模拟身边的对象和过程.虽然数学家基于他们的专业知识和对世界的理解来创建模型,而机器学习算法以某种隐蔽的不完全理解的方式描述世界, ...
- python 版 mldivide matlab 反除（左除）《数学建模算法与程序》Python笔记
今天在阅读数学建模的时候看到了差分那章 其中有一个用matlab求线性的代码,这里我贴出来 这里我送上 Python代码 In [39]: import numpy as np ...: from s ...
- 在数学建模中学MATLAB
为期三周的数学建模国赛培训昨天正式结束了,还是有一定的收获的,尤其是在MATLAB的使用上. 1. 一些MATLAB的基础性东西: 元胞数组的使用:http://blog.csdn.net/z1137 ...
- 数学建模--matlab基础知识
虽然python也能做数据分析,不过参加数学建模,咱还是用专业的 1. Matlab-入门篇:Hello world! 程序员入门第一式: disp(‘hello world!’) 2. 基本运算 先 ...
- 余胜威《MATLAB数学建模经典案例实战》2015年版
内容介绍 本书全面.系统地讲解了数学建模的知识.书中结合历年全国大学生数学建模竞赛试题,采用案例与算法程序相结合的方法,循序渐进,逐步引导读者深入挖掘实际问题背后的数学问题及求解方法.在本书案例的分析 ...
- Matlab 多个版本的安装包下载、安装和激活教程 + 多套数学建模视频教程
目录 1. 关键词 1.1. 说明 2. 下载地址 2.1. OneDrive高速云盘 2.1.1. 多版本的安装包 2.1.2. 多套数学建模的视频教程 2.2. 百度云 3. 安装教程 1. 关键 ...
- 数学建模学习笔记 | matlab基本命令及用法
前言 数学建模对matlab水平的要求 了解matlab的基本用法,如常用命令.脚本结构.矩阵的基本操作.绘图等: 熟悉matlab的程序结构,能创建和引用函数: 熟悉常见模型的求解算法和套路: 自主 ...
随机推荐
- Flask17 Flask_Script插件的使用
1 什么是Flask_Script 可以对flask应用做一些统一的操作 flask_script官网:点击前往 2 安装flask_script pip install -i https://pyp ...
- Luogu 2900 [USACO08MAR]土地征用Land Acquisition
斜率优化dp. 首先发现如果存在$x$和$y$使得$len(x) \geq len(y)$并且$wid(x) \geq wid(y)$,那么$y$直接不考虑就好了,因为在买$x$的时候就把$y$顺便带 ...
- 内核启动流程2-C语言部分的最后一个函数init_post()
最后分析最终调用用户空间init进程的函数init_post(). static noinline int init_post(void)这是一个非_init函数.强制让它为非内联函数,以防gcc让它 ...
- hdu1070
#include <stdio.h>#include <string.h>struct milk{ char brand[128]; int price; int volume ...
- Torando 入门
1. 前言 Tornado 是使用 Python 编写的一个强大的.可拓展性的 Web 服务器/框架.与其他主流 Web 服务器框架有着明显区别:Tornado 支持异步非阻塞框架.同时它处理速度非常 ...
- iphone5 A1429国行IOS8.4.1 越狱 完美使用电信3G
国航 8.4.1,越狱,使用电信3G,能打电话,能发短信,正常上网使.(越狱降级方式,请参照本文末端连接) 有好多人说不管怎么试都是无服务,请查看 设置-蜂窝移动网络-漫游里面,必须只有 语音漫游 这 ...
- Sharding-JDBC 使用入门和基本配置
一.什么是Sharding-JDBC Sharding-JDBC定位为轻量级Java框架,在Java的JDBC层提供的额外服务.它使用客户端直连数据库,以jar包形式提供服务,无需额外部署和依赖,可理 ...
- scrapy-redis源码解读之发送POST请求
1 引言 这段时间在研究美团爬虫,用的是scrapy-redis分布式爬虫框架,奈何scrapy-redis与scrapy框架不同,默认只发送GET请求,换句话说,不能直接发送POST请求,而美团的数 ...
- [CentOS7] 设置开机启动方式(图形界面或命令行)
由于CenOS之前一直都是通过修改inittab文件来修改开机启动模式,于是 通过 vim /etc/inittab 打开inittab来查看 如上所示,CentOS 7由于使用systemd而不是i ...
- elasticsearch shard 和 replica
(1)index包含多个shard(2)每个shard都是一个最小工作单元,承载部分数据,lucene实例,完整的建立索引和处理请求的能力(3)增减节点时,shard会自动在nodes中负载均衡(4) ...