原煤平均密度的Matlab编程算法.pdf

《原煤平均密度的Matlab编程算法.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、煤炭加工与综合利用No．5，2011COALPROCESSING＆COMPREHENSIVEUTILIZATION35原煤平均密度的Matlab编程算法吴文华，范肖南(安徽理工大学材料科学与工程学院，安徽淮南232001)摘要：介绍了利用Matlab软件编程计算原煤浮沉试验数据密度级平均密度的具体过程，在实际的运用中，可以此为模板，减少实际运算中的繁琐重复操作，有效提高工作效率。关键词：Matlab；平均密度；选煤中图分类号：TQ533文献标识码：A文章编号：1005．8397(2011)05-0035-03——第i个密度级的浮物累计产率，％；

2、——一第(i一1)个密度级的浮物累计产率，％；1——密度的倒数，可视为累计浮物产P率y的函数。1=，，)(2)P1式(2)是密度的倒数与浮物累计产率l，的P函数关系，可以采用多项式函数拟合，函数形式较为简单，在数据点的拟合精度也较高，但选择不同函数次数时，在数据点之间可能存在震荡问题，造成其函数曲线形态有较大差异，使得平均密度计算结果差异较大。笔者采用复合双曲正切函数对密度曲线数据进行拟合处理，拟合精度较高。定义密度曲线拟合函数如下：Y=P)(3)将函数(3)代入平均密度计算公式(1)中，式(1)可以表达为：—～YP=(4)f厂(P)40P—在

3、求平均密度之前，首先要确定浮物累计产i—一㈩L1率和密度之间的函数关系f(P)，而且-厂(P)的表达式较为复杂，其积分的求解手工计算困难，利用Matlab软件编程求解可以较好的加以解决。收稿日期：2011-05-302实例分析作者简介：吴文华(1987一)，男，安徽池州人，安徽理工表1列出了潘三选煤厂6—3mm粒级浅槽人大学材料科学与工程学院2009级在读硕士研究生，研究方向：选矿工艺。料的密度组成，以此数据为例来求对应密度级的36煤炭加工与综合利用2011年第5期密度／kg·L<1．31．3～1．41．4～1．51．5～1．61．6～1．8>

4、1．8，J、计要确定表达式的形式，首先要求出b～bdisp(pjmd(3))值，然后代入式(5)对密度求导，再代入式disp(1．6—1．8的平均密度是：々(4)即可求得对应密度级的平均密度P。disp(pjmd(4))首先定义一个函数pjz(用来调用函数)，其中在Matlab应用软件中的操作为：打开eunlsum(A)用来求浮物累计产率。语句”myfuncMatlab，按照File—New—M—File的顺序打开M文件编辑器，在其中输入以下源程序。=inline(100：l：(beta(1)+beta(2)t+beta(3)：l=functi

5、on[]=pjz；tanh(beta(4)$(t—beta(5))))~beta；t’；”和symsX；语句”beta=nlinfit(t，Y，myfune，[00l11])；”A=input(请依次输入入料量：々；是用内联函数来求复合双曲正切函数的参数b1～B=[1．3141．51．61．8]；b5。digits函数用来规定运算精度，digits(5)规定C：elmsum(A)；了运算精度是5位有效数字；vpa函数用来控制运t=B；Y=C；算精度；int函数用来计算表达式关于某一自变量myfunc=inline("100(beta(1)+be

6、ta(2)t+beta的积分，dif函数用来求表达式对某一自变量的微(3)tanh(beta(4)(t—beta(5))))beta；t，；分，disp函数是屏幕输出函数。将编好的程序文beta：nlinfit(t，Y，myfunc，[00111]；件保存并命名为pjz．m文件。bl：beta(1)；b2=beta(2)；b3=beta(3)；Io4：在Matlab的CommandWindow窗口中输入beta(4)；b5=beta(5)；D=[b1b2b3b4h5]；pjz，然后按下回车键，程序会提示“请依次输入a=B(1，1)；b=B(1，

7、2)；c=B(1，3)；d=B(1，入料量：”，此时输入表1中第二行占本级的量，4)；e=B(1，5)；一直输人到1．6—1．8的占本级的量，用中括号f=C(1，1)；g=C(1，2)；h=C(1，3)；1=C(1，“[]”将它们括起来，最后按下回车键，得到程4)；m=C(1，5)；序的运行结果如下。digits(5)；n=vpa(int(dif(100(b1+b2{x+b3>>pjztanh(b4(X—b5))))／x，a，b))；请依次输入入料量：[10．8723．9113．824．97digits(5)；O=vpa(int(dif(100

8、(b1+b2X+b3％4．66]tanh(b4(x—b5))))／x，b，c))；Warning：Explicitintegra]couldnotbe