Field II 超声声场仿真(2): 接收灵敏度(原理&代码)

系列二：利用Field II仿真计算接收灵敏度

系列一关注的是发射声场，接收跟发射类似，只不过用灵敏度receive sensitivity代替声场强度transmit pressure进行衡量。

1. 定义发射和接收孔径

接收孔径receive aperture与发射孔径emit aperture定义一致。

emit_aperture = xdc_focused_array (N_elements, element_width, element_height, kerf, Rfocus, N_sub_x, N_sub_y, focus);
receive_aperture = xdc_focused_array (N_elements, element_width, element_height, kerf, Rfocus, N_sub_x, N_sub_y, focus);

2. 设置换能器脉冲响应和激励脉冲

发射和接收孔径的脉冲响应一致。

t_ir = -2/f0:1/fs:2/f0;
Bw = 0.6;
impulse_response=gauspuls(t_ir,f0,Bw);
set_sampling(fs);
xdc_impulse (emit_aperture, impulse_response);
xdc_impulse (receive_aperture, impulse_response);

但不同的是，接收孔径不需要定义激励脉冲，只需要对发射孔径定义激励脉冲。

ex_periods = 1.5;
t_ex=(0:1/fs:ex_periods/f0);
excitation=square(2*pi*f0*t_ex);
xdc_excitation (emit_aperture, excitation);

3. 计算xz平面的接收灵敏度

使用`calc_hp'函数计算接收孔径xz平面的灵敏度。

[simData, startTime] = calc_hp(receive_aperture, measurement_points);

receive sensitivity

代码请加QQ：2971319104