海岸动力学实验报告专业年级:港航2班姓名:薛辉学号:1014020135二〇一三年五月中国南京2目录实验一:波浪数据采集与波高统计实验⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯3一、实验目的⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯3二、实验要求⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯3三、实验过程⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯3四、结果分析⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯4五、实验结论⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯8实验二:波压力测量实验⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯9一、实验目的⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯9二、实验要求⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯9三、实验水文要素⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯10四、实验仪器⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯10五、实验过程⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯10六、结果分析⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯11七、实验结论⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯133实验一:波浪数据采集与波高统计实验一、试验目的了解波浪中规则波及不规则波的区别,波浪模型实验的一般方法,规则波波高、周期、不规则波高的统计方法。二、试验要求1、规则波及不规则波的测量与特征值统计。2、明确实验的目的,掌握实验原理,掌握基本仪器的使用,包括波浪数据采集系统及水槽造波机的使用方法,通过自己设计出不同波长波高的规则波和不规则波,参与造波数据及数据采集的全过程,了解波浪物理模型实验的最基本方法,正确处理实验数据,能通过处理采样数据文件系统各种累积频率波高,发现规律,得出实验结论。分析实验误差,提出减少误差方法3、试验报告的编写,要求报告能准确的反映试验目的、方法、过程及结论。三、试验过程试验中共设置四根波高传感器,四个同学为一组,每人采用其中一根传感器的数据计算波高,规则波采样时间为20s,不规则波采样时间为80s左右。规则波试验结果主要统计平均波高。波峰减波谷即为波高,将采集到的所有波高进行算术平均,得到规则波的平均波高。不规则波试4验结果主要统计有效波高。波峰减波谷即为波高,将采集到的所有波高进行排序,取前1/3大波进行算术平均,得到不规则波的有效波高。四、结果分析:本次实验使用fortran90语言编写计算程序。5程序截图(上部为运行程序,下部为运行结果)规则波的具体代码:programmainimplicitnonereal::m,e,f,tinteger::i,j,k,preal,dimension(2000)::areal,dimension(20)::b,creal,dimension(9)::hopen(1,file='1.txt')doi=1,2000read(1,*)a(i)enddoclose(1)b(1)=1j=1doi=1,1999if(a(i)*a(i+1)<0)thenj=j+1b(j)=iendifenddodoi=1,19doj=b(i)+1,b(i+1)-1if(abs(a(j))>abs(a(j-1)).and.abs(a(j))>abs(a(j+1)))thenc(i)=a(j)endifenddoenddoe=0doi=1,9h(i)=c(2*i)-c(2*i-1)e=e+h(i)enddot=e/9doi=1,9print'(1x,f7.3)',h(i)enddoprint*,'规则波的平均波高为:'print'(1x,f7.3)',tendprogram6规则波运行结果如下:波高依次为:15.67516.0379.0300.0005.76014.2009.94213.2615.102规则波的平均波高:9.884结果修正:运行结果中出现0.000,所以要对结果进行修正,修正后的规则波的平均波高为:11.126规则波的波形图不规则波的程序及结果截屏如下7不规则波的程序代码:programmainimplicitnonereal,dimension(8000)::areal,dimension(114)::b,creal,dimension(56)::hinteger::i,jreal::e,topen(2,file='2.txt')8doi=1,8000read(2,*)a(i)enddoclose(2)b(1)=1j=1doi=1,7999if(a(i)*a(i+1)<0)thenj=j+1b(j)=iendifenddodoi=1,113doj=b(i)+1,b(i+1)-1if(abs(a(j))>abs(a(j-1)).and.abs(a(j))>abs(a(j+1)))thenc(i)=a(j)endifenddoenddodoi=1,56h(i)=c(2*i-1)-c(2*i)enddoprint*,'不规则波的波高分别为:'doi=1,14print'(1x,4f8.3)',h(4*i-3),h(4*i-2),h(4*i-1),h(4*i)enddodoi=1,55doj=i+1,56if(h(i)