Fortran6创新.pptVIP

第六章 Fortran数据结构 * * * * 第六章 Fortran的数据结构 在计算机语言中，用数据类型来表示不同 的数据结构。高级语言一般有以下三类数据结 构：基本类型、构造类型和指针类型。 Fortran77提供基本类型和构造类型。 Fortran90提供指针类型。 Fortran数据的基本类型有整型、实型、 双精度型、复型和逻辑型、字符型。 一、双精度型数据 双精度数在计算机系统中以8个字节来存放，它 提供15~17位有效数字。 第六章 Fortran的数据结构 双精度常数要以指数形式表示，并以d代替e。 如： 1.1e2 1.1d2 0.0 0.0d0 1.0 1.0d0 双精度型数据，要先定义后使用。 double precision amount, income implicit double precision (r, x-z) real*8 imax, volt 求： 前200项乘积。 double precision i,ai,y y=1d0 i=1d0 10 if(i.le.100d0) then ai=4.0*i*i/(4.0*i*i-1) y=y*ai i=i+1 goto 10 end if write(*,*) y=,y end 二、复型数据 复型数据分实部虚部，使用8个字节来存放，必 须先定义后使用。 第六章 Fortran的数据结构 complex pose, flux implicit complex (r, x-z) 复型常数表示为： 1+2.5i (1.0, 2.5) -5.4i (0.0, -5.4) 7.8 (7.8, 0.0) 如果实部或虚部是表达式，则需通过cmplx函数转 换成复型数。 complex cordxy, trans trans=(3.0, 6.5) cordxy=cmplx(3.0*x, 6.5*y) cmplx(3.0)则 转成（3.0, 0.0) complex x1,x2 read(*,*) a,b,c d=b*b-4.0*a*c part1=-b/(2.0*a) part2=sqrt(abs(d))/(2.0*a) if(d.ge.0.0) then x1=cmplx(part1+part2, 0.0) x2=cmplx(part1-part2, 0.0) else x1=cmplx(part1,part2) x2=cmplx(part1,-part2) end if write(*,*) x1=,x1 write(*,*) x2=,x2 end 已知交流电路：U=220v, r0=10Ω, L0=10-3 H, r1=100 Ω, C1=100μF, r2=50 Ω, L2=10-2 H, C2=200 μF。 求：(1)等效电阻Z；(2)电流I, I1, I2 implicit real(l) implicit complex(u,i,z) parameter (omega=314.1592) u=(220.0,0.0) write(*,*) enter r0,l0,r1,c1,r2,l2,c2 read(*,*) r0,l0,r1,c1,r2,l2,c2 z0=cmplx(r0,omega*l0) z1=cmplx(r1,-1.0/(omega*c1)) z2=cmplx(r2,omega*l2-1.0/(omega*c2)) z12=(z1*z2)/(z1+z2) z=z0+z12 i=u/z i1=i*z12/z1 i2=i*z12/z2 write(*,*) z=,z write(*,*) i=,i write(*,*) i1=,i1 write(*,*) i2=,i2 end 如果x为复数，ABS(x)意为求模。 ATAN2(AIMAG(x), REAL(x))函数 的意思是： 可用来求相位。 线段AB：A(1, 1)，B(4.5, 4.5)； 求：AB长度及黄金分割点C坐标。 complex a,b,c real l read(*,*) a,b c=a+0.618*(b-a) l=abs(a-b) write(*,*) l=,l write(*,*) c=,c end 不同数值类型数据间运算遵循由低级向高级转化原 则：（低）整型?实型?双精度型?复型（高）。 三、字符型数据 第六章 Fortran