機(jī)械優(yōu)化設(shè)計(jì)上機(jī)報(bào)告-資料下載頁
2025-05-14 01:24本頁面
【正文】 []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = ceq = []c = * ceq = []c = * ceq = []c = * ceq = []c = * ceq = []Local minimum possible. Constraints satisfied.fmincon stopped because the predicted change in the objective functionis less than the default value of the function tolerance and constraints are satisfied to within the default value of the constraint tolerance.stopping criteria detailsActive inequalities (to within = 1e06): lower upper ineqlin ineqnonlin 1 2x = fval = （3）結(jié)果分析：由迭代程序得x=[ 。 。 ],fval =,與理論結(jié)果相等。 4《平面連桿機(jī)構(gòu)中再現(xiàn)已知運(yùn)動(dòng)規(guī)律的優(yōu)化設(shè)計(jì)》上機(jī)實(shí)踐報(bào)告寫出所針對此實(shí)際的優(yōu)化問題，所確定的優(yōu)化目標(biāo)函數(shù)和約束條件。 設(shè)計(jì)一曲柄搖桿機(jī)構(gòu)（如圖9所示），要求曲柄從轉(zhuǎn)到時(shí)，搖桿的轉(zhuǎn)角最佳再現(xiàn)已知的運(yùn)動(dòng)規(guī)律： ， 且已知,為極位角，其傳動(dòng)角允許在范圍內(nèi)變化。MATLAB程序清單：function f=objfun(x)n=30。L1=1。L2=5。fx=0。fa0=acos(((L1+x(1))^2x(2)^2+L2^2)/(2*(L1+x(1))*L2))。 %199。177。179。245。202。188。189。199。pu0=acos(((L1+x(1))^2x(2)^2L2^2)/(2*(x(2))*L2))。 %210。161。184。203。179。245。202。188。189。199。for i=1:n fai=fa0+*pi*i/n。 pui=pu0+2*(faifa0)^2/(3*pi)。 ri=sqrt(L1^2+L2^22*L1*L2*cos(fai))。 alfai=acos((ri^2+x(2)^2x(1)^2)/(2*ri*x(2)))。 betai=acos((ri^2+L2^2L1^2)/(2*ri*L2))。 if fai0amp。amp。fai=pi psi=pialfaibetai。 elseif faipiamp。amp。fai=2*pi psi=pialfai+betai。 end fx=fx+(puipsi)^2。endf=fx。function [c,ceq]=nonlinconstr(x)L1=1。L2=5。m=45*pi/180。n=135*pi/180。c(1)=x(1)^2+x(2)^2(L1L2)^22*x(1)*x(2)*cos(m)。%238。208。161。180。171。182。175。189。199。212。188。202。248。c(2)=x(1)^2x(2)^2+(L1+L2)^2+2*x(1)*x(2)*cos(m)。%238。180。243。180。171。182。175。189。199。212。188。202。248。ceq=[]。 options=optimset(39。LargeScale39。,39。on39。,39。TolFun39。,1e6)。 [x,fval]=fmincon(@objfun,x0,a,b,[],[],lb,ub,@nonlinconstr,options)（2）運(yùn)行結(jié)果：x0=[。]。lb=[1。1]。ub=[10。10]。a=[1 1 。1 1 。1 1 ]。b=[6。4。4]。options=optimset(39。LargeScale39。,39。on39。,39。TolFun39。,1e6)。 [x,fval]=fmincon(@objfun,x0,a,b,[],[],lb,ub,@nonlinconstr,options)Warning: The default trustregionreflective algorithm does not solve problems with the constraints you havespecified. FMINCON will use the activeset algorithm instead. For information on applicable algorithms, seeChoosing the Algorithm in the documentation. In fmincon at 486 No feasible solution found.fmincon stopped because the size of the current search direction is less thantwice the default value of the step size tolerance but constraints are not satisfied to within the default value of the constraint tolerance.stopping criteria detailsx = fval = 46 / 4
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1