程序如下:
4单元均匀线阵自适应波束形成图
clear
clc
format long
v=1
M=4
N=1000%%%%%%%快拍数
f0=21*10^3%%%%%%%%%%%信号与干扰的频率
f1=11*10^3
f2=15*10^3
omiga0=2*pi*f0%%%%%%%信号与干扰的角频率
omiga1=2*pi*f1
omiga2=2*pi*f2
sita0=08%信号方向
sita1=04%干扰方向1
sita2=21%干扰方向2
for t=1:N %%%%%%%%%%%%信号
adt(t)=sin(omiga0*t/(N*f0))
a1t(t)=sin(omiga1*t/(N*f1))
a2t(t)=sin(omiga2*t/(N*f2))
end
for i=1:M%%%%%%%%%%%%信号的导向矢量:线阵的银旦形式
ad(i,1)=exp(j*(i-1)*pi*sin(sita0))
a1(i,1)=exp(j*(i-1)*pi*sin(sita1))
a2(i,1)=exp(j*(i-1)*pi*sin(sita2))
end
R=zeros(M,M)
for t=1:N
x=adt(t)*ad+a1t(t)*a1+a2t(t)*a2%阵列对信号的完整响应
R=R+x*x'%信号的协方差矩阵
end
R=R/N%%%%%%%%%协方差矩阵,所有快拍数的平均
miu=1/(ad'*inv(R)*ad)%%%%%%这个貌似是LMS算法的公式,具体我记不太清,这里是求最优权值,根据这个公式求出,然后加权
w=miu*inv(R)*ad
%%%%%%形成波束%%%%%%%%%%%%%%%%%%%
for sita=0:pi/100:pi
for i=1:M
x_(i,1)=exp(j*(i-1)*pi*sin(sita))
end
y(1,v)=w'*x_%%%%%%%对信号进行加权,消除干扰
v=v+1
end
y_max=max(y(:))%%%%%%%%%%%%%%%归一化
y_1=y/y_max
y_db=20*log(y_1)
sita=0:pi/100:pi
plot(sita,y)
Xlabel(‘sitaa’)
Ylabel(‘天线增益db’)
4单元均匀线阵自适应波束形成
目标
clear
clc
format long
v=1
M=4阵元数
N=100
f0=21*10^3
omiga0=2*pi*f0
sita0=06%信号方向
for t=1:N
adt(t)=sin(omiga0*t/(N*f0))
end
for i=1:M
ad(i,1)=exp(j*(i-1)*pi*sin(sita0))
end
R=zeros(4,4)
r=zeros(4,1)
for t=1:N
x=adt(t)*ad
R=R+x*x'
end
R=R/N
miu=1/(ad'*inv(R)*ad)
w=miu*inv(R)*ad
for sita=0:pi/100:pi/2
for i=1:M
a(i,1)=exp(j*(i-1)*pi*sin(sita))
end
y(1,v)=w'斗贺*a
v=v+1
end
sita=0:pi/100:pi/2
plot(sita,y)
xlabel('sita')
ylabel('天线增益’)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%我的程序%%%%%%%%%%%%%%%
function jieshousignal
%期望信空搏派号数:1个
%干扰信号数:4个
%信噪比已知
%干燥比已知
%方位角已知
clc
clear all
close all
%//参数设置===========================================
1=0
2=0
3=0
% for rrr=1:16000
signal_num=1 %signal number
noise_num=5 %interference number
R0=06 %圆的半径
SP=2000 %Sample number
N=8 %阵元数
snr=-10%Signal-to-Noise
sir1=10 %Signal-to-Interference one
sir2=10 %Signal-to-Interference two
sir3=10 %Signal-to-Interf
sir4=10
sir5=10
%//================noise Power-to-signal Power====================
factor_noise_1=10^(-sir1/10)
factor_noise_2=10^(-sir2/10)
factor_noise_3=10^(-sir3/10)
factor_noise_4=10^(-sir4/10)
factor_noise_5=10^(-sir5/10)
factor_noise_targe=10^(-snr/10)
% //======================== ===============
d1=85*pi/180%%干扰1的方位角
d2=100*pi/180%干扰2的方位角
d3=147*pi/180%干扰3的方位角
d4=200*pi/180%干扰4的方位角
d5=250*pi/180%干扰5的方位角
d6=150*pi/180%目标的方位角
e1=15*pi/180%%干扰1的俯仰角
e2=25*pi/180%干扰2的俯仰角
e3=85*pi/180%干扰3的俯仰角
e4=50*pi/180%干扰4的俯仰角
e5=70*pi/180%干扰5的俯仰角
e6=85*pi/180%目标的俯仰角
% //====================目标信号==========================
t=1:1:SP
fc=2e7
Ts=1/(3e10)
S0=5*cos(2*pi*fc*t*Ts)%目标信号
for kk=1:N
phi_n(kk)=2*pi*(kk-1)/N
end
%//==================== *** 纵矢量==========================================
A=[conj(exp(j*2*pi*R0*cos(d6-phi_n)*sin(e6)))conj(exp(j*2*pi*R0*cos(d1-phi_n)*sin(e1)))conj(exp(j*2*pi*R0*cos(d2-phi_n)*sin(e2)))conj(exp(j*2*pi*R0*cos(d3-phi_n)*sin(e3)))conj(exp(j*2*pi*R0*cos(d4-phi_n)*sin(e4)))conj(exp(j*2*pi*R0*cos(d5-phi_n)*sin(e5)))]'
A1=[conj(exp(j*2*pi*R0*cos(d1-phi_n)*sin(e1)))conj(exp(j*2*pi*R0*cos(d2-phi_n)*sin(e2)))conj(exp(j*2*pi*R0*cos(d3-phi_n)*sin(e3)))conj(exp(j*2*pi*R0*cos(d4-phi_n)*sin(e4)))conj(exp(j*2*pi*R0*cos(d5-phi_n)*sin(e5)))]'
% //==========================================================Power of the interference
% // depending on the signal power and SIR
Ps1=0
Ps2=0
Ps3=0
Ps4=0
Ps5=0
S1=zeros(1,SP)
S2=zeros(1,SP)
S3=zeros(1,SP)
S4=zeros(1,SP)
S5=zeros(1,SP)
Ps0=S0*S0'/SP % signal power
Ps1=Ps0*factor_noise_1
Ps2=Ps0*factor_noise_2
Ps3=Ps0*factor_noise_3
Ps4=Ps0*factor_noise_4
Ps5=Ps0*factor_noise_5
% //==========================干扰信号的随机包络=========================
S1=normrnd(0,sqrt(Ps1/2),1,SP)+j*normrnd(0,sqrt(Ps1/2),1,SP)
S2=normrnd(0,sqrt(Ps2/2),1,SP)+j*normrnd(0,sqrt(Ps2/2),1,SP)
S3=normrnd(0,sqrt(Ps3/2),1,SP)+j*normrnd(0,sqrt(Ps3/2),1,SP)
S4=normrnd(0,sqrt(Ps4/2),1,SP)+j*normrnd(0,sqrt(Ps4/2),1,SP)
S5=normrnd(0,sqrt(Ps5/2),1,SP)+j*normrnd(0,sqrt(Ps5/2),1,SP)
%//
S=[S0S1S2S3S4S5]
SS1=[S1S2S3S4S5]
X=A*S%信号加干扰
XX2=A1*SS1%接收到的干扰
Pw_noise=sqrt(Ps0*factor_noise_targe/2)
a1=randn(N,SP)
a2=randn(N,SP)
a1=a1/norm(a1)
a2=a2/norm(a2)
W=Pw_noise*(a1+sqrt(-1)*a2)
X=X+W
% //--------------------------SMI算法----------------------------------------
Rd=X*S0'/SP
R=X*X'/(SP*1)
Wc_SMI=pinv(R)*Rd/(Rd'*pinv(R)*Rd)%权向量
Wc_SMI=Wc_SMI/norm(Wc_SMI)
Y_SMI=Wc_SMI'*X %SMI算法恢复出来的信号
%//-------------------------------------GS算法------------------
m=1
for i=1:400:2000
X2(:,m)=XX2(:,i)
m=m+1
end
a=zeros(1,8)
phi_n=zeros(1,8)
phi=0:pi/180:2*pi
theta=0:pi/180:pi/2
for kk=1:8
a(kk)=1
phi_n(kk)=2*pi*(kk-1)/8
end
x1=zeros(1,8)
x2=zeros(1,8)
x3=zeros(1,8)
x4=zeros(1,8)
x5=zeros(1,8)
x1=X2(:,1)'
x2=X2(:,2)'
x3=X2(:,3)'
x4=X2(:,4)'
x5=X2(:,5)'
Z1=x1
Z1_inner_product=Z1*conj(Z1)
Z1_mode=sqrt(sum(Z1_inner_product))
Y1=Z1/Z1_mode
Inner_product=sum(x2*conj(Y1))
Z2=x2-Inner_product*Y1
Z2_inner_product=sum(Z2*conj(Z2))
Z2_mode=sqrt(Z2_inner_product)
Y2=Z2/Z2_mode
Inner_product1=sum(x3*conj(Y1))
Inner_product2=sum(x3*conj(Y2))
Z3=x3-Inner_product1*Y1-Inner_product2*Y2
Z3_inner_product=sum(Z3*conj(Z3))
Z3_mode=sqrt(Z3_inner_product)
Y3=Z3/Z3_mode
Inner_product1_0=sum(x4*conj(Y1))
Inner_product2_0=sum(x4*conj(Y2))
Inner_product3_0=sum(x4*conj(Y3))
Z4=x4-Inner_product1_0*Y1-Inner_product2_0*Y2-Inner_product3_0*Y3
Z4_inner_product=sum(Z4*conj(Z4))
Z4_mode=sqrt(Z4_inner_product)
Y4=Z4/Z4_mode
Inner_product1_1=sum(x5*conj(Y1))
Inner_product2_1=sum(x5*conj(Y2))
Inner_product3_1=sum(x5*conj(Y3))
Inner_product4_1=sum(x5*conj(Y4))
Z5=x5-Inner_product1_1*Y1-Inner_product2_1*Y2-Inner_product3_1*Y3-Inner_product4_1*Y4
Z5_inner_product=sum(Z5*conj(Z5))
Z5_mode=sqrt(Z5_inner_product)
Y5=Z5/Z5_mode
%Y1
%Y2
%Y3
%Y4
%Y5
w0=zeros(1,8)
w=zeros(1,8)
for mm=1:8
w0(mm)=exp(-j*2*pi*R0*cos(d6-phi_n(mm))*sin(e6))
end
dd1=sum(w0*conj(Y1))*Y1
dd2=sum(w0*conj(Y2))*Y2
dd3=sum(w0*conj(Y3))*Y3
dd4=sum(w0*conj(Y4))*Y4
dd5=sum(w0*conj(Y5))*Y5
w=w0-dd1-dd2-dd3-dd4-dd5
Wc_GS=w
Wc_GS=Wc_GS/(norm(Wc_GS))
Y_GS=Wc_GS*X %GS算法恢复出来的图像
%//----------------------------------MMSE算法-----------------------
Rd=X*S0'/SP
R=X*X'/(SP*1)
Wc_MMSE=pinv(R)*Rd
Wc_MMSE=Wc_MMSE/norm(Wc_MMSE)
Y_MMSE=Wc_MMSE'*X %MMSE算法恢复出来的信号
S0=S0/norm(S0)
Y_GS=Y_GS/norm(Y_GS)
Y_SMI=Y_SMI/norm(Y_SMI)
Y_MMSE=Y_MMSE/norm(Y_MMSE)
% figure(1)
% plot(real(S0))
% title('原始信号')
% xlabel('采样快拍数')
% ylabel('信号幅度')
% figure(2)
% plot(real(Y_SMI))
% title('运用SMI算法处理出的信号')
% xlabel('采样快拍数')
% ylabel('信号幅度')
% figure(3)
% plot(real(Y_GS))
% title('运用G-S算法处理出的信号')
% xlabel('采样快拍数')
% ylabel('信号幅度')
% figure(4)
% plot(real(Y_MMSE))
% for i=1:SP
% ss(i)=abs(S0(i)-Y_SMI(i))^2
% end
% q_1=mean(ss)
% for i=1:SP
% ss1(i)=abs(S0(i)-Y_GS(i))^2
% end
% q_2=mean(ss1)
% for i=1:SP
% ss2(i)=abs(S0(i)-Y_MMSE(i))^2
% end
% q_3=mean(ss2)
%
% 1=1+q_1
% 2=2+q_2
% 3=3+q_3
% end
% 1/16000
% 2/16000
% 3/16000
phi=0:pi/180:2*pi
theta=0:pi/180:pi/2
%
% % //------------------------ 形成波束-----------------------------------------
F_mmse=zeros(91,361)
F_smi=zeros(91,361)
F_gs=zeros(91,361)
for mm=1:91
for nn=1:361
p1=sin(theta(mm))
p2=cos(phi(nn))
p3=sin(phi(nn))
q1=sin(e6)
q2=cos(d6)
q3=sin(d6)
for hh=1:8
w1=cos(phi_n(hh))
w2=sin(phi_n(hh))
zz1=q2*w1+q3*w2
zz2=p2*w1+p3*w2
zz=zz2*p1-zz1*q1
F_mmse(mm,nn)= F_mmse(mm,nn)+conj(Wc_MMSE(hh))*(exp(j*2*pi*R0*(zz2*p1)))
F_smi(mm,nn)=F_smi(mm,nn)+conj(Wc_SMI(hh))*(exp(j*2*pi*R0*(zz2*p1)))
F_gs(mm,nn)=F_gs(mm,nn)+conj((Wc_GS(hh))')*(exp(j*2*pi*R0*(zz2*p1)))
end
end
end
F_MMSE=abs(F_mmse)
F_SMI=abs(F_smi)
F_GS=abs(F_gs)
figure(5)
mesh(20*log10(F_MMSE))
figure(6)
mesh(20*log10(F_SMI))
title('SMI算法波束形成图')
xlabel('方位角')
ylabel('俯仰角')
zlabel('幅度/dB')
figure(7)
mesh(20*log10(F_GS))
title('G-S算法波束形成图')
xlabel('方位角')
ylabel('俯仰角')
zlabel('幅度/dB')
余平 刘方兰 肖波
第一作者简介:余平,男,高级工程师,1993年毕业于长春地质学院仪器系电子仪器及测量技术专业,现主要从事多波束技术应用与海洋地质调查技术管理工作。
(广州海洋地质调查局 广州 510760)
摘要 换能器阵元的不同排列组合决定其指向性,波束形成是多波束测量的关键技术。文中通过数学计算总结了不同换能器阵进行波束形成的工作原理,并介绍了利用二维DFT进行频域波束形成的一般方法。最后结合现役多波束测深系统,简单解释说明不同系统所采用的波束形成技术。
关键词 多波束 阵元 指向性 波束形成 测深
1 前言
我国自20世纪90年代初以来,为满足近海航道、大洋调查和国家经济专署区及大陆架勘测的需要,陆续从欧美等国家引进了大量的多波束测深系统(见表1),这些多波束测深系统涵盖了深水、中深水和浅水等不同海域,我国多波束技术应用迎来了第一个高峰期。
进入21世纪后,随着旧多波束测深系统的老化以及多波束新技术的推出,多波束测深系统的更新换代已经展开,高精度、高覆盖、高波束数的多波束系统在一些专项中开始应用。在多波束测深系统的实际使用中,从事多波束测量的技术人员针对不同多波束测深系统所存在的问题进行了大量的研究工作,并出版了多波束技术专著,撰写了大量的论文。在这些应用型的研究成果中,关于多波束测深系统工作原理的关键技术——波束形成技术,要么是一个简单的比喻,要么是笼统大概的说明。本文试图在总结不同形式的波束形成原理的基础上,结纳磨合实际应用,阐述不同系和茄此统波束形成的模式,从而进一步理解多波束测深系统的工作原理。
2 波束形成原理
所谓波束形成是指将一定几何形状(直线、圆柱、弧形等)排列的多元基阵各阵元输出经过处理(例如加权、时延、求和等)形成空间指向性的方法(田坦等,2000)。波束形成也是将一个多元阵经适当处理使其对某些空间方向的声波具有所需响应的方法。波束形成的方法有很多,特别是在实际应用中,随着微电子技术、计算技术的快速发展,数字信号处理技术使时域、频域下的波束形成方法相互贯穿。
表1 我国目前已安装并使用的多波束测深系统(2004年前)Table1 Multibeam sound system has been installed and used in China(Before 2004)
2.1 波束形成一般原理
波束形成技术来自于基阵具有方向性的原理(蒋楠祥,2000)。设一个由N个无方向性阵元组成的接收换能器阵(如图1)。各阵元位于空间唤迅点(xn,yn,zn)处,将所有阵元的信号相加得到输出,就形成了基阵的自然指向性。此时,若有一远场平面波入射到这一基阵上,它的输出幅度将随平面入射角的变化而变化。
当信号源在不同方向时,由于各阵接收信号与基准信号的相位差不同,因而形成的和输出的幅度不同,即阵的响应不同。
如果上述阵是一N元线阵,阵元间距为d,各阵元接收灵敏度相同,平面波入射方向为θ(如图2)。各阵元输出信号为:
F0(t)=Acos(ωt)(1)
南海地质研究.2005
……
图2 线阵几何形状
Fig.2 Geometry shape of line array transducer
南海地质研究.2005
其中A为信号幅度;ω为信号角频率;φ为相邻阵元接收信号间的相位差,Re为取实部,有:
南海地质研究.2005
所以,阵的输出为:
南海地质研究.2005
因为:
南海地质研究.2005
则:
南海地质研究.2005
所以:
南海地质研究.2005
上式两边同时除以NA进行归一化处理,得:
南海地质研究.2005
R(θ)表明,一个多元阵输出幅度大小随信号入射角而变化。一般而言,对于一个任意的阵形,无论声波从哪一个方向入射,均不可能形成同相相加或得到最大输出,只有直线阵或空间平面阵才会在阵的法线方向形成同相相加,得到最大输出。然而,任意阵形的阵经过适当的处理,可在预定方向形成同相相加,得到最大输出,这就是波束形成的一般原理。
2.2 直线阵相移波束形成
在前面讨论的基础上,直线阵相移波束形成的根本目的是:在相邻阵元之间插入相移β,则直线阵的求和输出为:
南海地质研究.2005
归一化阵输出幅度变为:
南海地质研究.2005
所以主波束方向满足:
φ-β=0
即:
南海地质研究.2005
所以:
南海地质研究.2005
或:
南海地质研究.2005
上式表明:在阵元间插入不同的相移β,可以控制主波束位于不同的方向,这种在阵元之间插入相移使主波束方向控制于不同方位的方法称为相移波束形成。在窄带(主动声呐)应用中,一般常用相移波束形成方法。
2.3 直线阵时延波束形成
在直线阵相移波束形成的讨论中,有:
南海地质研究.2005
因为:
β=2πfτ
所以:
南海地质研究.2005
上式表明:在阵元间插入不同的时延τ,可以控制主波束位于不同的方向,这种在阵元之间插入时延使主波束方向控制于不同方位的方法称为时延波束形成。在宽带(被动声呐)应用中,一般常用时延波束形成方法。
2.4 圆阵波束形成
圆形阵的阵元一般均匀分布在圆周上。由于圆阵是几何上关于原点对称的,因而没有方向性。无自然的指向性波束,必须对阵元信号进行延迟或相移才能形成方向性,即使其补偿成一个等效的线阵。简单的实现方法是电子开关波束形成方法,这种方法利用电子开关进行控制,将一组延迟线接入不同阵元,以形成不同方位的波束。
以16元圆阵为例说明。假定只用圆弧上的七个阵元形成波束(如图3),如果目标信号从正前方来,为了形成同相相加,必须将各阵元信号延迟补偿到图中所示的直线(蓝色)上。设两相邻阵元所在圆弧的圆心角为α0,则各阵元所需的相应延迟为:
南海地质研究.2005
τ1=τ7=0(15)
南海地质研究.2005
2.5 弧形阵波束形成
弧形阵的波束形成是圆阵波束形成的一种特殊情况,分布在弧形阵上阵元最终必须投影到一个等效的线阵中。如以时延来完成指向性的控制,各阵元的时延算法与“圆阵波束形成”的例子相同。
2.6 频域波束形成
从前面讨论中可知,一个波束形成器可对空间某方位的信号有响应,而抑制其它方位的信号,因此,波束形成实际上是一种空间滤波过程。根据线性系统理论,波束形成也是一种卷积运算,因而可用频域的乘积实现。所以波束也可以在频域内形成,这就是频域波束形成。频域波束形成常采用离散傅里叶变换(DFT),可以用数字信号处理中的快速傅里叶变换(FFT)加以实现,因此频域波束形成比时间域波束形成运算量要小(曹洪泽等,2002)。
设均匀间隔直线阵有N个阵元,间距为d。对阵元i的输出信号xi(t)进行采样,取L点作DFT运算,即:
南海地质研究.2005
其中i为阵元号,k为谱线号,l为时间序号。因此Xi(k)表示第i号阵元接收的时间序列的谱。
其次,对同一序号k的谱线作空间傅里叶变换,将Xi(k)重排为Xk(i),进行下列运算:
南海地质研究.2005
其中m为波束号;wi为阵元的幅度权值;Yk(m)代表k号频率分量的第m号波束输出。这就是利用二维DFT实现频域波束形成的方法。
3 结论
综上所述,换能器的指向性是波束形成原理的基础。目前我国现役的多波束测深系统主要包括SeaBea m系列、Elac Botto mChart系列、EM系列、SeaBat系列和Atlas DS系列等[4],由于各系统生产厂家和工作水深范围不同,多波束系统采用的换能器、发射频率不同,因此,不同系统采用的波束形成方法也不尽相同。
Sea Bea m 2112深水多波束测深系统发射频率12 KHz,发射器和水听器独立安装,其中发射器14个模块,水听器8个模块共80个通道。水听器是4个模块一组共两组呈“V”型安装,换能器是典型的“米勒十字交叉”(Mill s Cr oss)安装模式。即便如此,波束形成原理符合直线阵相移波束形成原理。1998年8月,厂家根据合同对系统进行升级,在仅更换DSP 板的情况下,使系统的波束数从121个升级为151个,应该是运用了高级数字信号处理器完成的直线阵相移波束形成下的数字内插波束形成技术(移位边带波束形成)。EM120深水多波束测深系统的发射接收器也是独立安装,属于线性的“米勒十字交叉”结构阵,其基本的波束形成原理也是符合直线相移波束形成原理,由于其波束数已大大提高,应该还综合有频域波束形成技术。
EM950(或EM1002)中深水多波束测深系统发射频率95kHz,发射器和水听器二合一安装,波束数120个。换能器是一个半径为45cm的半圆弧形阵,作为一个高发射频率的主动声呐系统,采用的是弧形阵时延和相移波束形成技术的综合。EM3000浅水多波束测深系统发射频率300kHz,波束数120个,换能器是一个圆形阵(李家彪等,周兴华等,1999),采用技术与EM950类似。
SeaBat系列多波束系统在国内主要以浅水多波束测深系统为主,浅水多波束系统的换能器一般都是采用发射器和水听器二合一安装方式。SeaBat8101多波束测深系统的发生频率240kHz,波束数101个。换能器是一个直径为32cm的圆形阵,采用的波束形成方式与EM系列的类似。
Atlas Fansweep系列是利用侧扫声呐技术计算多个水深数据的多波束测深系统,与真正多波束测深系统比较起来技术指标相对落后。由于厂家产品开发战略转变的原因,深水多波束系统在近两年才推出。Atlas DS系列多波束系统在国内还没有用户,据称其新一代多波束系统采用了Chirp技术,接收波束数将超过300个,因此其波束形成技术应该主要以频域波束形成技术为主。
参考文献
曹洪泽,李蕾等.2002.一种基于FFT 波束形成的BDI 算法分析研究.海洋技术,21(2),55~59
蒋楠祥.2000.换能器与基阵.哈尔滨:哈尔滨工程大学出版社,50~75
李家彪等.1999.多波束勘测原理技术与方法.北京:海洋出版社,6~9
田坦,刘国枝,孙大军.2000.声呐技术.哈尔滨:哈尔滨工程大学出版社,63~120
周兴华,刘忠臣,傅命左等.1999.多波束海底地形勘测技术规程.8~14
Multibeam Pivotal Technology——Beam Forming
Yu Ping Liu Fanglan Xiao Bo
(Guangzhou Marine Geological Survey,Guangzhou,510760)
for different type of transducer,and introduces a universal way of frequency domain beam forming by using 2⁃dimension DFT.Finally,the author simply explains the different beam forming technology which the multibeam have in use.Abstract:Different arranged transducer deter mines the directional property of a transducer array.Multibeampivotal technology——the basis of Beamformingis howto control the directional property of transducer.This article summarizes the theory of beamfor ming with mathematics operation
Key words:Multibeam Transducer Directional Property Beam Forming Sound
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)