- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验二分层介质中的声线传播路径研究.doc
实验二分层介质中的声线传播路径研究
一、 实验目的
1、 了解海水介质的垂直分层特性;
2、 掌握Snell折射定律和声线弯曲规律;
3、 掌握基本的声线传播轨迹的搜索方法;
二、 实验仪器
计算机
三、实验原理
在本课程的研究中,一般只考虑海水介质的垂玄分层的特性,叩声速(折射率)不随水平方向变化, 仅是海水深度的断数。若令x、y代表水平坐标,z代表垂直坐标,则在分层介质中:
c(x, y, Z)= c(z) ( 1 )
z) = n⑵ (2)
除了超远程肖传播问题外,海水介质的分层模型是对实际海洋介质的-种接近的理想模型。本课程中 所获得的结果,都是在分层介质模型下得到的。
3.1 Snell折射定律和声线弯曲
已知射线声学所遵循的Snell定律为:
cos 7 cos二常数
cos 7 cos
二常数
(3)
a为声传播方向与水平坐标ox的夹角,称为掠射角;C为该处声速。々和Co为声线出射处的对应值。 若和声线的垂直分布c⑵已知,可以按照Snell定律求出海洋中任意深度处声线的掠射角,也就确定了任 意深度处声波传播方向。不同的起始掠射角他求出不同的Q值,对应不同的声线轨迹。负梯度下声速随深 度増加而下降,掠射角a随深度増加而増加,声线弯向海底;正梯度下声速随深度增加而增加,掠射角a 随深度増加而减小,声线弯向海而。
3.2声线轨迹
二维平面问题的声线曲率的表达式为:
d0 sind de cos a de
= = ?—— (4)
ds c dz c dz
(4)式中,〃为声线入射角,即声线长度微元ds与垂直轴线z的夹角;c为微元出处的声速。由Snell定 律可知,当起始值给定时(血和co),比値cosa/c等于常数。对于恒定声速梯度而言c=co(l+az), a为和对 声速梯度(a为常数),则血/能=cg=常数,即恒定声速梯度情况下声线Illi率到处相等,轨迹是圆弧。
3.3声线轨迹方程
对于恒定声速梯度c=c°(l+az), 67=l/c0-dc/dzmo若声线在海而以掠射角闪)=0 (0=n/2)出射,该声线曲 率半径等于:
ds
1
d0
进-?步可写出以《o=O掠射的声线轨迹方程:
对丁?声源在海面以任意掠射角cq出射的声线,同样可给出其声线轨迹方程:
WC0SQ]丿(7)
WC0SQ]丿
(7)
3.4声线经过的水平距离
若声源位于(圮习)处,接收点位于(x,z)处,声速按c=c⑵分布,可以从下列积分求出声线经过的水平距 离:
(8)
根据Snell定律可导出tancr = COS? Q| ]/cOSQ| ,其中〃⑵=c(zj/c⑵。因而水平距离:
(9)
(9)
对于恒定声速梯度情况,轨迹是一圆弧,玄接从声线轨迹图来求水平距离x更为方便。
3.5声线经过的水平距离
声线因深度而异,声线经过微元山距离所需要的时间dt=ds/c,声线从习声度传播到z深度所需要的时 间等于:
根据 Snell 定律,Csin a — 2— COS2 ax,其I1 n=c(Z])/ c(z)= c/ c,贝II: n
(II)1『 A72(z)dz c\ 1 ^772(Z)-COS2 ax
(II)
c sin ci
(II )式为计算声线传播时间的一般表示式。在恒定声速梯度下,从Snell定律可求得dz = — - da,
g cos ax
da代入公式(10)中得到:
da
(12)
◎ COS6Z
⑷和。分别为深度习和z处的声线掠射角。对(12)式积分可得:
. 1 + sina
. 1 + sina
In
I 1 - sina 丿
(a tC\
tan
-In
辽4;
/ 、
g
tan
鱼+兰
i
2 4,
丿」
(13)
四、实验内容
画出仿真条件下,从声源出发到接收点处三条主要声线的传播规律,三条声线分别肓达
声信号、海而一次反射信号、海底一次反射信号。 仿真条件1:
20iiL1 Out40m
20iiL
1 Out
40m
TOC \o 1-5 \h \z V lkm
仿真条件2:
lOui 20in *
I-. 30 川
(? …40(il
50iii
五、实验步骤
mk: function y=mk(xl,yl,x2,y2z x)
k=(y2-yl)/(x2-xl);
y=k* (x-x2)+y2;
End
mp: function [mz n]=mp (az b)
syms xy;
eql=(y-b)/(x-a)-100;
eq2=sqrt (xA2+ (y-50) A2) -sqrt (「2+ (b-50)八2);
s=solve(eql,eq2z x,y);
m=double (s.x);
n=double (s.y);
end clear;clc;
axis ( [0 1000 0 50]);
文档评论(0)