数学建模作业 2副本.docx

1、数学建模作业 2 副本矿物加工工程数学模型课程大作业学院（部）： 材料科学与工程学院 专业班级： 矿物加工11-2 学 号： 2011302978 学生姓名： 授课教师： 范肖南 2014年 5月8日学生姓名： 学号： 2011302978 专业班级： 矿物加工工程11-2班课程大作业题目： 潘三3-0.5mm粒级浅槽入料和产品的密度分析 授课教师评语：成绩： 授课教师： 年 月 日安徽理工大学课程大作业成绩评定表课题：潘三3-0.5mm粒级浅槽入料和产品的密度分析（11矿加2班耿宴作业）课题数据： 密度kg/L入料精煤矸石占本级(%)灰分(%)占本级(%)灰分(%)占本级(%)灰分(%)-1

2、.316.31 5.38 16.58 4.21 0.32 5.13 1.31.421.82 11.16 44.82 9.57 0.81 10.84 1.41.511.99 20.04 19.43 19.13 0.81 21.88 1.51.64.80 28.30 8.03 28.65 0.81 31.86 1.61.85.04 39.48 6.99 39.96 3.76 47.46 +1.840.05 83.65 4.15 68.49 93.50 84.44 小计100.00 42.56 100.00 16.64 100.00 81.27 课题要求：1.绘制原煤可选性曲线； 2.计算产品产率；

3、3.计算分配率；4.建立分配曲线模型；5.求实际原煤；6.绘制实际原煤与精煤产率、精煤灰分及分选密度的关系曲线；7.求=20对应的精煤产率、精煤灰分和分选密度；8.要求有完整的封面、目录、参考文献，步骤清晰、计算准确、图形完整、曲线平滑。利用Excel电子表格软件进行原煤可选性曲线数据计算十分方便，为选矿工作者广泛采用。采用Excel进行数据处理结果可以很容易地转换成图表，用图形的方式来观察数值的变化趋势和数据之间的关系，更加直观、易理解。在参考范肖南老师的论文采用Excel绘制可选性曲线的基础上，具有较高的精度和参考价值。一、绘制原煤可选性曲线绘制原煤可选性曲线所需原始数据如下：密度kg/L

4、入料浮物累计沉物累计分选密度kg/L占本级(%)灰分(%)重量，%灰分，%重量，%灰分，%密度级产率，%-1.316.31 5.38 16.31 5.38 100.00 42.56 1.30 38.13 1.31.421.82 11.16 38.13 8.69 83.69 49.81 1.40 33.81 1.41.511.99 20.04 50.12 11.40 61.87 63.44 1.50 16.79 1.51.64.80 28.30 54.92 12.88 49.88 73.87 1.60 6.68 1.61.85.04 39.48 59.95 15.11 45.08 78.72 1

5、.70 5.04 +1.840.05 83.65 100.00 42.56 40.05 83.65 1.80 63.11 小计100.0042.56注：原煤可选性曲线共5条，浮物累计曲线、基元灰分曲线、沉物累计曲线、密度曲线、以及密度+/-0.1含量曲线。1.1向图表中添加曲线单击工作表上的“图表导向”按钮，选取“散点图”中的“平滑曲线散点图”，单击完成按钮，获得嵌入数据工作表形式的图表。自动生成的图坐标需要进行调整，用鼠标右键双击Y轴，出现“坐标轴格式”对话框，打开“刻度”对话框，选中“数值次序反转”选项，为避免X轴在Y坐标转置后处于图形上方，选中“数值（X）轴置于最大值”项。用鼠标右键单击

6、图表中的任意空白处，出现图表选项框，打开数据源对话框，选择系列对话框，单击添加按钮，在名称栏中输入欲添加的曲线的名称，单击“X值”输入框的右端，出现“数据源-X值”输入框，用鼠标拖拽的方法选择浮物累计灰分值单元格区域，该单元格区域地址自动输入到“数据源-X值”输入框中，单击输入框的右端，回到系列对话框，浮物累计灰分单元格区域已经输入了。用同样的方法输入浮物累计产率所在单元格区域的地址。与上述步骤相同，分别添加其他几条曲线。下表为各曲线绘图基础数据：，基元灰分，浮物，沉物,密度,邻近物密度灰分，产率，灰分，产率，灰分，产率，密度,kg/l产率，密度,kg/l产率，505042.56 0.00 1

7、.383.69 1.361.87 5.38 8.15 5.38 16.31 49.81 16.31 1.461.87 1.466.19 11.16 27.22 8.69 38.13 63.44 38.13 1.549.88 1.583.21 20.04 44.12 11.40 50.12 73.87 50.12 1.645.08 1.693.32 28.30 52.52 12.88 54.92 78.72 54.92 1.840.05 1.794.96 39.48 57.43 15.11 59.95 83.65 59.95 1.836.89 83.65 79.98 42.56 100.00 9

8、5100.00 951001.2数据的细化与曲线绘制1.2.1数据细化对潘三数据可选性曲线的细化，主要针对浮物曲线、基元灰分特性曲线、沉物曲线、密度曲线进行了细化，采用了符合双曲正切数学模型，公式为“y=100*(b1*tanh(b2*(x-b3)+b4*x+b5)”,进行规划求解，然后对各项数据进行细化，利用规划求解出的b1、b2、b3、b4、b5再次进行计算，得到的具体数据如下表。规划求解主要步骤：1. 找到“规划求解”选项。用鼠标单击菜单栏“工具”，选择“加载宏”，之后再选中“规划求解”；2. 之后任选五个值，例如全取“1”，分别对应b1至b5,然后在对应Excel单元格中输入公式，其中

9、应该注意的是对b1至b5的绝对引用，然后就可以在求出一个值后直接用鼠标拖拽即可求出其余的值；3. 刚开始求出的数值偏差比较大，需要连续进行规划求解，直至所求出的值与原数值相差较小为止，即Q值尽量小，其中可以借助设置“工具”下拉菜单中的“选项”，选择“重新计算”，设置“最多迭代次数”为“10000”；4. 在进行细化时，注意变量的间隔要适当，求出适当数量的数据即可，此时需要规划求解最后得到的最合适的b1到b5的数值，还有上述提到的公式，四组数据的细化原理相同。，基元灰分复合双曲正切，基元灰分细化，浮物复合双曲正切，浮物细化灰分，产率，灰分，产率，灰分，产率，灰分，产率，5050.00 50505

10、.388.15 8.17 817.25 5.38 16.31 16.31 834.12 11.1627.22 27.15 1126.69 8.69 38.13 38.11 1148.76 20.0444.12 44.47 1434.36 11.40 50.12 50.22 1457.59 28.352.52 51.72 1740.15 12.88 54.92 54.78 1763.51 39.4857.43 58.03 2044.42 15.11 59.95 60.00 2068.28 83.6579.98 79.89 2347.61 42.56 100.00 100.00 2372.67 9

11、510085.48 2650.11 Q=0.03 2676.91 Q=1.12 2952.18 2981.11 3254.00 3285.29 参数估值3555.67 参数估值3589.47 b1=0.07 3857.26 b1=0.09 3893.65 b2=0.32 4158.81 b2=0.32 4197.83 b3=0.19 4460.32 b3=0.36 4299.22 b4=5.76 4761.82 b4=5.40 42.56 100.00 b5=0.00 5063.31 b5=0.01 5364.79 5666.27 5967.75 6269.23 6570.71 6872.18

12、 7173.66 7475.14 7776.62 8078.09 8379.57 8681.05 8982.53 9284.00 9585.48 9886.96 10087.94 ，沉物复合双曲正切,密度双曲正切,密度细化灰分，产率，密度,kg/l产率，密度,kg/l产率，42.56 0.00 3.96 49.81 16.31 14.33 1.383.69 83.71 1.383.71 63.44 38.13 33.85 1.461.87 61.81 1.3278.41 73.87 50.12 48.78 1.549.88 49.96 1.3473.58 78.72 54.92 55.72 1

13、.645.08 45.04 1.3669.22 83.65 59.95 62.79 1.840.05 40.05 1.3865.30 95100Q=0.01 1.461.81 Q=48.40 1.4258.72 1.4456.02 参数估值参数估值1.4653.67 b1=-0.28 b1=11.72 1.4851.66 b2=1.00 b2=1.23 1.549.96 b3=1.00 b3=1.67 1.5248.54 b4=0.01 b4=-14.53 1.5447.37 b5=-0.28 b5=24.73 1.5646.42 1.5845.65 1.645.04 1.6244.56 1.

14、6444.16 1.6643.82 1.6843.49 1.743.16 1.7242.77 1.7442.30 1.7641.71 1.7840.97 1.840.05 1.2.2 曲线绘制可选性曲线的5条曲线全部位于同一图表中，需要考虑到次坐标轴的创建，参照范肖南老师的论文，步骤如下：1.首先创建第二条Y轴，用鼠标双击要按照此坐标轴绘制的曲线，出现快捷菜单，单击“数据系列格式”项，打开“数据系列格式”对话框，在“坐标轴”对话框中选择“次坐标轴”，按“确定”即可出现第二条Y轴；2.创建第二条X轴，用鼠标右键双击图表中的空白处，出现快捷菜单，选择并打开“图表选项”，单击“坐标轴”表，选中“数值

15、X值”选项，单击确定之后，对应曲线出现变化，用鼠标双击第二条Y轴，打开“坐标轴格式对话框”，单击“刻度表”，选择“数值X值置于最大值”选项，并将其刻度设置为1.3至2.2，单击“数值次序反转”，取消“X轴置于最大值”选项，单击确定，即可完成次坐标轴的建立。3.可选性曲线的进一步完善。原煤可选性曲线绘制：二、计算产品产率主要运用格式法计算产品产率，公式为linest(f-t,c-t,m-t,0,0)。在进行计算产率之前需要求出f-t、c-t、m-t，这三项数据在下表中：密度kg/L入料精煤矸石f-t，c-t，占本级(%)灰分(%)占本级(%)灰分(%)占本级(%)灰分(%)-1.316.31 5

16、.38 16.58 4.21 0.32 5.13 15.98 16.26 1.31.421.82 11.16 44.82 9.57 0.81 10.84 21.02 44.01 1.41.511.99 20.04 19.43 19.13 0.81 21.88 11.18 18.62 1.51.64.80 28.30 8.03 28.65 0.81 31.86 3.99 7.23 1.61.85.04 39.48 6.99 39.96 3.76 47.46 1.27 3.23 +1.840.05 83.65 4.15 68.49 93.50 84.44 -53.45 -89.35 小计100.0

17、0 42.56 100.00 16.64 100.00 81.27 产品产率计算结果：精煤 矸石产率0.59 0.41三、分配率的计算潘三3-0.5mm粒级浅槽入料和产品的密度数据需要计算分配率E1和E2，结合下表，其中E1=（K5+L5）/M5*100；E2=L6/（L6+K6）*100；J、K、L三列数据即精煤、中煤、矸石三列数据是用上一题题目中的产品产率的计算结果求得，公式也相当简单，分别是各自原始数据“占本级产率”与产品产率的乘积，算的结果也在下表中;其余各项数据均可以非常容易的从原数据中直接获得或者经过简单的求和计算得到。分配率计算表格如下：密度kg/L入料精煤矸石精煤矸石计算入料，

18、占本级(%)灰分(%)占本级(%)灰分(%)占本级(%)灰分(%)占入料产率（%）占入料产率（%）-1.316.31 5.38 16.58 4.21 0.32 5.13 9.71 0.13 9.84 1.31.421.82 11.16 44.82 9.57 0.81 10.84 26.24 0.33 26.57 1.41.511.99 20.04 19.43 19.13 0.81 21.88 11.37 0.33 11.71 1.51.64.80 28.30 8.03 28.65 0.81 31.86 4.70 0.33 5.04 1.61.85.04 39.48 6.99 39.96 3.7

19、6 47.46 4.09 1.56 5.65 +1.840.05 83.65 4.15 68.49 93.50 84.44 2.43 38.77 41.19 小计100.00 42.56 100.00 16.64 100.00 81.27 四、建立分配曲线模型首先需要进行计算入料的方差以及分配率表格的整理，前五列数据均可以从上面计算结果得到，“计算入料”等于“产品占入料”中的“精煤”和“矸石”的和，“偏差”是“计算入料”与“原煤”的差值，“偏差平方和”是“偏差”的平方，“分配率”是“矸石”除以“计算入料”再乘以100即可。1.4.1计算入料的方差及分配率表：密度kg/L平均密度入料精煤矸石计算

20、入料,偏差偏差平方和分配率占本级(%)占入料产率（%）占入料产率（%）-1.31.25 16.31 9.71 0.13 9.84 -6.47 41.83 1.36 1.31.41.35 21.82 26.24 0.33 26.57 4.75 22.54 1.26 1.41.51.45 11.99 11.37 0.33 11.71 -0.28 0.08 2.85 1.51.61.55 4.80 4.70 0.33 5.04 0.24 0.06 6.64 1.61.81.70 5.04 4.09 1.56 5.65 0.62 0.38 27.58 +1.82.20 40.05 2.43 38.77

21、 41.19 1.14 1.31 94.11 用复合双曲正切模型进行数据细化处理，表格如下（其中红字是分配率为50时对应的密度）：模型有：反正切函数、双曲正切模型、复合双曲正切模型，公式分别为：1. y=100*(b1-atan(b2*(x-b3)/(b1-b4)2. y=100*(b1+b2tanh(b3*(x-b4)3. y=100*(b1+b2*x+b3*tanh(b4*(x-b5)其中，公式1中的b1和b4在任意取时不能相等。平均密度分配率双曲正切复合双曲反正切1.25 1.36 1.27 1.22 0.78 1.35 1.26 1.48 1.77 1.71 1.45 2.85 2.6

22、6 1.62 3.26 1.55 6.64 6.70 8.48 6.34 1.70 27.58 27.58 26.43 27.59 2.20 94.11 94.11 94.28 94.11 Q=0.10 Q=6.55 Q=0.788192b1=0.48 b1=-1.83 b1=-1.48 b2=5.77 b2=0.23 b2=17.79 b3=1.78 b3=-13.80 b3=1.74 b4=-0.03 b4=1.36 b4=1.63 b5=0.52 b5=1.36 双曲正切细化平均密度分配率1.25 1.27 1.82 59.00 1.28 1.28 1.85 66.01 1.311.33

23、 1.88 72.48 1.34 1.43 1.91 77.89 1.37 1.60 1.94 82.24 1.41.87 1.97 85.62 1.43 2.28 2.00 88.18 1.46 2.89 2.03 90.08 1.493.76 2.06 91.45 1.52 4.99 2.09 92.44 1.55 6.70 2.12 93.13 1.589.04 2.14 93.61 1.61 12.16 2.17 93.93 1.64 16.22 2.20 94.13 1.6721.35 2.23 94.26 1.70 27.58 2.26 94.33 1.73 34.79 2.29 94.35 1.7642.73 分选密度1.79 1.79 50.00 1.4.2分配率曲线经处理后如下：1.4.3 密度曲线平移数据处理以及绘图