本科第二组的文章：自动化车床管理

石伟、向中辉、喻欢摘要本文建立的是自动化车床单刀具连续加工零件工序定期检查和刀具更换的感谢阅读法，确定了刀具寿命服从600,196.62)的正态分布。为了使总的期望损失达到感谢阅读精品文档放心下载模型。针对问题一：我们将检查间隔和道具更换策略的问题确定为单个零件期望感谢阅读感谢阅读隔的分布函数F(t)。然后列出以单个零件的期望损失为目标函数，关于检测间谢谢阅读t隔和刀具定期更换间隔为变量的目标函数方程建立了一个单目标的期望值模型。

最后，利用计算机采用穷举搜索法求得模型一的最优解为每生产一次，定期更换刀具间隔u为，相应的单个零件的最小费用C为4.595元。精品文档放心下载谢谢阅读针对问题二：我们采用分摊法建立了单个零件的效益函数。首先，求出在感谢阅读举法求得检查间隔n为，换刀间距u为为问题二的最好检查间隔和换刀策略，相应的单个零件的最小费用C为元。谢谢阅读针对问题三：在问题二的建模基础上，我们增加每次检查零件的次数来减感谢阅读内平均损失除以周期长度的单目标优化模型。得到的结果为检查间隔n为

29，换刀间距u为为最优解，相应的单个零件的最小费用C为5.238元。感谢阅读最后，针对问题一的结果，对部分参数进行了灵敏度分析。感谢阅读关键词：卡方拟合正态分布检验穷举法单目标期望感谢阅读1.问题重述1.1问题背景量适用于各种专业技术规X能的X的好坏，便于对机器进行修理，减少不必要的损失，以提高经济效益。感谢阅读1.2问题提出其它故障仅占的,100零件数如附表1。现计划在刀具加工一定件数后定期更换新刀具。精品文档放心下载、故障时产出的零件损失费用f=200元/件；

、进行检查的费用t=10元/次；、发现故障进行调节使恢复正常的平均费用d=3000元/次(包括刀具费)；

、未发现故障时更换一把新刀具的费用k=1000元/次。谢谢阅读精品文档放心下载1.3需要解决的问题,②如果该工序正常时产出的零件不全是合格品，有为不合格品；而工序感谢阅读故障时产出的零件有为合格品，60%为不合格品。工序正常而误认有故障仃机产生的损失费用为1500策略。感谢阅读③在②的情况,可否改进检查方式获得更高的效益。

表1100次刀具故障记录(完成的零件数)精品文档放心下载3625425844338154529827425655936534876084281153593552781388824862859649515628771960610837473358634570416606484654339246687790724512496468544558765763217310精品文档放心下载精品文档放心下载精品文档放心下载感谢阅读谢谢阅读精品文档放心下载谢谢阅读谢谢阅读精品文档放心下载感谢阅读2.问题分析2.1问题一的分析失费用有刀具在未发生故障而采取预防性更换的损失费用和刀具发生故障而带设备的工作，到了定期更换刀具的时期，即使设备未出现故障，也要更换刀具。

件出现故障的几率相同且相互独立，对于刀具更换周期u来说，因为它是s,n的感谢阅读目标函数为：谢谢阅读=周期的区间，再利用有穷列举法得出最优化组合，则认为该组合为问题一的解。谢谢阅读2.2问题二的分析刀具正常产出的产品是合格品的为求出目标函数。谢谢阅读中图12.3问题三的分析查次数作考虑，可以对零件进行检查两次，也可以检查三次，实际上，由于时间论检查周期的长短都要检查n个零件，但是实际上，随着检查周期的不断缩短，

检查的零件数n间隔，采取不等距检查的策略来降低损失费用，即在稳定期内检查时间间隔大，着换刀周期递减的。感谢阅读3.模型假设与符号说明3.1模型的假设（）假定在生产任一零件时出现故障的机会均相同；

（）假设工序出现故障是完全随机的；（）假设题目所给的数据真

实可靠，且不考虑其他的费用；（）假设工序在发生故障后，无法

判断是刀具故障还是非刀具故障；（）假设非刀具故障在任一零件

处发生的概率均相同；（）假设从换上新刀到下一次工序故障或预

防性换刀构成一个周期；感谢阅读3.2符号的说明符号u预防性换刀的固定间隔(刀具更换周期)u刀具更换周期的长度的平均值n检查的固定间隔s预防性换刀对的间隔中检查的次数h一个周期内合格零件件数件f故障时产出的每件零件损失费用f元/精品文档放心下载t每次进行检验的费用t元/次d每次发现故障进行调节使恢复正常的平均费用（包含刀具费）d3000元/次感谢阅读k每次未发现故障时更换一把新刀具的费用k1000元/次谢谢阅读40%b刀具出现故障后生产合格零件的概率b谢谢阅读a刀具没有出现故障时生产合格零件数a98%g工序正常而误认为有故障停机产生的损失费用g次谢谢阅读T一个周期内总的故障间隔C每个零件的平均费用p工序中每一零件处非刀具故障率C0一个周期内平均费用C1在一周期内未出现故障所产生的费用C2一周期内出现故障所产生的费用T0出现非刀具故障在一周期内所产生的零件数T1出现非刀具故障时在一周期内所生产的零件数f(u)u的概率密度函数（）间隔期（）刀具故障发生次数的数学期望感谢阅读F(u)uq工序故障是由刀具故障引起的的概率Cd刀具未故障引起的损失Cf由刀具发生故障引起的损失4.数据处理与假设检验由附表1的谢谢阅读略了解这些数据的分布情况，我们先根据所给数据画直方图：精品文档放心下载252015数次的围范定指1050200600800图2刀具故障时已生产的零件个数的频数直方图次刀具故障记

录的完成零件数观察研究及用Excel服从正态分布x~N(,2)根据记录数据求出了感谢阅读谢谢阅读111(x)2xx600,(xx)f(x)e谢谢阅读22222iii1i1发现这次刀具故障时完成的零件数近似服从x~N(600,196.622)精品文档放心下载Excel检验和卡方优度检验，经软件处理得到如下表格谢谢阅读表2雅克-贝拉检验和卡方拟合优度检验结果感谢阅读精品文档放心下载9精品文档放心下载p值p值谢谢阅读雅克-贝拉检验综合考虑了正态分布的宽度和偏斜对分布的影响，由结果知统计量为0.813878402且=0.665684663>显著性水平0.05，即假设正确。而卡方拟合优度检验得结果=0.980290368>显著性水平0.05,进一步证明了假设的正确性。谢谢阅读综上检验我们接受原假设，即表1N(600,196.62)的正态分布。感谢阅读y0.9970.990.980.950.900.75率0.50概0.250.100.050.020.010.00310020030040050060070080090010001100图3正态性分布检验图所以，根据以上分析，我们可以得到一下表达式：刀具寿命的概率密度函数为f(x)12e所以刀具寿命的发布函数为F(x)12x222x2xe2205.问题一的解答针对问题一我们建立了模型一。5.1模型一的建立5.1.1确定目标函数首先求非刀具故障和刀具故障总的分布函数：谢谢阅读设T服从几何分布其分布函数为：111F(T)p(1p)m1,(0p1),E)精品文档放心下载1p1m1由题意知，刀具损坏故障占工序故障的，而其它故障占，近似求得：精品文档放心下载E)1求得1E(T)p0p所以总的故障间隔T,T),T感谢阅读12F(t)1(1FF))ttt精品文档放心下载01tt111p(1p)k1精品文档放心下载m1t11(1p)t然后求一个周期内平均损失费用对于模型一，可分为两种情况进行考虑第一种情况：设每检查n个零件检查一次，第s次检查换刀具。若第s次检查零件仍然合格，则前面生产的零件都合格，即工序正常，这时的费用为：感谢阅读(stk)PTuC1第二种情况：设每检查n个零件检查一次，第m次检查零件不合格，则设故谢谢阅读障出现在第(m-1)n+i个，这时费用记为C。C2nTs=dfi1P(m1)nTP(mTm1…………第s次检查预防性换刀一个周期图4模型一求解图最后求刀具更换周期：对于刀具更换周期u来说，因为它是s,n的函数，所以也是随机变量，求其谢谢阅读数学期望值，即得刀具更换的平均周期：uE(s,n)udF(u)snf(u)dusn感谢阅读0snsnF(sn)F(u)dusn1F(sn)sn感谢阅读0sn1F(u)du0所以问题二的目标函数为CC1uC25.1.2综上所述，得到问题一的单目标最优化模型minCCC12uC(stk)PT

u1nTs=dfC感谢阅读i1P(mTP(mT2m15.2模型一的求解对于模型一，我们用matlab编程进行穷举法进行搜索，由刀具寿命呈正态分布

600u应

在600s和nX围.谢谢阅读我们判定s和nu的取值以5为步长计算费用。感谢阅读一个检查间隔内完成的在一个换刀周期内完成单个零件的损失费用精品文档放心下载零件件数n，单位（件）的零件件数uC单位（元）精品文档放心下载4.5955.3模型一结果分析我们对模型一中部分结果用如下表格表示：由表3可以看出u和n的取值都会制约C的取值，当u和n和18时，目标函数的值最小，而当u和n的取值越接近这个值时，目标函

数的值越接近这个最小值，所以我们认为=4.595为模型一的最优解。精品文档放心下载表3问题一的数据分析un194.6214.6104.6014.6024.6174.6014.5994.5984.6114.5944.5974.5964.6074.5974.5954.5964.6024.5984.5964.5984.6044.5994.5994.6024.6114.6024.6034.607精品文档放心下载感谢阅读谢谢阅读精品文档放心下载精品文档放心下载精品文档放心下载感谢阅读结果说明：当工序设计换刀周期为生产360个零件，检查周期为个零件时，所获得的效益最好，即生产单位合格零件的消耗费用最少。感谢阅读6.问题二的解答针对问题二，我们建立了模型二。6.1模型二的建立6.1.1确立目标函数一个换刀周期内换刀次数的期望值。精品文档放心下载u感谢阅读时，设m为第一次故障发生在第i个零件时.iE(u)1E(uiu1

i1f(t)dt(uii0然后对可能会发生的不同情况，在模型一的基础上建立两个过渡模型：感谢阅读过渡模型一：假设工序处于正常状态下，次品率为，而当工序不正常时，生产的零件全精品文档放心下载感谢阅读检查费、换零件损失费。则每个零件的平均费用为：Cu1d[kE(u)unutfE(u)[E)](1a)(gf)]/uqn2n过渡模型二：假设工序处于正常状态时产品全为正品，而当工序处于不正常状态下，次品率为，有则有平均每个零件的费用包含四个部分：预防刀具费用、故障调整费、检查费、换零件损失费，所以感谢阅读Cu2[kdE(u)(1b)fE(u)n(1b)]/uqn2综合过渡模型一和过渡模型二，问题二的模型为：CC1C2C3其中C1故障调节恢复费用和检查费用：C2C3dE(u)(1b)Ct1+n为坏零件损失费：uCf2E(u)(1b)a)/un为误停机造成的损失费：uCg[E(u)(1b)](13a)/u

n所以，模型二的目标函数为minCCCC1236.1.2综上所述，得到问题二的最优化模型minCCCC123感谢阅读dE(u)(1b)Ct+n1CufE(u)(1b)a)/u谢谢阅读n2uCg[E(u)(1b)](1a)/u精品文档放心下载3n6.2模型二的求解精品文档放心下载致分布X围，用Matlab进行计算，最终求得的主要结果如下表（详细结果见附精品文档放心下载一个检查间隔内完成的在一个换刀周期内完成单个零件的损失费用零件件数n，单位（件）的零件件数u感谢阅读C单位（元）9.2158=2727件零感谢阅读u=297297精品文档放心下载用9.2158元。6.3模型二的结果分析针对模型二，我们对部分结果进行分析（见表u固定时，目标函数C的

会发生变化，说明检查间隔制约着最小费用。同理，当n固定时，目标函数C的

取值会随着u的改变而发生变化，说明换刀周期也对目标函数产生影响。所以，

只有当u和n=27,=297

时，目标函数C的取值最小为9.2158元。精品文档放心下载表4模型二的数据分析un289.23369.22959.22879.23079.22309.22029.22079.22419.21719.21589.21769.22249.21609.21619.21939.22569.21959.22109.22589.23359.22749.23059.23679.24609.23989.24449.25229.2630精品文档放心下载谢谢阅读谢谢阅读谢谢阅读感谢阅读精品文档放心下载谢谢阅读7.问题三的解答针对问题三，我们建立了模型三。7.1模型三的建立7.1.1确立目标函数长度，这样，每个零件的损失费用为：总的损失费用除以周期的平均长度(在问

题三中，由于我们只求解检查两次的情况，所以x的取值为2)。感谢阅读刀具未发生故障而采取预防性更换的损失费用为：xutxutxut

C(k)a2(h)fdga)2h)fddg]a(1a)nnn刀具发生故障引起的损失费用为：xutxutxut

C(k2((uh)fdb)2(uh)ffdb)nnn刀具更换的平均周期为：uE(s,n)udF(u)snf(u)dusn感谢阅读0snF(sn)F(u)dusn1F(sn)sn感谢阅读0sn1F(u)du0一个周期内的平均损失为：CF(u)所以问题三的目标函数为：C0fC单个零件的平均损失费用minC0ud(1F(u))7.1.2综上所述，得到问题三的最优化模型CCfF(u)Cdu(1F(u))7.2模型三的求解对问题三建立的模型，我们用Matlab进行积分运算，我们的计算方法同模就换刀，所以我们判定s和n必定落在（～u的取值以5为步长计算费用，确定最小费用的Xu的取值精度进一步提高，得到精品文档放心下载一个检查间隔内完成的在一个换刀周期内完成单个零件的损失费用零件件数n，单位（件）C单位（元）谢谢阅读5.23829件零感谢阅读u=319用=5.238元。谢谢阅读7.3模型三的结果分析对于模型三的结果，我们采用和模型二同样的分析方法。分别固定u和,感谢阅读此，我们认为我们求得的换刀周期=319，检查间隔费用=5.238为最优解。谢谢阅读表5模型三的数据分析un5.23905.23865.23865.23855.23815.23825.23805.23785.23805.23785.23775.23795.23775.23765.23805.23775.23785.23825.23795.23815.23865.23825.23855.23915.23875.23915.2398精品文档放心下载谢谢阅读谢谢阅读谢谢阅读精品文档放心下载精品文档放心下载精品文档放心下载谢谢阅读精品文档放心下载感谢阅读感谢阅读由问题二中的=9.2158元减小为5.238n由问题二中的27变为感谢阅读。这样采取模型三工序效益进一步提高，因此模型三有明显的优势。谢谢阅读．灵敏度分析为了方便，我们仅对问题一进行灵敏度分析。8.1=100,其它参数不变我们首先对参数（其他费用参数可以做同样处理）进行分析，由题目中的减小到100，其它参数不变，再对问题一进行求解，得到的结果如下表：

表6参数f谢谢阅读un4.30884.30214.29794.29434.29534.29004.28534.28214.28374.27884.27474.27284.27484.26954.26584.26444.26674.26194.25964.25774.26044.25614.25424.25284.25574.25284.25054.24964.25364.25024.24844.25214.24924.24804.24904.25214.24984.24994.2506感谢阅读谢谢阅读感谢阅读谢谢阅读精品文档放心下载感谢阅读感谢阅读谢谢阅读感谢阅读谢谢阅读4.25374.25194.25254.25374.25674.25634.25664.25834.26204.26124.26214.2653精品文档放心下载谢谢阅读感谢阅读与问题一的结果相比，f减小时，n增大，f减小后，产生的坏零件的损失

费减小，从而检查间隔也由模型一中的18增大到，一个换刀周期内每个零件

的平均费用也减小到了4.2484元。谢谢阅读8.2=1/1140,其它参数不变然后对非刀具故障率p进行考虑，将非刀具故障率调整到1/1140.再对问题一进行求解，得到的结果如下表：感谢阅读表7参数p的灵敏度分析un145.25845.25225.24585.25465.24355.23825.24645.23645.23795.24555.23655.23325.24005.23215.22995.23605.23385.23895.23315.23315.23735.23275.23375.24275.23345.24125.24345.24105.24425.24535.24405.2483感谢阅读精品文档放心下载精品文档放心下载精品文档放心下载谢谢阅读感谢阅读精品文档放心下载精品文档放心下载谢谢阅读精品文档放心下载感谢阅读n减少，每个零件的平均成本降低，这是由于p生到了5.2291元。感谢阅读．模型优缺点9.1模型的优点优点一：本文建模思想易于理解,模型可操作性强,有广泛的应用价值。

,得出的结论对实际的刀具生产

和零件加工的过程都具有比较实用的价值。感谢阅读感谢阅读100谢谢阅读障分布函数服从正态分布N(,)模型比较容易理解，也能反映实谢谢阅读际生产中的问题。谢谢阅读模型利用到的数据准确，可靠，最优。感谢阅读的最优解，对以后工序长期生产有指导性的价值。9.2模型的缺点带来较大的影响。感谢阅读谢谢阅读型的求解带来的不利因素。10.模型的改进与推广模型的改进数进行分析，分别求出目标函数。然后进行比较，得到最优解。然后针对检查次

让检查间隔是非固定的,由于刀具寿命服从正态分布，我们考虑可以在开始检查

化了，但可以在我们的模型基础上进一步完善，在实际生活中也是这样处理的。这样的改进也可以对目标函数做进一步优化。谢谢阅读样对工厂的效益会进一步增强。从而使该种方案得到优化。精品文档放心下载模型的推广本解的基础，再次通过一定的假设和判断建立了一个随机优化模型，并采用了matlab程序编写主要模型函数，得出了一组比较合理的最终结果，在机械零件于８的零件来说，本文的指导意义不大，具有一定的局限性。另一方面，由于我多个行业领域，例如各种机械零件更换，汽车轮胎更换领域等。谢谢阅读11.参考文献感谢阅读[2]宋来忠，王志明,《数学建模与实验》,：科学,2005感谢阅读精品文档放心下载[4]盛骤,谢式千潘承毅,《概率论与数理统计》第三版,，谢谢阅读12.附录附录一：正态分布检验结果数据F(Y)1-F(Y)安德森-达令检验i45910.0043422670.00245866136220.007320420.009396784数据感谢阅读62430.0132988430.026023591零假设是正态分布谢谢阅读54240.0257177130.033560904显著性水平0.05精品文档放心下载50950.0359032270.035903227数据个数10058460.0518230150.048664463平均值600感谢阅读43370.0586274610.073608922标准偏差196.6291695谢谢阅读74880.0701254330.09135495981590.0921928740.093885748中间计算谢谢阅读5050.1092094750.10088684AD统计量0.22273681谢谢阅读6120.1130626690.107318656调整了的0.224457452精品文档放心下载4520.1294428370.114040946p值0.823045761精品文档放心下载4340.1404785760.12730982谢谢阅读9820.1569741980.137102784检验结果谢谢阅读6400.1746953890.166950562接受零假设精品文档放心下载7420.1908568640.178651982精品文档放心下载5650.1978531640.18673289感谢阅读7060.1992707930.1922438475930.218250770.2006945226800.2258193480.202124339精品文档放心下载9260.2366605550.203560235精品文档放心下载6530.251010120.215264749精品文档放心下载1640.2591766390.2258193484870.260826270.2350944947340.2776144360.2477820536080.2827518920.26414166谢谢阅读4280.2984324150.27932184511530.303745820.2949138155930.3144971140.3073113718440.3217538290.310895229精品文档放心下载5270.3272415390.329079182精品文档放心下载5520.3290791820.3420564245130.3327668490.3476768627810.3552229770.3685646854740.3628262970.3820667593880.3762614930.391800783感谢阅读8240.3781935460.39375596精品文档放心下载5380.3840078280.401603018精品文档放心下载8620.3878986620.409489956感谢阅读6590.4035711330.419400017谢谢阅读7750.4094899560.423378777谢谢阅读8590.4154294810.431359386谢谢阅读7550.4273653660.443382019精品文档放心下载6490.4293614720.451426871精品文档放心下载6970.4393681150.457473814谢谢阅读5150.4534413370.475668222感谢阅读6280.4615106740.479719674谢谢阅读9540.4675732780.481746212精品文档放心下载7710.4858006510.483773222感谢阅读6090.4858006510.487828447感谢阅读4020.5121715530.514199349精品文档放心下载9600.5162267780.514199349谢谢阅读8850.5182537880.532426722感谢阅读6100.5202803260.538489326感谢阅读2920.5243317780.546558663精品文档放心下载8370.5425261860.560631885精品文档放心下载4730.5485731290.570638528感谢阅读6770.5566179810.572634634感谢阅读3580.5686406140.584570519感谢阅读6380.5766212230.590510044谢谢阅读6990.5805999830.596428867精品文档放心下载6340.5905100440.612101338谢谢阅读5550.5983969820.615992172谢谢阅读5700.606244040.621806454谢谢阅读840.6081992170.623738507精品文档放心下载4160.6179332410.637173703感谢阅读6060.6314353150.644777023感谢阅读1062680.6523231380.667233151感谢阅读4840.6579435760.670920818感谢阅读1200.6709208180.672758461精品文档放心下载4470.6891047710.678246171谢谢阅读6540.6926886290.685502886谢谢阅读5640.7050861850.69625418精品文档放心下载3390.7206781550.701567585谢谢阅读2800.735858340.717248108感谢阅读2460.7522179470.722385564感谢阅读6870.7649055060.73917373精品文档放心下载5390.7741806520.740823361谢谢阅读7900.7847352510.74898988精品文档放心下载5810.7964397650.7633394456210.7978756610.774180652精品文档放心下载7240.7993054780.78174923谢谢阅读5310.8077561530.8007292075120.813267110.8021468365770.8213480180.809143136精品文档放心下载4960.8330494380.8253046114680.8628972160.843025802精品文档放心下载4990.872690180.8595214245440.8859590540.8705571636450.8926813440.886937331精品文档放心下载7640.899113160.890790525精品文档放心下载5580.9061142520.9078071263780.9086450410.929874567精品文档放心下载7650.9263910780.941372539精品文档放心下载6660.9513355370.948176985精品文档放心下载7630.9640967730.9640967732170.9664390960.9742822877150.9739764090.9867011573100.9906032160.99267958精品文档放心下载8511000.9975413390.995657733精品文档放心下载[83.5,1153.5]等分为个小区间，小区间的长度记为,=（1153.5-83.5）

/10=107，称为组距，小区间的端点称为组限。数出落在每个小区间内的数

据的频数fif如下表所示：谢谢阅读附录二：各个区间数据频率图频数fi频率f3477321检验结果显示：①布尔变量，表示不拒绝零假设，说明提出的假设“均

值为是合理的；②95%的置信区间为［560.9845639.0155］，它完

全包括均值600，且精度很高；③sig的值为，远超过0.5，不能拒绝零假设。精品文档放心下载谢谢阅读附录三：问题一的结果un4.8330952854.8155162414.801920614.791679433

4.8011967244.7843162174.7714173674.761871491

4.7719600334.7557962824.7436124284.734780068

4.746152794.7298667024.7184162074.711173176

4.721975924.7073009064.6966022614.689252138

4.7002075934.6863048074.6763766214.669795466

4.681631114.667661394.6585223694.652728887

4.6644605324.6521584144.6438273454.638840319

4.6494869874.6380134544.630509194.626347467

4.6366473394.6260212934.6193627284.616045202

4.6267403634.616980734.6103279994.607873587

4.6179915884.6091173034.6042069144.602634541

4.6112048924.6032348784.5992269564.598555535

4.6071892514.5992826644.5961961374.596444587

4.6041722374.5980674194.595921074.59624802

4.6029663144.5978222524.5966348454.598779356

4.6043857684.5993600984.5991502844.602270853

4.6066637714.6034970954.6034219034.607537061

4.6106169714.6084666674.6102675194.615395659

4.6170651934.6150869274.6179205464.624079899

4.6242421514.6233159754.6272004114.63440902

4.6329690684.6339765724.6380658934.646341593

4.6440703624.6453015814.6513409494.660701345精品文档放心下载谢谢阅读精品文档放心下载谢谢阅读精品文档放心下载谢谢阅读感谢阅读精品文档放心下载精品文档放心下载谢谢阅读精品文档放心下载谢谢阅读精品文档放心下载感谢阅读感谢阅读谢谢阅读谢谢阅读谢谢阅读精品文档放心下载精品文档放心下载感谢阅读感谢阅读感谢阅读附录四：模型一所用的源程序clearclci=0;h=0;p=10.*ones(400,30);

f=200;t=10;d=300;k=1000;

q=1/11400;foru=100:5:800

i=i+1;h=0;forn=10:30

h=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;

谢谢阅读C0=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;谢谢阅读r=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.m精品文档放心下载functiony=myfun2(t)y=t.*exp(-(t-600).^2/(2*195.63^2))./(sqrt(2*pi)*195.63);

Myfun3.m谢谢阅读functiony=myfun3(u)y=(4194304*08498^(1/2)*pi^(1/2)*(erf((131072*08498^(1/2)*u)/54249-

(78643200*08498^(1/2))/54249)+

erf((78643200*08498^(1/2))/54249)))/80293;感谢阅读附录五：问题二的结果un9.36499.35479.34779.32779.31869.31289.29619.28829.28369.27019.26349.25999.24939.24399.24179.23369.22959.22879.22309.22029.22079.21719.21589.2176精品文档放心下载精品文档放心下载谢谢阅读精品文档放心下载精品文档放心下载感谢阅读感谢阅读精品文档放心下载9.21609.21619.21939.21959.22109.22589.22749.23059.23679.23989.24449.25229.25649.26269.27209.27729.28519.29619.30229.31179.32459.33129.34259.35709.36429.37729.39359.40109.41599.43409.44189.45859.47859.48639.50509.52699.53459.55529.57919.58649.60919.63509.64199.66679.69469.7256

9.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,9999999);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);谢谢阅读he=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;精品文档放心下载r2=min(min(q));

[x2,y2]=find(q==r2)

u2=(x2-1).*5+1n2=10+(y2-1)disp('平均每个零件的费用：')disp(r2)附录七：问题三的结果nu5.2526869955.2514062725.2504820025.2498932185.2506541235.2494561485.2486147185.2481088615.2487692195.2476550095.2468974375.2464755265.2470321885.2460027715.2453300865.2449931495.2454429465.2444993625.2439126035.2436616825.2440014235.24314472