专利名称：：一种海洋环境与声场不确实性的表征和传递的快速计算方法
技术领域：
：本发明涉及声学领域，具体涉及一种海洋环境与声场不确实性的表征和传递的快速计算方法。
背景技术：
：由于海洋环境参量例如海深、声速剖面、海底地声特性等不确实，即使发射确定信号，海洋中接收点声压场也表现为随机起伏场在空间和时间上都具有不确实性和变异性。而先进的声纳信号处理都基于传播模型对海洋中传播的信号作预测，再与接收的信号作拷贝相关处理。由于处理中使用确定的传播模型，即输入一组特定的环境参量，得到确定的预测声压场。这样的传播模型不能反映空时变异和不确实的实际声场，往往导致声纳信号处理的性能不稳定。因而，寻求海洋环境和海洋声场的不确实性表征和传递规律是研究减弱、消除和进一步利用不确实性对声纳处理性能影响的基础。数学上，不确实性定义为误差的概率密度函数(probabilitydensityfunctionPDF)。误差是指实际场与估计/测量场的差。在基于模型的预测问题中，初始条件误差、模型及边界条件的误差都会影响预测的精度，引起预测声场的不确实性，是初始误差和模型集成过程引入的误差的综合影响。在海洋环境中，对于近、中距离，通常我们认为声传播模型的精度已足够高，对预测声场不确实性的贡献可以忽略，因而认为引起预测声场不确实性的主要因素是作为初始条件和边界条件的海洋环境知识的不确实性。环境知识不确实性主要来fi于水体，内波、潮、锋、涡、流及其动力特性引起的声速梯度变化。边界不确实性主要来自于界面，表面波，海底粗糙度、地形变化、海底构成等。本发明中，不确实性者，指不完备的环境知识及所导致的声场的随机性，它们的定量度量用一种正交随机多项式展开(PCE)所形成的随机谱作表征。本发明土要针对假定声场预测模型是精确的，由于海洋环境知识不完备所引起的不确实性导致声场不确实性的问题。不确实声场建模将海洋环境及在其中传播的声场两部分合在一起考虑，将不确实性作为一个必要的组成嵌入，并提供它对声场影响的概率度量。波方程是单变量(声压场)和单参量(声速场)方程，环境部分不确实，通过声速场耦合到声场部分，建立一组耦合差分方程，既描述声场的传播，又描述所关联的不确实性的传播。
发明内容本发明所要解决的技术问题之一是环境参量不确实性和声场不确实性的表征方法。用于驱动声传播模型的环境参量包括水体声速梯度剖面，海深，沉积层声速梯度、沉积层密度、沉积层衰减、沉积层厚度，基底声速、基底密度、基底衰减，用一个向量表示。不确实环境参量由一个确知的向量及其邻域构成，邻域的大小决定不确实性程度。同样一个基阵采集的随机声场也由一个确知的声场向量与邻域构成。不确实性的领域由多项式混沌展开作随机谱表征。本发明所要解决的技术问题之二是环境参量的不确实性如何传递到声场不确实性。声传播模型，例如射线模型、抛物方程模型、简正波模型等，能够解决由确定的环境参量如何产生确定的声场问题，而不能直接解决不确实性由环境参量到声场的传递，因而也无法直接给出用于不确实信号检测所需的不确实驾驶向量空间。将不确实性嵌入到传播模型产生随机的传播模型既描述声场的传播也实现不确实性从环境参量到声场的传递。为解决以上技术问题，本发明是提出以下技术方案实现的-本发明所述的这种海洋环境与声场不确实性的表征与传递的实现方法，步骤如下1.海洋环境不确实性和声场不确实性的表征(1)海洋环境不确实性表征用向量b表征环境参量(由水中声速梯度，海深，沉积层声速梯度，沉积层密度，沉积层衰减，沉积层厚度，基底声速，基底密度，基底衰减构成)。用b。表示我们利用测量或反演方法获得的己知环境参量，作为不确实环境参量的均值，不确实环境参量采用参量不确实性域"S(b。)来表征，它表示了环境参量落在bo附近的一个领域之中，用概率密度函数p(b)表征不确实环境参量的分布。S(b。)表征了环境参量不确实性的范围。下面对不确实环境参量作多项式混沌展开随机谱表征。歩骤一将不确实性环境参量空间(b。)表示成随机变量b;步骤二将不确实环境参量b分解为确实性部分和随机扰动部分b(r，z)=b。+Sb(r,z)(1)其中b。为确定性部分，5b为随机扰动部分，r，z分别表示距离和深度。步骤三根据环境参量b的概率概率密度分布p(b)，对随机扰动部分作多项式混沌展开-5b(r,Z;e)=|>s(r,Z)As(,,(2)式中As(^0))是对应环境参量b概率密度分布的随机多项式混沌基函数，9为随机参量空间。例如，b是高斯分布，As(《(9))为Gauss-Hermite多项式；b是均匀分布，As(《(0))为Legendre多项式；b是Y分布，A^(e))为Lague汀e多项式，b是(3分布，As(《(G))为Jacohi多项式，b是泊松分布，As(《(e))为Charlier多项式，等等。ocs(r,z)是随机扰动环境参量投影到随机多项式混沌基函数得到的确定性系数。步骤四对无穷多阶多项式混沌展开作有限截取对不同分布的环境参量利用对应的多项式混沌基函数展开，展开的阶数是有限的。对于高斯分布分布的b，只要一阶就够了，其中ot。(r，z)等于0，ot"r,z)二cT，式中cj是随机扰动环境参量的标准离差。具有类似的特性，其它分布的环境参量也可以利用投影方法得到有限阶(阶数为S)的确定性系数。步骤五多项式混沌项数确定随机环境参量b由多个元素构成，b:^[b,^…bn]，根据最小均方误差准则，n维随机变量b所必需的多项式项数为(n+S)!/n!S!。下面的过程我们针对11=1时的情况进行讨论。其中的环境参量选取水层的声速梯度，令b(r,z,9)=c(r,z,9)=c。(z)+5c(r，z,9)(3)(2)声场不确实性表征①确实的环境参量b。驱动声传播模型产生确实的声场p。；②不确实环境参呈空间(b。)对应不确实声场空间吖p。J;③用概率密度函数p(p)表征不确实声场的分布；④^fp。J表征了声场不确实性的范围。同样，对不确实声场也作多项式混沌展开随机谱表征。步骤一不确实声场作多项式混沌展开p(r,Z;e)=|>q(r,Z)Aq(,(4)步骤二对无穷多阶多项式混沌展开作有限截取根据最小均方误差准则，对无穷多阶多项式混沌展开作有限阶Q的截取p(r，z;e)-lx(r,z)Aq(柳)(5)q=0式中p(r,z;0)为预测的随机卢压场。2.环境参量不确实性到声场不确实性的传递方法将环境参呈与预测声场的多项式随机谱表征与确实性波方程相结合，得到嵌入不确实性的随机波方程，实现环境参量不确实性到声场不确实性的传递。具体步骤如下--(1)确定性波方程到随机波方程步骤一基于简正波传播模型的确定性波方程<formula>formulaseeoriginaldocumentpage10</formula>式中V表示拉普拉氏算子，k是波数，k=c0/c(z),co为角频率。步骤二随机扰动声速作二项式展开对随机扰动声速作二项式展开，并保留一阶项，则不确实声速的波数可近似为，<formula>formulaseeoriginaldocumentpage10</formula>(6)<formula>formulaseeoriginaldocumentpage10</formula>(7)步骤三确定性波方程推广到随机波方程将近似的波数代入到确定性波方程，将确定性波方程推广到随机波方程<formula>formulaseeoriginaldocumentpage10</formula>(8)将式(3)和式(4)代入(8)式，并对等式两边同时乘以八,g(9))，然后通过系总平均将各项投影到基空间，可得<formula>formulaseeoriginaldocumentpage10</formula><formula>formulaseeoriginaldocumentpage10</formula>将混沌基函数的正交性应用子上述方程，得到一组声场不确实性系数的耦合差分方程，<formula>formulaseeoriginaldocumentpage10</formula>以高斯分布随机声速为例进行下面公式的推导。5c(r，z,e)=ocflHefl+aiHei+.=ct《，g是均值为0，方差为1的高斯随机变量。将5c代入(10)式给出第1个不确实系数的耦合方程8、(r,z)3、(r,z)2,、2k2<jf,、/TTeTT\n,"、二2+丄2+k、(r,z)—~^;p^丄八(r,z)(He《He,〉=0(11)利用Hermite基函数的回归特性和正交性，计算耦合项的系总平均，得到系数^的偏微分回归方程<formula>formulaseeoriginaldocumentpage10</formula>式中D三K2，Es2K2a/c。。用向量-矩阵形式，上述偏微分方程可化简为~f+~^"+Ay=0式中，A为环境参量不确实性到声场不确实性的耦合方程。A是对称的三角阵。(2)求解随机波方程得到不确实声场步骤一耦合矩阵A特征值分解q=l式中G二[(g3,(gJ,…]为耦合矩阵A的特征向量，GG—'-I，(13)式可以表示为令f二G—、，上式可以简化为(13)(14)(15)(16)5r5z式中，Q二G爿AG为对角阵，对角阵元素为特征值、，^，…。利用特征值和特征向量分解技术，得到第1个确定系数的波方程44+^=0(17)步骤二确定系数求解等式(17)与标准波方程(6)有相同的形式。冈而求解过程与标准的波方程求解过程是一样的。f,的解可以表示为式中、m、M^(z)分别是水平波数和模深度函数，由简正波传播模型求解得到。a』表示矩阵A的特征向量的逆的第一列中的第w个元素。在声速梯度不确实情况下，模深度函数没有受到不确实声速的影响，但声速梯度不确实导致了水平波数的不确实，krmkz2m，其中k^是深度波数。根据yi与夂的关系式，系数^是f,的线性组合<formula>formulaseeoriginaldocumentpage12</formula>步骤三Q个多项式混沌基函数的确定性系数加权累加得到不确实性声场最后，在远场近似下的不确实性声场可表示为:<formula>formulaseeoriginaldocumentpage12</formula>若SSP是确实的，则上式中只有y。是不等于0的，收敛到声速确实时的声场简正模表达式^<formula>formulaseeoriginaldocumentpage12</formula>步骤四不确实性声场的均值和相关函数由多项式混沌展开的不确实声场是随机声场，其均值可表示为<formula>formulaseeoriginaldocumentpage12</formula>相关函数为<formula>formulaseeoriginaldocumentpage12</formula>本发明能带来以下有益效果(1)对海洋环境不确实性和声场不确实性进行了定量化表征，分别表示为海洋参量不确实邻域空间(b。)和声场不确实性邻域空间^fp。J，由概率密度函数p(b)和p(p)定量化描述，并将不确实性部分表征为多项式混沌随机谱，为减弱、消除和利用不确实性的信号处理所需的驾驶向量邻域空间^(a。)提供输入。(2)不确实性海洋环境参量和不确实性的声场的多项式混沌随机谱表示与简正波声传播模型相结合，得到嵌入不确实性的随机波方程，提供了不确实性从海洋环境参量到声场的耦合和传递。通过这种传递方法比通常采用的MonteCarlo方法可以提高计算效率，节省计算时间。利用MonteCarlo方法获得(21)和(22)的结果，需耍进行至少3000次传播模型的运行，而利用多项式混沌展开方法，只需耍运行Q次传播模型，而Q《3000。图1海洋环境与声场不确实性的表征和传递结构框2(a)确实声场与高斯分布不确实声速N(O，2.52)导致的100Hz不确实声场的MC结果与PCE结果比较:确实声场图2(b)确实声场与高斯分布不确实声速N(O,2.52)导致的100Hz不确实声场的MC结果与PCE结果比较:不确实声场MC仿真图2(c)确实声场与高斯分布不确实声速N(O,2.52)导致的100Hz不确实声场的MC结果与PCE结果比较:不确实声场PCE仿真图3(a)高斯分布不确实声速cN(0，2.52)100Hz声场的均值、方差和归一化自相关均值图3(b)高斯分布不确实声速cN(0，2.52)100Hz声场的均值、方差和归一化自相关方差图3(c)高斯分布不确实声速cN(0，2.52)100Hz声场的均值、方差和归一化自相关归一化自相关图4(a)确实声场与高斯分布不确实声速N(O，2.。导致的200Hz不确实声场的MC结果与PCE结果比较:确实声场图4(b)确实声场与高斯分布不确实声速N(O，2.5"导致的200Hz不确实声场的MC结果与PCE结果比较:不确实声场MC仿真图4(c)确实声场与高斯分布不确实声速N(O,2.52)导致的200Hz不确实声场的MC结果与PCE结果比较:不确实声场PCE仿真图5(a)高斯分布不确实声速cN(0，2.52)200Hz声场的均值、方差和归一化自相关均值图5(b)髙斯分布不确实声速cN(0，2.52)200Hz声场的均值、方差和归一化自相关方差图5(c)高斯分布不确实声速cN(0,2.52)200Hz声场的均值、方差和归一化自相关归一化自相关图6(a)确实声场与高斯分布不确实声速N(0,50导致的100Hz不确实声场的MC结果与PCE结果比较确实声场图6(b)确实声场与高斯分布不确实声速N(0，5"导致的100Hz不确实声场的MC结果与PCE结果比较不确实声场MC仿真图6(c)确实声场与高斯分布不确实声速N(O，52)导致的100Hz不确实声场的MC结果与PCE结果比较不确实声场PCE仿真图7(a)高斯分布不确实声速cN(0,52)100Hz声场的均值、方差和归一化自相关均值图7(b)高斯分布不确实声速cN(0，52)100Hz声场的均值、方差和归一化自相关方差图7(c)高斯分布不确实声速cN(0,52)100Hz声场的均值、方差和归一化自相关归一化自相关图8(a)均匀分布不确实声速cN(0，2.52)100Hz声场的均值、方差和归一化自相关均值图8(b)均匀分布不确实声速cN(0，2.52)100Hz声场的均值、方差和归一化自相关方差图8(c)均匀分布不确实声速cN(0，2.52)100Hz声场的均值、方差和归一化自相关归一化自相关图9(a)均匀分布不确实声速cN(0,2.52)200Hz声场的均值、方差和归一化自相关均值图9(b)均匀分布不确实声速cN(0，2.52)200Hz声场的均值、方差和归一化自相关方差图9(c)均匀分布不确实声速cN(0，2.52)200Hz声场的均值、方差和归一化自相关归一化自相关图IO(a)均匀分布不确实声速cN(0，52)100Hz声场的均值、方差和归一化自相关均值图IO(b)均匀分布不确实声速cN(0，52)100Hz声场的均值、方差和归一化自相关方差图IO(c)均匀分布不确实声速cN(0,52)100Hz声场的均值、方差和归一化自相关归一化自相关具体实施例方式下面结合具体实施例子和附图对本发明做进一步的描述图l是本发明的海洋环境与声场不确实性的表征和传递结构框图。从图中可以看出本发明的实施过程对接收阵接收的阵数据，首先进行放大、滤波，然后进行A/D转换，将接收的模拟信号进行数字采样。对采样后的离散随机声压场与预测随机声压场作拷贝相关处理，处理后的结果送给显示，给出距离和深度两维的能量分布。最后，根据处理结果进行有无目标的判决和目标距离/深度的估计。其中预测随机声压场由本发明提供不确实海洋环境参量经多项式混沌展开后与确定性波方程相结合，得到将不确实性由海洋环境参量耦合到声场的随机波方程，对随机波方程的耦合矩阵作特征值分解，并作适当变换，得到一组确定性系数方程，确定性的解与多项式混沌线性组合得到不确实声场，不确实声场作为不确实驾驶向量的输入。实施例子采用的浅海环境参量为Pekeris波导，深度200米，垂直布阵40元，阵元间距5米，声源位丁.距离5km，深度50m。共实现3000次MonteCarlo仿真。实施例子一不确实声速剖面服从高斯分布当不确实声速剖面服从高斯分布时，式(13)中的矩阵A为对称三角阵，其中对角线给出确定性部分的信息，而非对角线给出了不确实环境参量与不确实声场相耦合的信息。例如取不确实声压场的五阶多项式混沌展开，则A为6X6的矩阵，<formula>formulaseeoriginaldocumentpage14</formula>(23)图2给出了f=100Hz时确实声场与不确实性声速服从高斯分布的垂直阵采集的随机声场的实部与虚部，与MonteCarlo仿真结果的比较。(a)声速为1500m/s时垂直阵采集的声场(实部与虚部)，(b)不确实声速服从高斯分布N(O，2.52)时MonteCarlo仿真结果，(c)不确实声速服从高斯分布N(0，2.52)时IOO阶PCE展开的仿真结果。图3给出了图2中接收水听器在50m深度处不确实声场的均值、方差和归一化自相关。图4给出了f=200Hz时确实声场与不确实性声速服从高斯分布N(0，2.52)时的垂直阵采集的随机声场的实部与虚部，与MonteCarlo仿真结果的比较，Q=300。图5给出了图4中接收水听器在5(k深度处不确实声场的均值、方差和归一化自相关。图6给出了f=100Hz时确实声场与不确实性声速服从高斯分布N(0，52)时垂直阵采集的随机声场的实部与虚部，与MonteCarlo仿真结果的比较，Q=200。图7给出了图2中接收水听器在50m深度处不确实声场的均值、方差和归一化自相关。表1给出了f=100Hz当不确实声速服从不同方差高斯分布时所需的最小多项式混沌阶数。实施例子二不确实声速剖面服从均匀分布式(13)给出的随机波方程不仅适合高斯分布，也适合其它分布，如均匀分布，所不同的只是A的取值。不确实声速均匀分布时，取不确实声压场的五阶多项式混沌展开，A仍为6X6的矩阵，<formula>formulaseeoriginaldocumentpage15</formula>(24)<formula>formulaseeoriginaldocumentpage15</formula>图8给出了f=100Hz、不确实性声速服从均匀分布U(O，2.52)、接收水听器在50m深度处时不确实声场的均值、方差和归一化自相关，Q=25。图9给出了f=200Hz、不确实性声速服从均匀分布U(0,2.5:)、接收水听器在50m深度处时不确实声场的均值、方差和归一化自相关，Q=50。图IO给出了f=100Hz、不确实性声速服从均匀分布U(O，52)、接收水听器在50m深度处时不确实声场的均值、方差和归一化自相关，Q=50。表2给出了f=100Hz不确实声场服从不同方差均匀分布时所需的最小多项式混沌阶数。表3给出了当均匀分布不确实声场标准离差为5m/s时，不同频率所需的最小多项式混沌展开阶数。从以上处理结果可以看山，无论从不确实声场的实部和虚部，或是不确实声场的均值、方差和归一化自相关，多项式混沌展开的不确实声场与MonteCarlo仿真结果都是一致的。而多项式混沌展开方法所需的计算量要较MonteCarlo方法小得多，在上述示例中，不确实声速取值大到标准离差5m/s时，在频率f=200Hz时，所需的多项式混沌阶数为300,计算量是3000次MonteCarlo方法的1/10。如果环境参量由l增大到n,MonteCarlo方法的计算量是3000X2n，而多项式混沌方法则是(n+Q)!/n!Q!，计算的有效性更明显。同时本发明提供的方法还提供了不确实声速是如何传递到不确实声场的内在机理。表1不确实声速服从高斯分布时不确实声场所需的最小多项式混沌阶数<table>tableseeoriginaldocumentpage16</column></row><table>以上对本发明的描述不具有限制性，如果本领域的普通技术人员受其启示，在不脱离本发明权利耍求的保护的情况，作出本发明的其它结构变形和实施方式，均属于本发明的保护范围。权利要求1、一种海洋环境与声场不确实性的表征和传递的快速计算方法，其特征在于一、海洋环境不确实性和声场不确实性的表征(1)海洋环境不确实性表征用向量b表征环境参量，用b0表示我们利用测量或反演方法获得的已知环境参量，作为不确实环境参量的均值，不确实环境参量采用参量不确实性域id="icf0001"file="A2009101005950002C1.tif"wi="10"he="4"top="66"left="109"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>来表征，它表示了环境参量落在b0附近的一个领域之中，用概率密度函数p(b)表征不确实环境参量的分布，id="icf0002"file="A2009101005950002C2.tif"wi="10"he="4"top="77"left="133"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>表征了环境参量不确实性的范围，下面对不确实环境参量作多项式混沌展开随机谱表征；步骤一将不确实性环境参量空间id="icf0003"file="A2009101005950002C3.tif"wi="10"he="4"top="97"left="87"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>表示成随机变量b；步骤二将不确实环境参量b分解为确实性部分和随机扰动部分b(r，z)＝b0+δb(r，z)(1)其中b0为确定性部分，δb为随机扰动部分，r，z分别表示距离和深度。步骤三根据环境参量b的概率概率密度分布p(b)，对随机扰动部分作多项式混沌展开<mathsid="math0001"num="0001"><math><![CDATA[<mrow><mi>δb</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>s</mi><mo>=</mo><mn>0</mn></mrow><mo>∞</mo></munderover><msub><mi>α</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>Λ</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>ξ</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中Λs(ξ(θ))是对应环境参量b概率密度分布的随机多项式混沌基函数，θ为随机参量空间αs(r，z)是随机扰动环境参量投影到随机多项式混沌基函数得到的确定性系数。步骤四对无穷多阶多项式混沌展开作有限截取对不同分布的环境参量利用对应的多项式混沌基函数展开，展开的阶数是有限的；步骤五多项式混沌项数确定随机环境参量b由多个元素构成，b＝[b1b2...bn]，根据最小均方误差准则，n维随机变量b所必需的多项式项数为(n+S)！/n！S！；(2)声场不确实性表征①确实的环境参量b0驱动声传播模型产生确实的声场p0；②不确实环境参量空间id="icf0005"file="A2009101005950002C5.tif"wi="10"he="4"top="261"left="71"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>对应不确实声场空间id="icf0006"file="A2009101005950002C6.tif"wi="11"he="4"top="260"left="116"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>；③用概率密度函数p(p)表征不确实声场的分布；④id="icf0007"file="A2009101005950003C1.tif"wi="10"he="4"top="38"left="35"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>表征了声场不确实性的范围；同样，对不确实声场也作多项式混沌展开随机谱表征；步骤一不确实声场作多项式混沌展开<mathsid="math0002"num="0002"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mo>∞</mo></munderover><msub><mi>γ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>Λ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>ξ</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math></maths>步骤二对无穷多阶多项式混沌展开作有限截取根据最小均方误差准则，对无穷多阶多项式混沌展开作有限阶Q的截取<mathsid="math0003"num="0003"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi><mo>)</mo></mrow><mo>≅</mo><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>Q</mi><mo></mo></mrow></munderover><msub><mi>γ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>Λ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>ξ</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中p(r，z；θ)为预测的随机声压场；二、环境参量不确实性到声场不确实性的传递方法将环境参量与预测声场的多项式随机谱表征与确实性波方程相结合，得到嵌入不确实性的随机波方程，实现环境参量不确实性到声场不确实性的传递，具体步骤如下(1)确定性波方程到随机波方程步骤一基于简正波传播模型的确定性波方程<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msup><mo>▿</mo><mn>2</mn></msup><mi>p</mi><mo>+</mo><msup><mi>κ</mi><mn>2</mn></msup><mi>p</mi><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中id="icf0011"file="A2009101005950003C5.tif"wi="2"he="3"top="179"left="35"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>表示拉普拉氏算子，κ是波数，κ＝ω/c(z)，ω为角频率；步骤二随机扰动声速作二项式展开对随机扰动声速作二项式展开，并保留一阶项，则不确实声速的波数可近似为，<mathsid="math0005"num="0005"><math><![CDATA[<mrow><msup><mover><mi>κ</mi><mo>~</mo></mover><mn>2</mn></msup><mo>=</mo><mfrac><msup><mi>ω</mi><mn>2</mn></msup><mrow><msubsup><mi>c</mi><mn>0</mn><mn>2</mn></msubsup><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>δc</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>≈</mo><msup><mi>κ</mi><mn>2</mn></msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mn>2</mn><mi>δc</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi></mrow></mrow><mrow><msub><mi>c</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math></maths>步骤三确定性波方程推广到随机波方程将近似的波数代入到确定性波方程，将确定性波方程推广到随机波方程<mathsid="math0006"num="0006"><math><![CDATA[<mrow><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>p</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>p</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><msup><mi>κ</mi><mn>2</mn></msup><mi>p</mi><mo>-</mo><msup><mrow><mn>2</mn><mi>κ</mi></mrow><mn>2</mn></msup><mfrac><mrow><mi>δc</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>c</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mi>p</mi><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math></maths>将式(3)和式(4)代入(8)式，并对等式两边同时乘以Λ1(ξ(θ))，然后通过系总平均将各项投影到基空间，可得<mathsid="math0007"num="0007"><math><![CDATA[<mrow><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mi>γ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo><</mo><msub><mi>Λ</mi><mi>q</mi></msub><msub><mi>Λ</mi><mn>1</mn></msub><mo>></mo><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mi>γ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo><</mo><msub><mi>Λ</mi><mi>q</mi></msub><msub><mi>Λ</mi><mn>1</mn></msub><mo>></mo><mo>+</mo><msup><mi>κ</mi><mn>2</mn></msup><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>γ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo><</mo><msub><mi>Λ</mi><mi>q</mi></msub><msub><mi>Λ</mi><mn>1</mn></msub><mo>></mo></mrow>]]></math>id="icf0014"file="A2009101005950004C1.tif"wi="110"he="11"top="28"left="48"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>(9)<mathsid="math0008"num="0008"><math><![CDATA[<mrow><mo>-</mo><mfrac><msup><mrow><mn>2</mn><mi>κ</mi></mrow><mn>2</mn></msup><mrow><msub><mi>c</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo><</mo><msubsup><mi>Λ</mi><mn>1</mn><mn>2</mn></msubsup><mo>></mo></mrow></mfrac><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><munderover><mi>Σ</mi><mrow><mi>s</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>α</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>γ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo><</mo><msub><mi>Λ</mi><mi>q</mi></msub><msub><mi>Λ</mi><mi>s</mi></msub><msub><mi>Λ</mi><mn>1</mn></msub><mo>></mo><mo>=</mo><mn>0</mn></mrow>]]></math>id="icf0015"file="A2009101005950004C2.tif"wi="79"he="12"top="41"left="57"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>将混沌基函数的正交性应用于上述方程，得到一组声场不确实性系数的耦合差分方程，<mathsid="math0009"num="0009"><math><![CDATA[<mrow><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><mfrac><mrow><msubsup><mrow><mo>∂</mo><mi>γ</mi></mrow><mn>1</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mi>γ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><msup><mi>κ</mi><mn>2</mn></msup><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>γ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><mfrac><msup><mrow><mn>2</mn><mi>κ</mi></mrow><mn>2</mn></msup><mrow><msub><mi>c</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo><</mo><msubsup><mi>Λ</mi><mn>1</mn><mn>2</mn></msubsup><mo>></mo></mrow></mfrac><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><munderover><mi>Σ</mi><mrow><mi>s</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>α</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>γ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo><</mo><msub><mi>Λ</mi><mi>q</mi></msub><msub><mi>Λ</mi><mi>s</mi></msub><msub><mi>Λ</mi><mn>1</mn></msub><mo>></mo><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math></maths>以高斯分布随机声速进行下面公式的推导，δc(r，z，θ)＝α0He0+α1He1+…＝σξ，ξ是均值为0，方差为1的高斯随机变量。将δc代入(10)式给出第1个不确实系数的耦合方程<mathsid="math0010"num="0010"><math><![CDATA[<mrow><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mi>γ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mi>γ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><msup><mi>κ</mi><mn>2</mn></msup><msub><mi>γ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msup><mrow><mn>2</mn><mi>κ</mi></mrow><mn>2</mn></msup><mi>σ</mi></mrow><mrow><msub><mi>c</mi><mn>0</mn></msub><mo><</mo><msubsup><mi>He</mi><mn>1</mn><mn>2</mn></msubsup><mo>></mo></mrow></mfrac><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mo>∞</mo></munderover><msub><mi>γ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo><</mo><msub><mi>He</mi><mi>q</mi></msub><mi>ξ</mi><msub><mi>He</mi><mn>1</mn></msub><mo>></mo><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math></maths>利用Hermite基函数的回归特性和正交性，计算耦合项的系总平均，得到系数γ1的偏微分回归方程<mathsid="math0011"num="0011"><math><![CDATA[<mrow><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mi>γ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mi>γ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><msub><mi>Dγ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><mi>E</mi><msqrt><mn>1</mn></msqrt><msub><mi>γ</mi><mrow><mn>1</mn><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><mi>E</mi><msqrt><mn>1</mn><mo>+</mo><mn>1</mn></msqrt><msub><mi>γ</mi><mrow><mn>1</mn><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中D≡κ2，E≡2κ2σ/c0。用向量-矩阵形式，上述偏微分方程可化简为<mathsid="math0012"num="0012"><math><![CDATA[<mrow><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>γ</mi></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>r</mi></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mi>Aγ</mi><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中，A为环境参量不确实性到声场不确实性的耦合方程，A是对称的三角阵；(2)求解随机波方程得到不确实声场步骤一耦合矩阵A特征值分解<mathsid="math0013"num="0013"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>1</mn></mrow><mi>Q</mi></munderover><msub><mi>λ</mi><mi>q</mi></msub><msub><mi>g</mi><mi>q</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中G＝[{g1}，{g2}，…]为耦合矩阵A的特征向量，GG-1＝I，(13)式可以表示为<mathsid="math0014"num="0014"><math><![CDATA[<mrow><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msup><mi>G</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>γ</mi></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msup><mi>G</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>γ</mi></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><msup><mi>G</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>AG</mi><msup><mi>G</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>γ</mi><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math></maths>令<mathsid="math0015"num="0015"><math><![CDATA[<mrow><mover><mi>γ</mi><mo>^</mo></mover><mo>=</mo><msup><mi>G</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>γ</mi><mo>,</mo></mrow>]]></math>id="icf0022"file="A2009101005950004C9.tif"wi="17"he="4"top="246"left="23"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>上式可以简化为<mathsid="math0016"num="0016"><math><![CDATA[<mrow><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mover><mi>γ</mi><mo>^</mo></mover></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mover><mi>γ</mi><mo>^</mo></mover></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mi>Ω</mi><mover><mi>γ</mi><mo>^</mo></mover><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中，Ω＝G-1AG为对角阵，对角阵元素为特征值λ1，λ2，…，利用特征值和特征向量分解技术，得到第1个确定系数的波方程<mathsid="math0017"num="0017"><math><![CDATA[<mrow><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mover><mi>γ</mi><mo>^</mo></mover><mn>1</mn></msub></mrow><msup><mrow><mo>∂</mo><mi>r</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mover><mi>γ</mi><mo>^</mo></mover><mn>1</mn></msub></mrow><msup><mrow><mo>∂</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><msub><mi>λ</mi><mn>1</mn></msub><msub><mover><mi>γ</mi><mo>^</mo></mover><mn>1</mn></msub><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math></maths>步骤二确定系数求解等式(17)与标准波方程(6)有相同的形式，因而求解过程与标准的波方程求解过程是一样的。id="icf0025"file="A2009101005950005C2.tif"wi="3"he="5"top="70"left="178"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的解表示为<mathsid="math0018"num="0018"><math><![CDATA[<mrow><msub><mover><mi>γ</mi><mo>^</mo></mover><mi>u</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>i</mi><msqrt><mn>8</mn><mi>πr</mi></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mfrac><mi>π</mi><mn>4</mn></mfrac></mrow></msup><msub><mi>a</mi><mrow><mi>u</mi><mn>0</mn></mrow></msub><munderover><mi>Σ</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mfrac><mrow><msub><mi>ψ</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>s</mi></msub><mo>)</mo></mrow><msub><mi>ψ</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow><msqrt><msub><mi>k</mi><mi>rm</mi></msub></msqrt></mfrac><msup><mi>e</mi><mrow><msub><mi>ik</mi><mi>rm</mi></msub><mi>r</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中krm、ψm(z)分别是水平波数和模深度函数，由简正波传播模型求解得到，au0表示矩阵A的特征向量的逆的第一列中的第u个元素，在声速梯度不确实情况下，模深度函数没有受到不确实声速的影响，但声速梯度不确实导致了水平波数的不确实，<mathsid="math0019"num="0019"><math><![CDATA[<mrow><msub><mi>k</mi><mi>rm</mi></msub><mo>=</mo><msqrt><msub><mi>λ</mi><mn>1</mn></msub><mo>-</mo><msubsup><mi>k</mi><mi>zm</mi><mn>2</mn></msubsup></msqrt><mo>,</mo></mrow>]]></math>id="icf0027"file="A2009101005950005C4.tif"wi="29"he="7"top="123"left="94"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>其中kzm是深度波数，根据γ1与id="icf0028"file="A2009101005950005C5.tif"wi="2"he="4"top="125"left="179"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的关系式，系数γ1是id="icf0029"file="A2009101005950005C6.tif"wi="2"he="4"top="136"left="50"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的线性组合<mathsid="math0020"num="0020"><math><![CDATA[<mrow><msub><mi>γ</mi><mn>1</mn></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>u</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>g</mi><mi>lu</mi></msub><msub><mover><mi>γ</mi><mo>^</mo></mover><mi>u</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow>]]></math></maths>步骤三Q个多项式混沌基函数的确定性系数加权累加得到不确实性声场最后，在远场近似下的不确实性声场可表示为<mathsid="math0021"num="0021"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><mo>[</mo><munderover><mi>Σ</mi><mrow><mi>u</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>g</mi><mi>qu</mi></msub><msub><mover><mi>γ</mi><mo>^</mo></mover><mi>u</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>]</mo><msub><mi>Λ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>ξ</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math></maths><mathsid="math0022"num="0022"><math><![CDATA[<mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><mo>[</mo><munderover><mi>Σ</mi><mrow><mi>u</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>g</mi><mi>qu</mi></msub><mrow><mo>(</mo><mfrac><mi>i</mi><msqrt><mn>8</mn><mi>πr</mi></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mfrac><mi>π</mi><mn>4</mn></mfrac></mrow></msup><msub><mi>a</mi><mrow><mi>u</mi><mn>0</mn></mrow></msub><munderover><mi>Σ</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mfrac><mrow><msub><mi>ψ</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>s</mi></msub><mo>)</mo></mrow><msub><mi>ψ</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow><mroot><mrow><msub><mi>λ</mi><mi>u</mi></msub><mo>-</mo><msubsup><mi>k</mi><mi>zm</mi><mn>2</mn></msubsup></mrow><mn>4</mn></mroot></mfrac><msup><mi>e</mi><mrow><mi>i</mi><msqrt><msub><mi>λ</mi><mi>u</mi></msub><mo>-</mo><msubsup><mi>k</mi><mi>zm</mi><mn>2</mn></msubsup></msqrt><mi>r</mi></mrow></msup><mo>)</mo></mrow><mo>]</mo><msub><mi>Λ</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>ξ</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow>]]></math></maths>若SSP是确实的，则上式中只有γ0是不等于0的，收敛到声速确实时的声场简正模表达式；步骤四不确实性声场的均值和相关函数由多项式混沌展开的不确实声场是随机声场，其均值可表示为<mathsid="math0023"num="0023"><math><![CDATA[<mrow><mo><</mo><mi>p</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>θ</mi><mo>)</mo></mrow><mo>></mo><mo>=</mo><msub><mover><mi>γ</mi><mo>^</mo></mover><mn>0</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>21</mn><mo>)</mo></mrow></mrow>]]></math></maths>相关函数为<mathsid="math0024"num="0024"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>p</mi><mi>p</mi></mrow></msub><mrow><mo>(</mo><msub><mi>r</mi><mn>1</mn></msub><mo>,</mo><msub><mi>z</mi><mn>1</mn></msub><mo>;</mo><msub><mi>r</mi><mn>2</mn></msub><mo>,</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>q</mi><mo>=</mo><mn>1</mn></mrow><mi>Q</mi></munderover><msub><mover><mi>γ</mi><mo>^</mo></mover><mi>q</mi></msub><mo>*</mo><mrow><mo>(</mo><msub><mi>r</mi><mn>1</mn></msub><mo>,</mo><msub><mi>z</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mover><mi>γ</mi><mo>^</mo></mover><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>r</mi><mn>2</mn></msub><mo>,</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>22</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math></maths>2、根据权利要求1所述的海洋环境与声场不确实性的表征和传递的快速计算方法，其特征在于b是高斯分布，As(《(e))为Gauss-Hemiite多项式；b是均匀分布，As(g(e))为Legendre多项式；b是Y分布，As(《(e))为Lagueire多项式，b是卩分布，八8(《(0))为Jacobi多项式，b是泊松分布，As(《(e))为Charlier多项式。3、根据权利要求2所述的海洋环境与声场不确实性的表征和传递的快速计算方法，其特征在于对于高斯分布分布的b，只要一阶就够了，其中a。(r，z)等于0，a"r,z)-cr，式中cr是随机扰动环境参量的标准离差，具有类似的特性，其它分布的环境参量也可以利用投影方法得到有限阶(阶数为S)的确定性系数。全文摘要本发明涉及声学领域，具体涉及一种海洋环境与声场不确实性的表征和传递的快速计算方法，解决的问题是不确实环境参量和不确实声压场的随机谱表征，以及不确实性如何从环境参量传递到声压场。本发明对不确实环境参量和不确实声压场分别建立以已知量为中心的领域空间，并对不确实量作概率密度函数描述，根据概率密度函数作不确实量的多项式混沌展开随机谱表征，进一步将多项式混沌随机谱表征与确定波方程相结合，导出嵌入不确实性的随机波方程，获得不确实性系数的偏微分方程组，最后得到由不确实性系数加权的多项式混沌基函数线性叠加构成的不确实声压场。本发明有益的效果得到的不确实声压场与MonteCarlo得到的结果完全一致，其计算速度提高了10倍以上。文档编号G01S7/52GK101604019SQ20091010059公开日2009年12月16日申请日期2009年7月13日优先权日2009年7月13日发明者献祝,赵航芳,邹丽娜申请人:中国船舶重工集团公司第七一五研究所