系列二:利用Field II仿真計算接收靈敏度
系列一關注的是發(fā)射聲場,接收跟發(fā)射類似,只不過用靈敏度receive sensitivity代替聲場強度transmit pressure進行衡量。
1. 定義發(fā)射和接收孔徑
接收孔徑receive aperture與發(fā)射孔徑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. 設置換能器脈沖響應和激勵脈沖
發(fā)射和接收孔徑的脈沖響應一致。
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);
但不同的是,接收孔徑不需要定義激勵脈沖,只需要對發(fā)射孔徑定義激勵脈沖。
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