BELLHOP 关于Actup冲激响应的绘制

很多同学在BELLHOP Actup原指南中找不到 BELLHOPdata.m 文件,同时存在程序代码的缺失等问题，在利用Actup时无法很好的工作。这里提供了本人编写的Actup使用说明书点击打开链接，以及相关代码点击打开链接。

这里在Matlab代码中再次详细说明BELLHOP 关于Actup冲激响应的绘制。

（1）首先利用Actup得到幅度和相位.arr文件。详细配置方式请参考本人编写的使用说明书。

（2）利用工具箱自带的read_arrivals_asc(filename)函数（可以通过搜索的方式找到该函数）获取幅度、相位、角度等数据。如下列代码所示。

clear
clc
filename = 'Test_00100.arr';
Minimum_range=100  %（接收水听器的水平方向上接收范围最小值，m）----R
Maximum_range=1000 %（接收水听器的水平方向上接收范围最大值，m）---RB [ amp1, delay, SrcAngle, RcvrAngle, NumTopBnc, NumBotBnc, narrmat, Pos ]... = read_arrivals_asc(  filename );

（3）根据需要绘制单位冲激响应。注意这里amp和delay的维度。这里由于在水平方向设置了50个接收机，因此这里的amp和delay是2维数据。图中选取最大的m值，表示选取第50号接收机的幅度和时延。如果设置水平和垂直方向均有接收机，则会使得amp和delay是3维数据，这里就需要根据实际情况取自己设定位置的接收机的幅值和时延。（这里可以参考本人编写的说明书第8章内容）

%%单位冲激响应
[m,n]=size(amp1);
amp = abs(amp1); %取模
x = delay(m,:); %获取第50个接收机的时延和幅值
y = amp(m,:);
figure(1)
stem(x,y)
grid on
xlabel('相对时延/s')
ylabel('幅度')
title('单位冲激响应')

（4）绘制归一化冲激响应。

%%归一化冲激响应
Amp_Delay = [x;y];
Amp_Delay(:,all(Amp_Delay==0,1))=[]; %去掉0值
Amp_Delay=sortrows(Amp_Delay',1);  %按照时延从小到大排序
normDelay = Amp_Delay(:,1)-Amp_Delay(1,1);%归一化时延
normAmp = Amp_Delay(:,2)/Amp_Delay(1,2);%归一化幅度
figure(2)
stem(normDelay,normAmp,'^')
grid on
xlabel('相对时延/s')
ylabel('归一化幅度')
title('归一化冲激响应')

（5）绘制距离时延图。

figure(3)
mum=1:m;
ReRange = Minimum_range+(Maximum_range-Minimum_range)/m*mum;
for i=1:min(narrmat)
plot(delay(:,i),ReRange,'o')
hold on
end
hold off
grid on
colorbar
xlabel('时延(sec)')
ylabel('Range(m)')
title(filename)

BELLHOP 关于Actup冲激响应的绘制相关推荐

matlab08a调节字体大小,MATLAB低通滤波器的设计代码
本次试验中设计了不同窗函数的低通滤波器的比较环节.各个窗函数设计低通滤波器的思路是相似的,只是其中的窗函数是不同的. 单位冲激响应及其幅频响应及其代码单位冲激响应的绘制思路比较简单,就是将一条直线用 ...
梅尔倒谱系数MFCC由浅入深
MFCC梅尔倒谱系数(Mel-scale Frequency Cepstral Coefficients) 在语音识别(Speech Recognition)和话者识别(Speaker Recogni ...
【Matlab】海底声学模拟（Bellhop）以及滤波器的设计
一.设计要求某单波束测深仪最大测量水深为300米,请根据<水声学原理>和<数字信号处理>相关知识,仿真设计该单波束测深仪的数字信号处理系统(包括模拟滤波器参数.采样频率.量化 ...
Bellhop 从入门到上手
文章目录前言一.Bellhop 简介二.Bellhop 结构 1.输入文件 2.输出文件三.Bellhop 环境文件 1.OPTIONS1 2.OPTIONS2 3.OPTIONS3 4.OP ...
模拟滤波器的单位冲激响应+单位阶跃响应+斜坡响应+抛物线响应matlab实现(转载+整理)
有人会纠结: 单位冲激响应和单位脉冲响应到底是不是一个东西? 参考zhihu回答: 该回答已被删除 - 知乎自动控制原理第四版(胡寿松)-P77 ------------------------- ...
matlab 冲激响应不变法,matlab实验七冲激响应不变法IIR数字滤波器设计.doc
(数字信号处理)实验报告实验名称实验七冲激响应不变法IIR数字滤波器设计实验时间年月日专业班级学号姓名成绩教师评语: 一.实验目的 1.掌握构成一个频率响应与给定的滤波特 ...
IIR滤波器设计之冲激响应不变法与双线性变换法
前言 IIR(Infinite Impulse Response)数字滤波器是指无限长冲激响应数字滤波器,因其具有良好的幅频响应特性,被广泛应用于通信.控制.生物医学.振动分析.雷达和声呐等领域. I ...
matlab冲激响应不变发,matlab实验七冲激响应不变法IIR数字滤波器设计
matlab实验七冲激响应不变法IIR数字滤波器设计 (数字信号处理)实验报告实验名称实验七冲激响应不变法IIR数字滤波器设计实验时间年月日专业班级学号姓名成绩教师评语: 一 ...
MATLAB：零状态响应(lsim(连续)；filter(离散))、冲激响应(impulse或impz)和阶跃响应(step)、卷积(conv)
例1:求系统y''(t)+2y'(t)+100y(t)=10x(t)的零状态响应,已知x(t)=sin(2pt)u(t).[连续时间系统零状态响应] 例2:求系统y''(t)+2y'(t)+100y( ...

BELLHOP 关于Actup冲激响应的绘制

BELLHOP 关于Actup冲激响应的绘制相关推荐

最新文章

热门文章