Fortran Coder

查看: 11166|回复: 3
打印 上一主题 下一主题

[求助] 关于复杂函数,,如何用过程写

[复制链接]

40

帖子

8

主题

0

精华

熟手

F 币
91 元
贡献
105 点
跳转到指定楼层
楼主
发表于 2014-5-22 15:25:37 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
[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的子程序,子程序要怎么写,希望知道下语句,谢谢各位
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

2033

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1641 元
贡献
709 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

沙发
发表于 2014-5-22 15:55:52 | 只看该作者
使用 function 会更方便。
请阅读教科书,关于子程序与函数的章节

66

帖子

5

主题

2

精华

版主

院士级水师

F 币
481 元
贡献
273 点

管理勋章帅哥勋章爱心勋章规矩勋章

QQ
板凳
发表于 2014-6-9 22:08:27 | 只看该作者
围观大神恢复,顺便学点经验
科研穷三代,读博毁一生

136

帖子

3

主题

0

精华

版主

F 币
1964 元
贡献
1677 点

帅哥勋章管理勋章爱心勋章新人勋章热心勋章元老勋章

地板
发表于 2014-6-10 09:05:33 | 只看该作者
岸边的鱼 发表于 2014-6-9 22:08
围观大神恢复,顺便学点经验

楼上在灌水么?
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-24 00:59

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表