ECG signals’complex network recognition for health and Myocardial Infarction based on multiplex visibility graph

doi:10.7498/aps.70.20211656

Abstract

FullText HTML

1. 引　言
2. 浅海全波动混响理论
2.1 浅海海底粗糙界面散射
2.2 浅海海底粗糙界面混响
3. 海底反射系数的三参数模型
4. 海底粗糙界面反向散射模型
4.1 角度特性
4.2 强度特性
5. 结　论

Cited By

Abstract

The Visibility Graph algorithm has been proved to be a simple and efficient method to transform time series into complex network and has been widely used in time series analysis because it can inherit the dynamic characteristics of original time series in topological structure. Now, visibility graph analysis of univariate time series has driven to mature, however, most complex systems in real world are multi-dimensional, it is difficult for univariate analysis to describe the global characteristics when applied to multi-dimensional series. In this paper, a novel method for multivariate time series analysis is proposed. For patients with myocardial infarction and healthy subjects, the 12-lead electrocardiogram signals of each individual are considered as a multivariate time series, which is transformed into a multiplex visibility graph through visibility graph algorithm and then mapped to fully connected complex network. Each node of the network corresponds to a lead and the inter-layer mutual information between visibility graph of two leads represents the weight of edges. Due to the fully connected network of different groups show identical topological structure, the dynamic characteristics of different individuals cannot be uniquely represented. Therefore, we reconstruct the fully connected network according to inter-layer mutual information, when the value of inter-layer mutual information is less than the threshold we set, edge corresponding to the inter-layer mutual information is deleted. We extract average weighted degree and average weighted clustering coefficient of reconstructed networks so as to recognize the 12-lead ECG signals of healthy subjects and myocardial infarction patients. Moreover, multiscale weighted distribution entropy is also introduced to analyze relationship between the length of original time series and final recognition results. Due to higher average weighted degree and average weighted clustering coefficient of healthy subjects, their reconstructed networks show a more regular structure, higher complexity and connectivity, and could be distinguished from patients with myocardial infarction, whose reconstructed networks are sparser. Experiment results show that the identification accuracy of both parameters, average weighted degree and average weighted clustering coefficient, reach 93.3%, which can distinguish the 12-lead ECG signals of healthy people and patients with myocardial infarction, and realize the automatic detection of myocardial infarction.

Keywords:

PACS:

03.67.Dd	(Quantum cryptography and communication security)
03.65.Ud	(Entanglement and quantum nonlocality)
03.67.Hk	(Quantum communication)

FullText HTML

1. 引　言

海底作为浅海混响的主要散射源, 浅海海底界面混响是很多学者探究的课题. 浅海混响根据其形成过程可以分为传播和散射两个过程. 入射声传播和散射声传播过程同属信道声传播问题, 研究较为完善. 而散射过程则比较复杂, 散射强度的测量也较难实现, 尤其是小掠射角的反向散射的测量, 所以针对散射的研究相对比较受限. 在混响研究伊始, 各国学者均采用经验的散射模型描述海底对声场的散射作用. 由于其形式较为简单, 运算速度较快, 在主动声纳预报系统中有很大的优势^[1-4], 至今仍被延续使用. 但是该类型的散射模型也存在着很大的局限性, 很难分析海底散射机理. 最初的经验散射模型是借鉴光学理论中提出的半无限自由空间中的Lambert散射定律给出的海底散射模型, 是海底散射强度与平面波掠射角之间的关系. 而在波导环境中, 由于频散效应的存在, 平面波散射理论不再适用于声场处理, 遂采用简正波理论, 避免了频散产生的多途现象. 然而简正模态的掠射角不同于平面波的掠射角, 因此对波导环境中的声散射问题不能够直接采用经验的Lambert散射模型. 近年来, 越来越多的学者开始从海底散射的物理机理着手建立散射模型. 国际上关于物理散射模型的建模方法包括有限元方法^[5]、Kirchhoff近似方法^[6]、微扰近似方法^[7-10]. 国内关于物理散射模型的研究最具代表性的是中国科学院声学研究所的高天赋和尚尔昌. Shang等^[11]在Gao^[12]和Tang^[13]的研究基础上提出了全波动混响理论, 并对该模型进行了一系列的发展. 2001年, Gao等^[14]提出了依据混响数据反演海底反向散射矩阵的方法. Wu等^[15,16]将海底反射模型引入到浅海混响模型中简化了混响衰减特性与海底沉积层参数之间的关系, 为海底地声参数的反演提供了一种新方法.

本文在尚尔昌和吴金荣的研究基础上, 对反向散射模型进行进一步分析. 通过引入海底反射系数的相移参数, 描述海底对声场的散射作用, 分析海底反向散射模型的角度特性和强度特性. 在实际应用允许的误差范围内, 进一步分析海底粗糙界面反向散射模型的角度特性和强度特性, 并与经验散射模型进行对比, 说明两者之间的异同.

2. 浅海全波动混响理论

在浅海波导环境中, 海底粗糙界面是海底混响的主要散射源. 在图1所示的浅海波导环境中, 平坦海底的平均海深为 ${H_0}$ , 任意水平位置 $r$ 的海底粗糙界面相对于平均海深的起伏高度为 $\eta (r)$ , 其量值远小于 ${H_0}$ , 均值为0, 方差为 $\sigma _\eta ^{\rm{2}}$ . 海底粗糙界面对声场的散射强度远小于入射声场的强度, 即满足弱散射的条件. 海底为半无限均匀介质, 其声速、密度和声吸收系数分别为 ${c_b},{\rho _b}\,{\text{和}}\,\alpha _b^{(w)}$ (声吸收系数的单位是 ${{{\rm{dB}}}/{\rm{\lambda }}}$ ). 海面为平坦的自由边界, 不考虑其对声场的散射作用.

$图 1 浅海波导环境\r\nFig. 1. Shallow water waveguide environment.$

图 1 浅海波导环境

Fig. 1. Shallow water waveguide environment.

下载: 全尺寸图片幻灯片

2.1 浅海海底粗糙界面散射

全波动混响理论的基础思想是将海底粗糙界面对声场的散射作用看作是“二次声源”向外辐射声能量, 从而将海底粗糙界面的散射作为“声源”向波导环境中辐射声能量^[11]. 该“声源”的声源强度与入射声强度成正比关系. 在水平均匀的波导环境中, 单位简谐点源入射声场满足波动方程为

${\nabla ^2}{u_{{\rm{w,b}}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}} + k_{{\rm{w,b}}}^{\rm{2}}{u_{{\rm{w,b}}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}} = 4{\rm{{\text π} \delta (}}{{{R}}_1} - {{{R}}_0}{\rm{)}},$

(1)

海底粗糙界面的“二次声源”满足

$\begin{split} & {\left. {\frac{{\partial {u_{\rm{w}}}({ R_1},{ R_0})}}{{\partial z}} - \frac{{\partial {u_{\rm{b}}}({ R_1},{ R_0})}}{{\partial z}}} \right|_{z = {H_0}}}\\ = & V({R_1}){G^{\rm{i}}}({ R_1},{ R_0}), \end{split}$

(2)

$\begin{split} & {\left. {{\rho _{\rm{w}}}{u_{\rm{w}}}({ R_1},{ R_0}) - {\rho _{\rm{b}}}{u_{\rm{b}}}({ R_1},{ R_0})} \right|_{z = {H_0}}}\\ = & p({ R_1}){G^{\rm{i}}}({ R_1},{ R_0}), \end{split}$

(3)

${G^{\rm{i}}}({{{R}}_1},{{{R}}_0}) = \sqrt {\frac{{2{\rm{\pi }}i}}{{{k_{\rm{w}}}r}}} \sum\limits_{m = 1}^M {{\phi _m}({z_0}){\phi _m}{\rm{(}}{H_0}{\rm{)}}{{\rm{e}}^{{\rm{i}}{k_m}r - {\beta _m}r}}} ,$

(4)

$V({{{R}}_1}) = \left(k_{\rm{w}}^{\rm{2}} - \frac{{k_{\rm{b}}^{\rm{2}}}}{\kappa }\right)\eta {\rm{(}}{r_1}{\rm{)}} + \left(1 - \frac{1}{\kappa }\right){\nabla _ \bot } \cdot {\rm{(}}\eta ({r_1}{\rm{)}}{\nabla _ \bot }{\rm{)}},$

(5)

$p{\rm{(}}{{{R}}_1}{\rm{)}} = {\rm{(}}{\rho _{\rm{w}}} - {\rho _{\rm{b}}}{\rm{)}}\eta {\rm{(}}{r_1}{\rm{)}}\frac{\partial }{{\partial z}},$

(6)

${k_{\rm{w}}} = \frac{{2{{\text π}}f}}{{{c_{\rm{H}}}}},\;\;\;\;\;\;{k_{\rm{b}}} = \frac{{2{{\text π}}f}}{{{c_{\rm{b}}}}},\;\;\;\;\;\;\kappa = \frac{{{\rho _{\rm{b}}}}}{{{\rho _{\rm{w}}}}},$

(7)

其中, ${{{R}}_0} = {\rm{(}}{r_{\rm{s}}},{z_0}{\rm{)}}$ , ${{{R}}_1} = ({r_1},{H_0})$ 分别是声源和散射点的空间位置, 两者之间的水平距离差为 $r = |{r_{\rm{s}}} - {r_1}|$ ; $u_{\rm{w}}, u_{\rm{b}}$ 分别是水介质和沉积层介质中“二次源”的辐射声场(下脚标w和b分别代表水介质和沉积层介质); ${k_{{\rm{w,b}}}}$ 是水介质和沉积层介质中的波数, 与声源频率 $f$ 和介质声速有关; ${c_{\rm{H}}} =$ ${c_{\rm{w}}}{\rm{(}}{H_0}{\rm{)}}$ 是海深处的声速, 沉积层介质和水介质的密度比用 $\kappa$ 表示; ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ 是波导环境中初级声场的格林函数. 根据简正波理论的思想, 在浅海波导环境中, ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ 可以表示为简正模态 ${\phi _m}({H_0})$ (实际上是 ${\phi _m}(z)$ , $z$ 是接收深度, 在这里的接收深度为海深 ${H_0}$ )叠加的形式, 如(4)式所示, k_m和 ${\beta _m}$ 分别是简正模态本征值的实部和虚部. $V({{{R}}_1})$ 和 $p({{{R}}_1})$ 则分别是质点振速算子和压力算子, 作用到初级声场的格林函数 ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ , 得到海底粗糙界面的“二次速度源”的声源强度 $V{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ 和“二次压力源”的声源强度 $p{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ . 在瑞利参数较小的情况下, Bass给出了(5)式和(6)式的近似表达式. ${\nabla _ \bot }$ 为水平梯度算子,

${\nabla _ \bot } = \frac{\partial }{{\partial r}}.$

(8)

根据格林定理, 水平均匀波导环境中, “二次速度源”在接收点 ${{R}} = {\rm{(}}{r_{\rm{r}}},z{\rm{)}}$ 处的辐射声场可以通过格林函数 $G_{{\rm{w,b}}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$ 给出形式解. 辐射声场的格林函数满足

$\begin{aligned} & {\nabla ^2}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} + k_{\rm{w}}^{\rm{2}}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = 0\;\;\;z < {H_0}\\ & {\nabla ^2}G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} + k_{\rm{b}}^2G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = 0\;\;\;z > {H_0} \end{aligned},$

(9)

同时满足边界条件

$\begin{gathered} {\left. {\frac{{\partial G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}}}{{\partial z}} - \frac{{\partial G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}}}{{\partial z}}} \right|_{z = {H_0}}} = {\rm{\delta (}}{{R}} - {{{R}}_1}{\rm{)}} \\ {\left. {\rho G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} - {\rho _{\rm{b}}}G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}} \right|_{z = {H_0}}} = 0 \\ \end{gathered} .$

(10)

在接收点 “二次速度源”的辐射声场满足形式解,

$u_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = \int {G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}V{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)d}}{{{R}}_1}} .$

(11)

“二次压力源”的辐射声场同样有(11)式的形式解,

$u_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = \int {G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}p{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)d}}{{{R}}_1}} ,$

(12)

其中, $G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$ 是“二次压力源”的辐射声场的格林函数. 根据边界处法向振速与压力之间的关系, 不难得到 $G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$ 与 $G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$ 满足

$G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = - \frac{1}{{{\rho _{\rm{w}}}}}\frac{\partial }{{\partial z}}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}.$

(13)

叠加“二次速度源”和“二次压力源”的辐射声场, 得到海底粗糙界面向波导环境中辐射的总声场

${u_{\rm{w}}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = u_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} + u_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}}.$

(14)

将(11)—(13)式代入(14)式, 得到海底粗糙界面散射声场的一般形式解

$\begin{split} & {u_{\rm{w}}}( R,{ R_1},{ R_0})\\ = & \int [G_{\rm{w}}^{\rm{V}}( R,{ R_1})V({ R_1})\\ & - \frac{1}{{{\rho _{\rm{w}}}}}\frac{\partial }{{\partial z}}G_{\rm{w}}^{\rm{V}}( R,{ R_1})p({ R_1})]{G^{\rm{i}}}({ R_1},{ R_0}){\rm{d}}{ R_1}. \end{split}$

(15)

实际上, 入射声场的格林函数、散射场的格林函数以及初级声场的格林函数有一致的形式解, 忽略他们之间的差异,

${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{R}}{\rm{)}} = G_{\rm{w}}^{\rm{V}}{\rm{(}}{{{R}}_1},{{R}}{\rm{)}} = G{\rm{(}}{{{R}}_1},{{R}}{)},$

(16)

所以对(15)式进行简化得到

$\begin{split} & {u_{\rm{w}}}( R,{ R_1},{ R_0})\\ = & \int_S {[G( R,{ R_1})V({ R_1})} \\ & - \frac{1}{{{\rho _{\rm{w}}}}}\frac{{\partial G( R,{ R_1})}}{{\partial z}}p({ R_1})]G({ R_1},{ R_0}){\rm{d}}{ R_1}. \end{split}$

(17)

对于本地混响, 声源和接收均在同一水平位置. 为了方便计算, 通常令其位于过源点的垂直轴线上, 即 ${r_{\rm{s}}} = {r_{\rm{r}}} = 0$ . 将(4)—(6)式代入(17)式, 并通过分部积分法得到简谐点源的海底粗糙界面的反向散射声场, 将其描述为简正模态的叠加形式^[11],

$\begin{split} & {u_{\rm{w}}}(r,z)\\ = & \displaystyle\frac{{2\pi i}}{{{k_{\rm{w}}}r}}\sum\limits_m^M {\sum\limits_n^M {{\phi _m}({z_0}){\phi _m}({H_0})} } \cdot C_{mn}^\eta K_{mn}^\eta ({k_m},{k_n})\\ & \times {\phi _n}({H_0}){\phi _n}(z){{\rm{e}}^{ - ({\beta _m} + {\beta _n})r}} \end{split}$

(18)

其中

$C_{mn}^\eta = k_{\rm{w}}^{\rm{2}} - \frac{{k_{\rm{b}}^2}}{\kappa } + {\rm{(}}1 - \frac{1}{\kappa }{\rm{)}}{k_m}{k_n} + \frac{{1 - \kappa }}{{{\kappa ^2}}}{\gamma _m}{\gamma _n},$

(19)

$K_{mn}^\eta {\rm{(}}{k_m},k{}_n{\rm{)}} = \int {\eta {\rm{(}}r'{\rm{)}}{{\rm{e}}^{{\rm{i(}}{k_m} + {k_n}{\rm{)}}r'}}{\rm{d}}r'} .$

(20)

2.2 浅海海底粗糙界面混响

实际上(18)式是平坦波导环境中海底粗糙界面散射场的稳定解形式. 对于声源脉冲信号s(t), 其功率谱 $F\left( \omega \right)$ 满足Fourie变换

$F{\rm{(}}\omega {\rm{)}} = \int_{ - \infty }^{ + \infty } {s{\rm{(}}t{\rm{)}}{{\rm{e}}^{ - {\rm{i}}\omega t}}{\rm{d}}t},\tag{21a}$

(21a)

$s{\rm{(}}t{\rm{)}} = \frac{1}{{2{\rm{\pi }}}}\int_{ - \infty }^{ + \infty } {F{\rm{(}}\omega {\rm{)}}{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}\omega }.\tag{21b}$

(21b)

根据Fourie变换的性质, 得到入射信号为s(t)的海底散射声场为

$u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{(}}r,z{\rm{)}} = \frac{1}{{2{\rm{\pi }}}}\int_{ - \infty }^{ + \infty } {F{\rm{(}}\omega {\rm{)}}{u_{\rm{w}}}{\rm{(}}r,z{\rm{)}}{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}\omega } .$

(22)

将声源信号的功率谱函数 $F\left( \omega \right)$ 和(18)式代入(22)式, 同时取水平波束在中心频率处的两阶泰勒级数展开, 进而得到脉宽为 $\tau$ 声源信号s(t)的反向散射声场,

$\begin{split} & u_{\rm{w}}^{s(t)}(r,z)\\ = & \frac{{2\pi i}}{{{k_{\rm{w}}}r}}s(t - {t_0})\sum\limits_m^M {\sum\limits_n^M {{\phi _m}({z_0}){\phi _m}({H_0}) } } \\ & \times C_{mn}^\eta K_{mn}^\eta \cdot {\phi _n}({H_0}){\phi _n}(z){{\rm{e}}^{ - ({\beta _m} + {\beta _n})r}} \end{split}$

(23)

其中, t₀是信号的传播时间, 与散射环的水平距离满足 ${t_0} \approx {{2r}/{{c_{\rm{H}}}}}$ .

混响是能够同一时刻被接收、来自各方向散射回波的叠加. 对于海底粗糙界面混响而言, 同一时刻接收到的散射回波来自于宽度为 $\Delta r = {{{c_{\rm{H}}}\tau }/2}$ 的散射圆环, 其内径为 ${r_{{\rm{in}}}} = r - \displaystyle\frac{{{c_{\rm{H}}}\tau }}{4}$ , 外径为 ${r_{{\rm{ex}}}} = r +$ $\displaystyle\frac{{{c_{\rm{H}}}\tau }}{4}$ . 对(23)式计算强度, 并在散射环面积上进行积分得到海底粗糙界面混响平均强度

${I_\eta }{\rm{(}}r,z{\rm{)}} = \int_A {\langle u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^{*}\rangle } {\rm{d}}A,$

(24)

其中, ${\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^*$ 是 $u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}$ 的复共轭.

假设散射环面积内的海底界面起伏高度满足各向同性的原则, 对(24)式进一步简化为

${I_\eta }{\rm{(}}r,z{\rm{)}} = 2{\rm{\pi }}r\int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\langle u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^{*}\rangle {\rm{d}}r} .$

(25)

将(23)式代入(25)式, 并假设脉宽远小于混响时间t, 将声源强度近似为

${E_0} = {s^2}\left( {t - {t_0}} \right) \cdot \tau .$

(26)

在这里, 只考虑其非相干特性, 得到短脉冲的非相干混响平均强度

$\begin{split} {I_\eta }(r,z) = & {E_0} \cdot {\left(\frac{{2{\text{π}}}}{{{k_{\rm{w}}}r}}\right)^2}\\ & \times {\text{π}} rc\sum\limits_m^M {\sum\limits_n^M {\phi _m^2({z_0})\phi _m^2({H_0}) \cdot S_{mn}^\eta } } \\ & \times \phi _n^2({H_0})\phi _n^2(z){{\rm{e}}^{ - 2({\beta _m} + {\beta _n})r}}, \end{split}$

(27)

$S_{mn}^\eta = {\left[ {C_{mn}^\eta } \right]^2}{S^\eta }\left( {2{k_0}} \right),$

(28)

$\begin{split} & {S^\eta }(2{k_0})\\ = & \int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\langle \eta (r')\eta (r'')\rangle {{\rm{e}}^{{\rm{i}}({k_m} + {k_n})r' - {\rm{i}}({k_m} + {k_n})r''}}} } \\ & \times {\rm{d}}r''{\rm{d}}r'\; = \;\sigma _\eta ^{\rm{2}}{P^\eta }(2{k_0}), \end{split}$

(29)

其中, ${S^\eta }(2{k_0})$ 是海底界面的粗糙度谱, 与海底沉积层的地声参数无关, 只是声源频率的函数; $S_{mn}^\eta$ 是海底粗糙界面引起的入射模态和散射模态之间的耦合系数, 反映海底粗糙界面的反向散射能力, 由海底界面粗糙度谱和地声参数决定. 定义

$\Theta _{mn}^\eta = \phi _m^2({H_0})S_{mn}^\eta \phi _m^2({H_0})$

(30)

是区别于经验散射模型的物理散射模型. 它结合了波导环境特性和简正波理论的思想, 建立了受格林函数约束的散射模型, 明确了海底地声参数以及海底粗糙界面与海底反向散射的定量关系. 因此, 地声参数的准确度直接影响到海底反向散射函数的精确度. 通常情况下, 很难直接获取大面积的地声数据, 通过反演得到的地声参数又受地声模型的影响较大, 而且反演的未知参数较多, 导致结果存在很大的不确定性. 从反射系数的角度思考, 海底对声场的影响主要体现在海底反射系数. 海底反射系数的三参数模型与海底地声模型无关, 且参数较少, 利用海底反射系数代替地声参数作为输入参数, 能很大程度上简化海底粗糙界面的反向散射模型.

3. 海底反射系数的三参数模型

尚尔昌在1979年提出了海底反射系数的三参数模型^[17], 随后给出了海底小掠射角的反射系数幅值和相移与掠射角之间的量化关系, 该参数与地声模型无关. 海底反射系数的幅值用Q参数表示

$\ln \left| {V\left( \theta \right)} \right| = - Q\theta ,\;\;\theta \to 0,$

(31)

海底反射系数的相移用P参数表示

$\arg \left[ {V\left( \theta \right)} \right] = -{\text{π}} + P/\theta ,$

(32)

$\theta$ 是海底界面处入射声波的掠射角; Q参数反映海底声吸收效果, 体现声能量的衰减; P参数则描述各阶简正模态在海深处的能量, 即控制各阶简正模态的水平波数. 海底反射系数的相移参数P与地声参数之间的转换关系满足^[18,19]

${P^{{\rm{Low}}}} = \frac{{2\kappa }}{{\sqrt {1 - {\upsilon ^{ - 2}}} }},\;\;{\rm{ }}\theta < \theta ^{\rm{*}},$

(33)

${{{P^{{\rm{Crit}}}} = {\text{π}}}/{\theta ^{*}}},\;\;{\rm{ }}\theta \to \theta ^{\rm{*}},$

(34)

其中, $\upsilon = {{{c_{\rm{b}}}}/{{c_{\rm{H}}}}}$ 表示沉积层介质和水介质的声速比; $\kappa$ 是密度比, 如(7)式所示; $\theta ^{*}$ 是临界掠射角.

(33)式和(34)式分别给出了不同海底掠射角时海底反射系数的相移参数的不同形式. ${P^{{\rm{Low}}}}$ 是甚小掠射角处的海底反射系数相移参数, 是海底沉积层和水介质的声速比和密度比的函数; ${P^{Crit}}$ 是临界角附近海底反射系数的相移参数, 是临界角的函数.

海底反射系数的幅值参数Q与地声参数之间的转换关系满足^[18,19]

$Q = 0.036 \cdot \frac{{{\upsilon ^{ - 1}}\kappa }}{{{{\left[ {1 - {\upsilon ^{ - 2}}} \right]}^{3/2}}}} \cdot \alpha _b^{\left( w \right)},$

(35)

其中, $\alpha _b^{\left( w \right)}$ 是海底沉积层的声吸收系数, 单位是 $dB/\lambda$ .

根据Snell折射定律, $\upsilon$ 与临界掠射角之间满足

${\rm{cos(}}\theta ^{*}{\rm{)}} = {{{c_{\rm{H}}}}/{{c_{\rm{b}}}}} = {\upsilon ^{ - 1}}.$

(36)

将(36)式代入(34)式, 临界角附近的海底反射系数的相移参数 ${P^{Crit}}$ 与声速之间满足

$\upsilon = {\rm{co}}{{\rm{s}}^{ - 1}}{\rm{(}}\frac{{\rm{\pi }}}{{{P^{{\rm{Crit}}}}}}{\rm{)}}.$

(37)

联合(37)式和(33)式得到甚小掠射角的海底反射系数的相移参数 ${P^{Low}}$ 与密度之间满足关系式

$\kappa = \frac{{{P^{{\rm{Low}}}}}}{2}\sin \left(\frac{{\text{π}}}{{{P^{{\rm{Crit}}}}}}\right).$

(38)

(37)式和(38)式分别是海底反射系数的参数与等效的均匀半无限海底介质的地声参数之间的转换关系.

4. 海底粗糙界面反向散射模型

常用的经验散射模型如(39)式所示:

${\Theta _E} = {\mu _E}{\sin ^l}\theta {\sin ^k}\varphi ,$

(39)

其中, ${\mu _{\rm{E}}}$ 是经验散射系数, 反映反向散射强度; $\theta$ 是入射掠射角; $\varphi$ 是散射掠射角. 当 $l = k = 1$ , 是经验的Lambert散射模型. 该类型的散射模型是可分离的经验海底反向散射模型, 不可分离的经验散射模型(参见附录)也被用来描述海底反向散射.

(30)式的海底反向散射模型是根据波导环境建立的散射模型, 明确了入射简正模态、散射简正模态以及海底地声参数对反向散射的强度特性和角度特性的影响. 浅海波导环境中简正模态是受格林函数严格约束的稳定解, 与海底沉积层的介质参数和水文环境有关. 在掠射角小于临界角的条件下, 忽略各阶简正模态的水平波数(掠射角)之间的差异, 即 ${k_m} \approx {k_n} \approx {k_{\rm{w}}}$ , 则(19)式可以近似描述为与模态无关的量.

$\begin{split} C_{mn}^\eta = & \left[1 \!- \!\frac{{k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^2\kappa }} \!+\! \left(1 - \frac{1}{\kappa }\right)\frac{{{k_m}{k_n}}}{{k_{\rm{w}}^2}} + \frac{{1 - \kappa }}{{{\kappa ^2}}}\frac{{{\gamma _m}{\gamma _n}}}{{k_{\rm{w}}^2}}\right]k_{\rm{w}}^2 \\ & \approx \;\left[1 - \frac{{k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^{\rm{2}}\kappa }} + \left(1 - \frac{1}{\kappa }\right) + \frac{{1 - \kappa }}{{{\kappa ^2}}}\frac{{k_{\rm{w}}^2 - k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^{\rm{2}}}}\right]k_{\rm{w}}^{\rm{2}} \\ & = \;\left[2 - 2{\kappa ^{ - 1}} + {\kappa ^{ - 2}} - {\left(\kappa \upsilon \right)^{ - 2}}\right]k_{\rm{w}}^2. \end{split}$

(40)

将(37)式和(38)式代入(40)式, 用海底反射系数的相移参数代替地声参数, 得到

$C_{mn}^\eta = k_{\rm{w}}^2 \cdot \varsigma {\rm{(}}P{\rm{)}},$

(41)

其中

$\varsigma {\rm{(}}P{\rm{)}} \approx 2 - \frac{4}{{P{\rm{sin(}}{{\text{π}}/P}{\rm{)}}}} + \frac{4}{{{P^2}}}.$

(42)

在这里忽略 ${P^{Low}}$ 和 ${P^{Crit}}$ 之间的差异, 文献[18]中也说明了这种近似的合理性.

同时, 海底反射系数的相移参数 $P$ 决定了各阶简正模态在海底粗糙界面平均深度处的能量, 即

${\phi _m}({H_0}) = \sqrt {\frac{2}{{{H_0}}}} \sin \left(\frac{{P{\theta _m}}}{2}\right).$

(43)

将(41)式和(43)式代入(30)式得到远距离(小掠射角)条件下的海底反向散射核函数,

${\Theta _P} = \frac{4}{{H_0^2}}S(2{k_0}){[k_w^2 \cdot \varsigma (P)]^2} \cdot {\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right),$

(44)

其中 $\theta = {\theta _m}$ 为入射掠射角, $\varphi = {\theta _n}$ 为散射掠射角.

采用海底反射系数的相移参数 $P$ 描述海底反射系数, 明确了海底对声散射的物理机理, 不同于经验反向散射模型.

4.1 角度特性

在远距离条件下, (44)式表明, 海底反向散射的角度特性受海底反射系数的相移参数 $P$ 的影响, 不同于经验散射模型. 在临界角附近,

$\mathop {\lim }\limits_{\theta \to \theta *} \frac{{P\theta }}{2} = \frac{{\text{π}}}{2}.$

(45)

以 $\displaystyle\frac{{P\theta }}{2}$ 为宗量, 不能对 $\sin \left(\displaystyle\frac{{P\theta }}{2}\right)$ 进行小角度近似, 所以在临界角附近, (44)式所描述的反向散射模型随掠射角并不满足线性变化的关系.

以 ${{P\theta }/2} = 1$ 为分界点, 当 ${{P\theta }/2} < 1$ 时, 为甚小掠射角范围, 小角度近似引起的海底反向散射强度小于3 dB, 在实际应用中可忽略不计; 当 ${{P\theta }/2} > 1$ 时, 为临界角附近, 小角度近似使海底反向散射强度出现大于3 dB的误差, 实际应用中不可忽略. 所以, 在随后的分析中, 对反向散射模型的强度特性和角度特性进行分段处理.

当 ${{P\theta }/2} > 1$ , 即 ${2/P} < \theta < \theta *$ 时, (44)式不能进行小角度近似, 其随角度的变化关系受P参数加权.

$\begin{split}\Theta _P^{Crit} = & \mu _P^{Crit}{\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right) \\ & \propto {\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right),\end{split}$

(46)

其中

$\mu _P^{Crit} = \frac{4}{{H_0^2}}S(2{k_0}){[k_0^2 \cdot \varsigma (P)]^2}.$

(47)

当 ${{P\theta }/2} < 1$ , 即 $\theta < {2/P}$ 时, 对(44)式的角度项进行小角度近似,

$\sin \left(\frac{{P\theta }}{2}\right) \approx \frac{{P\theta }}{2}.$

(48)

将(48)式代入(44)式得到

$\Theta _P^{Low} = \mu _P^{Low} \cdot {\theta ^2}{\varphi ^2} \propto {\theta ^2}{\varphi ^2},$

(49)

其中

$\mu _P^{Low} = \frac{{{P^4}}}{{4H_0^2}}S(2{k_0}){[k_0^2 \cdot \varsigma (P)]^2}.$

(50)

对(39)式的经验散射模型进行小角度近似得到

${\Theta _E} \approx {\mu _E}{\theta ^l}{\varphi ^k} \propto {\theta ^l}{\varphi ^k}.$

(51)

当 $l = k = 2$ 时, 比较(51)式和(49)式具有相同的角度关系. 所以当掠射角 $\theta < {2/P}$ 时, 基于物理散射机理的反向散射模型与 $l = k = 2$ 的可分离经验散射模型具有一致的角度特性.

4.2 强度特性

通过4.1节的分析, 掠射角不同时, 海底反向散射的角度特性不同, 同时也反映了海底反向散射系数的差异.

(47)式指出, 在临界角附近, 海底反向散射系数由海底反射系数的相移参数P、海底界面粗糙度谱, 海深以及声源信号的频率决定. 忽略临界角附近以及甚小掠射角的P参数之间的差异^[18], 则(42)式中等式右边第三项相比前两项为小量, 可以忽略不计. 第二项中采用小角度近似,

${\rm{sin}}\frac{{\text{π}}}{P} \approx \frac{{\text{π}}}{P},$

(52)

使得(42)式中的第二项近似与P参数无关, 所以 $\varsigma (P) \approx 0.73$ . 图2也表明了 $\varsigma (P)$ 随P参数的变化不是很明显. 所以, 在临界角附近的海底反向散射系数近似与海底沉积层介质无关, 只是海底界面粗糙度谱、海深以及声源频率的函数.

$\mu _P^{Crit} \approx \frac{{2.11}}{{H_0^2}}S(2{k_0})k_0^4.$

(53)

$图 2 $\zeta $随P参数的变化\r\nFig. 2. $\zeta $ varied with P parameters near critical angle.$

图 2

$\zeta$

随P参数的变化

Fig. 2.

$\zeta$

varied with P parameters near critical angle.

下载: 全尺寸图片幻灯片

对于甚小掠射角的海底反向散射系数满足(50)式. 结合 $\varsigma (P)$ 近似为常数的结果, (50)式可以近似为与 ${P^4}$ 成正比, 即

$\mu _P^{Low} \approx {0.13}\frac{{{P^4}}}{{H_0^2}}S(2{k_0})k_0^4 \propto {P^4}.$

(54)

图3以第一类海底为例说明了这种近似结果的合理性. 图中给出了甚小掠射角的海底反向散射系数随 $P$ 参数的变化曲线(仿真参数: 海深为50 m, 声源频率为600 Hz, 海底粗糙界面的标准差为 $0.1\;{\rm{m}}$ , 相关尺度为 ${\rm{10}}\;{\rm{m}}$ , 计算得到Goff-Jordan谱^[20]为－32.8 dB)以及对其进行 $P$ 参数的四次方拟合结果, 两者的变化趋势基本一致. 所以, 在甚小掠射角的条件下, 海底反向散射系数与海底反射系数的相移参数 $P$ 密切相关, 即与海底沉积层的声速比和密度比有关, 与声吸收系数无关.

$图 3 海底甚小掠射角反向散射系数对P参数的依赖\r\nFig. 3. Relationship between P and bottom backscattering coefficient at very low grazing angle.$

图 3 海底甚小掠射角反向散射系数对P参数的依赖

Fig. 3. Relationship between P and bottom backscattering coefficient at very low grazing angle.

下载: 全尺寸图片幻灯片

另外, 从混响平均强度的角度分析, 同样说明了这种近似的合理性. 海底界面粗糙度谱函数采用图3的仿真参数, 海底是第一类均匀半无限介质, 声速为 $18{\rm{36}}\;{\rm{m/s}}$ , 密度为 ${\rm{2}}{\rm{.03}}\;{\rm{g/c}}{{\rm{m}}^{\rm{3}}}$ , 声吸收系数为 $0.88\;{\rm{dB/}}\lambda$ , 海深为50 m, 水介质为1500 m/s的等声速剖面的均匀水体. 声源频率为600 Hz, 发射深度为10 m, 接收深度为30 m. 计算水平散射距离从5 km到50 km的混响平均强度衰减曲线. 图4中实线是根据(30)式仿真的非近似混响平均强度衰减曲线; “○”是根据(46)式和(53)式仿真的临界角的近似结果, 与非近似的混响平均强度衰减曲线在13 km以前完全重合; 随着水平距离的增加, 两者相差逐渐增大; “ $\bullet$ ”是根据(49)式和(54)式仿真的甚小掠射角的近似结果, 在20 km以前, 与非近似结果之间的差异随水平距离的增大而减小, 在20 km以后两者基本吻合. 在散射距离较近时, 掠射角相对较大, 临界角附近近似的海底反向散射模型与真实结果更接近; 在散射距离较远时, 掠射角相对较小, 甚小掠射角的近似结果与真实结果更接近.

$图 4 海底反向散射模型的比较\r\nFig. 4. Compare with bottom backscattering model with different approximate.$

图 4 海底反向散射模型的比较

Fig. 4. Compare with bottom backscattering model with different approximate.

下载: 全尺寸图片幻灯片

5. 结　论

在浅海混响平均强度模型中, 经验散射模型在混响特性分析中存在明显的局限性, 而现有的物理散射模型受地声模型的影响较大. 本文以全波动混响理论的物理散射模型为基础, 结合海底反射系数的三参数模型, 将海底反射系数的相移参数等效代替地声参数, 描述海底对声场的散射作用, 简化了海底反向散射模型. 通过理论分析, 明确了海底掠射角以 ${2/P}$ 为分界点, 海底反向散射的角度特性和强度特性对海底反射系数的相移参数存在不同的依赖关系. 在掠射角满足 $\theta < {2/P}$ 时, 海底反向散射的角度特性近似描述为与可分离经验散射模型的角度特性一致, 与 $P$ 参数无关; 而其散射系数则近似描述为与P⁴线性增强. 当海底掠射角满足 ${2/P} < \theta < \theta ^{*}$ 时, 海底反向散射的角度特性是受 $P$ 参数加权的, 即受海底地声参数的影响; 而其散射系数则近似为与 $P$ 参数无关. 所以不同掠射角范围, 海底对反向散射声场的强度特性和角度特性的贡献不同. 掠射角较大时海底对反向散射声场的影响主要体现在其角度特性; 掠射角非常小时, 海底的影响主要体现在其强度特性.

附录A　经验反向散射模型

经验反向散射模型有多种形式, 主要体现在其角度特性之间的差异. 最常用的海底散射模型是可分离的散射模型,

$\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\theta {\rm{si}}{{\rm{n}}^k}\varphi, \tag{A1}$

(A1)

其中, ${\mu _{\rm{E}}}$ 是海底反向散射系数; $l,k$ 取不同的值对应不同的散射模型. 常见的Lambert散射模型则是(A1)式中 $l = k = 1$ 时的散射模型. 该模型最初是Mackenzie^[1]在处理深海海底混响模型时由光学的散射原理引入, 给出经验的海底反向散射模型, 随后根据实验数据与理论结果对比得到 ${\mu _{\rm{E}}}$ 的分贝值大约在－27 dB左右, 也就是说 ${\mu _{\rm{E}}} \approx {10^{ - 2.7}}$ . 当 $l = 1,k = 0$ 表示散射源强正比于 $\mu {\rm{sin}}\theta$ 的均匀散射; 当 $l = 0,k = 0$ 表示与角度无关的均匀散射; 同时也存在 $l = k = 2$ 的反向散射模型.

在海底界面大尺度不均匀波导环境中, 小掠射角和基尔霍夫近似的条件下的海底散射模型可以表示为

$\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\left(\frac{{\theta + \varphi }}{2}\right) \;\;l = 0,\frac{1}{2},1,\frac{3}{2},....\tag{A2}$

(A2)

在海底界面小尺度不均匀波导环境中, 小掠射角和基尔霍夫近似的条件下的海底散射模型可以表示为

$\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\left[{\rm cos}^{- 1}\left(\frac{{{\rm{cos}}\theta + {\rm{cos}}\varphi }}{2}\right)\right]\;\;l = 0,\frac{1}{2},1,\frac{3}{2},....\tag{A3}$

(A3)

另外一种不可分离的散射模型, 同样被用来描述海底反向散射过程

$\Theta = {\mu _{\rm{E}}}\frac{{{\rm{sin}}\theta {\rm{sin}}\varphi }}{{{\rm{sin}}\theta + {\rm{sin}}\varphi }}.\tag{A4}$

(A4)

通常无论是可分离的经验反向散射强度还是不可分离的经验反向散射强度, 它们的共同点在于海底反向散射强度特性与角度特性相互分离, 相互独立.

Cited By

期刊类型引用(2)

1.	赵宁，江英华，周贤韬. 基于单光子的高效量子安全直接通信方案. 物理学报. 2022(15): 36-41 . 百度学术
2.	张晨涛，石小涛，朱文新，朱金龙，郝向英，金锐博. 利用域排列算法设计铌酸锂晶体实现3μm中红外波段频域纯态单光子源. 物理学报. 2022(20): 81-87 . 百度学术

其他类型引用(1)

[1]	Ma Zhi-Yi, Yang Xiao-Dong, He Ai-Jun, Ma Lu, Wang Jun. Complex network recognition of electrocardiograph signals in health and myocardial infarction patients based on multiplex visibility graph. Acta Physica Sinica, doi: 10.7498/aps.71.20211656
[2]	Ma Jin-Long, Zhang Jun-Feng, Zhang Dong-Wen, Zhang Hong-Bin. Quantifying complex network traffic capacity based on communicability sequence entropy. Acta Physica Sinica, doi: 10.7498/aps.70.20201300
[3]	Wang Li-Na, Cheng Yuan-Yuan, Zang Chen-Rui. A symbolized time series network based on seasonal-trend-loess method. Acta Physica Sinica, doi: 10.7498/aps.68.20190794
[4]	Chen Dan, Shi Dan-Dan, Pan Gui-Jun. Correlation between the electrical transport performance and the communicability sequence entropy in complex networks. Acta Physica Sinica, doi: 10.7498/aps.68.20190230
[5]	Xing Xue, Yu De-Xin, Tian Xiu-Juan, Wang Shi-Guang. Analysis of multi-state traffic flow time series properties using visibility graph. Acta Physica Sinica, doi: 10.7498/aps.66.230501
[6]	Wu Xi-Ping, Yang Hong-Yu, Han Song-Chen. Analysis on network properties of multivariate mixed air traffic management technical support system based on complex network theory. Acta Physica Sinica, doi: 10.7498/aps.65.140203
[7]	Hou Feng-Zhen, Dai Jia-Fei, Liu Xin-Feng, Huang Xiao-Lin. Phase synchrony in the cerebral infarction electroencephalogram based on the degree of network-links. Acta Physica Sinica, doi: 10.7498/aps.63.040506
[8]	Li Yu-Shan, Lü Ling, Liu Ye, Liu Shuo, Yan Bing-Bing, Chang Huan, Zhou Jia-Nan. Spatiotemporal chaos synchronization of complex networks by Backstepping design. Acta Physica Sinica, doi: 10.7498/aps.62.020513
[9]	Gao Zhong-Ke, Hu Li-Dan, Zhou Ting-Ting, Jin Ning-De. Limited penetrable visibility graph from two-phase flow for investigating flow pattern dynamics. Acta Physica Sinica, doi: 10.7498/aps.62.110507
[10]	Xing Chang-Ming, Liu Fang-Ai, Xu Ru-Zhi. Exact solution for mean trapping time of random walk on a scale-free Koch network. Acta Physica Sinica, doi: 10.7498/aps.61.200503
[11]	Gao Xiang-Yun, An Hai-Zhong, Fang Wei. Research on fluctuation of bivariate correlation of time series based on complex networks theory. Acta Physica Sinica, doi: 10.7498/aps.61.098902
[12]	Gao Zhong-Ke, Jin Ning-De, Yang Dan, Zhai Lu-Sheng, Du Meng. Complex networks from multivariate time series for characterizing nonlinear dynamics of two-phase flow patterns. Acta Physica Sinica, doi: 10.7498/aps.61.120510
[13]	Zhou Ting-Ting, Jin Ning-De, Gao Zhong-Ke, Luo Yue-Bin. Limited penetrable visibility graph for establishing complex network from time series. Acta Physica Sinica, doi: 10.7498/aps.61.030506
[14]	Cui Ai-Xiang, Fu Yan, Shang Ming-Sheng, Chen Duan-Bing, Zhou Tao. Emergence of local structures in complex network:common neighborhood drives the network evolution. Acta Physica Sinica, doi: 10.7498/aps.60.038901
[15]	Li Tao, Pei Wen-Jiang, Wang Shao-Ping. Optimal traffic routing strategy on scale-free complex networks. Acta Physica Sinica, doi: 10.7498/aps.58.5903
[16]	Lü Ling, Zhang Chao. Chaos synchronization of a complex network with different nodes. Acta Physica Sinica, doi: 10.7498/aps.58.1462
[17]	Wang Dan, Yu Hao, Jing Yuan-Wei, Jiang Nan, Zhang Si-Ying. Study on the congestion in complex network based on traffic awareness algorithm. Acta Physica Sinica, doi: 10.7498/aps.58.6802
[18]	Song Qing-Song, Feng Zu-Ren, Li Ren-Hou. Multiple clusters echo state network for chaotic time series prediction. Acta Physica Sinica, doi: 10.7498/aps.58.5057
[19]	Xu Dan, Li Xiang, Wang Xiao-Fan. An investigation on local area control of virus spreading in complex networks. Acta Physica Sinica, doi: 10.7498/aps.56.1313
[20]	Li Ji, Wang Bing-Hong, Jiang Pin-Qun, Zhou Tao, Wang Wen-Xu. Growing complex network model with acceleratingly increasing number of nodes. Acta Physica Sinica, doi: 10.7498/aps.55.4051

期刊类型引用(2)

1.	赵宁，江英华，周贤韬. 基于单光子的高效量子安全直接通信方案. 物理学报. 2022(15): 36-41 . 百度学术
2.	张晨涛，石小涛，朱文新，朱金龙，郝向英，金锐博. 利用域排列算法设计铌酸锂晶体实现3μm中红外波段频域纯态单光子源. 物理学报. 2022(20): 81-87 . 百度学术

其他类型引用(1)

Metrics

Abstract views: 3329
PDF Downloads: 33
Cited By: 3

姓名
邮箱
手机号码
标题
留言内容
验证码

Search

Article

留言板