0编程日记.docx

2017年1月15日星期日 17:571、自写bessel函数，bes(m,x)real function bes(m,x) real::x integer::m,n,k integer fact n=1 bes=0 do k=0,n bes=bes+(-1)**k*(x/2)**(m+2*k)/(fact(k)*fact(m+k)) !print*,(-1)**k*(x/2)**(m+2*k) !print*,fact(k)*fact(m+k) !print*,(-1)**k*(x/2)**(m+2*k)/(fact(k)*fact(m+k)) !print*,bes enddo end note：m表示阶数，x必须为实数（即，bes(0,3.0)）。2、自写阶乘函数，fact(m)integer function fact(m) integer::m fact=1 if(m==0) then fact=1 else do i=1,m fact=fact*i enddo endifend note：m表示阶乘的变量3、系统bessel函数besjN(x)，besyN(x)，（N=0,1,2……）分别表示N阶Bessel函数和N阶诺依曼函数。同样，注意参数x为实数，即：besj0(3.0)√besj0(3)×4、hankel函数complex function han1(m,x) integer::m integer::pi=3.1415 real::x real::bes integer::fact complex,parameter::i=(0.,1.) if (m==0) then han1=besj0(x)+i*besy0(x) else han1=besj1(x)+i*besy1(x) endif endnote：m表示第几阶hankel函数，此处只取m=0或1。注意给函数赋值时，不能带参量。即：han1=……√han1(m,x)=……×2017年1月23日星期一 19:111、程序参数（1）流场参数：ubar（，平均流流速）rbar（，平均流密度）（2）计算观测点的位置参数（球坐标）：theta（，观测点极角，又叫仰角）r（R，观测点径向坐标）方位角为0图中z轴即为喷流中的x轴（极轴），因此方位角为0（3）计算声源点的位置参数（柱坐标）：x（）ry（）fai（）（4）湍流参数：ls（，长度尺度）taos（，时间尺度）ks（，湍动能）epis（，湍动能消散率）（5）经验系数：ecl（，长度尺度经验系数）ect（，时间尺度经验系数）ea（，第三个经验系数）Tam模型：（4）有名常量：c0（，无穷远处声速）2017年2月17日星期五 20:581、三角函数的角度单位是“弧度”，而不是“度”sin（90*3.14159/180）2、word《0编程version 1.0》 2017年2月18日星期六 10:071、数组（1）数组特性：数组名、维数（即下标的个数，1~7）、下标取值范围（上界、下界）、类型（由数组名类型决定）（2）数组的说明：integer s(1:10),a(1:4,1:6),b(3,4)dimension a(2,3),b(1:4,1:3)integer a,breal,dimension(2,3)::a,b(3,4)（3）a（2,3）是2行3列的二维表数据，存储是按列连续存放的数据块。a(1,1)a(2,1)a(1,2)a(2,2)a(1,3)a(2,3)a(3,3)（4）数组的输入输出：real,dimension(2,3)::aread*,ado i=1,2 do j=1,3 read*,a(i,j) enddoenddoread(*,*) ((a(i,j),j=1,3),i=1,2)注意，一个read语句可以分多行输入，但是每个read语句必须从新的一行开始读数。（5）给数组赋初值【1】用data命令real c,b,ddata c,b,d/2,3,4/dimension k(2,3)data ((k(i,j),j=1,3),i=1,2)/1,2,3,4,5,6/或：data k/1,4,2,5,3,6/dimension k(10)data k/10*-1/注意，data属于说明语句，但它可以放在end语句之前的任意行，其作用相当于赋值。【2】用数组赋值符赋初值real m(6)m=(/1,2,3,4,5,6/)real m(9) m=(/1,2,(k,k=1,9,2),1,2*3/)real a(5)a=20！所有的a(1)~a(5)都为20.000000 2017年2月19日星期日 13:271、数组运算允许把整个数组或数组片段作为一个单独的对象进行运算。例如当数组A,B,C有相同大