环评报告中水、气等的计算方法统计
环评报告中水、气等的计算方法统计
◇生活水量统计
◇生活用水量：城镇居民生活用水量按150 l/人〃d计算；
农村居民生活用水量按60 l/人〃d计算；
办公楼、商场生活用水量按50 l/人〃d计，
无具体人数时按50 l/50m2〃d（具体酌情考虑）
餐饮：500 l/100m2〃d
◇排水量：排水量按生活用水量80%计（有些朋友建议说按60％－75％计算）。◇废气
◇用煤量计算：
①有建筑面积时，一个采暖期（5个月）40Kg/m2计；
②无建筑面积时，常压茶浴炉：无烟煤用煤量按每吨锅炉5t/月；（上限，可酌情减量）
卧式热水、蒸汽锅炉：每吨锅炉30t/月，360t/年；
◇常用煤质分析
砟子煤煤质成份分析表
矿井名称地点灰分（％）含水分（％）含硫分（％）低位发热量大卡/kg
磁窑堡灵武12.559.990.88
5500
无烟煤煤质成份分析表
矿井名称地点灰份（％）含水分（％）含硫分（％）
低位发热
大卡/k 汝箕沟石嘴山6-8 2.00.257500
◇锅炉烟气计算： 1 Kg煤产生8.5m3烟气（经验）；
链条炉1 Kg煤产生10.35m3烟气；
◇源强计算公式：
SO2：Q=G×2×0.8×S×(1－ηso2)；
式中：G－耗煤量；S－煤中含硫量；
ηso2－脱硫效率，麻石水浴、水膜按40%计，石灰石－石膏法按90%。
烟尘：Q=G×A×ηA×(1－η)
式中：G－耗煤量；
A－煤中灰分；
ηA－煤中飞灰占灰分的百分比；其值与燃烧方式有关，见下表；
η－除尘效率；麻石水浴按95%，陶瓷多管90%，旋风80%，布袋99%;
煤中飞灰占灰分的百分比（环境统计手册）
炉型ηA炉型ηA 手烧炉15～25沸腾炉40～50链条炉15～25煤粉炉75～85往复炉20振动炉排20～40抛煤机炉25～40
燃煤锅炉烟尘初始排放浓度限值（GB13721—2001）
锅炉类别燃煤收到基灰分
（%）
烟尘初始排放浓度(mg/m
Ⅰ时段Ⅱ时
层然锅炉
自然通风锅炉
（＜0.7MW，1t/h＝
/15012
其他锅炉
（≤2.8MW，4t/h）
Aar≤25%180016
Aar＞25%200018其他锅炉
（＞2.8MW，4t/h）
Aar≤25%200018
Aar＞25%220020
沸腾锅炉
循环流化床锅炉/15 000150其他沸腾锅炉/20 000180抛煤机锅炉/500050
注：锅炉烟尘初始排放浓度达到表中所列初始浓度，锅炉为合格锅炉，否则，为
不合格锅炉。
（出自锅炉大气污染物排放标准）
◇锅炉补水：
一般锅炉，正常运行时，有热交换器(大型锅炉房)，10t锅炉补水250 Kg/d，
蒸汽锅炉：10t热水锅炉每天补水量2t；
水处理水量：锅炉水力冲渣冷却水量：
◇锅炉灰渣排放量计算：
无烟煤按耗煤量20%；砟子煤按耗煤量15%；
◇天然气用量计算：（自供热办胡科长）
有建筑面积时，一个采暖期（5个月）12～18m3/m2计（节能住宅区下限，非节
能住宅取上限）；一般按15m3/m2计。
天然气主要成份参数表
成份CH4C2H6C3H8CO2H2O 指标（%）95.950.910.14 3.0062ppm
◇生活垃圾：
城镇居民生活垃圾产生量按1Kg/人〃d计算；
农村居民生活垃圾产生量按0.8Kg/人〃d计算；
办公楼、商场生活垃圾量按0.5Kg/人〃d计，无具体人数时按0.5Kg /50m2〃d；
餐饮垃圾：10Kg/100m2〃d
◇噪声：（自环境上岗培训教材）
点声源计算公式：△L=20 lg r1/r2;(△L－衰减量)
式中：r－点声源至受声点的距离。
点声源距离衰减值表
距离（米）△L dB（A）距离（米）△L dB（A）距离（米）△L 5144032100
10205034200
1523.560353004 20267037400
25288038500
3029.59039
线声源计算公式：
△L=10 lg1/2πrl；(△L－衰减量)
r－线声源至受声点的距离，m；
l线声源的长度，m。
当r/l＜1/10时，可视为无限长线声源，此时，在距离线声源r1～r2处的衰减
值为：
△L=10 lg r1/r2
当r2=2r1时，线声源传播距离增加一倍，衰减值3dB（A）。
◇绿化
绿化投资：按25元/m2计
绿化用水：按0.5 t/m2计。
绿化浇洒用水定额为1～3L/m2˙d，道路浇洒用水定额为2～3L/m2˙d。《建筑
给水排水设计规范》（GB50015－2003）
《计算方法》课内实验报告
《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 2011年10月26日
一、实验题目： 数值积分 二、实验目的： 1．熟悉matlab 编写及运行数值计算程序的方法。 2．进一步理解数值积分的基础理论。 3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。 三、实验内容： 1．分别用复合梯形求积公式及复合辛普森求积公式计算积分xdx x ln 10 ? , 要求计算精度达到410-,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分dx x x ?+3 021,使误差不超过510-. 3．用3=n 的高斯-勒让德公式计算积分?3 1 sin x e x ,给出计算结果. 4．用辛普森公式(取2==M N ) 计算二重积分.5 .00 5 .00 dydx e x y ? ? - 四、实验结果： 1.（1）复合梯形法: 将区间[a,b]划分为n 等份，分点n k n a b h kh a x k ,2,1,0,,=-=+=在每个区间[1,+k k x x ](k=0,1,2,···n-1)上采用梯形公式，则得 )()]()([2)()(1 11 1 f R x f x f h dx x f dx x f I n n k k k b a n k x x k k ++===∑?∑? -=+-=+ 故)]()(2)([21 1 b f x f a f h T n k k n ++=∑-=称为复合梯形公式 计算步长和划分的区间 Eps=1E-4 h1=sqrt(Eps/abs(-(1-0)/12*1/(2+1))) h1 =0.0600 N1=ceil(1/h1) N1 =17 用复合梯形需要计算17个结点。 复合梯形： function T=trap(f,a,b,n) h=(b-a)/n;
数值分析实验报告1
实验一误差分析 实验1.1（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 其中ε（1.1）和（1.221,,,a a 的输出b ”和“poly ε。 （1（2 （3）写成展 关于α solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。 实验过程： 程序： a=poly(1:20); rr=roots(a); forn=2:21 n form=1:9 ess=10^(-6-m);
ve=zeros(1,21); ve(n)=ess; r=roots(a+ve); -6-m s=max(abs(r-rr)) end end 利用符号函数：（思考题一）a=poly(1:20); y=poly2sym(a); rr=solve(y) n
很容易的得出对一个多次的代数多项式的其中某一项进行很小的扰动，对其多项式的根会有一定的扰动的，所以对于这类病态问题可以借助于MATLAB来进行问题的分析。 学号：06450210 姓名：万轩 实验二插值法
太原理工大学数值计算方法实验报告
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日
y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"
心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。
污染物排放j计算方法
工业污染物排放统计方法 一、工业污染物估算常用方法 工业企业环境统计工作中对废气、废水和固体废物及所含污染物产生量、排放量的计算通常采用三种方法，即实测法、物料衡算法和产排污系数法。 1、实测法 实测法是通过监测手段或国家有关部门认定的连续计量设施，测量废气、废水的流速、流量和废气、废水中污染物的浓度，用环保部门认可的测量数据来计算各种污染物的产生量和排放总量的统计计算方法。 G=KC i Q 式中：G——污染物产生量或排放量； Q——介质流量； C i——介质中i污染物浓度； K——单位换算系数。 浓度和流量的单位不一致时，单位换算系数K取不同的值。废水中污染物的浓度单位常取mg/L，系数K取10-3；废气中污染物的浓度一般取mg/L，系数K取10-6。 实测法的基础数据主要来自于环境监测站。监测数据是通过科学、合理地采集样品、分析样品而获得的。监测采集的样品是对监测的环境要素的总体而言，如采集的样品缺乏代表性，尽管测试分析很准确，不具备代表性的数据也毫无意义。 因受现有监测技术和监测条件的约束，实测法有一定的局限性。这主要是目前除了重点污染源有比较准确的监测数据外，其他多数非重点污染源不能得到有效的监测；而且很多重点污染源还未实现连续监测，监测结果的代表性有待提高。 例某炼油厂年排废水2万t，废水中废油浓度C油为500mg/L，COD浓度C COD为300mg/L，水未处理直接排放。计算该厂废油和COD的年排放量。 解：G油＝K C油Q ＝10-6×500×2×104 ＝10（t） G COD＝K C COD Q ＝10-6×300×2×104 ＝6（t） 例某冶炼厂排气筒截面0.4m2，排气平均流速12.5m/s，实测所排废气中SO2平均浓度12mg/m3，粉尘浓度8mg/L计算该排气筒每小时SO2和粉尘的排放量。 解：每小时废气流量Q=12.5×0.4×3600 = 1.8×104（m3/h） 每小时SO2排放量Gso2 = 10—6×12×1.8×104 = 0.216（kg/h) 每小时粉尘排放量G粉尘= 10—6×8×1.8×104 = 0.144（（kg/h) 2、物料衡算法 物料衡算法是指根据物质质量守恒原理，对生产过程中使用的物料变化情况进行定量分析的一种方法。即： 投入物料量总和＝产出物料量总和 ＝主副产品和回收及综合利用的物质量总和＋排出系统外的废物质量这里的排出系统外的废物质量包括可控制与不可控制生产性废物及工艺过程的泄漏等物料流失。
水环境容量计算方法
水环境容量计算方法 中国环境规划院李云生 2004.5 ?基本涵义 ?计算模型 ?计算步骤 ?校核方法 第一部分水环境容量的基本涵义 容量涵义 技术指南中的概念定义 ?在给定水域范围和水文条件，规定排污方式和水质目标的前提下，单位时间内该水域最大允许纳污量，称作水环境容量。 ?从上述定义可知，水环境容量主要决定于三个要素：水资源量、水环境功能区划和排污方式。 要素之一：水资源量 ?从某种意义上讲，水资源量是水环境容量基础； ?为了确保用水安全，水环境容量计算采用的是较高保证率的水文设计条件； ?并不是所有的水资源量都用来计算环境容量。 要素之二：水环境功能区 ?水环境功能区划体现人们对水环境质量的需求，反映了人们对水资源的态度：开发、利用或保护。 ?已划分水环境功能区的水域，要从时间、空间两个方面规范功能区达标标准； ?未划分水环境功能区的水域可不进行容量计算；若考虑计算，按较高功能标准进行（II类）。 要素之三：排污方式 ?排污口沿河（或其他水体）位置布设，对河流整体水环境容量影响较大； ?排污口排放方式（岸边或中心，浅水或深水），对局部的污染物稀释混合影响很大； ? ? 第二部分水环境容量的计算模型 ?1、流域概化模型 ?2、水动力学模型 ?3、污染源概化模型 ?4、水质模型 1、流域概化 ?将天然水域（河流、湖泊水库）概化成计算水域，例如天然河道可概化成顺直河道，复杂的河道地形可进行简化处理，非稳态水流可简化为稳态水流等。水域概化的结果，就是能够利用简单的数学模型来描述水质变化规律。同时，支流、排污口、取水口等影响水环境的因素也要进行相应概化。若排污口距离较近，可把多个排污口简化成集中的排污口。 2、水动力学模型 ?最枯月设计条件
c 计算器实验报告
简单计算器 姓名: 周吉祥 实验目的：模仿日常生活中所用的计算器，自行设计一个简单的计算器程序，实现简单的计算功能。 实验内容： (1)体系设计： 程序是一个简单的计算器，能正确输入数据，能实现加、减、乘、除等算术运算，运算结果能正确显示，可以清楚数据等。 (2)设计思路： 1)先在Visual C++ 6.0中建立一个MFC工程文件，名为 calculator. 2)在对话框中添加适当的编辑框、按钮、静态文件、复选框和单 选框 3)设计按钮，并修改其相应的ID与Caption. 4)选择和设置各控件的单击鼠标事件。 5)为编辑框添加double类型的关联变量m_edit1. 6)在calculatorDlg.h中添加math.h头文件，然后添加public成 员。 7)打开calculatorDlg.cpp文件，在构造函数中，进行成员初始 化和完善各控件的响应函数代码。 (3)程序清单：
●添加的public成员： double tempvalue; //存储中间变量 double result; //存储显示结果的值 int sort; //判断后面是何种运算：1.加法2.减法3. 乘法 4.除法 int append; //判断后面是否添加数字 ●成员初始化： CCalculatorDlg::CCalculatorDlg(CWnd* pParent /*=NULL*/) : CDialog(CCalculatorDlg::IDD, pParent) { //{{AFX_DATA_INIT(CCalculatorDlg) m_edit1 = 0.0; //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME); tempvalue=0; result=0; sort=0; append=0; }
SPC常用计算方法
SPC常用计算方法 SPC基础知识及常用计算方法 SPC基础知识 一、 SPC定义： 1、 SPC——统计制程管制：是指一套自制程中去搜集资料，并加以统计分析，从分析中去发气掘制程的异常，立即采取修正行动，使制程恢复正常的方法。 也就是说：品质不应再依赖进料及出货的抽样检验，而应该采取在生产过程中，认良好的管理方法，未获得良好的品质。 2、良好品质，必须做到下面几点： ①变异性低 ②耐用度 ③吸引力 ④合理的价格 3、变异的来源：大概来自5个方面： ①机器②材料③方法④环境⑤作业人员 应先从机器，材料方法，环境找变异，最后考虑人。 4、 SPC不是一个观念，而是要行动的 步骤一、确立制程流程——首先制程程序要明确，依据制程程序给制造流程图，并依据流程图订定工程品质管理表。 步骤二、决定管制项目——如果把所有对品质有影响的项目不论大小，轻重缓急一律列入或把客户不很重视的特性一并管制时，徒增管制成本浪费资料且得不赏失，反之如果重要的项目未加以管制时，则不能满足设计者，后工程及客户的需求，则先去管制的意义。 步骤三、实施标准化——欲求制程管制首先即得要求制程安定，例如：在风浪很大的船上比赛乒乓球，试部能否确定谁技高一筹，帮制程作业的安定是最重要的先决条件，所以对于制程上影响产品口质的重要原因，应先建立作业标准，并透过教育训练使作业能经标准进行。 步骤四、制程能力调查——为了设计、生产、销售客户满意且愿意购买的产品，制造该产品的制程能力务必符合客户的要求。因此制程的能力不足时，必顺进行制程能力的改善，而且在制程能力充足后还必须能继续，所以在品质管理的系统中制程能力的掌握很重要。 步骤五、管制图运用——SPC的一个基本工具就是管制图，而管制图又分计量值管制图与计数值管制图。 步骤六、问题分析解决——制程能力调查与管制图是可筛提供问题的原因系由遇原因或非机遇原因所造成，但无法告知你确切的原因为何及如何解决决问题？解决问题？而问题的解决技巧，在于依据事实找出造成变异的确切原因，并提此对策加以改善，及如何防止再发生。 步骤七、制程之继续管制——经过前6个步骤，人制程能力符合客户的要求，且管制图上的点未出管制界限时，则可将此管制界限沿有作为制程之继续管制，但当制程条件如有变动时，如机器，材料，方法等产生异动时，则须回到步骤三，不可沿原先之管制界限。 SPC的应用步骤其流程图如下： Ca制程准确度 Cp制程精密度 Cpk制程能力指数 二、管制图的运用 管制图的种类又依数值资料是计量值或计数值者，划分为二大类即计量值管制图与计数值管制图，计量值管制图不但只告诉你制程有问题了，还可以告诉你制程在什么地方出了问题，是中心值产生了问题还是变异量产生了问题。而在计量值管制图应用不便或应用时，则可采用计数值
环境容量
1．面积法 游人容量的计算公式为： 瞬时容量＝空间面积/单位规模指标 日容量＝瞬时容量×日周转率 年容量＝日容量×年可游天数 计算结果见下表： （1）按风景名胜区各区分类面积计算 东湖风景名胜区游人容量计算表一 东湖风景名胜区游人容量计算表二 2．线路法 到规划期末（2020年），东湖风景名胜区的游览性道路总面积约238240平方米，按人均占有道路面积10平方米计，计算结果见下表： 按游览道路总面积计算： 东湖风景名胜区游人容量计算三
分析并满足该地区的生态允许标准、游览心理标准、功能技术标准等因素而确定。并应符合下列规定： 1.生态允许标准应符合表3.5.1的规定； 2.游人容量应由一次性游人容量、日游人容量、年游人容量三个层次表示。 (1) 一次性游人容量(亦称瞬时容量)，单位以“人/次”表示； (2)游人容量，单位以“人次/日”表示； (3)游人容量，单位以“人次/年”表示。 3.游人容量的计算方法宜分别采用：线路法、卡口法、面积法、综合平衡法，并将计算结果填入表3.5.1.1： 表3.5.1.1 游人容量计算一览表(1) 游览用地名称(2) 计算面积(m2) (3) 计算指标(m2/人) (4) 一次性容量(人/次) (5) 日周转率(次) (6) 日游人容量(人次/日) (7) 备注 4.游人容量计算宜采用下列指标：(1)线路法：以每个游人所占平均道路面积计，5-10m2/人。(2)面积法：以每个游人所占平均游览面积计。其中：主景景点：50-100m2/人(景点面积)；一般景点：100-100m2/人(景点面积)；浴场海域：10-20m2/人(海拔0~-2以内水面)；浴场沙滩：5-10m/人(海拔0~+2m以内沙滩)。
数值计算实验报告
(此文档为word格式，下载后您可任意编辑修改！) 2012级6班###(学号)计算机数值方法 实验报告成绩册 姓名:宋元台 学号： 成绩：
数值计算方法与算法实验报告 学期： 2014 至 2015 第 1 学期 2014年 12月1日课程名称: 数值计算方法与算法专业:信息与计算科学班级 12级5班 实验编号： 1实验项目Neton插值多项式指导教师：孙峪怀 姓名：宋元台学号：实验成绩： 一、实验目的及要求 实验目的： 掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。 实验要求： 1. 给出Newton插值算法 2. 用C语言实现算法 二、实验内容 三、实验步骤(该部分不够填写.请填写附页)
1.算法分析： 下面用伪码描述Newton插值多项式的算法： Step1 输入插值节点数n，插值点序列{x(i),f(i)},i=1,2,……,n,要计算的插值点x. Step2 形成差商表 for i=0 to n for j=n to i f(j)=((f(j)-f(j-1)(x(j)-x(j-1-i)); Step3 置初始值temp=1,newton=f(0) Step4 for i=1 to n temp=(x-x(i-1))*temp*由temp(k)=(x-x(k-1))*temp(k-1)形成 (x-x(0).....(x-x(i-1)* Newton=newton+temp*f(i); Step5 输出f（x）的近似数值newton(x)=newton. 2.用C语言实现算法的程序代码 #include
计算方法实验报告格式
计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格
形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 1．理解数值计算稳定性的概念． 2．了解数值计算方法的必要性． 3．体会数值计算的收敛性与收敛速度． 三、实验内容 计算dx x x I n n ? += 1 10 ，1,2,,10n = ． 四、算法描述 由 dx x x I n n ? += 1 10 ，知 dx x x I n n ?+=--101110，则
数值分析实验报告1
实验一 误差分析 实验（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（）中19x 的系数作一个小的扰动。我们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。 实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b =
的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve = ))20:1((ve poly roots + 上述简单的Matlab 程序便得到（）的全部根，程序中的“ess ”即是（）中的ε。 实验要求： （1）选择充分小的ess ，反复进行上述实验，记录结果的变化并分析它们。 如果扰动项的系数ε很小，我们自然感觉（）和（）的解应当相差很小。计算中你有什么出乎意料的发现表明有些解关于如此的扰动敏感性如何 （2）将方程（）中的扰动项改成18x ε或其它形式，实验中又有怎样的现象 出现 （3）（选作部分）请从理论上分析产生这一问题的根源。注意我们可以将 方程（）写成展开的形式， ) 3.1(0 ),(1920=+-= x x x p αα 同时将方程的解x 看成是系数α的函数，考察方程的某个解关于α的扰动是否敏感，与研究它关于α的导数的大小有何关系为什么你发现了什么现象，哪些根关于α的变化更敏感 思考题一：（上述实验的改进） 在上述实验中我们会发现用roots 函数求解多项式方程的精度不高，为此你可以考虑用符号函数solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。
环境统计6
第六次作业 1、分发统一的含铜0.100mg/L的样品到六个实验室，各实验室5次测定值如表1，试比较不同实验室之间是否存在显著性差异？ 表1 6个实验室测定结果比较 实验室铜测定值（mg/L） 10.0980.0990.0980.1000.099 20.0990.1010.0990.0980.097 30.1010.1010.1010.1010.102 40.1000.1000.0970.0970.095 50.0980.0940.1020.1000.100 60.0980.0940.0980.0980.098解：单因素方差分析 (1)H0: 不同实验室之间不存在差异 H A:. 不同实验室之间存在差异。 (2)确定显著水平α=0.05 (3)计算 进行F 检验，P=0.017 由此判断组间差异极显著 为了确定各个实验室之间的差异是否显著，需要进行多重比较。结果表明3与1、2、4、5、6差异显著。其余差异不显著。 2、用3种方法测定水中硫酸盐含量，结果如表2，问3种方法测定结果是否有显著性差别？
表2 3种方法测定水中硫酸盐含量 甲法乙法丙法 279229210 334274285 303310117 378 198 解：组内观测次数不相等的方差分析 (1)H0: 3种方法测定结果没有显著性差别 H A: 3种方法测定结果有显著性差别 (2)确定显著水平α=0.05 (3)计算 进行F 检验，P=0.212 组间差异无显著性差异 3、某地区通过大量饮用水源调查得知，压力井细菌总数合格率为63%，先抽查压水井水样80份，细菌总数合格的58份，合格率为72.5%。问这批抽查水样的合格率与大量调查的合格率有无显著性差别？ 解：进行适合性检验 (1)H0: 这批抽查水样的合格率与大量调查的合格率无差别。 H A:.这批抽查水样的合格率与大量调查的合格率有差别。 (2)确定显著水平α=0.05 (3)计算
计算方法实验报告 拟合
南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用matlab编制以函数为基的多项式最小二乘拟合程序，并用于对下列数据作三次多项式最小二乘拟合(取权函数wi=1) x -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 y -2.30 -1 -0.14 -0.25 0.61 1.03 1.75 2.75 4.42 6.94 给定直线方程为：y=1/4*x3+1/2*x2+x+1 三、实验结论： 最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。 一般地。当测量数据的散布图无明显的规律时，习惯上取n次代数多项式。 程序运行结果为： a = 0.9731 1.1023 0.4862 0.2238 即拟合的三次方程为：y=0.9731+1.1023x+0.4862*x2+0.2238*x3
-2.5 -2-1.5-1-0.5 00.51 1.52 2.5 -4-20246 81012 x 轴 y 轴 拟合图 离散点 y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3 结论： 一般情况下，拟合函数使得所有的残差为零是不可能的。由图形可以看出最小二乘解决了残差的正负相互抵消的问题，使得拟合函数更加密合实验数据。 优点：曲线拟合是使拟合函数和一系列的离散点与观测值的偏差平方和达到最小。 缺点：由于计算方法简单，若要保证数据的精确度，需要大量的数据代入计算。
环境影响评价 常用计算系数
环境影响评价必须掌握的方法 计算系数 烧一吨煤，产生1600×S%千克SO2，1万立方米废气；产生200千克烟尘。 烧一吨柴油，排放2000×S％千克SO2，1.2万立米废气；排放1千克烟尘。 烧一吨重油，排放2000×S％千克SO2，1.6万立米废气；排放2千克烟尘。 大电厂，烟尘治理好，去除率超98%，烧一吨煤，排放烟尘3-5千克。 普通企业，有治理设施的，烧一吨煤，排放烟尘10-15千克； 砖瓦生产，每万块产品排放40-80千克烟尘；12-18千克二氧化硫。 规模水泥厂，每吨水泥产品排放3-7千克粉尘；1千克二氧化硫。 乡镇小水泥厂，每吨水泥产品排放12-20千克粉尘；1千克二氧化硫。 【物料衡算公式】 1吨煤炭燃烧时产生的SO2量＝1600×S千克；S含硫率，一般0.6-1.5%。若燃煤的含硫率为1%，则烧1吨煤排放16公斤SO2 。 1吨燃油燃烧时产生的SO2量＝2000×S千克；S含硫率，一般重油 1.5-3%，柴油0.5-0.8%。若含硫率为2%，燃烧1吨油排放40公斤SO2 。 排污系数：燃烧一吨煤，排放0.9-1.2万标立方米燃烧废气，电厂可取小值，其他小厂可取大值。燃烧一吨油，排放1.2－1.6万标立方米废气，柴油取小值，重油取大值。 【城镇排水折算系数】 0.7~0.9，即用水量的70-90%。 【生活污水排放系数】采用本地区的实测系数。。 【生活污水中COD产生系数】60g/人.日。也可用本地区的实测系数。 【生活污水中氨氮产生系数】7g/人.日。也可用本地区的实测系数。使用系数进行计算时，人口数一般指城镇人口数；在外来较多的地区，可用常住人口数或加上外来人口数。 【生活及其他烟尘排放量】 按燃用民用型煤和原煤分别采用不同的系数计算： 民用型煤：每吨型煤排放1~2公斤烟尘 原煤：每吨原煤排放8~10公斤烟尘 【工业废气排放总量计算】 1.实测法 当废气排放量有实测值时，采用下式计算： Q年= Q时× B年/B时/10000 式中： Q年——全年废气排放量，万标m3/y； Q时——废气小时排放量，标m3/h；
环境统计主要计算方法
工业污染物排放统计方法 工业企业环境统计工作中对废气、废水和固体废物及所含污染物产生量、排放量的计算通常采用三种方法，即实测法、物料衡算法和产排污系数法。 1、实测法 实测法是通过监测手段或国家有关部门认定的连续计量设施，测量废气、废水的流速、流量和废气、废水中污染物的浓度，用环保部门认可的测量数据来计算各种污染物的产生量和排放总量的统计计算方法。 G=KC i Q 式中：G——污染物产生量或排放量； Q——介质流量； C i——介质中i污染物浓度； K——单位换算系数。 浓度和流量的单位不一致时，单位换算系数K取不同的值。废水中污染物的浓度单位常取mg/L，系数K取10-3；废气中污染物的浓度一般取mg/L，系数K取10-6。 实测法的基础数据主要来自于环境监测站。监测数据是通过科学、合理地采集样品、分析样品而获得的。监测采集的样品是对监测的环境要素的总体而言，如采集的样品缺乏代表性，尽管测试分析很准确，不具备代表性的数据也毫无意义。监测样品的代表性由以下环节来决定：（1）采样点的布设。应充分考虑采样点的代表性，满足概率随机性的要求，尽量减少主观误差。废水污染物的监测要求，一类污染物一律在各车间或车间处理设施排放口取样监测；二类污染物在企业各个废水排放口取样监测。 （2）采样时问和频率。应根据监测的目的及监测组分的时间变化而定。污染源的监测频率要求一年监测2～4次，每次间隔时间不得少于1个月；一般监测两次(在正常生产条件下)，上半年和下半年各监测一次。 （3）样品的完整性。数据的完整性取决于采集到的样品的完整性，只有对所有采样点采集到的全套样品进行监测分析，才能得到完整的监测数据。 （4）监测数据的可比性。要使监测数据具有可比性，常采用的办法是使用标准样品(又称标准物质)和国家认可的环境监测分析方法。使用国家级标准样品可以使监测结果在很大范围内准确可比，使用国家认可的环境监测分析方法可减少系统误差，增加监测数据之间的可比性。 因受现有监测技术和监测条件的约束，实测法有一定的局限性。这主要是目前除了重点污染源有比较准确的监测数据外，其他多数非重点污染源不能得到有效的监测；而且很多重点污染源还未实现连续监测，监测结果的代表性有待提高。 例某炼油厂年排废水2万t，废水中废油浓度C油为500mg/L，COD浓度C COD为300mg/L，水未处理直接排放。计算该厂废油和COD的年排放量。 解：G油＝K C油Q ＝10-6×500×2×104 ＝10（t） G COD＝K C COD Q ＝10-6×300×2×104 ＝6（t） 例某冶炼厂排气筒截面0.4m2，排气平均流速12.5m/s，实测所排废气中SO2平均浓度12mg/m3，
(完整版)哈工大-数值分析上机实验报告
实验报告一 题目：非线性方程求解 摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。 前言：（目的和意义） 掌握二分法与Newton法的基本原理和应用。 数学原理： 对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。 对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)<0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)<0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。重复运行计算，直至满足精度为止。这就是二分法的计算思想。
Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式 产生逼近解x*的迭代数列{x k}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。 程序设计： 本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x); y=-x*x-sin(x); 写成如上形式即可，下面给出主程序。 二分法源程序： clear %%%给定求解区间 b=1.5; a=0;
%%%误差 R=1; k=0;%迭代次数初值 while (R>5e-6) ; c=(a+b)/2; if f12(a)*f12(c)>0; a=c; else b=c; end R=b-a;%求出误差 k=k+1; end x=c%给出解 Newton法及改进的Newton法源程序：clear %%%% 输入函数 f=input('请输入需要求解函数>>','s') %%%求解f(x)的导数 df=diff(f);
环境统计指标体系编制说明.doc
“十二五”环境统计指标体系编制说明 2011年9月
一、背景情况 环境统计是环境管理的一项重要内容和基础工作。1980年11月，为了加强环境管理,掌握环境污染、治理情况和环境保护工作开展情况，为制定相关政策、编制规划和开展管理工作服务，国务院环境保护领导小组办公室在北京主持召开了全国第一次环境统计工作会议，针对我国县及县以上工业“三废”排放，及其治理情况和环保队伍自身建设、工作发展情况开展环境统计，这标志着我国环境统计工作的起步和环境统计报表制度的建立。 经过30多年的不断探索和完善，环境统计已经在为各级政府制定社会经济和环境保护方针政策，开展各项环境管理工作，向社会提供环境保护情况信息等方面都发挥了十分重要的作用。而且已经初步形成了较为系统、完整的部门统计体系。 现行的环境统计报表制度，包括环境统计综合报表、环境统计专业报表制度和环境统计年报工作技术要求等一系列文件。其中，环境统计综合报表和环境统计专业报表制度是核心内容。 环境统计综合年报制度的实施范围为有污染排放的工业企业、医院、城镇生活及其他排污单位、实施污染物集中处置的危险废物集中处置厂和城市污水处理厂等。环境统计综合年报制度由16张年报表、5张季报表、351项指标组成，主要统计内容为企业污染物排放和环境治理设施运行情况。 环境统计专业年报制度的实施范围是除核安全和外事（单独统
计）以外的所有环境保护业务部门的工作，包括环保系统能力建设，污染防治、生态保护、环境法制和监督执法，环境监测、科技标准、环境影响评价等环境管理工作情况。“十一五”环境统计专业报表制度由22张年报表、7张季报表、639项指标组成。 环境统计报表按年度每年初上报一次。其中综合报表由企业和相关单位填报，各级环保部门逐级汇总上报；专业报表由各级环保部门的主管业务部门统计填报并逐级上报。国家环境统计年报和公报在完成全国环境统计数据汇总之后，一般在六五世界环境日前发布。 随着环境管理工作的不断深入，尤其是“十一五”期间主要污染物总量减排工作的推动，环境统计数据质量引起各级领导和各相关方面更高的重视与关注，环境统计工作得到了明显的强化和提高。 二、修改的必要性 “十二五”期间我国将进一步推进主要污染物总量减排工作，各项环境管理工作也在不断深入和实行量化管理，对环境统计的要求也将越来越高，从而使环境统计面临着严峻的挑战和良好的发展机遇。 近年来，环境统计工作虽然不断进行积极的探索和必要的调整，但从总体上看，主要统计指标体系、技术体系和管理模式仍沿用上世纪九十年代甚至更早期的做法，而且由于长期以来环境统计能力建设和工作基础薄弱、方法相对落后，加之客观上存在统计对象种类多、差异大、变化快、数据获取难等具体困难，尤其是一部分企业和地方，为了应付各种检查和考核，人为调整数据，更使得数据质量受到严重影响，许多方面已不能适应当前环境保护工作的需要，改革势在必行。
地表水水环境容量计算方法回顾与展望_董飞
第25卷第3期 2014年5月水科学进展ADVANCES IN WATEＲSCIENCE Vol．25，No．3May ，2014 地表水水环境容量计算方法回顾与展望 董飞1，2，刘晓波1，2，彭文启1，2，吴文强 1，2（1．中国水利水电科学研究院水环境研究所，北京100038； 2.流域水循环模拟与调控国家重点实验室，北京100038） 摘要：为厘清中国地表水水环境容量计算方法演变历史，探讨计算方法发展趋势，在系统调研大量水环境容量研 究文献基础上，详细梳理水环境容量从概念引入到研究至今的过程，归纳出中国地表水水环境容量研究过程中产 生的五大类计算方法：公式法、模型试错法、系统最优化法（线性规划法和随机规划法）、概率稀释模型法和未确 知数学法。解析了各类方法的基本思路、产生过程及应用进展，评述了各类方法的优缺点及适用范围。通过与国 外水环境容量计算方法的比较，基于水环境系统复杂性及中国水资源管理特点与应用需求，认为中国应强化对概 率稀释模型法、未确知数学法及随机规划法等3种方法的研究和改进。 关键词：地表水；水环境容量；计算方法；概率稀释模型；系统最优化；未确知数学 中图分类号：TV131，X143；G353.11文献标志码：A 文章编号：1001- 6791（2014）03-0451-13收稿日期：2013- 10-11；网络出版时间：2014-04-10网络出版地址：http ：//https://www.wendangku.net/doc/2c19206781.html, /kcms /detail /32.1309.P.20140410.0950.010.html 基金项目：国家自然科学基金资助项目（51209230）；水体污染控制与治理科技重大专项（2013ZX07501- 004）作者简介：董飞（1983—），男，山东淄博人，博士研究生，主要从事流域容量总量控制理论与方法等研究。 E-mail ：dongfei99999@https://www.wendangku.net/doc/2c19206781.html, 通信作者：彭文启，E- mail ：pwq@https://www.wendangku.net/doc/2c19206781.html, 环境容量是环境科学的基本理论问题之一，是环境管理的重要实际应用问题之一［1］。水环境容量是环 境容量的重要组成部分，是容量总量技术体系的核心内容之一。随着中国水环境管理体系从浓度控制、目标 总量控制向容量总量控制的转变，实现流域水质目标管理 ［2］与水功能区限制纳污红线管理［3］，水环境容量理论及计算方法研究的重要性更加凸显。 早在20世纪70年代后期，随着环境容量概念的引入，中国学者即开始了对水环境容量的研究［4］。在经 过短时期的对水环境容量基本概念的强烈争论后，迅速实现从基本理论到实际应用，从定性研究到定量化计 算的转变［5］；同时注重吸收欧美等国的研究成果［6］。随着研究的不断深入，特别是水环境数学模型应用及 计算机技术的不断进步，逐渐形成了公式法 ［7］、系统最优化法［5］、概率稀释模型法［6］、模型试错法［8］等计算方法，盲数理论等不确定性数学方法也引入其中［9］。在地表水方面，水环境容量计算中所用的水环境数学模型从Streeter- Phelps 简单模型［5］发展到WASP 、Delft 3D 等大型综合模型软件［10］，计算区域从河段、河流发展到河口、湖库、河网、流域［11］，计算维数从一维发展到二维和三维［12］，计算条件从稳态发展到动 态［13］，所针对的污染物从易降解有机物、重金属发展到营养盐等［7］。近年来，常见关于水环境容量总体研究进展的文献 ［14-15］，然而未有专门系统论述水环境容量计算方法研究进展的文献；同时，文献中通常将中国水环境容量计算方法分为3类或4类 ［8，10］，笔者认为这难以对水环境容量计算方法作全面概括，本研究旨在弥补这一不足。以地表水水环境容量为重点，兼顾海洋水环境容量，大量调研中外文献，系统研究中国在地表水水环境容量计算方面从起步到当前的各种方法；同时对照欧美国家的计算方法，对中国地表水水环境容量计算方法进行重新归类。在解析各类计算方法研究及应用情况的基础上，对各类计算方法的优缺点及适用范围作了评述。在比较分析国内外计算方法特征的基础上，结合各类计算方法对复杂水环境系统的适应性及中国水资源管理特点对水环境容量计算的需求，对中国今后地表水水环境容量计算方法的发展趋势作了展望。DOI:10.14042/https://www.wendangku.net/doc/2c19206781.html,ki.32.1309.2014.03.020
数值分析实验报告1
实验一 误差分析 实验1.1（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（1.1）中19x 的系数作一个小的扰动。我们希望比较（1.1）和（1.2）根的差别，从而分析方程（1.1）的解对扰动的敏感性。 实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b = 的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve =