收稿日期:2001-01-13
作者简介:黄伟朵(1972)),女,浙江诸暨人,助教.
第13卷 第2期
浙江水利水电专科学校学报V ol.13 N o.2
2001年6月J.Zhejiang W at.Cons &Hy dr.Colleg e Jun.2001
基于参考面的快速计算可视域算法
黄伟朵
(浙江省水利水电专科学校水利工程系,浙江杭州 310016)
摘 要:与传统基于视线的可视域算法相比,利用参考面快速计算可视域算法不需进行DEM 内插,无冗余计算,计算时间与视点位置和可视域面积无关,受DEM 区域大小影响不大.关键词:可视域;DEM 应用;GI S 空间分析;可视性分析
中图分类号:P2 文献标识码:A 文章编号:1008-536X(2001)01-0034-02
An Algorithm for Computing Viewsheds based on R eference Planes
H UA N G Wei -duo
(Dept.of Water Conservancy Eng ineering,Zhejiang Water Conservancy and Hydr opower College,Hang zhou 310016,China)
Abstract:Compar ed with traditional sightline -based algorithms,the algor ithm for comput ing v iewsheds based on reference planes has no use for DEM interpllations,w ithout r edundance calculations.T he time for computing a viewshed has no thing to do w ith the lo cation of the view -piont and the area of viewshed,and it is slightly infuenced by the area of DEM S.Key words:v iew sheds;DEM applicat ions;GIS spat ila analysis;v isibility analysis 可视性分析作为一重要的G IS 空间分析,已广泛地应用于土木工程、环境规划、军事分析等领域.可视性分析包括两个基本因子,一个是两点间的通视性(Inter visibility),另一个是可视域(V iew shed),即对于给定的观察点所覆盖的区域.对于规则格网DEM ,计算可视域的一种简单方法就是沿着视线的方向,从视点开始到目标格网点,计算与视线相交的格网单元(边或面),判断相交的格网单元是否可视,从而确定视点与目标点之间是否可视.这种算法存在大量的冗余计算
[1]
.为提高效率,
Van 和Kreveld 提出了一基于/线扫描0的算法,Franklin 和R ay 提出利用辅助格网等算法,这些算法都是基于视线的DEM 格网单元,但须内插计算,计算量仍较大.
本文介绍的是一种完全脱离视线,不需DEM 内插计算,通过视点和目标点间的空间关系所形成的参考面来判断视点与所有目标点是否可视的算法,简称参考面法.该算法比任何基于视线的可视域算法都快[2].
1 基本原理
计算某视点s i ,j 的可视域时,该算法涉及到/参考面0的概念.此时,参考面定义为通过视点s i,j 和目标点d m,n 附近的两个或同行号或同列号的辅助格网点的平面P =p (m ,n).辅助格网具有与已知DEM 相同的点数,记为R={r m,n }(m =1,2,,M ;n =1,2,,N ),M 和N 为对应的DEM 行列数.
该算法基于这样的一个原则:视点s i ,j 与目标点d m,n 可见的前提是,它必须位于参考面p (m ,n )之内或之上,p (m ,n)由视点s i ,j 和在目标点d m,n /之前0(比d i,j 更接近s i,j )的两个相邻的辅助格网点形成.如果目标点d m,n 可见,则相应的辅助格网点r m,n 的值等于d m,n 的高程值,则可视矩阵点v m,n 值为真;否则r m,n 值等于正好使得从s i,j 到d m,n 可见的最小高程值Z,则将可视矩阵点v m,n 的值置为假.通过从视点向外计算,由辅助格网点和视点形成的参考面本质地定义了一个局部眼界,这可用于判断下一个DEM 格网点与视点s i,j 的可视性,这样反复计算,直到遇到DEM 的边界.最后通过对可视矩阵V 的统计和显示得到可视域的面积及分布.
2 具体步骤
一旦选定视点s i,j ,可将坐标原点移到s i,j ,再利用s i,j 与其邻近格网点的连线将DEM 区域分割为八条方向线和八个扇区.
(1)初始化辅助格网R 和可视矩阵V,将s i,j 及其八个邻近格网点的辅助格网点值置为对应D EM 格网点的高程值,相应可视矩阵点位置为真,其它辅助格网点值为零,相应可视矩阵点值为假,这样就假设这八个邻近格网点与s i,j 均可见.由于s i ,j 与邻近格网点之间无其它格网点遮挡,在DEM 精度范围内假设成立.
(2)处理八条方向线上的点,由于位于八条方向线上的目标点d m,n 的几何特征简单,形成判断其遮挡关系的参考面的两个相邻辅助格网点t1,t2重合为一点t1,此时参考面简化为参考
线,所以这些点的遮挡关系可由视点s i,j和在目标点d m,n/之前0的点t1m1,n1形成的参考线来判断.另外,判断扇区内目标点的遮挡关系时,需要利用方向线上的点的结果,因此程序首先从视点s i,j开始向外依次处理八条方向线上的所有点:
1)用下列公式计算s i,j与d m,n正好可见的最小高程值Z:
若m1-i X0,则Z=r i,j+
m-i
m1-i
(r m1,n1-r i,j);
否则,Z=r i,j+
n-i
n1-i
(r m1,n1-r i,j);
当目标点d m,n在西方向(W):点t1位于(m+1,n);
西北方向(N W):点t位于(m+1,n+1);
北方向(N):点t位于(m,n+1);
东北方向(N E):点t位于(m-1,n+1);
其它方向依此类推.
2)判断并赋值:
若d m,n>Z则视点s i,j与目标点d m,n可见,令r m,n=d m,n, v m,n=true;
否则s i,j与d m,n不可见,令r m,n=Z,v m,n=false;
3)将目标点d m,n向外移动,重复上述1)、2)步骤,如遇DEM边界则进入下一方向的循环,直到八个方向上的点均计算完毕,则进入213.
(3)处理八个扇区内的点,经过上面流程处理后,需处理的剩余点分别位于八个扇区内,这些点的遮挡关系可由视点s i,j 和在d m,n/之前0的两个或同行号或同列号的辅助格网点t1m1,n1和t2m2,n2形成的平面来判断,所以可分扇区地从视点s i,j 开始向外依次处理这些点:
1)计算s i,j与d m,n正好可见的最小高程值Z:
已知空间三点(x1,y1,z1),(x2,y2,z2),(x3,y3,z3)形成的空间平面方程为
x y z1
x1y1z11
x2y2z21
x3y3z31
=0
令x21=x2-x1;y21=y2-y1;z21=z2-z1;
x31=x3-x1;y31=y3-y1;z31=z3-z1;
有:z=z1-
(x-x1)(y21z31-y31z21)+(y-y1)(z21x31-z31x21)
x21y31-x31y21
所以,将下列值代入上式可计算从s i,j到d m,n正好可见的最小高程值Z:
x1=m1#dx;y1=n1#dy;z1=r m1,n1;
x2=m2#dx;y2=n2#dy;z2=r m2,n2;
x3=i#dx;y3=j#dy;z3=r i,j;
x=m#dx;y=n#dy;z=Z;
(dx和dy为DEM格网在X和Y方向的间距)
在此,关键在于正确指定形成参考平面P的两个辅助格网点t1m1,n1和t2m2,n2的坐标值:
当目标点d m,n位于W
-
N W扇区:t1位于(m+1,n),t2位于(m+1,n+1);
NW
-
N扇区:t1位于(m,n+1),t2位于(m+1,n+1);
N-N E扇区:t1位于(m-1,n+1),t2位于(m,n+1);
NE-E扇区:t1位于(m-1,n),t2位于(m-1,n+1);
其它扇区依此类推.
2)判断并赋值,若d m,n>Z则视点s i,j与目标点d m,n可见,此时r m,n=d m,n,v m,n=true;否则s i,j与d m,n不可见,此时r m,n =Z,v m,n=false
3)将目标点d m,n向远离视点s i,j的方向移动,并重复上述1)、2)步骤,如遇DEM边界则进入下一扇区的循环,直到八个扇区内的点均计算完毕,则结束可视域计算.
3算法比较
为了测试该算法的效率,用C语言编程实现该算法和最简单沿视线计算可视域算法,两者均以相同DEM、相同视点为参数在PC机上运行,除去输出输入时间,用这两种算法计算可视域所花CPU的时间进行比较可知.
用参考面法计算同一DEM、同一点、不同高度视点的可视域所花CPU时间为一常数,它与该点在DEM中的位置和可视域面积大小无关.而最简单沿视线法所花CPU时间却随视点在DEM中的位置和视点高度面变化,一般来说,计算位于DEM 中心的视点可视域所花CPU的时间要比其它点的时间少;对于同一视点,当视点高度增加,可视域面积增大,计算可视域花费的CPU时间增加增快.另外,DEM格网数增加对参考面算法效率的影响不大,而使最简单沿视线法运算效率降低很快.
4结语
基于参考面在格网DEM上计算可视域的新算法利用的是参考面而不是视线来计算可视域.与传统基于视线的可视域算法相比,由于不用进行视线相交格网的DEM内插计算,因此该算法简单、高效.另外,该算法还有一个有用的特性,它的运行时间与视点的位置和可视域面积大小无关.通过实验证明该算法结果与基于视线算法的结果一致.
参考文献:
[1]李志林,朱庆.数字高程模型[M],武汉:武汉测绘科技大学出版
社,2000.
[2]Jianjun Wang,Gray j.Robinson.Kevin White.Generating View sheds
without Using Sightines[J].Photogram metric Enginneering&Re-mote S ensing,2000,66(1):87~90.
[3]张祖勋,张剑清.数字摄影测量学[M],武汉:武汉测绘科技大学出
版社,1997.
[4]Peter F.Fisher.Extending the Applicaility of View sheds in Landscape
Planning[J].Photogrammetric Engineering&Remote Sensing,1996, 62(11):297~1302.
35
第13卷第2期黄伟朵.基于参考面的快速计算可视域算法
程建设其他费用参考计算方法 工程建设其他费用参考计算方法 工程建设其他费用参考计算方法 一、固定资产其他费用 (一) 建设管理费 1、以建设投资中的工程费用为基数乘以建设管理费率计算。 建设管理费= 工程费用×建设管理费费率 2、由于工程监理是受建设单位委托的工程建设技术服务,属建设管理范畴。如采用监理,建设单位部分管理工作量转移至监理单位。监理费应根据委托的监理工作范围和监理深度在监理合同中约定具体 3、如建设管理采用工程总承包方式,其总包管理费由建设单位与总包单位根据总包工作范围在合同中商定,从建设管理费中支出. 4、改扩建项目的建设管理费率应比新建项目适当降低。 5、建设项目按批准的设计文件所规定的内容组成,工业项目经负荷试车考核(引进国外设备项目按合同约定试车考核期满)或试运行期能够正常生产合格产品,非工业项目符合设计要求,能够正常使用时,应及时组织验收,移交生产或使用。凡已超过批准的试运行期,并已符合验收条件但未及时办理竣工验收手续的建设项目,视同项目已交付生产,其费用不得从基建投资中支付,所实现的收入作为生产经营收入,不再作为基建收入。 (二) 建设用地费 1、根据征用建设用地面积、临时用地面积,按省、市人民政府制定颁发的土地征用补偿费、安置补助费标准和耕地占用税、城镇土地使用税标准计算. 2、建设用地上的建(构)筑物如需迁建,其迁建补偿费应按迁建补偿协议计列或按新建同类工程造价计算。建设场地平整中的余物拆除清理费在“场地准备及临时设施费”中计算。 3、建设项目采用“长租短付”方式租用土地使用权,在建设期间支付的租地费用计入建设用地费,在生产经营期间支付的土地使用费应进入营运成本中核算。 (三) 可行性研究费 1、依据前期研究委托合同计列,或参照《国家计委关于印发〈建设项目前期工作咨询收费
《计算方法》期中复习试题 一、填空题: 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ,则用辛普生(辛卜生)公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案:2.367,0.25 2、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数为 , 拉格朗日插值多项式为 。 答案:-1, )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字; 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( ); 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 ); 6、计算方法主要研究( 截断 )误差和( 舍入 )误差; 7、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 ( 1 2+-n a b ); 8、已知f (1)=2,f (2)=3,f (4)=5.9,则二次Newton 插值多项式中x 2系数为( 0.15 ); 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f ),代数精 度为( 5 ); 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少,应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ,为了减少舍入误差,应将表达式 19992001-改写为 199920012 + 。
《数值计算方法》复习试题 一、填空题: 1、????? ?????----=410141014A ,则A 的LU 分解为 A ??? ?????????=? ?????????? ?。 答案: ?? ????????--??????????--=1556141501 4115401411A 2、已知3.1)3(,2.1)2(,0.1)1(===f f f ,则用辛普生(辛卜生)公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案:, 3、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数为 , 拉格朗日插值多项式为 。 答案:-1, )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字; 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( ); ( 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 ); 7、计算方法主要研究( 截断 )误差和( 舍入 )误差; 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 ( 1 2+-n a b ); 9、求解一阶常微分方程初值问题y '= f (x ,y ),y (x 0)=y 0的改进的欧拉公式为
( )] ,(),([2111+++++=n n n n n n y x f y x f h y y ); 10、已知f (1)=2,f (2)=3,f (4)=,则二次Newton 插值多项式中x 2系数为( ); 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f ),代数精 度为( 5 ); 12、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均 不为零)。 13、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少,应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ,为了减少舍入误差,应将表达式 19992001-改写为 199920012 + 。 14、 用二分法求方程01)(3 =-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间 为 ,1 ,进行两步后根的所在区间为 , 。 15、 、 16、 计算积分?1 5 .0d x x ,取4位有效数字。用梯形公式计算求得的近似值为 ,用辛卜 生公式计算求得的近似值为 ,梯形公式的代数精度为 1 ,辛卜生公式的代数精度为 3 。 17、 求解方程组?? ?=+=+042.01532121x x x x 的高斯—塞德尔迭代格式为 ?????-=-=+++20/3/)51()1(1)1(2)(2)1(1 k k k k x x x x ,该迭 代格式的迭代矩阵的谱半径)(M ρ= 121 。 18、 设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ,)(x f 的二次牛顿 插值多项式为 )1(716)(2-+=x x x x N 。 19、 求积公式 ?∑=≈b a k n k k x f A x x f )(d )(0 的代数精度以( 高斯型 )求积公式为最高,具 有( 12+n )次代数精度。
27×99 541×67-67×441 48×101-48 34×201 256×7-56×7 103×37 125×16 420÷35 76×23+24×23 103×23 25×(40+8)75×3×4
123×67-23×67 38×7+62×7 25×16×5 68×48+68×2 52×32+48×32 5×27+63×5 64×9-14×9 125×18 67+42+33+58 18×137-18×37 18×45+18×55 250×28
199×9+199 50×(60+8)49+49×49 304×22 12×(40-5)75×141-75×40 55×25+25×45 (30+4)×25 27×37+37×23 47+99+47 25×65+25×25 24×250
163×8+37×8 256×9-46×9 63×8+91×63+63 28×111-28×11 201×34 78×101-78 560÷16 373×9-73×9 2×46+46×1813×125×8 44×25 28×57+43×28
99×64+64 16×401 36×25 199×53+53 12+19×12 (30+2)×15 226×13-26×13 125×16 402×15 25×19 (30+8)×25 48×125
63+15×2 202×41 21+254+79+46 125×(8+16)202×13 13+13×49 304+297 25×124-24×25 41×99 56+56×49 15×301-15 250×9×4
建筑安装工程费用参考计算方法附件一: 建筑安装工程费用参考计算方法各组成部分参考计算公式如下: 一、直接费 (一)直接工程费 直接工程费=人工费+材料费+施工机械使用费 1、人工费 人工费=∑(工日消耗量×日工资单价) (1)基本工资 (2)工资性补贴
(3)生产工人辅助工资 (4)职工福利费 职工福利费(G4)=(G1+G2+G3)×福利费计提比例(%) (5)生产工人劳动保护费 2、材料费 材料费=∑(材料消耗量×材料基价)+检验试验费 (1)材料基价 材料基价=[(供应价格+运杂费)×(1+运输损耗率(%))]×(1+采购保管费率(%)) (2)检验试验费 检验试验费=∑(单位材料量检验试验费×材料消耗量) 3、施工机械使用费 施工机械使用费=∑(施工机械台班消耗量×机械台班单价) 机械台班单价 台班单价=台班折旧费+台班大修费+台班经常修理费+台班安拆费及场外运费+台班人工费+台班燃料动力费+台班养路费及车船使用税 (二)措施费
本处只列部分通用措施费项目的计算方法,专用措施费项目按施工组织设计计算。 1、环境保护 环境保护费=直接工程费(或其中人工费)×环境保护费费率(%) 2、文明施工 文明施工费=直接工程费(或其中人工费)×文明施工费费率(%) 3、临时设施费 临时设施费有以下三部分组成: (1)周转使用临建(如,活动房屋) (2)一次性使用临建(如,简易建筑) (3)其他临时设施(如,临时管线) 临时设施费=(周转使用临建费+一次性使用临建费)×(1+其他临时设施所占比例(%))其中: ①周转使用临建费
②一次性使用临建费 ③其他临时设施。 或:临时设施费=直接工程费(或其中人工费)×临时设施费费率(%) 4、夜间施工增加费
2002-2003 第一学期 一.计算及推导( 5*8) 1.已知 x* 3.141, x ,试确定 x * 近似 x 的有效数字位数。 * * * 0.100 * * * 2.有效数 x 1 3.105, x 2 0.001, x 3 1 x 2 3 ,试确定 x x 的相对误差限。 3.已知 f ( x) 0.5 x 3 0.1x 2 ,试计算差商 f 0,1,2,3 4.给出拟合三点 A (0,1), B (1,0) 和 C (1,1) 的直线方程。 5.推导中矩形求积公式 b (b a) f ( a b ) 1 f '' ( )(b a)3 f (x)dx a 2 24 b n f (x)dx A i f ( x i ) a 6.试证明插值型求积公式 i 0 的代数精确度至少是 n 次。 7.已知非线性方程 x f (x) 在区间 a, b 内有一实根,试写出该实根的牛顿迭代 公式。 8.用三角分解法求解线性方程组 1 2 1 x 1 0 2 2 3 x 2 3 1 3 0 x 3 2 二.给出下列函数值表 0.4 0.5 0.6 0.7 0.8 x i 0.38942 0.47943 0.56464 0.64422 0.71736 f ( x i ) 要用二次插值多项式计算 f (0.63891) 的近似值,试选择合适的插值节点进行计 算,并说明所选用节点依据。 (保留 5 位有效数字)(12 分) 三. 已知方程 x ln x 0 在 (0,1) 内有一实根 ( 1)给出求该实根的一个迭代公式,试之对任意的初始近似 x 0 (0,1) 迭代法都收 敛,并证明其收敛性。 ( 2) x 0 0.5 试用构造的迭代公式计算 的近似值 x n ,要求 x n x n 1 10 3 。 四. 设有方程组
计算方法考试题(一) 满分70分 一、选择题:(共3道小题,第1小题4分,第2、3小题3分,共10分) 1、将A 分解为U L D A --=,其中),,(2211nn a a a diag D =,若对角阵D 非奇异(即),1,0n i a ii =≠,则b Ax =化为b D x U L D x 1 1)(--++=(1) 若记b D f U L D B 111 1),(--=+= (2) 则方程组(1)的迭代形式可写作 ) 2,1,0(1 )(1)1( =+=+k f x B x k k (3) 则(2)、(3)称 【 】 (A)、雅可比迭代。(B)、高斯—塞德尔迭代 (C)、LU 分解 (D)、Cholesky 分解。 2、记*x x e k k -=,若0lim 1≠=+∞→c e e p k k k (其中p 为一正数)称序列}{k x 是 【 】 (A)、p 阶收敛; (B)、1阶收敛; (C)、矩阵的算子范数; (D)、p 阶条件数。 3、牛顿切线法的迭代公式为 【 】 (A)、 ) () (1k x f x f x x k k k '- =+ (B)、 )()())((111--+--- =k k k k k k k x f x f x x x f x x 1 )() ()1()()()(x x f x f x f k i k i k i ??+=+ (D)、 )() ()()1(k k k x f x x -=+ 二、填空题:(共2道小题,每个空格2分,共10分) 1、设0)0(f =,16)1(f =,46)2(f =,则一阶差商 ,二阶差商=]1,2,0[f ,)x (f 的二次牛顿 插值多项式为 2、 用二分法求方程 01x x )x (f 3 =-+=在区间]1,0[内的根,进行第一步后根所在的区间为 ,进行第二步后根所在的区间 为 。 三、计算题:(共7道小题,第1小题8分,其余每小题7分,共50分) 1、表中各*x 都是对准确值x 进行四舍五入得到的近似值。试分别指出试用抛物插值计算115的近似值,并估计截断误差。 3、确定系数101,,A A A -,使求积公式 ) ()0()()(101h f A f A h f A dx x f h h ++-≈? -- (1) 具有尽可能高的代数精度,并指出所得求积公式的代数精度。
附件1:工程建设其他费用参考计算方法 计算工程建设其他费应符合财政部、文物局《重点文物保护专项补助资金管理办法》(财教[2013]116号)、发改委《关于降低部分建设项目收费标准规行为等有关问题的通知》(发改价格[2011]534号),以及《基本建设财务管理规定》(财建[2002]394号)文件要求。 一、固定资产其他费用 (一)建设单位管理费。 建设单位管理费:是指建设单位从项目开工之日起至办理竣工财务决算之日止发生的管理性质的开支。包括:不在原单位发工资的工作人员工资、基本养老保险费、基本医疗保险费、失业保险费,办公费、差旅交通费、劳动保护费、工具用具使用费、固定资产使用费、零星购置费、招募生产工人费、技术图书资料费、印花税、业务招待费、施工现场津贴、竣工验收费和其他管理性质开支。 参照财政部关于印发《基本建设财务管理规定》的通知(财建[2002]394号)规定计算。 建设单位管理费总额控制数费率表 单位:万元
(二)青苗补偿费。 青苗补偿费是指征用土地时,农作物正处在生长阶段而未能收获,应给予土地承包者或土地使用者的经济补偿。 《土地管理法》规定,被征用土地,在拟定征地协议以前已种植的青苗和已有的地上附着物,也应当酌情给予补偿。但是,在征地方案协商签订以后抢种的青苗、抢建的地上附着物,一律不予补偿。被征用土地上的附着物和青苗的补偿标准,由省、自治区、直辖市规定。 以省为例:省物价局、省财政厅、省国土资源厅下发的《关于等三市调整征地地面附着物和青苗补偿标准的批复》(鲁价费发[2008]178号)计算市历城区100亩地的青苗补偿费如下: 市历城区100亩地青苗补偿费=100亩*900元/亩=90000元。 (三)研究试验费。 1. 按照研究试验容和要求进行编制。 研究试验费是指为项目提供和验证设计参数、数据、资料等所进行的必要的试验费用以及设计规定在项目实施中必须进行试验、验证所需费用。包括自行或委托其他部门研究试验所需人工费、材料费、试验设
二、设(2)0,(0)2,(2)8f f f -===,求 )(x p 使 )()(i i x f x p =,)2,1,0(=i ;又 设 M x f ≤''')( ,则估计余项 )()()(x p x f x r -= 的大小 。(15分) 三、设(0)1,(0.5)5,(1)6,(1.5)3,(2)2f f f f f =====,()k f M ≤(2,3,4)k =, (1)计算? 20 )(dx x f , (2)估计截断误差的大小(12分) 寂涯网络 https://www.wendangku.net/doc/a67427758.html, xx ~xx 学年第 1学期 《计算方法》课程试卷A 第 1 页 共 4 页
寂涯网络 https://www.wendangku.net/doc/a67427758.html, xx ~xx 学年第 1学期 《计算方法》课程试卷A 第 2 页 共 4 页 四、设方程012523=-+x x 在 [2,1] 内有实根α,试写出迭代公式 ,,2,1,0)(1 ==+k x x k k ? 使 {}α→k x ,并说明迭代公式的收敛性。 (10分) 五、设有线性方程组b Ax =,其中 ???? ? ??=?? ????????=582,3015515103531b A (1)求A LU =分解; (2) 求方程组的解 (3) 判断矩阵A 的正定性(14分)
寂涯网络 https://www.wendangku.net/doc/a67427758.html, xx ~xx 学年第 1学期 《计算方法》课程试卷A 第 3 页 共 4 页 六、设有线性方程组b Ax =,其中 1 442 12441A -?? ??=?? ???? , 试讨论Jacobi 迭代法和Gauss-Seidel 迭代法的收敛性。(14分) 七、设()i j n n A a ?=是n 阶实对称正定矩阵,A 经过一次高斯消元计算变为 ?? ? ???211A O T a , 其中T 为行向量,O 是零列向量,试证明2A 是对称正定矩阵(8分)
《计算方法》期末考试试题 一 选 择(每题3分,合计42分) 1. x* = 1.732050808,取x =1.7320,则x 具有 位有效数字。 A 、3 B 、4 C 、5 D 、6 2. 取7 3.13≈(三位有效数字),则 ≤-73.13 。 A 、30.510-? B 、20.510-? C 、10.510-? D 、0.5 3. 下面_ _不是数值计算应注意的问题。 A 、注意简化计算步骤,减少运算次数 B 、要避免相近两数相减 C 、要防止大数吃掉小数 D 、要尽量消灭误差 4. 对任意初始向量)0(x ?及常向量g ?,迭代过程g x B x k k ? ??+=+)() 1(收敛的充分必要条件是_ _。 A 、11< B B 、1<∞ B C 、1)(简便方法
整数简算·四则混合运算 【练习】 14.用简便方法计算下面各题. 3×999+3+99×8+8+2×9+9 125×128-125×27-125 ※(11×9+11)×(111×999+111)×(7×11×13-1001) (24×21×45)÷(15×4×7) (125×72×24)÷9÷8 111×111 1111×1111 999×999 9999×9999 15.利用数的分解法计算下面各题. (1)9+99+999+9999+99999 (2)2772÷28 (3)579999971÷29 (4)1986+331×594 (5)1111×58+6666×7 (6)99999×77778+33333×66666 (7)321×17+107×39+1070 (8)2999998+299997+29996+2995+294+23 16.用简便方法计算下列各题. (1)54+38+46 (2)37+44+56 (3)88+(37+22) (4)67+15+33 (5)375+342+658+625 (6)827+74+36+163
(7)428+267+(733+572) (8)536+(541+464)+469 (9)327+108(10)325+98 (11)872-48-272 (12)384-(184+36) (13)528-(138-72) (14)387-124 (15)564-387+187 (16)843+78-43 (17)274-87+26-13 (18)936-867-99+267 (19)813-(613-237) (20)537-(543-163)-57 (21)36×(468÷9) (22)58÷17×34 (23)48×5 (24)24×25 (25)56×125 (26)26×64×625 (27)84×(25×37) (28)68×36+36+31×36 (29)84×29-18×84-21×4 (30)72×(51÷12) (31)4321-1996+1998
建筑安装工程费用参考计算方法 一、各费用构成要素参考计算方法如下: (一)人工费 公式1: 人工费=刀(工日消耗量X 日工资单价) 生产工人平均月工资(计时、计件)?平均月(奖金?津贴补贴?特殊情况下支付的工资) 年平均每月 法定工作日 注:公式1主要适用于施工企业投标报价时自主确定人工费, 也是工程造价管理机构编 制计价定额确定定额人工单价或发布人工成本信息的参考依据。 公式2: 人工费=E (工程工日消耗量X 日工资单价) 日工资单价是指施工企业平均技术熟练程度的生产工人在每工作日 (国家法定工作时间 内)按规定从事施工作业应得的日工资总额。 工程造价管理机构确定日工资单价应通过市场调查、 根据工程项目的技术要求, 参考实 物工程量人工单价综合分析确定, 最低日工资单价不得低于工程所在地人力资源和社会保障 部门所发布的最低工资标准的:普工 1.3倍、一般技工2倍、高级技工3倍。 工程计价定额不可只列一个综合工日单价, 应根据工程项目技术要求和工种差别适当划 分多种日人工单价,确保各分部工程人工费的合理构成。 注:公式2适用于工程造价管理机构编制计价定额时确定定额人工费, 是施工企业投标 报价的参考依据。 (二) 材料费 1. 材料费 材料费=刀(材料消耗量X 材料单价) 材料单价=[(材料原价+运杂费)X 〔 1+运输损耗率(%)〕] X [1+采购保管费率(%)] 2. 工程设备费 工程设备费=E (工程设备量X 工程设备单价) 工程设备单价=(设备原价+运杂费)X [1+采购保管费率(% ] (三) 施工机具使用费 1. 施工机械使用费 施工机械使用费=刀(施工机械台班消耗量X 机械台班单价) 机械台班单价=台班折旧费+台班大修费+台班经常修理费+台班安拆费及场外运费 +台班 人工费+台班燃料动力费+台班车船税费 注:工程造价管理机构在确定计价定额中的施工机械使用费时, 附件3 : 日工资单价二 应根据《建筑施工机械
四年级数学用简便方法计算的几种类型 类型一:(注意:一定要括号外的数分别乘括号里的两个数,再把 积相加) (40+8)×25 125×(8+80)36×(100+50) 24×(2+10)86×(1000-2)15×(40-8) 类型二:(注意:两个积中相同的因数只能写一次) 36×34+36×66 75×23+25×23 63×43+57×63 93×6+93×4 325×113-325×13 28×18-8×28 类型三:(提示:把102看作100+1;81看作80+1,再用乘法分配 律)78×102 69×102 56×101 52×102 125×81 25×41 类型四:(提示:把99看作100-1;39看作40-1,再用乘法分配 律)31×99 42×98 29×99 85×98 125×79 25×39 类型五:(提示:把83看作83×1,再用乘法分配律) 83+83×99 56+56×99 99×99+99 75×101-75 125×81-125 91×31-91 四年级数学简便计算:方法归类
一、交换律(带符号搬家法) 当一个计算题只有同一级运算(只有乘除或只有加减运算)又没有括 号时,我们可以“带符号搬家”。适用于加法交换律和乘法交换律。256+78-56 450×9÷50 =256-56+78 =450÷50×9 =200+78 =9×9 =278 =81 二、结合律 (一)加括号法 1.当一个计算题只有加减运算又没有括号时,我们可以在加号后面直 接添括号,括到括号里的运算原来是加还是加,是减还是减。但是在 减号后面添括号时,括到括号里的运算,原来是加,现在就要变为减; 原来是减,现在就要变为加。(即在加减运算中添括号时,括号前是 加号,括号里不变号,括号前是减号,括号里要变号。) 例:345-67-33 789-133+33 =345-(67+33) =789-(133-33) =345-100 =789-100 =245 =689 2.当一个计算题只有乘除运算又没有括号时,我们可以在乘号后面直 接添括号,括到括号里的运算,原来是乘还是乘,是除还是除。但是 在除号后面添括号时,括到括号里的运算,原来是乘,现在就要变为
建设项目总投资费用参考计算方法 第一章建设项目总投资费用 第一条建设项目总投资费用 建设项目总投资=工程造价+增值税+资金筹措费+流动资金 第二条工程造价 工程造价=工程费用(不含税)+工程建设其他费用(不含税)+预备费(不含税) 第三条增值税 增值税应按工程费、工程建设其他费、预备费和资金筹措费分别计取。 第四条资金筹措费 1.自有资金额度应符合国家或行业有关规定。 2.建设期利息:根据不同资金来源及利率分别计算。 式中Q——建设期利息; Pj-1——建设期第(j-1)年末贷款累计金额与利息累计金额之和; Aj——建设期第j年贷款金额; i——贷款年利率; n——建设期年数。 3. 其他方式资金筹措费用按发生额度或相关规定计列。 第五条流动资金 流动资金的估算方法有扩大指标估算法和分项详细估算法两种。 (1)扩大指标估算法,此方法是参照同类企业的流动资金占营业收入、经营成本的比例或者是单位产量占用营运资金的数额估算流动资金,并按以下公式计算: 流动资金额=各种费用基数×相应的流动资金所占比例(或占营运资金的数额)
式中,各种费用基数是指年营业收入,年经营成本或年产量等。 (2)分项详细估算法,可简化计算,其公式如下: 流动资金=流动资产-流动负债 流动资产=应收账款+预付账款+存货+库存现金 流动负债=应付账款+预收账款 第二章工程费用 第六条工程费用 工程费用=建筑工程费+设备购置费+安装工程费 第七条建筑工程费 建筑工程费=直接费+间接费+利润 第八条直接费 直接费=人工费+材料费+施工机具使用费+其他直接费 第九条人工费 人工费=∑(工日消耗量×日工资单价) 日工资单价由工程造价管理机构通过市场调查、根据工程项目的技术要求,参考实物工程量人工单价综合分析确定。 第十条材料费 材料费=∑(材料消耗量×材料单价) 材料单价=[(材料原价+运杂费)×〔1+运输损耗率(%)〕]×[1+采购保管费率(%)] 第十一条施工机具使用费 施工机具使用费=施工机械使用费+施工仪器仪表使用费 1.施工机械使用费
计算方法模拟试题 一、 单项选择题(每小题3分,共15分) 1.近似值210450.0?的误差限为( )。 A . 0.5 B. 0.05 C . 0.005 D. 0.0005. 2. 求积公式)2(3 1 )1(34)0(31)(2 0f f f dx x f ++≈ ?的代数精确度为( )。 A. 1 B. 2 C. 3 D. 4 3. 若实方阵A 满足( )时,则存在唯一单位下三角阵L 和上三角阵R ,使LR A =。 A. 0det ≠A B. 某个0 det ≠k A C. )1,1(0det -=≠n k A k D. ),,1(0det n k A k =≠ 4.已知?? ?? ? ?????=531221112A ,则=∞A ( )。 A. 4 B. 5 C. 6 D 9 5.当实方阵A 满足)2(,221>>-=i i λλλλ,则乘幂法计算公式1e =( )。 A. 1+k x B. k k x x 11λ++ C. k x D. k k x x 11λ-+ 二、填空题(每小题3分,共15分) 1. 14159.3=π,具有4位有效数字的近似值为 。 2. 已知近似值21,x x ,则=-?)(21x x 。 3.已知1)(2-=x x f ,则差商=]3,2,1[f 。 4.雅可比法是求实对称阵 的一种变换方法。