[Fortran] 纯文本查看 复制代码 G(17)=45. * 3.14 / 180. - delt1
G(20)=0.85 - TAN(delt1) / ( TAN(beta) + TAN(delt1) )
PI_1=sin(alpha - theta) + f2 * sin(delt2) * cos(alpha - theta)
PI_2=cos(alpha - theta) - f2 * sin(delt2) * sin(alpha - theta) * cos(beta)
A1=PI_2 * COTAN(alpha) + PI_1 * cos(beta)
B1=( f2 * COTAN(alpha) * ( cos(delt2) + sin(delt2) * sin(alpha - theta) * sin(beta) ) - PI_1 * sin(beta) )
C1=( 1 / f1 ) * ( f2 * cos(beta) * ( cos(delt2) + sin(delt2) * sin(alpha - theta) * sin(beta) ) - PI_2 * sin(beta) )
delt1=ASIN( A1 / SQRT( A1**2 + B1**2) ) - Beta - ASIN( C1 / SQRT( A1**2 + B1**2) )
PRINT*,'delt1 = ', delt1
delt2=ATAN( sin(beta) * sin(alpha) / cos(theta) )
PRINT*,'delt2 = ', delt2
请教一下,我建立优化数学模型,其中约束条件里,有G(17)和G(20)两个类型的约束,其中用到delt1这个参数,但是其中delt1,又要一列式子计算得知,
我如何再约束力调用delt1的子程序,子程序要怎么写,希望知道下语句,谢谢各位 |