轉(zhuǎn)向梯形優(yōu)化設(shè)計(jì)matlab程序

1、精選優(yōu)質(zhì)文檔-----傾情為你奉上專心---專注---專業(yè)優(yōu)化計(jì)算 MATLAB 程序首先，將目標(biāo)函數(shù)寫成 M 文件，其程序語(yǔ)句如下；function f = fun (x)global K L thetamax alphafor i=1:61f = 0betae = atan(tan(alpha(i)/(1-(K/L)*tan(alpha(i))));A(i)=2*x(1).^2*sin(x(2)+alpha(i));B(i)=2*K

2、*x(1)-2*x(1).^2*cos(x(2)+alpha(i));C(i)=2*x(1).^2-4*x(1).^2*(cos(x(2)).^2+4*K*x(1)*cos(x(2))-2*K*x(1)*cos(x(2)+alpha(i));theta3(i)= 2*acot((A(i)+sqrt(A(i).^2+B(i).^2-C(i).^2))/(B(i)+C(i)));beta(i)=x(2)+theta3(i)-pi;if al

3、pha(i)=pi/18,alpha(i)<=pi/9;f(i)=abs(betaa(i)-betae3(i));elsef(i)=0.5*abs(beta(i)-betae3(i));global K L thetamax alphaK=inputL=inputthetamax=inputx0(1)=input精選優(yōu)質(zhì)文檔-----傾情為你奉上專心---專注---專業(yè)B( i) = 2* K* x( 1) - 2* ( x (

4、1) ) . ∧2* cos( x( 2) + alpha( i) ) ) ;C ( i ) = 2* ( x ( 1) ) . ∧ 2- 4 * ( x ( 1) ) . ∧ 2*( cos( x( 2) ) ) . ∧2+ 4* K * x ( 1) * cos( x ( 2) ) - 2* K* x ( 1) * cos( x ( 2) + alpha( i) ) ;theta3(i)= 2* acot ( A( i) + sqr