计算深海远场情况下的声源定位

计算深海远场情况下的声源定位，首先要编译海洋环境文件（包含阵元信息），利用Krakren程序计算本征函数和本征模态，再利用本程序的子程序计算拷贝向量，然后合成声源信号，最后根据匹配场的程序进行

clc;clear;close all;
%%%%这个匹配场为什么会有问题，感觉不理解，为什么那么多虚峰？
%%%参数    $$$$$$$很明显这是个错误的$$$$$声压公式不对
f=200;                           %声源频率
mi=1.03;                         %海水密度
zs1=20;                           %实际声源1深度
rs1=17000;                        %实际距离1
w=2*pi*f;                        %圆频率
i=sqrt（-1）;
Q=7;                             %所取简正波号数，因为只有前几号简正波起主要作用，省略一些应该对结果的影响不大
zy_ch=1;                        %实际接收水听器基阵的第一个阵元深度
zy_d=59;                        %实际接收水听器基阵的最后一个阵元深度
zy_jiange=1;                     %实际的阵元间隔，这和扫描的没有关系，扫描的是意义是依次假设声源在某个位置上
k_N=59;                          %阵元数目 
sm_r1=16000;                      %%%扫描距离的近点
sm_r2=20000;                     %%%扫描距离的远点
sm_chushi=1;                     %%%扫描深度的初始点
sm_sd=59;                        %%%扫描距离的最深点
sm_jiange=1;

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

     文件     236278  2011-07-24 08:55  匹配场单元定位\1.bmp

     文件     236278  2011-07-24 15:32  匹配场单元定位\3_4.bmp

     文件     236278  2011-07-24 15:32  匹配场单元定位\5_6.bmp

     文件     236278  2011-07-23 09:01  匹配场单元定位\bartlett_danyuan.bmp

     文件     236278  2011-07-24 08:54  匹配场单元定位\ck_xushi.bmp

     文件       2060  2014-04-17 08:24  匹配场单元定位\danyuandingwei_main.m

     文件     236278  2011-07-24 08:54  匹配场单元定位\jianzhengboxiangsutu.bmp

     文件        792  2011-07-23 08:53  匹配场单元定位\mode_fangzhenzhi.m

     文件     236278  2011-07-23 09:41  匹配场单元定位\mv_+10.bmp

     文件     236278  2011-07-24 08:50  匹配场单元定位\mv_+3.bmp

     文件     236278  2011-07-23 22:03  匹配场单元定位\mv_wuzao.bmp

     文件        542  2011-07-19 22:25  匹配场单元定位\pipeichang.m

     文件        548  2011-07-23 09:03  匹配场单元定位\pipeichang1.m

     文件       1042  2014-04-17 08:24  匹配场单元定位\rep_vector.m

     文件       1848  2011-07-22 15:54  匹配场单元定位\yellow_examk.txt

     文件     186648  2011-07-22 15:54  匹配场单元定位\yellow_examphi.txt

     目录          0  2014-04-17 08:26  匹配场单元定位

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

              2319982                    17