自动化车床管理相关知识

{管理信息化OA自动化}自动化车床管理相关知识石伟、向中辉、喻欢自动化车床管理摘要本文建立的是自动化车床单刀具连续加工零件工序定期检查和刀具更换的谢谢阅读谢谢阅读精品文档放心下载序得到最好的效益，我们针对三个不同的问题建立了三个最优化模型。精品文档放心下载针对问题一我们将检查间隔和道具更换策略的问题确定为单个零件期望损感谢阅读谢谢阅读感谢阅读谢谢阅读用计算机采用穷举搜索法求得模型一的最优解为每生产18个零件检查一次，定感谢阅读期更换刀具间隔u为360，相应的单个零件的最小费用C为4.595元。感谢阅读针对问题二我们采用分摊法建立了单个零件的效益函数。首先，求出在一精品文档放心下载精品文档放心下载精品文档放心下载精品文档放心下载法求得检查间隔n为27u为297为问题二的最好检查间隔和换刀策略，感谢阅读相应的单个零件的最小费用C为9.216元。针对问题三在问题二的建模基础上，我们增加每次检查零件的次数来减少感谢阅读精品文档放心下载谢谢阅读n为29精品文档放心下载间距u为319为最优解，相应的单个零件的最小费用C为5.238元。感谢阅读最后，针对问题一的结果，对部分参数进行了灵敏度分析。感谢阅读关键词：卡方拟合正态分布检验穷举法单目标期望灵敏度分析精品文档放心下载1.问题重述1.1问题背景感谢阅读谢谢阅读感谢阅读感谢阅读谢谢阅读的好坏，便于对机器进行修理，减少不必要的损失，以提高经济效益。精品文档放心下载1.2问题提出感谢阅读现故障，其中刀具损坏故障占95%,其它故障仅占5%。工序出现故障是完全随机谢谢阅读的,假定在生产任一零件时出现故障的机会均相同。工作人员通过检查零件来确谢谢阅读100谢谢阅读零件数如附表1。现计划在刀具加工一定件数后定期更换新刀具。精品文档放心下载已知生产工序的费用参数如下：1、故障时产出的零件损失费用f=200元/件；2、进行检查的费用t=10元/次；3、发现故障进行调节使恢复正常的平均费用d=3000元/次(包括刀具费)；精品文档放心下载4、未发现故障时更换一把新刀具的费用k=1000元/次。感谢阅读1.3需要解决的问题,感谢阅读谢谢阅读②如果该工序正常时产出的零件不全是合格品，有2%谢谢阅读障时产出的零件有40%为合格品，60%感谢阅读生的损失费用为1500元/感谢阅读③在②的情况,可否改进检查方式获得更高的效益。表1100次刀具故障记录(完成的零件数)459362624542509584433748815505谢谢阅读612452434982640742565706593680谢谢阅读9266531644877346084281153593844感谢阅读527552513781474388824538862659感谢阅读775859755649697515628954771609精品文档放心下载402960885610292837473677358638谢谢阅读69963455557084416606106484120谢谢阅读447654564339280246687539790581谢谢阅读621724531512577496468499544645精品文档放心下载764558378765666763217715310851谢谢阅读2.问题分析2.1问题一的分析感谢阅读谢谢阅读精品文档放心下载精品文档放心下载精品文档放心下载失费用有刀具在未发生故障而采取预防性更换的损失费用和刀具发生故障而带谢谢阅读谢谢阅读设备的工作，到了定期更换刀具的时期，即使设备未出现故障，也要更换刀具。精品文档放心下载精品文档放心下载感谢阅读谢谢阅读谢谢阅读感谢阅读感谢阅读精品文档放心下载件出现故障的几率相同且相互独立，对于刀具更换周期来说，因为它是的函数，谢谢阅读感谢阅读数为：感谢阅读周期的区间，再利用有穷列举法得出最优化组合，则认为该组合为问题一的解。精品文档放心下载2.2问题二的分析精品文档放心下载精品文档放心下载精品文档放心下载谢谢阅读刀具正常产出的产品是合格品的为98%，只是在刀具损坏产出的产品合格率变化谢谢阅读精品文档放心下载感谢阅读精品文档放心下载谢谢阅读求出目标函数。检查到2%中因检查到不合格

刀具正常检查到40%因没有检查到不

刀具故障图1问题二分析图2.3问题三的分析精品文档放心下载谢谢阅读谢谢阅读谢谢阅读谢谢阅读们考虑采用不相等间隔的检查周期，在问题（12精品文档放心下载论检查周期的长短都要检查n个零件，但是实际上，随着检查周期的不断缩短，谢谢阅读检查的零件数n感谢阅读谢谢阅读间隔，采取不等距检查的策略来降低损失费用，即在稳定期内检查时间间隔大，精品文档放心下载感谢阅读着换刀周期递减的。3.模型假设与符号说明3.1模型的假设（1）假定在生产任一零件时出现故障的机会均相同；（2）假设工序出现故障是完全随机的；（3）假设题目所给的数据真实可靠，且不考虑其他的费用；谢谢阅读（4）假设工序在发生故障后，无法判断是刀具故障还是非刀具故障；谢谢阅读（5）假设非刀具故障在任一零件处发生的概率均相同；（6）假设从换上新刀到下一次工序故障或预防性换刀构成一个周期；精品文档放心下载3.2符号的说明符号符号说明预防性换刀的固定间隔(刀具更换周期)刀具更换周期的长度的平均值检查的固定间隔预防性换刀对的间隔中检查的次数一个周期内合格零件件数故障时产出的每件零件损失费用元/件每次进行检验的费用元/次每次发现故障进行调节使恢复正常的平均费用（包含刀具费）元/次感谢阅读每次未发现故障时更换一把新刀具的费用元/次刀具出现故障后生产合格零件的概率刀具没有出现故障时生产合格零件数工序正常而误认为有故障停机产生的损失费用元/次一个周期内总的故障间隔每个零件的平均费用工序中每一零件处非刀具故障率一个周期内平均费用在一周期内未出现故障所产生的费用一周期内出现故障所产生的费用出现非刀具故障在一周期内所产生的零件数出现非刀具故障时在一周期内所生产的零件数的概率密度函数间隔期（0，u）刀具故障发生次数的数学期望的分布函数工序故障是由刀具故障引起的的概率刀具未故障引起的损失由刀具发生故障引起的损失4.数据处理与假设检验由附表1的100谢谢阅读略了解这些数据的分布情况，我们先根据所给数据画直方图：精品文档放心下载图2刀具故障时已生产的零件个数的频数直方图100次刀具故障记谢谢阅读录的完成零件数观察研究及用Excel精品文档放心下载服从正态分布，根据记录数据求出了，则有发现这100次刀具故障时完成的零件精品文档放心下载数近似服从的正态分布Excel—贝拉感谢阅读检验和卡方优度检验，经软件处理得到如下表格表2雅克-贝拉检验和卡方拟合优度检验结果雅克-贝拉检验卡方拟合优度检验数据100数据个数100零假设是正态分布平均值600显著性水平0.05标准偏差196.最小值84中间计算最大值1153数据个数100区间个数16雅平均值600区间宽度66.875克-标准偏差196.假设检验偏度-0.25零假设服从正态分布贝拉峰度0.49自由度9检验JB统计量0.02卡方统计量2.19p值0.63p值0.68综合检验结果接受零假设显著性水平0.05考虑结果接受零假设了正态分布的宽度和偏斜对分布的影响，由结果知JB统计量为0.02且P=0.63>显谢谢阅读著性水平0.05，即假设正确。而卡方拟合优度检验得结果P=0.68>显著性水平感谢阅读0.05,进一步证明了假设的正确性。综上检验我们接受原假设，即表1中的数据服从的正态分布。谢谢阅读图3正态性分布检验图所以，根据以上分析，我们可以得到一下表达式：刀具寿命的概率密度函数为所以刀具寿命的发布函数为5.问题一的解答针对问题一我们建立了模型一。5.1模型一的建立5.1.1确定目标函数首先求非刀具故障和刀具故障总的分布函数：谢谢阅读感谢阅读感谢阅读从几何分布其分布函数为：由题意知，刀具损坏故障占工序故障的95%，而其它故障占5%，近似求得：精品文档放心下载求得所以总的故障间隔的分布函数然后求一个周期内平均损失费用对于模型一，可分为两种情况进行考虑感谢阅读合格，则前面生产的零件都合格，即工序正常，这时的费用为：谢谢阅读精品文档放心下载第(m-1)n+i个，这时费用记为C2。图4模型一求解图最后求刀具更换周期：谢谢阅读即得刀具更换的平均周期：所以问题二的目标函数为5.1.2综上所述，得到问题一的单目标最优化模型5.2模型一的求解对于模型一，我们用matlab编程进行穷举法进行搜索，由刀具寿命呈正态分布谢谢阅读600个零件附近时出现故障的概率最大，所以刀具更换周期应精品文档放心下载在600附近，又因为和有函数关系，所以检查周期也应有一定的范围.因此，我感谢阅读们判定和必定落在（10~100u的取值以5为步长计算费用。得到精品文档放心下载的主要结果如下表（详细结果见附录三）：一个检查间隔内完成的在一个换刀周期内完成单个零件的损失费用精品文档放心下载零件件数，单位（件）单位（元）谢谢阅读183604.5955.3模型一结果分析我们对模型一中部分结果用如下表格表示：由表3可以看出u和n的取值都会制约Cu和n的取值分别为360和感谢阅读18u和n精品文档放心下载接近这个最小值，所以我们认为C=4.595为模型一的最优解。谢谢阅读表3问题一的数据分析un161718193454.6214.6104.6014.602谢谢阅读3504.6174.6014.5994.598谢谢阅读3554.6114.5944.5974.596谢谢阅读3604.6074.5974.5954.596感谢阅读3654.6024.5984.5964.598谢谢阅读3704.6044.5994.5994.602精品文档放心下载3754.6114.6024.6034.607精品文档放心下载结果说明：当工序设计换刀周期为生产360个零件，检查周期为18个零件谢谢阅读时，所获得的效益最好，即生产单位合格零件的消耗费用最少。精品文档放心下载6.问题二的解答针对问题二，我们建立了模型二。6.1模型二的建立6.1.1确立目标函数精品文档放心下载感谢阅读谢谢阅读一个换刀周期内换刀次数的期望值。12……个零件时，感谢阅读设为第一次故障发生在第个零件时.然后对可能会发生的不同情况，在模型一的基础上建立两个过渡模型：感谢阅读过渡模型一：假设工序处于正常状态下，次品率为2%，而当工序不正常时，生产的零件全谢谢阅读精品文档放心下载检查费、换零件损失费。则每个零件的平均费用为：过渡模型二：假设工序处于正常状态时产品全为正品，而当工序处于不正常状态下，次品精品文档放心下载率为60%，有则有平均每个零件的费用包含四个部分：预防刀具费用、故障调整谢谢阅读费、检查费、换零件损失费，所以综合过渡模型一和过渡模型二，问题二的模型为：其中为故障调节恢复费用和检查费用：为坏零件损失费：为误停机造成的损失费：所以，模型二的目标函数为6.1.2综上所述，得到问题二的最优化模型6.2模型二的求解谢谢阅读谢谢阅读谢谢阅读致分布范围，用Matlab进行计算，最终求得的主要结果如下表（详细结果见附谢谢阅读录五）：一个检查间隔内完成的在一个换刀周期内完成单个零件的损失费用感谢阅读零件件数，单位（件）单位（元）感谢阅读272979.2158n=27件，即生产27件零精品文档放心下载=297件，即每生产297精品文档放心下载C=9.2158元。6.3模型二的结果分析6当uC的会感谢阅读nC的取感谢阅读值会随着u谢谢阅读有当u和nn=27,u=297感谢阅读时，目标函数C的取值最小为9.2158元。表4模型二的数据分析un262728292959.23369.22959.22879.2307感谢阅读2969.22309.22029.22079.2241谢谢阅读2979.21719.21589.21769.2224精品文档放心下载2989.21609.21619.21939.2256谢谢阅读2999.21959.22109.22589.2335谢谢阅读3009.22749.23059.23679.2460谢谢阅读3019.23989.24449.25229.2630谢谢阅读7.问题三的解答针对问题三，我们建立了模型三。7.1模型三的建立7.1.1确立目标函数精品文档放心下载精品文档放心下载谢谢阅读误判的概率更小，为了计算简单，我们只检查两次，检查三次建模方法类似），精品文档放心下载精品文档放心下载精品文档放心下载长度，这样，每个零件的损失费用为：总的损失费用除以周期的平均长度(在问精品文档放心下载题三中，由于我们只求解检查两次的情况，所以的取值为2)。谢谢阅读刀具未发生故障而采取预防性更换的损失费用为：刀具发生故障引起的损失费用为：刀具更换的平均周期为：一个周期内的平均损失为：所以问题三的目标函数为：单个零件的平均损失费用7.1.2综上所述，得到问题三的最优化模型7.2模型三的求解对问题三建立的模型，我们用Matlab进行积分运算，我们的计算方法同模感谢阅读精品文档放心下载就换刀，所以我们判定和必定落在（10～100）之间，我们以先对u的取值以5精品文档放心下载u谢谢阅读型二的结果如下表（详细结果见附录七）：一个检查间隔内完成的在一个换刀周期内完成单个零件的损失费用精品文档放心下载零件件数，单位（件）的零件件数u单位（元）精品文档放心下载293195.238n=29件，即生产29件零精品文档放心下载=319件，即每生产319感谢阅读C=5.238元。7.3模型三的结果分析u和n,然精品文档放心下载感谢阅读感谢阅读我们认为我们求得的换刀周期u=319n=29感谢阅读C=5.238为最优解。表5模型三的数据分析un282930313155.23905.23865.23865.2388精品文档放心下载3165.23855.23815.23825.2386精品文档放心下载3175.23805.23785.23805.2384感谢阅读3185.23785.23775.23795.2385谢谢阅读3195.23775.23765.23805.2386精品文档放心下载3205.23775.23785.23825.2390感谢阅读3215.23795.23815.23865.2395感谢阅读3225.23825.23855.23915.2401感谢阅读3235.23875.23915.23985.2409精品文档放心下载感谢阅读感谢阅读由问题二中的C=9.2158元减小为C=5.238n由问题二中的27变为感谢阅读29。这样采取模型三工序效益进一步提高，因此模型三有明显的优势。精品文档放心下载8．灵敏度分析为了方便，我们仅对问题一进行灵敏度分析。8.1f=100,其它参数不变我们首先对参数f（其他费用参数可以做同样处理）进行分析，由题目中的200谢谢阅读减小到100，其它参数不变，再对问题一进行求解，得到的结果如下表：感谢阅读表6参数f的灵敏度分析un222324253454.30884.30214.29794.2943谢谢阅读3504.29534.29004.28534.2821谢谢阅读3554.28374.27884.27474.2728谢谢阅读3604.27484.26954.26584.2644感谢阅读3654.26674.26194.25964.2577谢谢阅读3704.26044.25614.25424.2528感谢阅读3754.25574.25284.25054.2496感谢阅读3804.25364.25024.24844.2490感谢阅读3854.25214.24924.24804.2490精品文档放心下载3904.25214.24984.24994.2506精品文档放心下载3954.25374.25194.25254.2537谢谢阅读4004.25674.25634.25664.2583感谢阅读4054.26204.26124.26214.2653感谢阅读fnf谢谢阅读减小，从而检查间隔也由模型一中的18增大到24，一个换刀周期内每个零件的感谢阅读平均费用也减小到了4.2484元。8.2p=1/1140,其它参数不变然后对非刀具故障率p进行考虑，将非刀具故障率调整到1/1140.再对问题谢谢阅读一进行求解，得到的结果如下表：表7参数p的灵敏度分析un131415163505.25845.25225.24585.2442精品文档放心下载3555.25465.24355.23825.2435谢谢阅读3605.24645.23645.23795.2384感谢阅读3655.24555.23655.23325.2347谢谢阅读3705.24005.23215.22995.2382感谢阅读3755.23605.22915.23385.2373谢谢阅读3805.23895.23315.23315.2377感谢阅读3855.23735.23275.23375.2451感谢阅读3905.24275.23345.24125.2480感谢阅读3955.24345.24105.24425.2520谢谢阅读4005.24535.24405.24835.2628谢谢阅读与问题一的结果相比，当非刀具的平均故障率增大时，n减少，每个零件的精品文档放心下载p感谢阅读精品文档放心下载到了5.2291元。9．模型优缺点9.1模型的优点优点一：本文建模思想易于理解,模型可操作性强,有广泛的应用价值。精品文档放心下载优点二：由对已往数据通过概率统计建立的模型,得出的结论对实际的刀具生产谢谢阅读和零件加工的过程都具有比较实用的价值。100次刀具故障记录的完成零件数观察研究及处理验证。刀具故感谢阅读障分布函数服从正态分布精品文档放心下载问题。精品文档放心下载模型利用到的数据准确，可靠，最优。感谢阅读的最优解，对以后工序长期生产有指导性的价值。9.2模型的缺点精品文档放心下载带来较大的影响。感谢阅读型的求解带来的不利因素。10.模型的改进与推广10.1模型的改进谢谢阅读谢谢阅读谢谢阅读谢谢阅读让检查间隔是非固定的,由于刀具寿命服从正态分布，我们考虑可以在开始检查感谢阅读600谢谢阅读化了，但可以在我们的模型基础上进一步完善，在实际生活中也是这样处理的。感谢阅读这样的改进也可以对目标函数做进一步优化。谢谢阅读精品文档放心下载精品文档放心下载感谢阅读样对工厂的效益会进一步增强。从而使该种方案得到优化。精品文档放心下载10.2模型的推广本文着重讨论了一个关于机械零件加工生产过称中的一个随机优化问题，即如何谢谢阅读精品文档放心下载感谢阅读解的基础，再次通过一定的假设和判断建立了一个随机优化模型，并采用了谢谢阅读matlab精品文档放心下载际加工生产中，具有比较重要的实际指导作用。对于实际一般车床加工过程中。精品文档放心下载精品文档放心下载精品文档放心下载感谢阅读感谢阅读多个行业领域，例如各种机械零件更换，汽车轮胎更换领域等。谢谢阅读11.参考文献[1]韩中耕，《数学建模竞赛-获奖论文精选与点评》，北京：科学出版社，谢谢阅读2006.[2]宋来忠，王志明,《数学建模与实验》,北京：科学出版社,2005精品文档放心下载[3]梁国业，廖键平，《数学建模》，北京，冶金工业出版社，2004谢谢阅读[4]盛骤,谢式千,潘承毅,,北京，高等教育出版社，感谢阅读200312.附录附录一：正态分布检验结果数据序号F(Yi)1-F(YN+1-i)安德森-达令检验精品文档放心下载45910.670.6136220.20.84数据62430.430.91零假设是正态分布54240.130.04显著性水平0.05感谢阅读50950.270.27数据个数10058460.150.63平均值60043370.610.22标准偏差196.74880.330.5981590.740.48中间计算505100.750.4AD统计量0.1精品文档放心下载612110.690.56调整了的AD0.52感谢阅读452120.370.46p值0.61434130.760.2982140.980.84检验结果640150.890.62接受零假设742160.640.82565170.640.9706180.930.47593190.70.22680200.480.39926210.550.35653220.20.49164230.390.48487240.70.94734250.360.53608260.920.6428270.150.451153280.20.15593290.140.71844300.290.29527310.390.82552320.820.24513330.490.62781340.770.85474350.970.59388360.930.83824370.460.6538380.280.18862390.620.56659400.330.17775410.560.77859420.810.86755430.660.19649440.720.71697450.150.14515460.370.22628470.740.74954480.780.12771490.510.22609500.510.47402510.530.49960520.780.49885530.880.22610540.260.26292550.780.63837560.860.85473570.290.28677580.810.34358590.140.19638600.230.44699610.830.67634620.440.38555630.820.72570640.40.5484650.170.07416660.410.03606670.150.231062680.380.51484690.760.18120700.180.61447710.710.71654720.290.86564730.850.8339740.550.85280750.40.08246760.470.64687770.060.3539780.520.61790790.510.8581800.650.45621810.610.52724820.780.3531830.530.07512840.10.36577850.180.36496860.380.11468870.160.02499880.80.24544890.540.63645900.440.31764910.60.25558920.520.26378930.410.67765940.780.39666950.370.85763960.730.73217970.960.87715980.090.57310990.160.88511000.390.33841153[84,1153]上，感谢阅读现取区间[83.5,1153.5]，它能覆盖区间[84,1153]所有的数据。我们将区间感谢阅读[83.5,1153.5]等分为10个小区间，小区间的长度记为,=（1153.5-83.5）精品文档放心下载/10=107，称为组距，小区间的端点称为组限。数出落在每个小区间内的数据的精品文档放心下载频数fi，算出频率fi/n(n=100,i=10)如下表所示：感谢阅读附录二：各个区间数据频率图组限频数fi频率fi/n累计频率83.5——190.530.030.03190.5——297.540.040.07297.5——404.570.070.14404.5——511.5160.16030511.5——618.5250.250.55618.5——725.5200.200.75725.5——832.5130.130.88832.5——939.570.070.95939.5——1046.530.030.981046.5——1153.520.021buel=0，表示不拒绝零假设，说明提出的假设“均值感谢阅读为600”是合理的；②95%的置信区间为［560..0155］，它完全包括均值600，精品文档放心下载且精度很高；③sig的值为1，远超过0.5，不能拒绝零假设。精品文档放心下载附录三：问题一的结果un161718192904.854.414.14.332954.244.174.674.913004.334.824.284.683054.94.024.074.763104.24.064.614.383154.934.074.214.663204.14.94.694.873254.324.144.454.193304.874.544.94.673354.394.934.284.023404.634.34.994.873454.884.034.144.413504.924.784.564.353554.514.644.374.873604.374.194.74.23654.144.524.454.563704.684.984.844.533754.714.954.034.613804.714.674.194.593854.934.274.464.993904.514.754.114.23954.684.724.934.934004.624.814.494.45附录四：模型一所用的源程序clearclci=0;h=0;p=10.*ones(400,30);f=200;t=10;d=300;k=1000;q=1/11400;foru=100:5:800i=i+1;h=0;forn=10:30h=h+1;Ca=(2100+100*n).*quad(@myfun1,0,u)+1000+10*u./n-10.*quad(@myfun3,0,u)./感谢阅读n;Ta=u-quad(@myfun3,0,u);T0=(1-(1-q).^u)/q;C0=0;forj=1:u-1C0=C0+(j.*t./n+d+f.*(n+1)/2).*(1-q).^(j-1).*q;精品文档放心下载endC0=C0+(floor(u/n).*t+k).*(1-q).^(u-1);谢谢阅读T=T0.*[1-(1-q).^u]+Ta.*(1-q).^u;精品文档放心下载C=C0.*[1-(1-q).^u]+Ca.*(1-q).^u;谢谢阅读L=C./T;p(i,h)=L;endendr=min(min(p));[x,y]=find(p==r)u=100+5.*(x-1)%检查的平均间隔n=10+y-1%预防性更换刀具的平均间隔disp('平均每个零件的费用：')disp(r)三个自定义函数：Myfun1.mfunctiony=myfun1(t)y=exp(-(t-600).^2/(2*195.6^2))./(sqrt(2*pi)*195.63);感谢阅读Myfun2.mfunctiony=myfun2(t)y=t.*exp(-(t-600).^2/(2*195.63^2))./(sqrt(2*pi)*195.63);感谢阅读Myfun3.mfunctiony=myfun3(u)y=(*^(1/2)*pi^(1/2)*(erf((131072*^(1/2)*u)/-(0*^(1/2))/)+erf((0*^(1/2谢谢阅读))/)))/;附录五：问题二的结果un262728292909.36499.35479.34779.3437精品文档放心下载2919.32779.31869.31289.3099感谢阅读2929.29619.28829.28369.2819谢谢阅读2939.27019.26349.25999.2595谢谢阅读2949.24939.24399.24179.2425谢谢阅读2959.23369.22959.22879.2307感谢阅读2969.22309.22029.22079.2241感谢阅读2979.21719.21589.21769.2224谢谢阅读2989.21609.21619.21939.2256感谢阅读2999.21959.22109.22589.2335精品文档放心下载3009.22749.23059.23679.2460谢谢阅读3019.23989.24449.25229.2630感谢阅读3029.25649.26269.27209.2844精品文档放心下载3039.27729.28519.29619.3102精品文档放心下载3049.30229.31179.32459.3402精品文档放心下载3059.33129.34259.35709.3744感谢阅读3069.36429.37729.39359.4127精品文档放心下载3079.40109.41599.43409.4551谢谢阅读3089.44189.45859.47859.5015感谢阅读3099.48639.50509.52699.5518谢谢阅读3109.53459.55529.57919.6059精品文档放心下载3119.58649.60919.63509.6639感谢阅读3129.64199.66679.69469.7256精品文档放心下载3139.70109.72799.75809.7910感谢阅读附录六：问题二所用的程序%%question2clearclcf=200;t=10;d=3000;k=1000;w=0.4;e=1500;v=0;p1=0.02;p2=0.4;i=0;H=zeros(1,);精品文档放心下载q=15000.*ones(600,60);foru=1:5:400i=i+1;j=0;forn=10:50j=j+1;forh=0:u-1v=v+1;H(1,v)=H(1,v)+(1+H(1,u-h))*quad(@myfun1,h,h+1);精品文档放心下载endhe=sum(H,2);L2=(k+d.*he./0.95*0.6+t.*u./n+f.*(u./n-he.*0.6)*0.02+f.*he.*n./2*0.6+精品文档放心下载e.*(u./n-he.*0.6).*0.02)./u;感谢阅读q(i,j)=L2-3;endendr2=min(mi