传统波束形成的算法实现「建议收藏」

传统波束形成的算法实现「建议收藏」最近学习了传统波束形成(CBF)的原理,尝试着写出识别一个单声源的波束形成程序。下面按照程序说明一下。1、初始化设置一些常数,例如抽样频率,所要计算的频率,时间步等。clearall;closeall;clc;%—————-初始化—————-%c=1500;%声速cfs=10000;%抽样频率fsT=0.1…

大家好,又见面了,我是你们的朋友全栈君。

最近学习了传统波束形成(CBF)的原理,尝试着写出识别一个单声源的波束形成程序。下面按照程序说明一下。

1、初始化

设置一些常数,例如抽样频率,所要计算的频率,时间步等。

clear all;
close all;
clc;
%---------初始化常量----------%
c = 334;      % 声速c
fs = 1000;    % 抽样频率fs
T = 0.1;       % ??
t = 0:1/fs:T;  % 时间 [0,0.1]
L = length(t); % 时间长度
f = 500;       % 感兴趣的频率
w = 2*pi*f;    % 角频率
k = w/c;       % 波数 k

2、阵元设置

下面是阵元部分,包括阵元个数,间距设置。

M = 18;        % 阵元个数

接下来是阵元坐标,声源位置,这里设置的是(12,10,12)点,同时设置阵元的高斯白噪声。

%---------各阵元坐标---------%
M = 18;        % 阵元个数
% Nmid = 12;      % 参考点
% d = 3;         % 阵元间距
% m = (0:1:M-1) 
yi = zeros(M,1); % 生成一个M*1维的零矩阵
zi = [0;3;6;9;12;15;18;21;24;12;12;12;12;12;12;12;12;12];
xi = [12;12;12;12;12;12;12;12;12;0;3;6;9;12;15;18;21;24];     
%xi = xi.'      % 列向量 m*d 阵元数*阵元间距


figure(1)
plot(xi,zi,'r*');
title('十字形麦克风阵列')
%--------- 声源位置----------%
x1 = 12;  
y1 = 10;
z1 = 12;  %声源位置 (12,10,12) x,z为水平面


x2 = 12;
y2 = 0;
z2 = 12;


Ric1 = sqrt((x1-xi).^2+(y1-yi).^2+(z1-zi).^2); % 声源到各阵元的距离
Ric2 = sqrt((x1-x2).^2+(y1-y2).^2+(z1-z2).^2);
Rn1 = Ric1 - Ric2;  %声源至各阵元与参考阵元的声程差矢量


s1 = cos(2*w*t);    % 参考阵元接收到的矢量
   
Am = 10^(-1);  % 振幅
n1 = Am * (randn(M,L)+j*randn(M,L));  % 各阵元高斯白噪声
p1 = zeros(M,L);

3、延迟求和

整个程序最关键的部分,延迟求和,同时得到各阵元接收的声压信号矩阵。以及协方差矩阵,这个还有疑问,要把论文读懂来理解。

%--------------------------各阵元的延迟求和--------------------------------%
for k1 = 1:M
    p1(k1,:) = Ric2/Ric1(k1)*s1.*exp(-j*w*Rn1(k1)/c);
    % 接收到的信号
end
p = p1+n1;  % 各阵元接收的声压信号矩阵
R = p*p'/L; % 接收数据的自协方差矩阵  A.'是一般转置,A'是共轭转置

4、扫描整个声源平面

我们设置步长为0.1,扫描范围是20×20的平面,双重for循环得到M*1矢量矩阵,最后得到交叉谱矩阵(cross spectrum matrix)。由DSP理论,这个就是声音的功率。

%-------扫描范围------%
step_x = 0.1;  % 步长设置为0.1
step_z = 0.1;
y = y1;
x = (9:step_x:15);  % 扫描范围 9-15
z = (9:step_z:15);


for k1=1:length(z)
    for k2=1:length(x)
        Ri = sqrt((x(k2)-xi).^2+(y-yi).^2+(z(k1)-zi).^2);
        Ri2 = sqrt((x(k2)-x2).^2+(y-y2).^2+(z(k1)-z2).^2);
        % 该扫描点到各阵元的聚焦距离矢量
        Rn = Ri-Ri2;
        % 扫描点到各阵元与参考阵元的程差矢量
        
        b = exp(-j*w*Rn/c); % 声压聚焦方向矢量
        Pcbf(k1,k2) = abs(b'*R*b); % CSM
    end
end

5、归一化处理

归一化处理的程序

%--------------------------------------归一化------------------------------%
for k1 = 1:length(z);
    pp(k1) = max(Pcbf(k1,:)); % Pcbf 的第k1行的最大元素的值
end

Pcbf = Pcbf/max(pp);  % 所有元素除以其最大值 归一化幅度

6、作图 

观察得到的结果

%-------------------------------作图展示-----------------------------------%
figure(2)
surf(x,z,Pcbf);
xlabel('x(m)'),ylabel('z(m)')
title('CBF三维单声源图')
colorbar


figure(3)
pcolor(x,z,Pcbf);
shading interp;
xlabel('x(m)');
ylabel('z(m)');
title('CBF单声源图')
colorbar

7、结果

十字型阵列最终得到的结果效果并不理想,没有达到一个点声源的理想结果。

传统波束形成的算法实现「建议收藏」

可以看出下图三维点声源的效果并不理想。

传统波束形成的算法实现「建议收藏」

传统波束形成的算法实现「建议收藏」

接下来是矩形阵列的结果:

传统波束形成的算法实现「建议收藏」

可以看出下图三维点声源的效果还可以。

传统波束形成的算法实现「建议收藏」

传统波束形成的算法实现「建议收藏」

论文中给出的结果,很明显的点声源分布。

传统波束形成的算法实现「建议收藏」

接下来是继续熟悉CBF算法,提高可扩展性,同时学习声源反卷积成像算法(DAMAS),下一篇进行算法总结。

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/148167.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • 软件测试用例常用七大方法

    软件测试用例常用七大方法第一:测试用例格式包括十大特点用例编号测试项测试标题用例属性重要级别:高中低预置条件测试输入操作步骤预期结果实际结果第二:等价类1,等价类定义2,等价类划分3,等价类划分规则4,进行等价类用例设计5,案例加以说明第三:边界值1,边界值的三点2,边界值应用场景3,边界值方法应用步骤第四:判定…

  • 前端语音信号处理

    前端语音信号处理1、语音活动检测语音活动检测(VoiceActivityDetection,VAD)用于检测出语音信号的起始位置,分离出语音段和非语音(静音或噪声)段。VAD算法大致分为三类:基于阈值的VAD、基于分类器的VAD和基于模型的VAD。基于阈值的VAD是通过提取时域(短时能量、短时过零率等)或频域(MFCC、谱熵等)特征,通过合理的设置门限,达到区分语音和非语音的目的;基于分类…

  • IBinder对象在进程间传递的形式(一)

    IBinder对象在进程间传递的形式(一)

  • linux 内网文件传输工具_Linux服务器局域网(内网)快速传输文件[通俗易懂]

    linux 内网文件传输工具_Linux服务器局域网(内网)快速传输文件[通俗易懂]最近新网一台服务器因为Thinkphp5.0.0~5.0.23版本远程代码执行漏洞导致网站被挂马,查杀病毒过程中为了不影响业务使用,临时开了一台新服务器,那么怎么快速在这两台服务器之间快速传送文件呢?在这里我们使用scp命令在两台内网服务器间快速传送文件,scp语法如下:scp–r远程用户名@IP地址:文件名本地目录比如我们要从内网IP为10.10.1.1的服务器中的www/wwwroo…

  • mysql数据库double类型_timestamp是什么数据类型

    mysql数据库double类型_timestamp是什么数据类型1、整型MySQL数据类型含义(有符号)tinyint(m)1个字节范围(-128~127)smallint(m)2个字节范围(-32768~32767)mediumint(m)3个字节范围(-8388608~8388607)int(m)4个字节范围(-2147483648~2147483647)bigint(m)8个字节范围(+-9.22*10的18次方)取值范围如果加了un…

  • C++并发实战19:lock free编程

    C++并发实战19:lock free编程涉及到并行/并发计算时,通常都会想到加锁,加锁可以保护共享的数据,不过也会存在一些问题:1.由于临界区无法并发运行,进入临界区就需要等待,加锁使得效率的降低。多核CPU也不能发挥全部马力2.在复杂的情况下,很容易造成死锁,并发进程、线程之间无止境的互相等待。3.在中断/信号处理函数中不能加锁,给并发处理带来困难。4.加锁影响实时性,等待时间不确定5.优先级反转,优先级

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号