系列二:利用Field II仿真計(jì)算接收靈敏度
系列一關(guān)注的是發(fā)射聲場(chǎng)歇拆,接收跟發(fā)射類(lèi)似吨岭,只不過(guò)用靈敏度receive sensitivity
代替聲場(chǎng)強(qiáng)度transmit pressure
進(jìn)行衡量甩鳄。
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. 設(shè)置換能器脈沖響應(yīng)和激勵(lì)脈沖
發(fā)射和接收孔徑的脈沖響應(yīng)一致在辆。
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);
但不同的是册踩,接收孔徑不需要定義激勵(lì)脈沖泳姐,只需要對(duì)發(fā)射孔徑定義激勵(lì)脈沖。
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. 計(jì)算xz平面的接收靈敏度
使用`calc_hp'函數(shù)計(jì)算接收孔徑xz平面的靈敏度暂吉。
[simData, startTime] = calc_hp(receive_aperture, measurement_points);
receive sensitivity
代碼請(qǐng)加QQ:2971319104