baumol wolfe matlab,MATLAB计算方法在选煤和动力配煤中的应用

随着市场经济的发展,人们对煤炭产品加工的规格质量及资源回收利用率要求越来越高,为了取得最大可能的经济效益和社会效益,各企业精打细算,各项生产统计和数据处理越来越繁重。以往,一些简单计算应用计算器即可完成,然后再用手工绘成图形。如煤质分析资料、浮沉试验资料等。而一些复杂计算,例如生产试验的统计分析、各种报表数据的规律研究以及新工艺的应用等等就要编写程序利用计算机计算。这时,不仅要对所采用的计算方法非常熟悉,而且对于所采用的算法语言(例如BASIC,FORTRAN,C语言等)还要有深刻了解,并且要有一定的编程技巧,编好程序后还要上机调试。这些工作,由于过程的复杂且要花费大量时间和精力往往使人望而却步。近年来,计算技术迅速发展,微机日益普及,同时出现了大量现代工程计算软件,如MAT-LAB,MATHEMATICA,MATHCAD等等。它们的共同特点是,集数值计算、符号解析运算、文字处理和图形功能为一体,算法效率高,数值稳定性好,使用非常方便。不必研究计算方法和算法语言以及编程等问题。现以MATLAB(矩阵实验室)为例,说明其在煤炭加工利用中的应用。1确定选煤厂利润模型某选煤厂处理量、总收入、总成本及总利润的统计结果列入表1中。试确定总利润与处理量间的关系模型。表1处理量、总收入、总成本及总利润表处理量Q/Mta-1总收入R/106元a-1总成本C/106元a-1总利润Y/106元a-1统计值计算值0.610.949.631.311.320.816.1610.825.345.331.020.5012.008.508.491.223.9613.2610.7010.701.426.5414.6911.8511.861.527.5015.5012.0012.001.628.2416.3911.8511.861.829.0618.4110.6510.622.029.0021.008.008.01首先,利用plot函数按表1中的数据将它们的关系绘至图1。由图1中总收入、总成本和总利润与处理量间的关系曲线形态可知,它可用二阶或三阶多项式来描述。现采用三阶多项式并列出计算式,即:x=[0.60.81.01.21.41.51.61.82.0];R=[10.9416.1620.5023.9626.5427.5028.2429.0629.00];C=[9.6310.8212.0013.2614.6915.5016.3918.4121.00];y=[1.315.348.5010.7011.8512.0011.8510.658.00];n=3图1总利润、总收入、总成本与处理量间的关系为了便于分析,我们首先对总收入和总成本进行拟合。利用MATLAB拟合函数:p=polyfit(x,R,n)计算得:R=-11x2+41.5x-10同理可算得:C=2.0727x3-5.3039x2+10.3815x+4.8566本例中,利润为收入与成本之差(不考虑其它因素),于是总利润的数学模型为:y=-2.0727x3-5.6961x2+31.1185x-14.8566由此模型计算出的利润值列入表1中。与实际统计值的对比可以发现,两者吻合得非常好。该模型利用无约束最优化解法(也可用其它方法,例如求导数后再求极值等)计算出,处理量为1.50Mt/a时总利润最大,达12106元/a。该模型由四项组成。第一项系数为负值,是促使利润下降的常数,该值越大对利润下降的影响也越大。由关系式R和C可以发现,该系数来自成本。当处理量增加时该项以三次方关系增大。它是由于处理量增大而引起生产技术、管理工作等附加投资的函