CN106950569B - 基于序贯回归方法的多阵元合成孔径聚焦波束形成方法 - Google Patents
基于序贯回归方法的多阵元合成孔径聚焦波束形成方法 Download PDFInfo
- Publication number
- CN106950569B CN106950569B CN201710076139.8A CN201710076139A CN106950569B CN 106950569 B CN106950569 B CN 106950569B CN 201710076139 A CN201710076139 A CN 201710076139A CN 106950569 B CN106950569 B CN 106950569B
- Authority
- CN
- China
- Prior art keywords
- submatrix
- signal
- array
- array element
- indicate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52003—Techniques for enhancing spatial resolution of targets
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法,该发明方法利用多阵元合成孔径聚焦波束形成和最小均方误差准则,通过期望信号与自适应系统的输出作差得到输出误差,并将这个误差引入到序贯回归方法中,再用序贯回归方法的输出作为线性传感器阵列各通道的权值,对多阵元合成孔径聚焦波束形成的子阵与子阵阵元分别进行动态幅度遍迹,有效减少了波束形成的主瓣宽度和旁瓣幅度,提高了成像质量。
Description
技术领域
本发明涉及合成孔径波束形成领域,特别是一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法。
背景技术
合成孔径聚焦成像相比于传统的波束形成技术,其优势在于合成孔径聚焦每次发射、接收只有一个阵元,大大简化了电路;同时可以实现发射和接收的双向聚焦(见文献:P.Karwat,Z.Klimonda,M.Lewandowski.A.Nowicki.Imaging quality of the classicalbeamforming,SAFT and plane wave imaging-Experimental results[A].2012 IEEEInternational Ultrasonics Symposium[C].New York:IEEE Trans,2012:1-4.)。但是其缺点也很明显,合成孔径聚焦采用单阵元换能器,但是单阵元换能器具有扫描频率低、成像信噪比差的缺点,同时,合成孔径聚焦的副瓣等级限制了它的性能(见文献:Mok KunJeong,Sung Jae Kwon.A novel side lobe estimation method in medical ultrasoundimaging systems[A].Ultrasonics Symposium(IUS),2015IEEE International[C].NewYork:IEEE Trans,2015:1-4.)。为进一步提高成像质量,一些改进的方法被提出来了。针对合成孔径聚焦技术信噪比低的特点,提出多阵元合成孔径聚焦的方法,提高了成像质量和分辨率(见文献:杜英华,祁欣.多阵元合成孔径聚焦波束形成研究[J].纳米技术与精密工程.2014,12(3):162-166)。C.J.Martín-Arguedas等提出一种加速多阵元合成孔径聚焦成像的方法,减少了计算量(见文献:C.J.Martín-Arguedas,O.Martínez-Graullera,D.Romero-Laorden,L.Gómez-Ullate.Method and architecture to accelerate multi-element syntheticaperture imaging[J].Digital Signal Processing.2013,23(4):1288-1295)。岳海林提出一些改进发射、接收模式的方法:如合成接收孔径(SRA)、合成发射孔径(STA)等(见文献:岳海林.超声合成聚焦信号噪声抑制与相位校正技术研究[D].四川:西南交通大学,2015:7-20)。Jonathan Reeg等提出空减影成像(Null SubtractionImaging,NSI)的方法,降低旁瓣、提高横向分辨率(见文献:Jonathan Reeg,MichaelL.Oelze.Improving lateral resolution in ultrasonic Imaging by utilizing nullsin the beam pattern[J].Ultrasonics Symposium(IUS),2015 IEEE International[C].New York:IEEE Trans,2015:1-4)。Y.Tasinkevych等提出对合成发射孔径进行变迹加权的方法,提高成像质量(见文献:Y.Tasinkevych,Z.Klimonda,M.Lewandowski,A.Nowicki,P.A.Lewin.Modified multi-element synthetic transmit aperture method forultrasound Imaging:A tissue phantom study[J].Ultrasonics.2013,53(2):570-57)。李荣兴引入空间平滑对角加载Capon自适应波束形成方法,实现真正的动态幅度变迹(见文献:李荣兴.波束形成自适应波束形成方法研究[D].哈尔滨:哈尔滨工业大学,2013:52-55)。刘广在李荣兴的研究基础上引入相干因数和平面波的概念,进一步提高了成像质量(见文献:刘广.合成孔径波束形成中波束合成方法研究[D].河北:燕山大学,2015:41-48)。
上述方法没有有效解决波束主瓣宽度宽、旁瓣幅度大的问题,影响了成像分辨率和对比度。
发明内容
本发明所要解决的技术问题是克服现有技术的不足而提供一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法,该方法利用自适应序贯回归方法代替传统的Hanning窗对接收信号进行幅度变迹,降低了主瓣宽度,压制了旁瓣幅度,提高了成像分辨率和对比度。
本发明为解决上述技术问题采用以下技术方案:
根据本发明提出的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法,包括如下步骤:
步骤1、将一个阵元个数为N的线性传感器阵列分割成长度为L的子阵,子阵的个数为M,发射、接收使用相同的阵列;
步骤2、对于空间中的一点P,激励超声阵列发射单一频率的超声平面波信号s(t),经过接收阵列得到接收信号其中,t表示连续的时间参量,r表示P点到阵列中心点的参考声程,c表示声速,θ表示P点与阵列中心点的连线和扫描线垂直平面的夹角,τn(θ)表示第n个子阵的发射延时,其中,1≤n≤M,τi,n(θ)表示第n个子阵中第i个阵元的接收延时,其中,1≤i≤L,不计参考声程r对波束形成的影响,则将p(t,r,θ)改写为p(t,θ)=s(t-τn(θ)-τi,n(θ)),p(t,θ)表示t时刻θ方位上的接收信号;
步骤3、将步骤2所述的接收信号p(t,θ)表示为p(t,θ)=s(t)s(-τn(θ))s(-τi,n(θ)),由于s(t)为单一频率的平面波信号,p(t,θ)只与θ有关,即p(t,θ)=p(θ)=s(-τn(θ))s(-τi,n(θ)),其中,s(-τn(θ))表示第n个子阵接收信号,s(-τi,n(θ))表示第n个子阵中第i个阵元的接收信号;
步骤4、将所述步骤3中s(-τn(θ))写成向量形式s(θ)=[s(-τ1(θ)),s(-τ2(θ)),…,s(-τM(θ))]T,其中,上标T表示转置,s(θ)表示阵列接收信号,考虑噪声noise的影响后得到的接收信号为y(θ)=[s(-τ1(θ)),s(-τ2(θ)),…,s(-τM(θ))]T+noise,其中,noise为一组与子阵个数一致的高斯噪声,对每个子阵设置一组权向量W1,利用幅度遍迹得到遍迹信号l1(θ),其中,W1=[W1,1,W1,2,…,W1,M];W1,n表示第n个子阵的1次迭代权值,将s(θ)与l1(θ)作差得序贯回归方法的约束误差ε1(θ)=s(θ)-l1(θ),利用最小方差准则与迭代方法获得子阵优化权向量Wk+1=[Wk+1,1,Wk+1,2,…,Wk+1,M],Wk+1,n表示第n个子阵的k+1次迭代权值,利用子阵优化权向量进行幅度遍迹得到子阵的第k+1迭代的遍迹信号lk+1(θ);
步骤5、将所述步骤3中s(-τi,n(θ))写成向量形式sn(θ)=[s(-τ1,n(θ)),s(-τ2,n(θ)),…,s(-τL,n(θ))],考虑noise′的影响后所得到的接收信号为yn(θ)=[s(-τ1,n(θ)),s(-τ2,n(θ)),…,s(-τL,n(θ))]+noise′,其中,sn(θ)表示第n个子阵接收信号,noise′为一组与子阵阵元个数一致的高斯噪声,利用最小均方误差准则、迭代方法以及幅度遍迹方法得到子阵阵元优化权向量Wk+1,n与第n个子阵的第k+1迭代的遍迹信号lk+1,n(θ),其中Wk+1,n=[Wk+1,1,n,Wk+1,2,n…,Wk+1,L,n],Wk+1,i,n表示第n个子阵中第i个阵元的k+1次迭代权值;
步骤6、利用叠加方法,由所述步骤4中lk+1(θ)与所述步骤5中lk+1,n得到基于序贯回归方法的多阵元合成孔径聚焦波束Gk+1(θ),Gk+1(θ)=∑lk+1(θ)∑lk+1,n(θ)。
作为本发明所述的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法进一步优化方案,所述步骤1中子阵的个数M=N-L+1。
作为本发明所述的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法进一步优化方案,所述步骤3中接收信号p(θ)的获取步骤如下:
步骤3-1、将单一频率平面波信号表示为其中,表示以自然常数为底的指数函数,j表示虚数单位,ω0为角频率;则将所述步骤3中p(t,θ)表示为
步骤3-2、步骤3-1中第n个子阵的发射延时第n个子阵中第i个阵元接收的接收延时其中,xn表示P点到第n个发射子阵中心的距离,xi,n表示P点到第n个子阵中第i个阵元的距离;由此得到步骤3中其中,k0表示波数,
步骤3-3、不计单一频率平面波信号对所述步骤3-2中接收信号p(t,θ)的影响,得
作为本发明所述的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法进一步优化方案,所述步骤4中的幅度遍迹方法为一种控制发射和接收声场分布的手段,幅度变迹方法通过收发通道幅度加权的方式,使中心阵元的信号幅度增强,而两旁阵元的信号幅度减弱;所述步骤4中l1(θ)=yT(θ)W1,lk+1(θ)=yT(θ)Wk+1。
作为本发明所述的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法进一步优化方案,所述步骤4中子阵优化权向量Wk+1获取步骤如下:
步骤A、将均方误差约束波束形成问题的均方误差表示为E(|εk(θ)|2)=E[(s(θ)-lk(θ))H(s(θ)-lk(θ))],其中,E表示数学期望,误差信号εk(θ)表示期望信号与实际信号的差,H表示共轭转置,且εk(θ)=s(θ)-lk(θ),lk(θ)=yT(θ)Wk,Wk表示第k次迭代权向量;E(|εk(θ)|2)对Wk求梯度为
步骤B、由最小均方误差方法获得其中,为的估计;μ控制自适应速度和稳定性的增益常数,当最小均方误差方法收敛时0<μ<1;λav为信号y(θ)自相关函数Rk特征值的平均;
步骤C、利用函数式经过k次迭代,得到进一步获得y(θ)的自相关函数的逆矩阵其中,α为遗忘因子且0<α<1,Qk表示估计自相关函数的临时参量,*表示共轭,表示Qk的逆矩阵;
步骤D、由步骤C所述的得
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)本发明采用多阵元合成孔径聚焦波束形成的方法,相比于合成孔径聚焦波束形成,解决了每次收发单阵元的缺陷,提高了信噪比;
(2)本发明中采用自适应序贯回归变迹方法代替传统Hanning窗幅度变迹的方法,具有动态幅度变迹、沿着最佳权值点迭代、收敛速度快的特点;
(3)本发明实施例表明,与采用Hanning窗幅度变迹的方法相比,波束图具有更好的减少主瓣宽度和压制旁瓣的效果。
附图说明
图1是本发明一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法原理图。
图2是本发明方法中MSAF的发射/接收模式。
图3是本发明方法中MSAF的扫描方式。
图4是SAFT波束图,(a)阵元间距(b)阵元间距
图5是本发明方法中MSAF聚焦波束图,(a)子阵阵元个数L=4,(b)子阵阵元个数L=16,(c)子阵阵元个数L=32。
图6是SER方法与LMS方法比较,其中,(a)为SER方法与LMS方法性能表面图,(b)为LMS方法最小均方误差图,(c)为迭代初值q0=10,SER方法最小均方误误差曲线,(d)为迭代初值q0=1,SER方法最小均方误误差曲线。
图7是Hanning窗变迹MSAF波束图。
图8是本发明方法波束图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法原理图,如图1所示;图1中s(t)为阵列发射信号,p(t,r,θ)为阵列接收信号,p(θ)为初步处理的信号,s(θ)表示接收信号,sn(θ)表示第n个子阵接收信号,y(θ)、yn(θ)分别为阵列与第n个子阵受高斯噪声影响后得到的接收信号,Wk、Wk,n为阵列与第n个子阵第k次迭代的权向量,lk(θ)、lk,n(θ)为阵列与第n个子阵第k次迭代的遍迹信号,误差信号εk(θ)为信号s(θ)与信号lk(θ)的差,误差信号εk,n(θ)为信号sn(θ)与信号lk,n(θ)的差,lk+1(θ)、lk+1,n(θ)分别为阵列和第n个子阵的第k+1迭代的遍迹信号,Gk+1(θ)为多阵元合成孔径聚焦波束。
图2为本发明方法中多阵元合成孔径聚焦(Multi-element Synthetic ApertureFocusing,MSAF)发射/接收模式,将一个阵元个数为N的线性传感器阵列分割成一系列长度为L的子阵,子阵的个数M=N-L+1,发射、接收使用相同的阵列。
图3为本发明方法中MSAF的扫描方式,P为空间任意一点;阵元中心点为参考阵元;r表示P点到阵列中心点的参考声程;rn为第n个子阵的发射声程;ri,n为第n个子阵中第i个阵元的接收声程;d为阵元间距。
发射信号为单一频率平面波信号,即
图1中,
不计参考距离r对波束形成的影响,式(2)中p(t,r,θ)写成
p(t,θ)=s(t-τn(θ)-τi,n(θ)) (3)
表示t时刻θ方位上的接收信号。
由于为单一频率平面波信号,p(t,θ)只与θ有关,将p(t,θ)写为
p(θ)=s(-τn(θ)-τi,n(θ)) (4)
将式(4)代入式(1),得
第n个子阵的发射延时为
第n个子阵中第i个阵元的接收延时τi,n(θ)为
将式(6)、(7)代入(5),得
式(8)中,k0表示波数。
图1中
均方误差约束波束形成问题表示为
E(|εk(θ)|2)=E[(s(θ)-lk(θ))H(s(θ)-lk(θ))] (11)
式中,εk(θ)=s(θ)-lk(θ),lk(θ)=yT(θ)Wk。
E(|εk(θ)|2)对Wk求梯度为为
用最小均方方法(LMS,Least mean square)表示权向量迭代公式为
式(12)代入式(13),得
利用函数式
式(15)化简为
式(15)中α,0<α<1为遗忘因子,起到加重现时刻分量的贡献及消退记忆的作用。
由式(16)估计自相关函数为
式(17)代入式(14)
式(18)中,的更新公式为
q0为迭代初值,I表示单位矩阵
最终的多阵元合成孔径聚焦波束为
实施例
以下实施例所用参数如下:
图3是MSAF的扫描方式,超声换能器个数N=32,阵元间距和其中,λ表示声波的波长;线阵中心位置为坐标原点;子阵阵元个数L=4、L=16和L=32,中心频率f0=3.5MHz;声速c=1540mm/μs,波束扫描角度θ=-90°~+90°,r表示参考距离,rn表示发射声程,ri,n表示接收声程;发射信号exp(j2πf0t);高斯噪声的均值为0,方差为1。
图4是合成孔径聚焦(Synthetic Aperture Focusing Technique,SAFT)波束图。图4中的(a)阵元间距图4中的(b)阵元间距图4表明,单阵元的合成孔径聚焦信噪比过低;而且当阵元间距时,会出现较大的栅瓣,成像时会出现伪影,影响成像质量;当阵元间距时,消除了旁瓣,然而阵元间距的减小,却增加了主瓣宽度。
图5为本发明方法中MSAF聚焦波束图焦波束图。图5中的(a)阵元间距子阵阵元个数L=4,图5中的(b)阵元间距子阵阵元个数L=16,图5中的(c)阵元间距子阵阵元个数L=32。与图4相比,多阵元合成孔径聚焦方法有效地压制了旁瓣幅度、提高了信噪比;当子阵阵元个数时,波束效果比较好。
图6为SER方法与最小均方误差(Least Mean Square,LMS)方法的比较图。图6中的(a)为SER方法与LMS方法性能表面图,图6中的(b)为LMS方法最小均方误差图,图6中的(c)为迭代初值q0=10,SER方法最小均方误差曲线,图6中的(d)为迭代初值q0=1,SER方法最小均方误差曲线。图6中的(a)中,SER方法的收敛过程是沿着直接指向最佳权值的方向收敛的,而LMS则是稍微偏离了一段距离,然后指向最佳权值收敛;从图6中的(b)、(c)、(d)可知,SER方法比LMS方法收敛速度快,且q0值越大,收敛速度越快。
图7为Hanning窗变迹MSAF波束图(HanningMSAF)。图中采用了Hanning窗变迹加权的方法,其中阵元间距子阵阵元个数L=4。与图4、5相比,压制了旁瓣,然而却增加了主瓣宽度。提高了成像对比度却降低了分辨率。
图8为本发明方法波束图。图中采用了SER变迹加权的方法,其中阵元间距子阵阵元个数L=4,迭代初值q0=10。与图7相比,成功实现了主瓣宽度和旁瓣幅度的双重压制、提高了成像分辨率和对比度。可见,将序贯回归方法应用到多阵元合成孔径聚焦中,能有效减少主瓣宽度和压制旁瓣幅度。
本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段,还包括由以上技术特征任意组合所组成的技术方案。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。
Claims (5)
1.一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法,其特征在于,包括如下步骤:
步骤1、将一个阵元个数为N的线性传感器阵列分割成长度为L的子阵,子阵的个数为M,发射、接收使用相同的阵列;
步骤2、对于空间中的一点P,激励超声阵列发射单一频率的超声平面波信号s(t),经过接收阵列得到接收信号其中,t表示连续的时间参量,r表示P点到阵列中心点的参考声程,c表示声速,θ表示P点与阵列中心点的连线和扫描线垂直平面的夹角,τn(θ)表示第n个子阵的发射延时,其中,1≤n≤M,τi,n(θ)表示第n个子阵中第i个阵元的接收延时,其中,1≤i≤L,不计参考声程r对波束形成的影响,则将p(t,r,θ)改写为p(t,θ)=s(t-τn(θ)-τi,n(θ)),p(t,θ)表示t时刻θ方位上的接收信号;
步骤3、将步骤2所述的接收信号p(t,θ)表示为p(t,θ)=s(t)s(-τn(θ))s(-τi,n(θ)),由于s(t)为单一频率的平面波信号,p(t,θ)只与θ有关,即p(t,θ)=p(θ)=s(-τn(θ))s(-τi,n(θ)),其中,s(-τn(θ))表示第n个子阵接收信号,s(-τi,n(θ))表示第n个子阵中第i个阵元的接收信号;
步骤4、将所述步骤3中s(-τn(θ))写成向量形式s(θ)=[s(-τ1(θ)),s(-τ2(θ)),…,s(-τM(θ))]T,其中,上标T表示转置,s(θ)表示阵列接收信号,考虑噪声noise的影响后得到的接收信号为y(θ)=[s(-τ1(θ)),s(-τ2(θ)),…,s(-τM(θ))]T+noise,其中,noise为一组与子阵个数一致的高斯噪声,对每个子阵设置一组权向量W1,利用幅度遍迹得到遍迹信号l1(θ),其中,W1=[W1,1,W1,2,…,W1,M];W1,n表示第n个子阵的1次迭代权值,将s(θ)与l1(θ)作差得序贯回归方法的约束误差ε1(θ)=s(θ)-l1(θ),利用最小方差准则与迭代方法获得子阵优化权向量Wk+1=[Wk+1,1,Wk+1,2,…,Wk+1,M],Wk+1,n表示第n个子阵的k+1次迭代权值,利用子阵优化权向量进行幅度遍迹得到子阵的第k+1次迭代的遍迹信号lk+1(θ);
步骤5、将所述步骤3中s(-τi,n(θ))写成向量形式sn(θ)=[s(-τ1,n(θ)),s(-τ2,n(θ)),…,s(-τL,n(θ))],考虑noise′的影响后所得到的接收信号为yn(θ)=[s(-τ1,n(θ)),s(-τ2,n(θ)),…,s(-τL,n(θ))]+noise′,其中,sn(θ)表示第n个子阵接收信号,noise′为一组与子阵阵元个数一致的高斯噪声,利用最小均方误差准则、迭代方法以及幅度遍迹方法得到子阵阵元优化权向量Wk+1,n与第n个子阵的第k+1次迭代的遍迹信号lk+1,n(θ),其中Wk+1,n=[Wk+1,1,n,Wk+1,2,n…,Wk+1,L,n],Wk+1,i,n表示第n个子阵中第i个阵元的k+1次迭代权值;
步骤6、利用叠加方法,由所述步骤4中lk+1(θ)与所述步骤5中lk+1,n(θ)得到基于序贯回归方法的多阵元合成孔径聚焦波束Gk+1(θ),Gk+1(θ)=∑lk+1(θ)∑lk+1,n(θ)。
2.根据权利要求1所述的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法,其特征在于,所述步骤1中子阵的个数M=N-L+1。
3.根据权利要求1所述的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法,其特征在于,所述步骤3中接收信号p(θ)的获取步骤如下:
步骤3-1、将单一频率平面波信号表示为其中,表示以自然常数为底的指数函数,j表示虚数单位,ω0为角频率;则将所述步骤3中p(t,θ)表示为
步骤3-2、步骤3-1中第n个子阵的发射延时第n个子阵中第i个阵元的接收延时其中,xn表示P点到第n个发射子阵中心的距离,xi,n表示P点到第n个子阵中第i个阵元的距离;由此得到步骤3中其中,k0表示波数,
步骤3-3、不计单一频率平面波信号对所述步骤3-2中接收信号p(t,θ)的影响,得
4.根据权利要求1所述的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法,其特征在于,所述步骤4中的幅度遍迹方法为一种控制发射和接收声场分布的手段,幅度变迹方法通过收发通道幅度加权的方式,使中心阵元的信号幅度增强,而两旁阵元的信号幅度减弱;所述步骤4中l1(θ)=yT(θ)W1,lk+1(θ)=yT(θ)Wk+1。
5.根据权利要求1所述的一种基于序贯回归方法的多阵元合成孔径聚焦波束形成方法,其特征在于,所述步骤4中子阵优化权向量Wk+1获取步骤如下:
步骤A、将最小方差准则约束波束形成问题的最小方差准则表示为E(|εk(θ)|2)=E[(s(θ)-lk(θ))H(s(θ)-lk(θ))],其中,E表示数学期望,误差信号εk(θ)表示期望信号与实际信号的差,H表示共轭转置,且εk(θ)=s(θ)-lk(θ),lk(θ)=yT(θ)Wk,Wk表示第k次迭代权向量;E(|εk(θ)|2)对Wk求梯度为
步骤B、由最小方差准则获得其中,为的估计;μ为控制自适应速度和稳定性的增益常数,当最小均方误差方法收敛时0<μ<1;λav为信号y(θ)自相关函数Rk特征值的平均;
步骤C、利用函数式经过k次迭代,得到进一步获得y(θ)的自相关函数的逆矩阵其中,α为遗忘因子且0<α<1,Qk表示估计自相关函数的临时参量,*表示共轭,表示Qk的逆矩阵;
步骤D、由步骤C所述的得
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710076139.8A CN106950569B (zh) | 2017-02-13 | 2017-02-13 | 基于序贯回归方法的多阵元合成孔径聚焦波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710076139.8A CN106950569B (zh) | 2017-02-13 | 2017-02-13 | 基于序贯回归方法的多阵元合成孔径聚焦波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106950569A CN106950569A (zh) | 2017-07-14 |
CN106950569B true CN106950569B (zh) | 2019-03-29 |
Family
ID=59466386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710076139.8A Expired - Fee Related CN106950569B (zh) | 2017-02-13 | 2017-02-13 | 基于序贯回归方法的多阵元合成孔径聚焦波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106950569B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108155957B (zh) * | 2017-10-25 | 2021-01-01 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 数字多波束系统发射零值在线监测方法 |
CN114047256B (zh) * | 2021-10-25 | 2023-10-20 | 扬州大学 | 基于动态阵元合成孔径聚焦的平板陶瓷膜缺陷超声成像方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010116153A1 (en) * | 2009-04-09 | 2010-10-14 | Ntnu Technology Transfer As | Optimal modal beamformer for sensor arrays |
CN102749622A (zh) * | 2012-07-03 | 2012-10-24 | 杭州边界电子技术有限公司 | 基于多波束测深的声速剖面及海底地形的联合反演方法 |
CN102895000A (zh) * | 2012-11-06 | 2013-01-30 | 重庆大学 | 一种基于自适应加权的双聚焦波束合成方法 |
EP2702426A1 (en) * | 2011-04-29 | 2014-03-05 | B-K Medical ApS | Harmonic ultrasound imaging using synthetic aperture sequential beamforming |
CN103728591A (zh) * | 2013-12-17 | 2014-04-16 | 河海大学 | 一种mimo雷达近场目标高效实波束方位聚焦方法 |
CN105223544A (zh) * | 2015-08-26 | 2016-01-06 | 南京信息工程大学 | 近场线性约束最小方差自适应加权频率不变波束形成方法 |
CN105681972A (zh) * | 2016-01-14 | 2016-06-15 | 南京信息工程大学 | 线性约束最小方差对角加载的稳健频率不变波束形成方法 |
CN105891835A (zh) * | 2016-06-16 | 2016-08-24 | 北京海卓同创科技有限公司 | 一种实时动态聚焦波束形成方法及系统 |
-
2017
- 2017-02-13 CN CN201710076139.8A patent/CN106950569B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010116153A1 (en) * | 2009-04-09 | 2010-10-14 | Ntnu Technology Transfer As | Optimal modal beamformer for sensor arrays |
EP2702426A1 (en) * | 2011-04-29 | 2014-03-05 | B-K Medical ApS | Harmonic ultrasound imaging using synthetic aperture sequential beamforming |
CN102749622A (zh) * | 2012-07-03 | 2012-10-24 | 杭州边界电子技术有限公司 | 基于多波束测深的声速剖面及海底地形的联合反演方法 |
CN102895000A (zh) * | 2012-11-06 | 2013-01-30 | 重庆大学 | 一种基于自适应加权的双聚焦波束合成方法 |
CN103728591A (zh) * | 2013-12-17 | 2014-04-16 | 河海大学 | 一种mimo雷达近场目标高效实波束方位聚焦方法 |
CN105223544A (zh) * | 2015-08-26 | 2016-01-06 | 南京信息工程大学 | 近场线性约束最小方差自适应加权频率不变波束形成方法 |
CN105681972A (zh) * | 2016-01-14 | 2016-06-15 | 南京信息工程大学 | 线性约束最小方差对角加载的稳健频率不变波束形成方法 |
CN105891835A (zh) * | 2016-06-16 | 2016-08-24 | 北京海卓同创科技有限公司 | 一种实时动态聚焦波束形成方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106950569A (zh) | 2017-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101190134B (zh) | 超声波诊断系统中的多波束发射和接收方法及其装置 | |
JP5973349B2 (ja) | 超音波システムのビーム形成方法及び超音波撮像システム | |
US8317712B2 (en) | Retrospective dynamic transmit focusing for spatial compounding | |
CN101683274B (zh) | 偏转凸阵数字扫描变换和凸阵空间复合的方法及成像系统 | |
US8905931B2 (en) | Subject information processing apparatus | |
JPH10507101A (ja) | 超音波ビーム生成器結像システムにおけるリアルタイムで同時の適応式焦点調節の方法と装置 | |
CN108836389B (zh) | 平面波相关点相干自适应波束合成成像方法 | |
JP2014133180A (ja) | コヒーレントでない遡及的で動的な送信の焦点合わせ | |
CN106950569B (zh) | 基于序贯回归方法的多阵元合成孔径聚焦波束形成方法 | |
CN102895000A (zh) | 一种基于自适应加权的双聚焦波束合成方法 | |
CN107102298B (zh) | 基于迭代互耦校正的雷达协方差矩阵重构波束形成方法 | |
CN103969651A (zh) | 自适应声学成像方法 | |
CN102435992A (zh) | 基于广义相干系数的合成聚焦成像方法 | |
CN109521401A (zh) | 一种合成孔径成像的快速波束形成方法 | |
US20160074016A1 (en) | Transmit beamforming apparatus, receive beamforming apparatus, ultrasonic probe having the same, and beamforming method | |
Hansen et al. | Compounding in synthetic aperture imaging | |
CN110687538A (zh) | 一种基于近场聚焦的超波束形成方法 | |
CN102846336B (zh) | 超声成像延时控制方法及装置 | |
US20180038947A1 (en) | Clutter suppression in ultrasonic imaging systems | |
CN105741236B (zh) | 超声系统图像广义消旁瓣方法 | |
CN105044706B (zh) | 一种自适应波束形成方法 | |
KR101911470B1 (ko) | 평면파를 이용하여 3차원 초음파 영상을 생성하는 장치 및 방법 | |
Kortbek et al. | P2b-1 synthetic aperture focusing applied to imaging using a rotating single element transducer | |
CN104188684B (zh) | 一种自适应医疗超声成像声速优化和信号修正方法及系统 | |
CN208969233U (zh) | 一种波束数可调多波束声呐 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190329 Termination date: 20220213 |