fortran指令大全
Fortran常用函数
1、RANDOM_NUMBERSyntax ['sintæks] n. 语法CALL RANDOM_NUMBER (harvest结果)Intrinsic Subroutine（固有子程序）：Returns a pseudorandom number greater than or equal to zero and less than one from the uniform distribution.返回大于或等于0且小于1，服从均匀分布的随机数2、RNNOA/ DRNNOA （Single/Double precision）Generate pseudorandom numbers from a standard normal distribution using an acceptance/rejection method.产生服从标准正态分布的随机数Usage（用法）CALL RNNOA (NR, R)Arguments（参数）NR— Number of random numbers to generate. (Input) 要产生随机数的个数R— Vector of length NR containing the random standard normal deviates. (Output)输出长度为NR，随机正态分布的向量Comments（注解）The routine RNSET can be used to initialize the seed of the random number generator. The routine RNOPT can be used to select the form of the generator.程序RNSET可以用来初始化随机数发生器的种子ExampleIn this example, RNNOA is used to generate five pseudorandom deviates from a standard normal distribution.INTEGER ISEED, NOUT, NRREAL R(5)EXTERNAL RNNOA, RNSET, UMACHCCALL UMACH (2, NOUT)NR = 5ISEED = 123457CALL RNSET (ISEED)CALL RNNOA (NR, R)WRITE (NOUT,99999) R99999 FORMAT (' Standard normal random deviates: ', 5F8.4)ENDOutputStandard normal random deviates: 2.0516 1.0833 0.0826 1.2777 -1.22603、RESHAPEIntrinsic Function（内部函数）Constructs an array of a specified shape from the elements of another array. 构造规定形式的数组Syntax（语法）result = RESHAPE (source, shape [ , pad][ , order])source(Input) Any type. Array whose elements will be taken in standard Fortran array order (see Remarks), and then placed into a new array.shape(Input) Integer. One-dimensional array that describes the shape of the output array created from elements of source. 描述输出数组的大小的一维数组，The elements of shape are the sizes of the dimensions of the reshaped array in order. If pad is omitted 省略, the total size specified by shape must be less than or equal to source.pad 可选参数(Optional; input) Same type as source. Must be an array. If there are not enough elements in source to fill the result array, elements of pad are added in standardFortran array order. If necessary, extra copies of pad are used to fill the array.order 可选参数(Optional; input) Integer. One-dimensional array. Must be the same length as shape.Permutes the order of dimensions in the result array. The value of order must be a permutation of (1, 2,...n) where n is the size of shape.Return Value（返回值）The result is an array the same data type and kind as source and a shape as defined in shape.ExamplesINTEGER AR1( 2, 5)REAL F(5,3,8)REAL C(8,3,5)AR1 = RESHAPE((/1,2,3,4,5,6/),(/2,5/),(/0,0/),(/2,1/))! returns 1 2 3 4 5! 6 0 0 0 0!! Change Fortran array order to C array orderC = RESHAPE(F, (/8,3,5/), ORDER = (/3, 2, 1/))END4、SUMIntrinsic Function（内部函数）Sums elements of an array or the elements along an optional dimension. The elements summed can be selected by an optional mask.将数组中的元素求和Syntax（语法）result = SUM (array [ , dim] [ , mask])array(Input) Integer, real, or complex. Array whose elements are to be summed.dim 可选参数(Optional; input) Integer. Dimension along which elements are summed.1 ≤dim≤n, where n is the number of dimensions in array.mask 可选参数(Optional; input) Logical. Must be same shape as array. If mask is specified, only elements in array that correspond to .TRUE. elements in mask are summed.Return Value（返回值）Same type and kind as array and equal to the sum of all elements in array or the sum of elements along dimension dim. If mask is specified, only elements that correspondto .TRUE. elements in mask are summed. Returns a scalar if dim is omitted or array is one-dimensional. Otherwise, returns an array one dimension smaller than array.ExamplesINTEGER array (2, 3), i, j(3)array = RESHAPE((/1, 2, 3, 4, 5, 6/), (/2, 3/))! array is 1 3 5! 2 4 6i = SUM((/ 1, 2, 3 /)) ! returns 6j = SUM(array, DIM = 1) ! returns [3 7 11]WRITE(*,*) i, jEND5、SEEDRun-Time Subroutine Changes the starting point of the pseudorandom number generator. 改变随机数发生器的起始点ModuleUSE MSFLIBSyntax（语法）CALL SEED (iseed)iseed(Input) INTEGER(4). Starting point for RANDOM.Remarks（注解）SEED uses iseed to establish the starting point of the pseudorandom number generator.A given seed always produces the same sequence of values from RANDOM.If SEED is not called before the first call to RANDOM, RANDOM always begins with a seed value of one. If a program must have a different pseudorandom sequence each time it runs, pass the constant RND$TIMESEED (defined in MSFLIB.F90) to the SEED routine before the first call to RANDOM.ExampleUSE MSFLIBREAL randCALL SEED(7531)CALL RANDOM(rand)6、RANDOMPurposeRun-Time Subroutine Returns a pseudorandom number greater than or equal to zero and less than one from the uniform distribution. 返回大于或等于0且小于1，服从均匀分布的随机数ModuleUSE MSFLIBSyntaxCALL RANDOM (ranval)ranval(Output) REAL(4). Pseudorandom number, 0 ≤ranval< 1, from the uniformdistribution.RemarksA given seed always produces the same sequence of values from RANDOM.If SEED is not called before the first call to RANDOM, RANDOM begins with a seed value of one. If a program must have a different pseudorandom sequence each time it runs, pass the constant RND$TIMESEED (defined in MSFLIB.F90) to SEED before the first call to RANDOM.All the random procedures (RANDOM, RAN, and RANDOM_NUMBER, and the PortLib functions DRAND, DRANDM, RAND, IRANDM, RAND, and RANDOM) use the same algorithms and thus return the same answers. They are all compatible and can be used interchangeably. (The algorithm used is a “Prime Modulus M Multiplicative Linear Congruential Generator,” a modified version of t he random number generator by Park and Miller in “Random Number Generators: Good Ones Are Hard to Find,” CACM, October 1988, Vol. 31, No. 10.)CompatibilityCONSOLE STANDARD GRAPHICS QUICKWIN GRAPHICS WINDOWS DLL LIBExampleUSE MSFLIBREAL(4) ranCALL SEED(1995)CALL RANDOM(ran)7、FFT2BCompute the inverse Fourier transform of a complex periodic two-dimensional array.计算二维复数数组的逆傅里叶变换Usage（用法）CALL FFT2B (NRCOEF, NCCOEF, COEF, LDCOEF, A, LDA)Arguments（参数）NRCOEF— The number of rows of COEF. (Input) 数组COEF的行数NCCOEF— The number of columns of COEF. (Input) 数组COEF的列数COEF—NRCOEF by NCCOEF complex array containing the Fourier coefficients to be transformed. (Input) NRCOEF行NCCOEF列数组LDCOEF— Leading dimension of COEF exactly as specified in the dimension statement of the calling program. (Input)A—NRCOEF by NCCOEF complex array containing the Inverse Fourier coefficients of COEF. (Output) NRCOEF行NCCOEF列复数数组，包含数组COEF的逆傅里叶系数LDA— Leading dimension of A exactly as specified in the dimension statement of the calling program. (Input)Comments（注解）1.Automatic workspace usage isFFT2B4 * (NRCOEF + NCCOEF) + 32 + 2 *MAX(NRCOEF, NCCOEF) units, orDFFT2B8 * (NRCOEF + NCCOEF ) + 64 + 4 *MAX(NRCOEF, NCCOEF) units.Workspace may be explicitly provided, if desired, by use of F2T2B/DF2T2B. The reference isCALL F2T2B (NRCOEF, NCCOEF, A, LDA, COEF, LDCOEF,WFF1, WFF2, CWK, CPY)The additional arguments are as follows:WFF1— Real array of length 4 *NRCOEF + 15 initialized by FFTCI. The initialization depends on NRCOEF. (Input)WFF2— Real array of length 4 *NCCOEF + 15 initialized by FFTCI. The initialization depends on NCCOEF. (Input)CWK— Complex array of length 1. (Workspace)CPY— Real array of length 2 *MAX(NRCOEF, NCCOEF). (Workspace)2.The routine FFT2B is most efficient when NRCOEF and NCCOEF are the product of small primes.3.The arrays COEF and A may be the same.4.If FFT2D/FFT2B is used repeatedly, with the same values for NRCOEF and NCCOEF, then use FFTCI to fill WFF1(N = NRCOEF) and WFF2(N = NCCOEF). Follow this with repeated calls to F2T2D/F2T2B. This is more efficient than repeated calls toFFT2D/FFT2B.AlgorithmThe routine FFT2B computes the inverse discrete complex Fourier transform of a complex two-dimensional array of size (NRCOEF = N) ⨯ (NCCOEF = M). The method used is a variant of the Cooley-Tukey algorithm , which is most efficient when N and M are both products of small prime factors. If N and M satisfy this condition, then the computational effort is proportional to N M log N M. This considerable savings has historically led people to refer to this algorithm as the "fast Fourier transform" or FFT.Specifically, given an N⨯M array c = COEF, FFT2B returns in aFurthermore, a vector of Euclidean norm S is mapped into a vector of normFinally, note that an unnormalized inverse is implemented in FFT2D. The routine FFT2B is based on the complex FFT in FFTPACK. The package FFTPACK was developed by Paul Swarztrauber at the National Center for Atmospheric Research.ExampleIn this example, we first compute the Fourier transform of the 5 ⨯ 4 arrayfor 1 ≤n≤ 5 and 1 ≤m≤ 4 using the IMSL routine FFT2D. The resultis then inverted by a call to FFT2B. Note that the result is an array a satisfying a = (5)(4)x = 20x. In general, FFT2B is an unnormalized inverse with expansion factor N M.INTEGER LDA, LDCOEF, M, N, NCA, NRACOMPLEX CMPLX, X(5,4), A(5,4), COEF(5,4)CHARACTER TITLE1*26, TITLE2*26, TITLE3*26INTRINSIC CMPLXEXTERNAL FFT2B, FFT2D, WRCRNCTITLE1 = 'The input matrix is below 'TITLE2 = 'After FFT2D 'TITLE3 = 'After FFT2B 'NRA = 5NCA = 4LDA = 5LDCOEF = 5C Fill X with initial dataDO 20 N=1, NRADO 10 M=1, NCAX(N,M) = CMPLX(FLOAT(N+5*M-5),0.0)10 CONTINUE20 CONTINUECCALL WRCRN (TITLE1, NRA, NCA, X, LDA, 0)CCALL FFT2D (NRA, NCA, X, LDA, COEF, LDCOEF)CCALL WRCRN (TITLE2, NRA, NCA, COEF, LDCOEF, 0)CCALL FFT2B (NRA, NCA, COEF, LDCOEF, A, LDA)CCALL WRCRN (TITLE3, NRA, NCA, A, LDA, 0)CENDOutputThe input matrix is below1 2 3 41 ( 1.00, 0.00) ( 6.00, 0.00) ( 11.00, 0.00) ( 16.00, 0.00)2 ( 2.00, 0.00) ( 7.00, 0.00) ( 12.00, 0.00) ( 17.00, 0.00)3 ( 3.00, 0.00) ( 8.00, 0.00) ( 13.00, 0.00) ( 18.00, 0.00)4 ( 4.00, 0.00) ( 9.00, 0.00) ( 14.00, 0.00) ( 19.00, 0.00)5 ( 5.00, 0.00) ( 10.00, 0.00) ( 15.00, 0.00) ( 20.00, 0.00) After FFT2D1 2 3 41 ( 210.0, 0.0) ( -50.0, 50.0) ( -50.0, 0.0) ( -50.0, -50.0)2 ( -10.0, 13.8) ( 0.0, 0.0) ( 0.0, 0.0) ( 0.0, 0.0)3 ( -10.0, 3.2) ( 0.0, 0.0) ( 0.0, 0.0) ( 0.0, 0.0)4 ( -10.0, -3.2) ( 0.0, 0.0) ( 0.0, 0.0) ( 0.0, 0.0)5 ( -10.0, -13.8) ( 0.0, 0.0) ( 0.0, 0.0) ( 0.0, 0.0) After FFT2B1 2 3 41 ( 20.0, 0.0) ( 120.0, 0.0) ( 220.0, 0.0) ( 320.0, 0.0)2 ( 40.0, 0.0) ( 140.0, 0.0) ( 240.0, 0.0) ( 340.0, 0.0)3 ( 60.0, 0.0) ( 160.0, 0.0) ( 260.0, 0.0) ( 360.0, 0.0)4 ( 80.0, 0.0) ( 180.0, 0.0) ( 280.0, 0.0) ( 380.0, 0.0)5 ( 100.0, 0.0) ( 200.0, 0.0) ( 300.0, 0.0) ( 400.0, 0.0)8、TIMEFPurposePortLib Function Returns the number of seconds since the first time it is called, or zero.ModuleUSE PORTLIBSyntaxresult=TIMEF ( )Return ValueREAL(8). Number of seconds that have elapsed since the first time TIMEF( ) was called. The first time called, TIMEF returns 0.0D0.CompatibilityCONSOLE STANDARD GRAPHICS QUICKWIN GRAPHICS WINDOWS DLL LIBExampleUSE PORTLIBINTEGER i, jREAL(8) elapsed_timeelapsed_time = TIMEF() DO i = 1, 100000j = j + 1END DOelapsed_time = TIMEF() PRINT *, elapsed_time END。
fortran常用程序
fortran常用程序----------------------- Page 1----------------------- Fortran 简单程序集杨洋1!大小写转换program ex39implicit nonecharacter*20 strinteger iprint*,'input the string:'read*,strdo i=1,len_trim(str)if(str(i:i)>='a'.and.str(i:i)<='z') str(i:i)=char(ichar(str(i:i))-32) end doprint*,'the inverted string:'print*,strend!将字符串转化为整数program ex104implicit noneexternal fcharacter*10 strinteger fprint*,'输入由数字组成的字符串:'read*,strprint*,'转化后的整数:'print*,f(str)endfunction f(str)implicit nonecharacter*(*) strinteger f,k,if=0k=len_trim(str)do i=1,kf=f+(ichar(str(i:i))-ichar('0'))*10**(k-i)end doend function!将十进制数转化为二进制数（用字符串保存）program ex1101 作者简介：杨洋，南京信息工程大学大气科学学院海洋科学系2008 级1----------------------- Page 2-----------------------program ex110implicit noneinteger acharacter*8::b=' 'print*,'输入一个十进制整数:'read*,ado while(a>0)if (mod(a,2)==1) thenb='1'//belseb='0'//bend ifa=a/2end doprint*,'对应的二进制数为:'print*,bend!将二进制数（用字符串保存）转化为十进制数program ex111implicit nonecharacter*8 ainteger::b=0,k,iprint*,'输入一个二进制数:'read*,ak=len_trim(a)do i=1,kb=b+(ichar(a(i:i))-ichar('0'))*2**(k-i)end doprint*,'对应的十进制数为:'print*,bend!统计大写、小写、数字及其他字符的个数program ex51implicit nonecharacter*20 strinteger::i,n1=0,n2=0,n3=0,n4=0print*,'输入字符串:'read*,strdo i=1,len_trim(str)select case(str(i:i))case('a':'z')n1=n1+12----------------------- Page 3-----------------------case('A':'Z')n2=n2+1case('0':'9')n3=n3+1case defaultn4=n4+1end selectend doprint*,’大写字母个数:’,n1print*,’小写字母个数:’,n2print*,’数字个数:’,n3print*,’其他字符个数:’,n4End!判断闰年program ex24implicit noneinteger yearprint*,'input a year:'read*,yeark=k+1work(k)=num(i)end ifend docall sort(work,k)print*,'the final array:'8----------------------- Page 9-----------------------print*,(work(i),i=1,k)endfunction prime(n)implicit noneinteger i,k,nlogical primek=sqrt(real(n))do i=2,kif (mod(n,i)==0) exitend doif (i>k) thenprime=.true.elseprime=.false.end ifend functionsubroutine sort(a,n)implicit noneinteger i,j,k,n,a(n),tdo i=1,n-1k=ido j=i+1,nif (a(j)>a(k)) k=jend doif (k/=i) thent=a(k)a(k)=a(i)a(i)=tend ifend doend subroutine! 求2～999 中同时满足下列条件的数：(a) 该数各位数字之和为奇数；(b) 该数是素数。
fortran语法手册
1F O R T R A N77四则运算符+ - */ ** (其中**表示乘方)在表达式中按优先级次序由低到高为：+或-→*或/→**→函数→()2 FORTRAN77变量类型隐含约定：I-N规则凡是以字母I，J，K，L，M，N六个字母开头的，即认为是整型变量，其它为实型变量。
如IMPLICIT REAL (I,J)三种定义的优先级别由低到高顺序为：I-N规则→IMPLICIT语句→类型说明语句，因此，在程序中IMPLICIT语句应放在类型说明语句之前。
数组的说明与使用使用I-N规则时用DIMENSION说明数组，也可在定义变量类型同时说明数组，说明格式为：数组名(下标下界,下标上界），也可省略下标下界，此时默认为1，例：DIMENSION IA(0:9),ND(80:99),W(3,2),NUM(-1:0),A(0:2,0:1,0:3)REAL IA(10),ND(80:99)使用隐含DO循环进行数组输入输出操作：例如WRITE(*,10) ('I=',I,'A=',A(I),I=1,10,2)10FORMAT(1X,5(A2,I2,1X,A2,I4))使用DATA语句给数组赋初值变量表中可出现变量名,数组名,数组元素名,隐含DO循环，但不许出现任何形式的表达式：例如DATA A,B,C/,,DATA A/,B/,C/DATA A,B,C/3*CHARACTER*6 CHN(10)DATA CHN/10*''/INTEGER NUM(1000)DATA (NUM(I),I=1,500)/500*0/,(NUM(I),I=501,1000)/500*1/3 FORTRAN77程序书写规则程序中的变量名，不分大小写；变量名称是以字母开头再加上1到5位字母或数字构成，即变更名字串中只有前6位有效；一行只能写一个语句；程序的第一个语句固定为PROGRAM 程序名称字符串某行的第１个字符至第５个字符位为标号区，只能书写语句标号或空着或注释内容；某行的第１个字符为Ｃ或＊号时，则表示该行为注释行，其后面的内容为注释内容；某行的第６个字符位为非空格和非０字符时，则该行为上一行的续行，一个语句最多可有19个续行；某行的第７至７２字符位为语句区，语句区内可以任加空格以求美观；某行的第７３至８０字符位为注释区，８０字符位以后不能有内容。
fortran语法手册
1 FORTRAN77四则运算符+ - * / ** (其中**表示乘方)在表达式中按优先级次序由低到高为： +或-→*或/→**→函数→()2 FORTRAN77变量类型2.1 隐含约定：I-N规则凡是以字母I，J，K，L，M，N六个字母开头的，即认为是整型变量，其它为实型变量。
2.2 用类型说明语句确定变量类型：可以改变I-N规则2.3 用IMPLICIT语句将某一字母开头的全部变量指定为所需类型如 IMPLICIT REAL (I,J)三种定义的优先级别由低到高顺序为：I-N规则→IMPLICIT语句→类型说明语句，因此，在程序中IMPLICIT语句应放在类型说明语句之前。
2.4 数组的说明与使用使用I-N规则时用DIMENSION说明数组，也可在定义变量类型同时说明数组，说明格式为：数组名(下标下界,下标上界），也可省略下标下界，此时默认为1，例：DIMENSION IA(0:9),ND(80:99),W(3,2),NUM(-1:0),A(0:2,0:1,0:3)REAL IA(10),ND(80:99)使用隐含DO循环进行数组输入输出操作：例如WRITE(*,10) ('I=',I,'A=',A(I),I=1,10,2)10FORMAT(1X,5(A2,I2,1X,A2,I4))2.5 使用DATA语句给数组赋初值变量表中可出现变量名,数组名,数组元素名,隐含DO循环，但不许出现任何形式的表达式：例如DATA A,B,C/-1.0,-1.0,-1.0/DATA A/-1.0/,B/-1.0/,C/-1.0/DATA A,B,C/3*-1.0/CHARACTER*6 CHN(10)DATA CHN/10*' '/INTEGER NUM(1000)DATA (NUM(I),I=1,500)/500*0/,(NUM(I),I=501,1000)/500*1/3 FORTRAN77程序书写规则程序中的变量名，不分大小写；变量名称是以字母开头再加上1到5位字母或数字构成，即变更名字串中只有前6位有效；一行只能写一个语句；程序的第一个语句固定为PROGRAM 程序名称字符串某行的第１个字符至第５个字符位为标号区，只能书写语句标号或空着或注释内容；某行的第１个字符为Ｃ或＊号时，则表示该行为注释行，其后面的内容为注释内容；某行的第６个字符位为非空格和非０字符时，则该行为上一行的续行，一个语句最多可有19个续行；某行的第７至７２字符位为语句区，语句区内可以任加空格以求美观；某行的第７３至８０字符位为注释区，８０字符位以后不能有内容。
Fortran语法（中级）
Fortran语法（中级）!-----------------------------------------------------------------------EXTERNAL属性和哑过程 (哑元为外部过程，即哑过程)指定EXTERNAL语句或属性说明实元实际上是外部过程类型定义语句：类型,EXTERNAL :: 外部函数名[,外部函数名]…或EXTERNAL语句：EXTERNAL [外部函数名][,⼦程序名][,块数据名]…哑元也可以是⼀个过程，这时作为哑元的过程称为哑过程。
（⾄少两层调⽤）例如：!-----------------------------------------------------------------------Programm mainReal x,yExternal Plus !外部过程名作实元，必须⽤External说明，或者具有External属性x=1.0 ; y=2.0Print,* Calculate(x,y,Plus) !调⽤Calculate函数，实元为外部过程PlusEnd Program mainReal Function Plus(a,b) !(第⼆层被调⽤的外部函数)Real, Intent(In) :: a,bPlus=a+bEnd Function PlusReal Function Calculate (x,y,func)Real, Intent(In) :: x,yReal, External func !类型定义语句, 说明哑元时⼀个外部过程, 也可以直接⽤External说明Calculate=func(x,y) !调⽤⾃定义的外部函数End Function Calculate!-----------------------------------------------------------------------或者将 Real, External func 改为接⼝程序：InterfaceReal Function Plus(a,b) !Plus被接⼝块说明为⼀个哑元，即⼀个哑过程Real, Intent(In) :: a,bEnd Function PlusEnd InterfaceINTENT属性 (过程的哑元说明)在类型定义语句中：类型,INTENT(意图说明符) :: 哑元名表或⽤INTENT语句： INTENT(意图说明符) :: 哑元名表意图说明符为以下字符串：IN 指明哑元仅⽤于向过程提供数据，过程的执⾏期间哑元不能被重定义或成为未定义的，相联合的实元可以是常数、变量、数组以及它们的算术表达式。
fortran语言并行计算程序
fortran语言并行计算程序Fortran是一种编程语言，用于科学和工程计算。
Fortran语言中可以使用并行计算来加速程序的执行。
以下是一个简单的Fortran并行计算程序的示例：```fortranprogram parallel_program! 使用OpenMP库启用并行计算use omp_libimplicit noneinteger :: iinteger, parameter :: N = 1000000real :: x(N), y(N), z(N)! 初始化数组x = 1.0y = 2.0! 使用并行循环计算!$omp parallel dodo i = 1, Nz(i) = x(i) + y(i)end do! 打印部分结果do i = 1, 10print *, z(i)end doend program parallel_program```在上面的示例中，使用OpenMP库来启用并行计算。
通过添加`use omp_lib`语句，可以在程序中使用OpenMP的并行化指令。
在循环计算的部分，使用了`!$omp parallel do`指令，它告诉编译器这个循环可以并行执行。
这样，循环中的每个迭代都可以在不同的处理器上并行执行。
在这个示例中，数组`x`和`y`被初始化为 1.0和2.0。
然后，使用并行循环计算数组`z`的每个元素，将`x`和`y`对应位置的元素相加。
通过循环打印了数组`z`的前10个元素。
请注意，实际的并行计算程序可能需要更复杂的逻辑和更多的并行化指令，具体取决于问题的性质和计算的需求。
上述示例只是一个简单的示例，用于说明Fortran中的并行计算的基本概念。
fortranwrite用法
fortranwrite用法Fortran中的fortranwrite是一个用于文件写入的子程序。
它提供了一种将数据以二进制格式写入文件的方法，并可以按照指定格式将数据转换为字符文本。
fortranwrite具有灵活的选项，可以根据不同的需求调整数据写入的方式。
fortranwrite的基本语法如下：`WRITE (unit, format) [input-list]`其中，unit指定了要写入的文件单元号，通常是一个整数或逻辑单元号。
这个单元号可以是事先由OPEN语句获得的逻辑单元号，也可以是内建的实体所定义的逻辑单元号（如*或者"stdout"表示控制台输出，或者trim文件的单元号）。
format是一个字符表达式，用于指定将数据转换为字符的格式。
格式指令是用于指定如何写入和读取数据的符号，也称为“格式编辑描述符”。
Fortran的格式编辑描述符包括格式宽度、字段宽度、编辑序列等。
input-list是需要写入文件的数据列表。
数据可以是常数、数组、变量或者复合结构。
input-list中的每个元素都按照format中的格式转换为字符并写入文件。
下面是使用fortranwrite的一个例子：fortranprogram write_exampleimplicit noneinteger :: i, jreal :: x(10), y(10)! 给数组x和y赋值do i = 1, 10x(i) = iy(i) = i2end do! 打开文件并写入数据open(unit=10, file='data.dat', status='replace', action='write') do i = 1, 10write(10, '(2F10.2)') x(i), y(i)end doclose(unit=10)! 从文件读取数据并打印open(unit=10, file='data.dat', status='old', action='read')do i = 1, 10read(10, '(2F10.2)') x(i), y(i)end doclose(unit=10)! 打印读取到的数据do i = 1, 10write(*, '(2F10.2)') x(i), y(i)end doend program write_example在上面的示例中，我们首先定义了两个数组x和y，然后使用do循环为它们赋值。
fortranwrite用法 -回复
fortranwrite用法-回复Fortran是一种编程语言，广泛应用于科学计算和数值分析领域。
它由IBM 公司于20世纪50年代初开发，旨在为科学家和工程师提供一种高级编程语言，以简化复杂的数学计算和数据处理。
Fortran具有强大的数值计算能力和高效的执行速度，因此在科学计算领域一直被广泛使用。
在Fortran中，fortranwrite函数是用于向文件中写入数据的一个子程序。
该函数可以将变量或数组的值写入一个已打开的文件。
fortranwrite函数的语法如下：CALL FORTRANWRITE(N, DATA)其中N是要写入文件的数据的数量，DATA是包含数据的变量或数组。
该函数将数据写入默认的输出文件，通常是标准输出设备（即终端窗口）。
如果要将数据写入特定的文件，可以使用另一个fortranopen函数来打开该文件，并将其作为fortranwrite函数的第一个参数。
下面我将一步一步回答关于fortranwrite函数的使用方法和示例。
第一步：声明变量和数组在使用fortranwrite函数之前，首先需要声明变量或数组。
例如，我们可以声明一个整数变量N和一个长度为10的实数数组X，如下所示：INTEGER :: NREAL :: X(10)第二步：打开文件接下来，我们需要打开一个文件以便写入数据。
可以使用fortranopen函数打开一个文件，并返回一个文件句柄（file handle）。
例如，我们可以打开一个名为"output.txt"的文件，如下所示：INTEGER :: FHCALL FORTRANOPEN(FH, "output.txt", "WRITE")在这个例子中，文件句柄存储在FH变量中，文件名是"output.txt"，打开模式为"WRITE"，表示我们只能写入文件但不能读取。
fortran语言语法
fortran语⾔语法FORTRAN是世界上最早出现的⾼级编程语⾔，是⼯程界最常⽤的编程语⾔，它在科学计算中（如航空航天、地质勘探、天⽓预报和建筑⼯程等领域）发挥着极其重要的作⽤。
经过40多年的发展，伴随着FORTRAN语⾔多次版本的更新及相应开发系统的出现，其功能不断完善，最新版本的开发系统⼏乎具备了VC、VB的所有特点，如图形界⾯编程、数据库等。
⽬前，⼯科院校开设的计算机编程语⾔课⾸选仍然是FORTRAN :<说实话，从科技发展的趋势来说这不是好事。
您可以设想⼀下，如果需要⽤鹅⽑笔抄写⼤量的古籍是什么感受！强烈建议阅读《发掘C#特性赋予科学计算项⽬以威⼒》1 FORTRAN77四则运算符+ - * / ** (其中**表⽰乘⽅)在表达式中按优先级次序由低到⾼为： +或-→*或/→**→函数→()2 FORTRAN77变量类型隐含约定：I-N规则凡是以字母I，J，K，L，M，N六个字母开头的，即认为是整型变量，其它为实型变量。
⽤类型说明语句确定变量类型：可以改变I-N规则⽤IMPLICIT语句将某⼀字母开头的全部变量指定为所需类型如IMPLICIT REAL (I,J)三种定义的优先级别由低到⾼顺序为：I-N规则→IMPLICIT语句→类型说明语句，因此，在程序中IMPLICIT语句应放在类型说明语句之前。
数组的说明与使⽤使⽤I-N规则时⽤DIMENSION说明数组，也可在定义变量类型同时说明数组，说明格式为：数组名(下标下界,下标上界），也可省略下标下界，此时默认为1，例：DIMENSION IA(0:9),ND(80:99),W(3,2),NUM(-1:0),A(0:2,0:1,0:3)REAL IA(10),ND(80:99)使⽤隐含DO循环进⾏数组输⼊输出操作：例如WRITE(*,10) ('I=',I,'A=',A(I),I=1,10,2)10FORMAT(1X,5(A2,I2,1X,A2,I4))使⽤DATA语句给数组赋初值变量表中可出现变量名,数组名,数组元素名,隐含DO循环，但不许出现任何形式的表达式：例如DATA A,B,C/,,DATA A/,B/,C/DATA A,B,C/3*CHARACTER*6 CHN(10)DATA CHN/10*''/INTEGER NUM(1000)DATA (NUM(I),I=1,500)/500*0/,(NUM(I),I=501,1000)/500*1/3 FORTRAN77程序书写规则·程序中的变量名，不分⼤⼩写；·变量名称是以字母开头再加上1到5位字母或数字构成，即变更名字串中只有前6位有效；·⼀⾏只能写⼀个语句；·程序的第⼀个语句固定为PROGRAM 程序名称字符串·某⾏的第１个字符⾄第５个字符位为标号区，只能书写语句标号或空着或注释内容；·某⾏的第１个字符为Ｃ或＊号时，则表⽰该⾏为注释⾏，其后⾯的内容为注释内容；·某⾏的第６个字符位为⾮空格和⾮０字符时，则该⾏为上⼀⾏的续⾏，⼀个语句最多可有19个续⾏；·某⾏的第７⾄７２字符位为语句区，语句区内可以任加空格以求美观；·某⾏的第７３⾄８０字符位为注释区，８０字符位以后不能有内容。
(完整版)fortran文件操作
fortran文件操作From: 《Fortran 95 程序设计》彭国伦1.文件读取的概念：读取：“顺序读取”和“直接读取”。
保存：“文本文件”和“二进制文件”。
二进制文件：读取快，因为保存格式和数据在内存中的保存方法一样，同时也节省空间。
－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－2. 文件的操作：－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－open的使用：使用open命令打开文件之后，就可以对文件来做输入输出。
example:program ex0901impicit noneopen(unit=10, file='hello.txt') ! 打开hello.txt文件, unit指定文件代码，file 指定文件名称。
write(10, *) "hello" !在代码为10的文件中写入hellostopend program ex0901open中有很多参数可以使用，详细如下：OPEN(UNIT=number, FILE='filename', FORM='...', STATUS='...', ACCESS='...', RECL=length, ERR=label, IOSTAT=iostat, BLANK='...', POSITION='...', ACTION=action, PAD='...', DELIM='...')UNIT='number'： number必须是一个正整数，它可以使用变量或是常量来赋值。
number最好避开1,2,5,6。
因为2,6是默认的输出位置，也就是屏幕。
1，5则是默认的输入位置，键盘。
FILE='filename'：指定要打开的文件名称，文件名要符合系统规定。
fortran基础知识1----循环
fortran基础知识1----循环do循环结构使得⼀个语句或⼀系列语句，以进⾏迭代，当⼀个给定的条件为真。
语法do循环的⼀般形式是：do var = start, stop [,step]! statement(s)…end do这⾥,循环变量var应该是⼀个整数start 是初始值stop 是最终值步骤step是递增，如果此被省略，则变量var以单位增加跳出循环Fortran中跳出循环的两个命令为exit和cycle，exit表⽰跳出离命令最近的循环，cycle表⽰结束此次循环并开始下⼀次循环。
pause：程序执⾏到pause 时，会暂停执⾏，直到⽤户按下Enter键才会继续执⾏。
continue这个命令没有实际的⽤途，它的功能就是继续向下执⾏程序，相当于python中的pass。
stop⽤来结束程序执⾏。
cycle命令可由略过循环的程序模块中，相当于python和C中的continue，直接跳过后⾯的所有代码，然后开始下⼀次循环。
exit的功能是可以直接跳出⼀个正在进⾏的循环，不论是do 循环还是do while 循环，相当于python和C中的break。
1 program sime23implicit none4 integer :: tn56 tn=27 print_tn : do8 print *, tn9 tn = tn+210if (tn>=200) exit print_tn11 end do print_tn12 pause1314 end1 program odd_num23implicit none4 integer :: num,num_odd5 print *, 'Input data number:'6 read *,num78 num_odd = 0910do11 num = num - 11213if (num<0) then14 exit15 elseif (mod(num,2)==0) then16 cycle17else18 num_odd = num_odd + 119 endif202122 print *, num23 enddo2425 print *, 'The number of odd number is',num_odd 2627 pause2829 end如果输⼊20，输出结果为：。
（完整）Fortran经典编程语言笔记（你值得拥有）
（完整）Fortran经典编程语言笔记（你值得拥有）FORTRAN笔记2014.10.07目录第七讲_FORTRAN的基本知识.ppt (2)FORTRAN语言程序设计初步 (2)FORTRAN源程序的书写格式（以77为例） (2)变量 (2)变量类型 (2)算术运算符和运算优先级 (3)赋值语句 (3)参数语句（PARAMETER语句） (3)END语句 (3)PAUSE语句 (3)逻辑运算和选择结构 (4)关系表达式 (4)FORTRAN中数组的定义及使用 (4)其他 (5)1. fortran语言定义CHARACTER*6 TTL(14,3),CNAM(400)是什么意思? (5)2. fortran里character*10 是什么意思 (5)3. Fortran中kind是什么函数？ (5)第七讲_FORTRAN的基本知识.pptFORTRAN语言程序设计初步FORTRAN是Formula Translation的缩写，意为“公式翻译”，它是为科学、工程问题或企事业管理中的那些能够用数学公式表达的问题而设计的，其数值计算的功能较强。
常用的是FORTRAN77和FORTRAN90两种标准。
1、一个程序由若干个程序单位组成。
主程序和每一个子程序分别是一个独立的程序单位。
2、每一个程序单位都是以“END”结束的。
3、一个程序单位包括若干行。
1）语句行。
由一个FORTRAN语句组成。
2）非语句行，即注释行。
4、FORTRAN程序中的语句可以没有标号，也可以有标号，根据需要而定。
标号的作用是标志一个语句以便被其他语句引用。
5、一个程序单位中各类语句的位置是有一定规定的。
6、FORTRAN源程序必须按一定的格式书写。
FORTRAN源程序的书写格式（以77为例）每一行有80列，分别如下：1、第1-5列为标号区。
一行中第一列为“C”或“*”，该行即被认为是注释行。
2、第6列为“续行标志区”，如果在一行的第6列上写一个非空格和非零的字符，则该行作为其上一行的续行。
Fortran函数大全(1)
Fortran函数大全(1)Fortran库参考SunStudio11SunMicrosystems,文件号码819-4757-102005年11月, 修订版A请将有关本文档的意见和建议提交至：/hwdocs/feedback版权所有#169;2005SunMicrosystems,Inc.,4150NetworkCircle,SantaCla ra,California95054,U.S.A.保留所有权利.美国政府权利－商业用途. 政府用户应遵循SunMicrosystems,Inc.的标准许可协议, 以及FAR （FederalAcquisitionRegulations, 即“联邦政府采购法规”）的适合条款及其补充条款. 必须依据许可证条款使用. 本发行版可能包含由第三方开辟的内容.本产品的某些部份可能是从BerkeleyBSD系统衍生出来的, 并获得了加利福尼亚大学的许可. UNIX是X/OpenCompany,Ltd.在美国和其他国家/地区独家许可的注册商标.Sun、SunMicrosystems、Sun徽标、Java和JavaHelp是SunMicrosystems,Inc.在美国和其他国家/地区的商标或者注册商标. 所有的SPARC商标的使用均已获得许可, 它们是SPARCInternational,Inc.在美国和其他国家/地区的商标或者注册商标. 标有SPARC商标的产品均基于由SunMicrosystems,Inc.开辟的体系结构.本服务手册所介绍的产品以及所包含的信息受美国出口控制法制约, 并应遵守其他国家/地区的进出口法律. 严禁将本产品直接或者间接地用于核设施、导弹、生化武器或者海上核设施, 也不能直接或者间接地出口给核设施、导弹、生化武器或者海上核设施的最终用户. 严禁出口或者转口到美国禁运的国家/地区以及美国禁止出口清单中所包含的实体, 包括但不限于被禁止的个人以及特殊指定的国家/地区的公民. 本文档按“原样”提供, 对于所有明示或者默示的条件、陈述和担保, 包括对适销性、适合性或者非侵权性的默示保证, 均不承担任何责任, 除非此免责声明的适合范围在法律上无效.目录阅读本书之前印刷约定xixiShell提示符支持的平台xiixiiixiii访问SunStudio软件和手册页访问编译器和工具文档访问相关的Solaris文档开辟者资源xviiixixxixxvixviii联系Sun技术支持Sun欢迎您提出意见1. Fortran库例程1.11.21.31–1数据类型注意事项64位环境1–21–1Fortran数学函数1.3.11.3.21.3.3单精度函数1–31–3双精度函数1–6四倍精度函数1–91–111–11iii1.4Fortran库例程参考1.4.1abort:终止和写入核心文件1.4.21.4.31.4.41.4.51.4.61.4.71.4.81.4.91.4.101.4.111.4. 121.4.131.4.141.4.151.4.161.4.171.4.181.4.191.4.201.4.21 1.4.221.4.231.4.241.4.251.4.261.4.271.4.281.4.29access:检查文件权限或者存在性1–11alarm:在指定的时间后调用子例程1–12bit:位函数：and、or、…、bit、setbit、…chdir：更改默认目录chmod:更改文件的模式1–131–161–171–18date：获取以字符串表示的当前日期dtime, etime：经过的执行时间exit：终止进程并设置状态1–201–231–23fdate：以ASCII字符串返回日期和时间flush：刷新逻辑单元的输出1–24fork：创建当前进程的副本1–25fseek,ftell：确定文件的位置以及重新确定文件的位置1–25fseeko64, ftello64：确定大文件的位置以及重新确定大文件的位置1–271–29getarg, iargc：获取命令行参数getc, fgetc：获取下一个字符1–301–32getcwd：获取当前工作目录的路径getenv：获取环境变量的值1–33getfd：获取外部单元编号的文件描述符getfilep：获取外部单元编号的文件指针getlog:获取用户的登录名getpid：获取进程ID1–341–341–361–361–36getuid, getgid：获取进程的用户ID或者组IDhostnm：获取当前主机的名称1–37idate：返回当前日期1–381–38ieee_flags, ieee_handler,sigfpe：IEEE算术index, rindex, lnblnk：子串的索引或者长度. 1–44inmax：返回最大正整数1–45itime：当前时间1–46ivFortran库参考2005年11月1.4.301.4.311.4.321.4.331.4.341.4.35kill：将信号发给进程1–471–47link, symlnk：链接到现有的文件loc：返回对象的地址1–491–491–50long, short：整型对象转换longjmp,isetjmp：返回至isetjmp设置的位置malloc, malloc64, realloc, free：分配/重新分配/释放内存1–521.4.36mvbits：挪移位字段1–561.4.37perror, gerror, ierrno：获取系统错误消息1–571.4.38putc, fputc：将字符写入逻辑单元1–581.4.39qsort, qsort64：对一维数组的元素进行排序1–601.4.40ran：生成一个介于0和1之间的随机数1–621.4.41rand, drand, irand：返回随机值1–631.4.42rename：重命名文件1–641.4.43secnds：获取以秒数表示的系统时间并减去参数1–651.4.44set_io_err_handler, get_io_err_handler：设置并获取I/O错误处理程序1–651.4.45sh：快速执行sh命令1–681.4.46signal：更改信号的操作1–691.4.47sleep：一段时间暂停执行1–701.4.48stat, lstat, fstat：获取文件状态1–701.4.49stat64, lstat64, fstat64：获取文件状态1–731.4.50system：执行系统命令1–731.4.51time, ctime, ltime, gmtime：获取系统时间1–741.4.52ttynam, isatty：获取终端端口的名称1–781.4.53unlink：删除文件1–791.4.54wait：等待进程终止1–802.Fortran95内部函数2–12.1标准Fortran95的通用内部函数2–1v目录2.1.12.1.22.1.32.1.42.1.52.1.62.1.72.1.82.1.92.1.102.1.1 12.1.122.1.132.1.142.1.152.1.162.1.172.1.182.1.192.1.202 .1.212.1.222.1.232.2参数存在查询函数数值函数数学函数字符函数。
mpif90编译参数
mpif90编译参数mpif90编译参数是用于编译Fortran程序的一组指令和选项。
在本文中，我们将探讨一些常用的mpif90编译参数，以及它们的作用和用法。
一、-o 参数-o 参数用于指定生成的可执行文件的名称。
例如，使用以下命令编译Fortran程序：mpif90 -o myprogram.exe myprogram.f90这将生成一个名为myprogram.exe的可执行文件。
二、-c 参数-c 参数用于将源代码文件编译为目标文件，而不进行链接操作。
例如，使用以下命令编译Fortran程序：mpif90 -c myprogram.f90这将生成一个名为myprogram.o的目标文件。
三、-I 参数-I 参数用于指定编译器在搜索头文件时要查找的目录。
例如，使用以下命令编译Fortran程序：mpif90 -I /path/to/include myprogram.f90这将告诉编译器在路径/path/to/include中搜索头文件。
四、-L 参数和-l 参数-L 参数用于指定链接器在搜索库文件时要查找的目录，而-l 参数用于指定要链接的库文件。
例如，使用以下命令编译Fortran程序：mpif90 -L /path/to/lib -l mylib myprogram.f90这将告诉链接器在路径/path/to/lib中搜索名为libmylib.so的库文件，并将其链接到可执行文件中。
五、-f 参数-f 参数用于指定编译器的优化级别和特定的编译选项。
例如，使用以下命令编译Fortran程序：mpif90 -fopenmp -O3 myprogram.f90这将启用OpenMP多线程支持，并使用最高级别的优化。
六、-D 参数-D 参数用于定义预处理器宏。
例如，使用以下命令编译Fortran程序：mpif90 -DDEBUG myprogram.f90这将在编译时定义一个名为DEBUG 的宏，可以在程序中使用条件编译。
第3章 最基本的FORTRAN语句
前一页 休息
2
算术赋值语句
算术赋值语句的一般格式：
xxxxx v=e 其中：v是变量名或数组元素。 e是算术表达式。
作用：将算术表达式e的结果赋给变量v
前一页 休息
3
参数说明语句
参数说明语句的一般格式：
parameter(参数名=常数,……) 作用：将常数赋给参数名。 例如：
parameter(pi=3.1415926,k=100,a0=214.0,imax=5) 程序中定义的参数名只能够在表达式中使用，
FORTRAN程序设计
第3章 最基本的FORTRAN语句
FORTRAN语句
FORTRAN语句分为可执行语句和不可语句。
可执行语句：计算机在执行的过程中是要执行 的。标准格式书写时，根据具体情况语句标号可以 写也可以不写。
不可执行语句：计算机在执行的过程中是不执 行的。标准格式书写时，语句标号一定不能写或语 句标号一定要写（仅指FORMAT语句）。
据宽度,d表示指数型小数部分的位数
前一页 休息
6
输入和输出
非数值型格式说明符：
X型格式符： nX n为十进制正整数，空n格
H型格式符：
nHh1h2…hn n为十进制正整数，原样输出h1h2 …hn 或‘h1h2…hn’ 原样输出h1h2 …hn
格式符的重复数：
n格式符 n为十进制正整数，表示格式符的重复数
不能对它赋值。
前一页 休息
4
输入和输出
一：带格式输出语句的一般形式
XXXXX write(设备号，格式语句标号）变量名表 语句标号 format(1x,格式说明符表） 变量名表:由若干个变量名或数组元素相互间用逗号分
割组成的内容。 格式说明符表：由若干个格式说明符相互间用逗号分
fortran语言中entry用法(一)
fortran语言中entry用法(一)Fortran语言中entry用法Fortran语言中的entry语句用于在程序中定义子程序的入口点。
以下是一些entry语句的用法和详细讲解：创建简单的入口点在Fortran中，您可以使用entry语句创建一个简单的入口点。
以下是一个示例：program mainimplicit noneinteger :: num, resultnum = 5result = 0call square(num, result)print *, "Square of", num, "is", resultstopcontainsentry square(x, y)integer, intent(in) :: xinteger, intent(out) :: yy = x * xreturnend entry squareend program main在上面的示例中，我们在contains块中定义了一个名为square 的子程序入口点。
square接受一个整数参数x，计算它的平方，并将结果存储在传入的整数参数y中。
注意，在entry语句中，我们使用了intent属性来指定参数的用途。
多个入口点Fortran中的entry语句还允许我们定义多个入口点，即多个子程序。
以下是一个示例：program mainimplicit noneinteger :: num1, num2, resultnum1 = 5num2 = 3result = 0call add(num1, num2, result)print *, "Result of adding", num1, "and", num2, "is", resultcall subtract(num1, num2, result)print *, "Result of subtracting", num2, "from", num1, "is", resultstopcontainsentry add(x, y, z)integer, intent(in) :: x, yinteger, intent(out) :: zz = x + yreturnend entry addentry subtract(x, y, z)integer, intent(in) :: x, yinteger, intent(out) :: zz = x - yreturnend entry subtractend program main在上面的示例中，我们定义了两个子程序入口点：add和subtract。
fortran内部命令
Categories of Intrinsic FunctionsCategory Subcategory DescriptionNumeric Computation Elemental functions that performtype conversions or simple numericoperations: ABS, AIMAG, AINT, AMAX0,AMIN0, ANINT, CEILING, CMPLX, CONJG,DBLE, DCMPLX, DFLOAT, DIM, DPROD,DREAL, FLOAT, FLOOR, IFIX, IMAG,INT, MAX, MAX1, MIN, MIN1, MOD,MODULO, NINT, QCMPLX, QEXT, QFLOAT,QREAL, RAN, REAL, SIGN, SNGL, ZEXTManipulation1Elemental functions that returnvalues related to the components ofthe model values associated with theactual value of the argument:EXPONENT, FRACTION, NEAREST,RRSPACING, SCALE, SET_EXPONENT,SPACINGInquiry1Functions that return scalar valuesfrom the models associated with thetype and kind parameters of theirarguments2: DIGITS, EPSILON, HUGE,ILEN, MAXEXPONENT, MINEXPONENT,PRECISION, RADIX, RANGE, SIZEOF,TINYTransformational Functions that perform vector andmatrix multiplication: DOT_PRODUCT,MATMULSystem Functions that return informationabout a process or processor:PROCESSORS_SHAPE, NWORKERS,NUMBER_OF_PROCESSORS, SECNDSKind type Functions that return kind typeparameters: KIND,SELECTED_INT_KIND,SELECTED_REAL_KINDMathematical Functions that perform mathematicaloperations: ACOS, ACOSD, ASIN,ASIND, ATAN, ATAND, ATAN2, ATAN2D,COS, COSD,COSH, COTAN, COTAND, EXP,LOG, LOG10, SIN, SIND, SINH, SQRT,TAN, TAND, TANHBit Manipulation Elemental functions that performsingle-bit processing, and logicaland shift operations; and allow bitsubfields to be referenced: AND,BTEST, IAND, IBCHNG, IBCLR, IBITS,IBSET, IEOR, IOR, ISHA, ISHC, ISHFT,ISHFTC, ISHL, LSHIFT, NOT, OR,RSHIFT, XORInquiry Function that lets you determineparameter s(the bit size) in the bitmodel4: BIT_SIZERepresentation Functions that return information onbit representation of integers:LEADZ, POPCNT, POPPAR, TRAILZCharacter Comparison Elemental functions that make alexical comparison of thecharacter-string arguments andreturn a default logical result:LGE, LGT, LLE, LLTConversion Elemental functions that takecharacter arguments and returninteger, ASCII, or charactervalues3: ACHAR, CHAR, IACHAR, ICHARString handling Functions that perform operations oncharacter strings, return lengths ofarguments, and search for certainarguments: ADJUSTL, ADJUSTR, INDEX,LEN_TRIM, REPEAT, SCAN, TRIM,VERIFYInquiry Function that returns length ofargument: LENArray Construction Functions that construct new arraysfrom the elements of existing array:MERGE, PACK, SPREAD, UNPACKInquiry Functions that let you determine ifan array argument is allocated, andreturn the size or shape of an array,and the lower and upper bounds ofsubscripts along each dimension:ALLOCATED, LBOUND, SHAPE, SIZE,UBOUNDLocation Transformational functions thatfind the geometric locations of themaximum and minimum values of anarray: MAXLOC, MINLOCManipulation Transformational functions thatshift an array, transpose an array,or change the shape of an array:CSHIFT, EOSHIFT, RESHAPE,TRANSPOSEReduction Transformational functions thatperform operations on arrays. Thefunctions "reduce" elements of awhole array to produce a scalarresult, or they can be applied to aspecific dimension of an array toproduce a result array with a rankreduced by one: ALL, ANY, COUNT,MAXVAL, MINVAL, PRODUCT Miscellaneous Functions that do the following:•Let you use assemblerinstructions in an executableprogram (ASM)•Check for pointer association(ASSOCIATED)•Check for end-of-file (EOF)•Return the class of afloating-point argument(FP_CLASS)•Count actual arguments passedto a routine (IARGCOUNT)•Return a pointer to an actualargument list for a routine(IARGPTR)•Return the INTEGER KIND thatwill hold an address(INT_PTR_KIND)•Test for Not-a-Number values(ISNAN)•Return the internal address ofa storage item (LOC)•Return a logical value of anargument (LOGICAL)•Allocate memory (MALLOC)•Return the upper 64 bits of a128-bit unsigned result(MULT_HIGH)•Return a disassociatedpointer (NULL)•Check for argument presence(PRESENT)•Convert a bit pattern(TRANSFER)1All of the numeric manipulation, and many of the numeric inquiry functions are defined by the model sets for integers and reals.2 The value of the argument does not have to be defined.3The Compaq Fortran processor character set is ASCII, so ACHAR = CHAR and IACHAR = ICHAR.4 For more information on bit functions, see Bit functions.。
fortran 函数总结
附录C 部分函数指令表(c). .(注解:本指令表只收集了部分常用指令, 有关全部指令请参照文档文件) + 加- 减* 矩阵乘数组乘 *.1. 通用指令^ 矩阵乘方数组乘方 ^.\ 反斜杠或左除在线帮助/ 斜杠或右除文档中关键词搜寻或.\ 数组除/.缺省变量名以及最新表达式的运算结果等号不等号从内存中清除变量和函数< 小于关闭> 大于退出<= 小于或等于把内存变量存入磁盘>= 大于或等于运行脚本文件逻辑与文件运行中的显示格式逻辑或显示版本逻辑非浮点运算溢出显示模式选择: 冒号列出工作内存中的变量名( ) 园括号文件编辑器[ ] 方括号变量类型{ } 花括号列出基本命令小数点 .设置数据输出格式, 逗号改变当前工作目录; 分号给出环境值注释号创建目录= 赋值符号显示当前工作目录' 引号执行表达式' 复数转置号转置号 '.最新表达式的运算结果2.运算符和特殊算符浮点误差容限, =2 -52≈2.22×10-16虚数单位= √(-1)正无穷大圆周率,π=3.14....3. 编程语言结构中止计算或循环终止最内循环同一起使用将控制转交给外层的或循环同一起使用同一起使用结束，，语句按规定次数重复执行语句条件执行语句可同一起使用暂停模式返回多个条件分支同一起使用不确定次数重复执行语句特定值计算函数特定值计算或多变量计算函数文件头定义全局变量检测变量是否为全局变量显示错误信息显示最近的错误信息按格式把数字转换为串显示警告信息4.基本数学函数反余弦反双曲余弦反余切反双曲余切反余割反双曲余割反正弦反双曲正弦反正切反双曲正切余弦双曲余弦余切双曲余切正弦双曲正弦正切双曲正切指数自然对数10 常用对数2 以2为底的对数平方根绝对值复数共轭复数虚部复数实部向上(正无穷大方向)取整向零方向取整向下(负无穷大方向)取整四舍五入取整符号函数降次排序误差函数补误差函数函数插值函数线性插值函数样条插值函数样条平滑函数样条函数方波函数符号函数将整数转换为双精度浮点数5.基本矩阵函数和操作单位阵全零矩阵全1 矩阵均匀分布随机阵生成随机矩阵线性等分向量对数等分向量矩阵对数运算矩阵元素累计乘矩阵元素累计和矩阵显示矩阵和文字内容确定向量的长度确定矩阵的维数创建对角阵或抽取对角向量找出非零元素1的下标矩阵变维90 矩阵逆时针旋转90度2 据全下标换算出单下标抽取下三角阵抽取上三角阵共轭矩阵伴随矩阵行列式的值矩阵或向量范数矩阵中非零元素个数清空向量或矩阵中的某个元素正交基矩阵秩矩阵迹矩阵条件数逆矩阵条件数矩阵的逆分解或高斯消元法伪逆分解变换求解线性方程方程矩阵特征多项式分解矩阵指数1 矩阵指数的逼近2 用泰勒级数求矩阵指数3 通过特征值和特征向量求矩阵指数计算一般矩阵函数矩阵对数矩阵平方根6. 特性值与奇异值矩阵特征值矩阵束特征值块矩阵, 广义特征向量正则化特征向量特征空间投影奇异值分解。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
附录C SCILAB 部分函数指令表(c)LIAMA. All rights reserved.(注解:本指令表只收集了部分常用指令, 有关全部指令请参照文档文件)+ 加- 减* 矩阵乘数组乘*.1. 通用指令^ 矩阵乘方数组乘方^.\ 反斜杠或左除help 在线帮助/ 斜杠或右除apropos 文档中关键词搜寻或.\ 数组除/.ans 缺省变量名以及最新表达式的运算结果== 等号~= 不等号clear 从内存中清除变量和函数< 小于exit 关闭SCILAB> 大于quit 退出SCILAB<= 小于或等于save 把内存变量存入磁盘>= 大于或等于exec 运行脚本文件&,and 逻辑与mode 文件运行中的显示格式|,or 逻辑或getversion 显示SCILAB 版本~,not 逻辑非ieee 浮点运算溢出显示模式选择: 冒号who 列出工作内存中的变量名( ) 园括号edit 文件编辑器[ ] 方括号type 变量类型{ } 花括号what 列出SCILAB 基本命令小数点 . format 设置数据输出格式, 逗号chdir 改变当前工作目录; 分号getenv 给出环境值// 注释号mkdir 创建目录= 赋值符号pwd 显示当前工作目录' 引号evstr 执行表达式' 复数转置号转置号'.ans 最新表达式的运算结果2.运算符和特殊算符%eps 浮点误差容限, =2-52≈2.22×10-16%i 虚数单位= √(-1)%inf 正无穷大%pi 圆周率,π=3.1415926535897....3. 编程语言结构abort 中止计算或循环break 终止最内循环case 同select 一起使用continue 将控制转交给外层的for或while循环else 同if一起使用elseif 同if一起使用end 结束for，while，if 语句for 按规定次数重复执行语句if 条件执行语句otherwise 可同switch 一起使用pause 暂停模式return 返回select 多个条件分支then 同if一起使用while 不确定次数重复执行语句eval 特定值计算feval 函数特定值计算或多变量计算function 函数文件头global 定义全局变量isglobal 检测变量是否为全局变量error 显示错误信息lasterror 显示最近的错误信息sprintf 按格式把数字转换为串warning 显示警告信息4.基本数学函数acos 反余弦acosh 反双曲余弦acot 反余切acoth 反双曲余切acsc 反余割acsch 反双曲余割asin 反正弦asinh 反双曲正弦atan 反正切atanh 反双曲正切cos 余弦cosh 双曲余弦cotg 余切coth 双曲余切sin 正弦sinh 双曲正弦tan 正切tanh 双曲正切exp 指数log 自然对数log10 常用对数log2 以2为底的对数sqrt 平方根abs 绝对值conj 复数共轭imag 复数虚部real 复数实部ceil 向上(正无穷大方向)取整fix 向零方向取整floor 向下(负无穷大方向)取整round 四舍五入取整sign 符号函数gsort 降次排序erf 误差函数erfc 补误差函数gamma gamma 函数interp 插值函数interpln 线性插值函数intsplin 样条插值函数smooth 样条平滑函数spline 样条函数quarewave 方波函数sign 符号函数double 将整数转换为双精度浮点数5.基本矩阵函数和操作eye 单位阵zeros 全零矩阵ones 全1 矩阵rand 均匀分布随机阵genmarkov 生成随机Markov 矩阵linspace 线性等分向量logspace 对数等分向量logm 矩阵对数运算cumprod 矩阵元素累计乘cumsum 矩阵元素累计和toeplitz Toeplitz 矩阵disp 显示矩阵和文字内容length 确定向量的长度size 确定矩阵的维数diag 创建对角阵或抽取对角向量find 找出非零元素1的下标matrix 矩阵变维rot90 矩阵逆时针旋转90度sub2ind 据全下标换算出单下标tril 抽取下三角阵triu 抽取上三角阵conj 共轭矩阵companion 伴随矩阵det 行列式的值norm 矩阵或向量范数nnz 矩阵中非零元素个数null 清空向量或矩阵中的某个元素orth 正交基rank 矩阵秩trace 矩阵迹cond 矩阵条件数rcond 逆矩阵条件数inv 矩阵的逆lu LU分解或高斯消元法pinv 伪逆qr QR分解givens Givens 变换linsolve 求解线性方程lyap Lyapunov 方程hess Hessenberg 矩阵poly 特征多项式schur Schur 分解expm 矩阵指数expm1 矩阵指数的Pade逼近expm2 用泰勒级数求矩阵指数expm3 通过特征值和特征向量求矩阵指数funm 计算一般矩阵函数logm 矩阵对数sqrtm 矩阵平方根6. 特性值与奇异值spec 矩阵特征值gspec 矩阵束特征值bdiag 块矩阵, 广义特征向量eigenmarkov 正则化Markov 特征向量pbig 特征空间投影svd 奇异值分解sva 奇异值分解近似7. 矩阵元素运算cumprod 元素累计积cumsum 元素累计和hist 统计频数直方图max 最大值mean 平均值median 中值min 最小值prod 元素积sort 由大到小排序std 标准差sum 元素和trapz 梯形数值积分corr 求相关系数或方差8. 稀疏矩阵运算sparse 稀疏矩阵(只存储非零元素)adj2sp 邻接矩阵转换为稀疏矩阵full 稀疏矩阵转换为全矩阵mtlb_sparse 将SCILAB 稀疏矩阵转换为MATLAB稀疏矩阵格式sp2adj 稀疏矩阵转换为邻接矩阵speye 稀疏矩阵方式单位阵sprand 稀疏矩阵方式随机矩阵spzeros 稀疏矩阵方式全零阵lufact 稀疏矩阵LU分解lusolve 稀疏矩阵方程求解spchol 稀疏矩阵Cholesky分解9. 输入输出函数diary 生成屏幕文本记录disp 变量显示file 文件管理input 用户键盘输入load 读已存的变量mclose 关闭文件mget 读二进制文件mgetl 按行读ASCII码文件mgetstr 读字符串中单个字mopen 打开文件mput 写二进制文件mfscanf 读ASCII 码文件print 将变量记录为文件read 读矩阵变量save 存变量为二进制文件strartup 启动文件write 按格式存文件xgetfile 对话方式获取文件路径x_dialog 建立Xwindow参数输入对话框Tk_Getvar 得到Tk文件变量Tk_EvalFile 执行Tk 文件10. 函数与函数库操作deff 在线定义函数edit 函数编辑器function 打开函数定义functions SCILAB 函数或对象genlib 在给定目录下建立所有文件的函数库get_function_path 读函数库的文件存储目录路径getd 读函数库中的全部文件getf 在文件中定义一个函数lib 函数库定义macro SCILAB函数或对象macrovar 输入变量个数newfun 输出变量个数11. 字符串操作code2str 将SCILAB数码转换为字符串convstr 字母大小转换emptystr 清空字符串grep 搜寻相同字符串part 字符提取str2code 将字符串转换为SCILAB数码string 字符串转换strings SCILAB 对象, 字符串strcat 连接字符strindex 字符串的字符位置搜寻strsubst 字符串中的字符替换12. 日期与时间date 日期getdate 读日期与时间timer CPU时间计时13. 二维图形函数plot2d 直角坐标下线性刻度曲线champ 2 维向量场champ1 由颜色箭头表示的2维向量场contour2d 等高线图errbar 曲线上增加误差范围框线条grayplot 应用颜色表示的表面xgrid 画坐标网格线histplot 统计频数直方图Matplot 散点图阵列14. 三维图形函数plot3d 三维表面plot3d1 用颜色或灰度表示的三维表面param3d 三维中单曲线param3d1 三维中多曲线contour 三维表面上的等高线图hist3d 三维表示的统计频数直方图geom3d 三维向二维上的投影15. 线条类图形xpoly 单线条或单多边形xpolys 多线条或多各多边形xrpoly 正多边形xsegs 非连接线段xfpoly 单个多边形内填充xfpolys 多个多边形内填充xrect 矩形xfrect 单个矩形内填充xrects 多个矩形内填充xarc 单个弧线段或弧园xarcs 多个弧线段或弧园xfarc 单个弧线段或弧园填充xfarcs 多个弧线段或弧园填充xarrows 多箭头16. 图形注释, 变换xstring 图形中字符xstringb 框内字符xtitle 图形标题xaxis 轴名标注plotframe 图形加框并画坐标网格线isoview 等尺寸比例显示(原图形窗口不改变)square 等尺寸比例显示(原图形窗口改变)xsetech 设置小窗口xchange 转换实数为图形象素坐标值subplot 设置多个子窗口17. 图形颜色及图形文字colormap 应用颜色图getcolor 交互式选择颜色图addcolor 增加新色于颜色图graycolormap 线性灰度图hotcolormap 热色(红到黄色)颜色图xset 图形显示方式设定xget 读当前图形显示方式设定getsymbol 交互式选择符号和尺寸18. 图形文件及图形文字xsave 将图形存储为文件xload 从磁盘中读出图形文件xbasimp 将图形按PS文件打印或存储为文件xs2fig 将图形生成Xfig 格式文件xbasc 取消图形窗及其相关内容xclear 清空图形窗driver 选择图形驱动器xinit 图形驱动器初始化xend 关闭图形xbasr 图形刷新replot 更改显示范围后的图形刷新xdel 关闭图形xname 改变当前图形窗名称19. 控制分析用图形bode 伯德图坐标gainplot 幅值图坐标(伯德图中的幅值图) nyquist 奈奎斯特图m_circle M－圆图chart 尼库拉斯图black Black－图evans 根轨迹图sgrid s 平面图plzr 零-极点图zgrid z 平面图20. 图形应用中的其它指令graphics 图形库指令表xclick 等待鼠标在图形上的点击输入locate 由鼠标点击读入图形中的多点位置坐标xgetmouse 由鼠标点击读入图形中的当前点位置坐标21. 系统与控制abcd 状态空间矩阵cont_mat 可控矩阵csim 线性系统时域响应dsimul 状态空间的离散时域响应feedback 反馈操作符flts 时域响应（离散、采样系统〕frep2tf 基于传递函数的频域响应freq 频域响应g_margin 幅值裕量imrep2ss 基于状态空间的脉冲响应lin 线性化操作lqe Kalman 滤波器lqg LQG补偿器lqr LQ补偿器ltitr 基于状态空间的离散时域响应obscont 基于观测器的控制器observer 观测器obsv_mat 观测矩阵p_margin 相位裕量phasemag 相位与幅值计算ppol 极点配置repfreq 频域响应ricc Riccati 方程rtitr 基于传递函数的离散时域响应sm2ss 系统矩阵到状态空间变换ss2ss 反馈连接的状态空间到状态空间变换ss2tf 状态空间到传递函数变换stabil 稳定性计算tf2ss 传递函数到状态空间变换time_id SISO系统最小方差辨识22. 鲁棒控制augment 被控对象增广操作bstap Hankel 矩阵近似ccontrg H∞控制器dhnorm 离散H∞范数h2norm H2 范数h_cl 闭环矩阵h_inf H∞控制器h_norm H∞范数hankelsv Hankel 矩阵奇异值leqr H∞控制器的LQ增益linf 无穷范数riccati Riccati 矩阵sensi 敏感函数23. 动态系统arma ARMA 模型arma2p 基于AR模型中获得多项式矩阵armac ARMAX 辨识arsimul ARMAX系统仿真noisegen 噪声信号发生器odedi 常微分方程仿真检测prbs_a 伪随机二进制序列发生器reglin 线性拟合24. 系统与控制实例artest Arnold 动态系统bifish 鱼群人口发展的离散时域模型boucle 具有观测器的动态系统相位图chaintest 生物链模型gpech 渔业模型fusee 登陆火箭问题lotest Lorennz 吸引子mine 采矿问题obscontl 可控可观系统portr3d 三维相位图portrait 二维相位图recur 双线性回归方程systems 动态系统tangent 动态系统的线性化tadinit 动态系统的交互初始化25. 非线性工具（优化与仿真〕bvode 边界值问题的常微分方程dasrt 隐式微分方程过零解dassl 代数微分方程datafit 基于测量数据的参数辨识derivative 导数计算fsolve 非线性函数过零解impl 线性微分方程int2d 二维定积分int3d 三维定积分intg 不定积分leastsq 非线性最小二乘法linpro 线性规划lmisolver 线性不等矩阵ode 常微分方程ode_discrete 离散常微分方程ode_root 常微分方程根解odedc 连续/离散常微分方程optim 非线性优化quapro 线性二次型规划semidef 半正定规划26. 多项式计算coeff 多项式系数coffg 多项式矩阵逆degree 多项式阶数denom 分母项derivat 有理矩阵求导determ 矩阵行列式值factors 因式分解hermit Hermit 型horner 多项式计算invr 有理矩阵逆lcm 最小公倍数ldiv 多项式矩阵长除numer 分子项pdiv 多项式矩阵除pol2des 多项式矩阵到表达式变换pol2str 多项式到字符串变换polfact 最小因式residu 余量roots 多项式根simp 多项式化简systmat 系统矩阵27. 信号处理%asn 椭圆积分%k Jacobi 完全椭圆积分%sn Jacobi 椭圆函数analpf 模拟量低通滤波器buttmag Butterworth 滤波器响应cepstrum 倒谱计算cheb1mag Chebyshev 一型响应cheb2mag Chebyshev 二型响应chepol Chebyshev 多项式convol 卷积corr 相关, 协方差cspect 谱估计(应用相关法)dft 离散富立叶变换fft 快速富立叶变换filter 滤波器建模fsfirlin FIR滤波器设计hank 协方差矩阵到Hankel矩阵变换hilb Hilbert 变换iir IIR数字滤波器intdec 信号采样率更改kalm Kalman 滤波器更新mese 最大熵谱估计mfft 多维快速富立叶变换mrfit 频率响应拟合phc Markov 过程srkf Kalman 滤波器平方根sskf 稳态Kalman 滤波器system 观测更新wfir 线性相位FIR滤波器weiener Weiener(维纳)滤波器window 对称窗函数yulewalk 最小二乘滤波器zpbutt Buthererworth 模拟滤波器zpch1 Chebyshev 模拟滤波器28. 音频信号analyze 音频信号频域图auread 读*.au 音频文件auwrite 写*.au 音频文件lin2mu 将线性信号转换为µ率码信号loadwave 取*.wav 音频文件mapsound 音频信号图示mu2lin 将µ率码信号转换为线性信号playsnd 音频信号播放savewave 存*.wav 音频文件wavread 读*.wav 音频文件wavwrite 写*.wav 音频文件29. 语言与数据转换工具ascii 字符串的ASCII码excel2sci 读ASCII 格式的Excel 文件fun2string 将SCILAB 函数生成ASCII 码mfile2sci 将MATLAB 的M 格式文件转换为SCI格式文件mtlb_load 取MATLAB第4版本文件中变量matlb_save 按MATLAB 第4 版本文件格式存变量pol2tex 将多项式转换为TeX格式sci2for 将SCILAB 函数转换为FORTRAN格式文件texprint 按TeX 格式输出SCILAB 对象translatepaths 将子目录下的所有MA TLAB 文件转换为SCI文件格式。