【正文】 。n+4(k)f(Xn+4(k)f(Xl(k)? Xh(k)163。189。n+3(k)????11 2() ()21[ ]?nk kj njfXfX?? ??? ???? ?? ??X*161。214。l(k)163。187。f(Xl(k)205。163。 214。185。Xh(k)163。189。n+3(k)199。243。209。185。203。245。181。227。 Xn+5(k)163。189。n+2(k)163。171。166。194。168。Xh(k)163。n+2(k)163。169。188。198。203。227。f(n5(k)197。208。182。207。209。185。203。245。181。227。202。199。241。212。218。191。201。208。211。242。196。218。yi(Xn+5(k)0 i=1,2161。163。161。m 166。194。163。189。161。193。166。202。213。203。245。181。165。180。191。208。206。Xj(k)163。189。l(k)163。171。p168。Xj(k)163。l(k)163。169。f(Xn+5(k)f(Xh(k)197。208。182。207。202。213。203。245。181。165。180。191。208。206。196。182。165。181。227。202。199。241。212。218。191。201。211。242。218。 yi(Xn+5(k)0i=1,2161。163。161。mXh(k)163。189。n+5(k) 166。194。163。189。161。193。166。k:=+1 205。188。42181。165。180。191。208。206。168。181。196。188。198。203。227。179。204。208。242。191。205。188。TFF T TT F FTFTFF TTFFTFT180。185。 編寫和調試程序的一些注意點編寫程序之前必須熟悉整個的所寫程序的數(shù)學理論和程序流程圖，對于這些東西要吃透，吃深，對整個的流程的原理問題也要理解，不能對著流程圖編程，因為在本次設計中我發(fā)現(xiàn)好多教材資料給出的流程圖也是有問題的，不理解的套，是肯定不可以的。在大的編程過程中涉及的變量是相當?shù)亩?，而且每個變量都代表一定的意思，故而在編寫程序中很容易忘記或寫錯變量名。所以對變量名的選擇要容易識別。在本次設計中目標函數(shù)和約束條件的表達式是非常的多，而且每個表達式的也非常的長，為了可讀性的要求，為了方便檢查錯誤，故而寫這些表達式的代碼時要注意把重復使用的表達式模塊利用函數(shù)或利用新的變量來代替，對利用 Matlab 擬合出的曲線則直接做成函數(shù)的形式，總之，每個表達式不要太長，不要太復雜。在調試程序中，有兩種錯誤，一種是語法錯誤，一種是邏輯錯誤。語法錯誤的更正要看我們自己語言的掌握情況，我們主要檢查與語言規(guī)則不符合的部分；邏輯錯誤要重要檢查數(shù)學模型的建立和程序流程圖的正確性；有時問題出現(xiàn)在變量賦值順序，語句先后順序的不正確。調試程序的常用方法是采用打斷點的方法，就是對每步進行檢查，利用輸出語句顯現(xiàn)每步的執(zhí)行結果，利用輸入語句來中斷程序的運行，以便檢查。程序的調試對于初學者來說是痛苦的，在本次畢業(yè)設計中程序調試工作占據(jù)了大部分的工作時間，有時我發(fā)現(xiàn)程序幾乎陷入一種無法調試的困境中，還在有同學和張彩麗老師的支持和協(xié)助解決一些問題，使得我不斷戰(zhàn)勝困難，堅持到最后調試成功。陜西科技大學畢業(yè)設計（論文）說明書 第 44 頁結 論1. 對減速器進行了一般設計，熟悉了目前企業(yè)中采用的設計方法，掌握了調研，查資料，使用工具書的技能；2. 在分析研究的基礎上，確定了目標函數(shù)，約束條件，建立了減速器的數(shù)學模型；3. 利用 Matlab 把數(shù)據(jù)圖表轉化為數(shù)學表達式；4. 利用優(yōu)化方法，編寫優(yōu)化計算機程序，并調試得出優(yōu)化結果； 5. 在輸入功率為 ，傳動比為 ,輸入端的轉速為 1000r/min 的條件下通過普通的設計得到的結果為：Xarray[0]=50。Xarray[1]=100。Xarray[2]=20。Xarray[3]=21。Xarray[4]=。Xarray[5]=。Xarray[6]=30。Xarray[7]=42。Xarray[8]=48。Xarray[9]=70。Xarray[10]=265。Xarray[11]=6。其體積為：10 7，其中數(shù)組 Xarray 分別依次對應變量 1 到 12 的設計結果值。優(yōu)化設計的運行的結果為：減速器的體積為funcperfect=+007各設計變量的取值為Xperfect[0]=Xperfect[1]=Xperfect[2]=Xperfect[3]=22Xperfect[4]=Xperfect[5]=Xperfect[6]=Xperfect[7]=Xperfect[8]=Xperfect[9]=Xperfect[10]=Xperfect[11]=其中 funcperfect 是優(yōu)化的體積結果，數(shù)組 Xperfect 分別依次對應變量 1 到 12 的優(yōu)化值。通過對比發(fā)現(xiàn)優(yōu)化前和優(yōu)化后體積相差 107，即優(yōu)化后體積在原基礎上減少了 %參考文獻[1] 程序設計[M]. 北京：[2] 程序設計與應用 [M].北京：[3] ++程序設計教程[M]. 北京：[4] 徐曉剛，高兆法， C++ 入門與提高[M].北京：[5] [M].西安：西北工業(yè)大學出版社 .2022[6] [M].北京：高等教育出版社 .1999[7] 濮良貴，[M]. 北京：高等教育出版社 .2022[8] 李德志，[M]. 西安：西北大學出版社 .1993[9] 孫桓，[M]. 北京：[10] 龔溎義，陳式樁， [M].北京：機械工業(yè)出[11] [M]. 北京：[12] [M]. 北京：[13] [M]. 北京：陜西科技大學畢業(yè)設計（論文）說明書 第 46 頁致 謝本次畢業(yè)設計是在張彩麗老師的悉心指導和耐心幫助下完成的，在整個設計過程中張老師付出了辛勤的勞動，在此我對她致以崇高的敬意。本次畢業(yè)設計也得到其它老師的指導和幫助，在此對在百忙之中給予我?guī)椭睦蠋?，李茜教授，陳樺教授，張云云教授致以誠摯的謝意；同時對我的好友王濤，鄭新剛，黃全明等的大力幫助表示感謝。最后感謝我的母校給我提供的求學機會，使我度過了四年寶貴的時光。這四年中，通過各位任課老師知識傳授，使我學到了不少的東西，這將使我在以后的學習和工作中處理問題的方法和經(jīng)驗更加豐富。附 錄：程序源代碼include include include double Kbada(double b)。double Kv(double v)。double Ysa(double x)。double Yfa(double x)。double func(double * Xarray,double I,double P,double N,double r)。double * PerfectDancunxing(double * Xarray,double I,double P,double N,double r)。double * yueshutiaojian( double * Xarray ,double I,double P,double N)。double yuanfunc(double * Xarray,double I)。const int n=12。//未知數(shù)的個數(shù)void main(){double I=。//傳動比double N=1000。//轉速double P=。//輸入功率double r=。//初始懲罰因子double eb=1。//未知數(shù)的中止準則double eb1=。//函數(shù)的中止準則double C=。//遞減因子int k。//循環(huán)次數(shù)int i。//列數(shù)循環(huán)標志double Xperfect[n]。//最終優(yōu)化參數(shù)返回值double * min。//調用單純形優(yōu)化設計方法的返回值double Xmin[n]。//單純形的優(yōu)化的參數(shù)值double funcmin。//單純形法的優(yōu)化函數(shù)值結果double funcperfect。//最終優(yōu)化體積結果double Xarray[n]。//初始可行點////////////////////////////////////////給出初始可行點Xarray[0]=50。Xarray[1]=100。Xarray[2]=20。Xarray[3]=21。陜西科技大學畢業(yè)設計（論文）說明書 第 48 頁Xarray[4]=。Xarray[5]=。Xarray[6]=30。Xarray[7]=42。Xarray[8]=48。Xarray[9]=70。Xarray[10]=265。Xarray[11]=6。coutfunc(Xarray,I,P,N,r)=endlfunc(Xarray,I ,P,N,r)endl。coutyuanfunc(Xarray,I)=endlyuanfunc(Xarray,I)endl。double * gg。gg=yueshutiaojian( Xarray , I, P, N)。coutyueshutiaojian=endl。for(int j=0。j29。j++){coutgg[j]=gg[j] 。}double che。coutendlpleaseinput any key!endl。cinche。///////////////////////////////////判斷初始點是否在可行域內(nèi)double * g。if ((g=new double[29])==NULL){coutcan39。t allocate more memory for g,terminate.\n。 exit(1)。}g= yueshutiaojian( Xarray , I, P, N)。for(i=0。i29。i++){if(g[i]0){cerrthe ith beginning nummber is not right。please check it!endl。}}k=1。//////////////////////////////////////////////////////////////////單純形法求解最優(yōu)點step1:min=PerfectDancunxing( Xarray, I, P, N, r)。funcmin=min[n]。for(i=0。in。i++){Xmin[i]=min[i]。}coutfunc(Xmin, I,P, N, r)=func(Xmin, I,P, N, r)endl。coutfunc(Xarray, I,P, N, r)=func(Xarray, I,P, N, r)endl。coutyuanfunc(Xmin,I)=yuanfunc(Xmin,I)endl。coutyuanfunc(Xarray,I)=yuanfunc(Xarray,I)endl。cout當 k＝k時:endl。coutfuncmin=funcminendl。for(i=0。in。i++){coutXmin[i]=Xmin[i]endl。}////////////////////////////////////////////////////////////////檢驗迭代中止準則double Mozhi=0。for(i=0。in。i++){Mozhi+=pow((Xmin[i]Xarray[i]),2)。}double MMozhi=sqrt(Mozhi)。cout精度fabs(yuanfunc(Xmin,I)yuanfunc(Xarray,I)) MMozhiendl。if((fabs(yuanfunc(Xmin,I)yuanfunc(Xarray,I))=eb)amp。amp。(MMozhi=eb1)){for(i=0。in。i++)Xperfect[i]=Xmin[i]。funcperfect=yuanfunc(Xmin,I)。goto end。}陜西科技大學畢業(yè)設計（論文）說明書 第 50 頁else{for(i=0。in。i++)//上一步的優(yōu)化結果為下一步的初始點{Xarray[i]=Xmin[i]。}r=C*r。cout總的第k次循環(huán)次數(shù)endl。k++。goto step1。}///////////////////////////////////////////////////////////結束標志end:cout優(yōu)化體積值endl。cout funcperfect=funcperfectendl。cout各個優(yōu)化變量的優(yōu)化值為endl。for(i=0。in。i++){coutXperfect[i]=Xperfect[i]endl。}cout最終精度值fabs(yuanfunc(Xmin,I)yuanfunc(Xarray,I)) MMozhiendl。cout總的循環(huán)次數(shù) k=kendl。return。}////////////////////////////////////////////////////////////單純形優(yōu)化子函數(shù)double * PerfectDancunxing (double * Xarray,double I,double P,double N,double r){double h=。//構成單純形定點的步長double eb=1。//中止準則double erfa=1。//發(fā)射系數(shù)double bada=。//壓縮系數(shù)do