射线声学bellhop程序

本资源主要讲解射线声学的matlab计算程序BELLHO的详细使用说明，压缩包的文档介绍了程序中各参数的 设置选择，并给出了深海和浅海典型声信道的计算结果以供参考

 % arr
clear all;
close all;
clc;

dc=load（‘munk.txt‘）;
D=dc（end1）;
R=30;    %最远距离km
N=length（dc（:1））;
F=300;    %频率
k=0;
for f=F
    k=k+1
    fname=[‘Test‘num2str（f）‘.env‘];
    fp=fopen（fname‘wt‘）;
    fprintf（fp‘‘‘haha‘‘\n‘）;
    fprintf（fp‘%f\n‘f）;
    fprintf（fp‘%d\n‘1）;
    fprintf（fp‘‘‘SVM‘‘\n‘）;
    tf=f/1000;
%     alphaw=（3.3e-3+0.11*tf^2/（1+tf^2）+44*tf^2/（4100+tf^2）+3e-4*tf^2）/1000;
    alphaw=0;
    fprintf（fp‘%d %f %f\n‘00.0D）;
    SSP=[dczeros（N1）1.024*ones（N1）alphaw*ones（N1）zeros（N1）];
    fprintf（fp‘%f %f %f %f %f %f\n‘SSP‘）;
    fprintf（fp‘%s %f\n‘‘A‘0.0）;
%     alphas=0.2*tf^0.6/（1546/f）;
    alphas=0;
    fprintf（fp‘%f %f %f %f %f %f\n‘[D  1700.000  0.000000  1.800  alphas  0.000000]）;
    fprintf（fp‘%d\n‘1）;       % 声源数目
    fprintf（fp‘%f\n‘1000.0）;  % 声源深度
    fprintf（fp‘%d\n‘1）;       % 接收点数目
    fprintf（fp‘%f\n‘1000.0）;  % 接收点深度
    fprintf（fp‘%d\n‘1）;       % r方向取点数目
    fprintf（fp‘%f  /\n‘R）;    % r方向取点距离
    fprintf（fp‘‘‘AG‘‘\n‘）;
    fprintf（fp‘%d\n‘100000）;
    fprintf（fp‘%f %f /\n‘[-89.000000 89.000000]）;
    fprintf（fp‘%f %f %f/\n‘[0.0 D+D*0.01 R+R*0.01]）;
    fclose（fp）;

end
bellhop（fname（1:end-4））;
filename=[fname（1:end-4）‘.arr‘];
plotarr（ filename 11 1 ）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1339  2017-05-09 13:42  OASES软件的使用\arr.asv

     文件       1340  2017-05-09 14:18  OASES软件的使用\arr.m

     文件     876049  2010-09-01 17:56  OASES软件的使用\bellhop.exe

     文件        212  2010-01-26 10:01  OASES软件的使用\bellhop.m

     文件       1574  2017-05-08 18:44  OASES软件的使用\chuanbosunshi.asv

     文件       1579  2018-05-29 10:12  OASES软件的使用\chuanbosunshi.m

     文件       4803  2017-04-12 14:57  OASES软件的使用\deepsurface.txt

     文件       1622  2017-03-15 19:31  OASES软件的使用\envcalculation1.asv

     文件       1848  2017-04-28 08:33  OASES软件的使用\envcalculation1.m

     文件       1801  2017-03-10 09:05  OASES软件的使用\envcalculation2_bty.m

     文件       2502  2017-06-09 08:48  OASES软件的使用\munk.txt

     文件       1658  2017-05-08 18:43  OASES软件的使用\plotarr.m

     文件       3393  2015-03-13 00:44  OASES软件的使用\plotray.m

     文件       1434  2017-05-08 17:01  OASES软件的使用\ray.asv

     文件       1414  2017-06-09 09:17  OASES软件的使用\ray.m

     文件       2588  2009-04-21 18:25  OASES软件的使用\read_arrivals_asc.m

     文件       2817  2010-06-13 16:33  OASES软件的使用\read_shd_bin.m

     文件       2113  2017-05-07 18:09  OASES软件的使用\ssp.asv

     文件       2265  2017-05-11 20:57  OASES软件的使用\ssp.m

     文件       2478  2017-05-09 14:17  OASES软件的使用\Test1000.arr

     文件       6386  2017-05-09 14:14  OASES软件的使用\Test1000.env

     文件       6967  2017-05-09 14:17  OASES软件的使用\Test1000.prt

     文件        358  2017-06-09 09:06  OASES软件的使用\Test300.arr

     文件       6406  2017-06-09 10:36  OASES软件的使用\Test300.env

     文件       4857  2017-06-09 10:36  OASES软件的使用\Test300.prt

     文件    1864457  2017-06-09 10:36  OASES软件的使用\Test300.ray

     文件   24095496  2017-05-11 20:54  OASES软件的使用\Test300.shd

     文件     457195  2015-03-14 10:14  OASES软件的使用\第3讲 OASES声学软件使用.docx

     目录          0  2018-05-17 17:09  OASES软件的使用

----------- ---------  ---------- -----  ----

............此处省略2个文件信息