系列二:利用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