clc

clear

syms AF1;

%**************************************************************************

%       AF阵列因子

%**************************************************************************

f = 40000;                   %超声波传感器的中心频率为40kHz

omega = 2*pi*f;              %ω角速度

rho_air = 1.205;             %20摄氏度时的空气密度

c0 = 344;                    %20摄氏度时，超声波在空气中的传播速度

p_min = 6.32;                %测试声压

r0 = 0.3;                    %测试点距离

ua = 0.7826;                 %振动速度幅值

k = 730.6029;                %k=w/c0波数

lamda = 0.0086;              %λ=0.0086米

ra=0.0055;                  %超声传感器半径

slipNum=201;

wphi1=60;

theta=-pi/2:0.01:pi/2;      %仰角的范围变化

theta2=-wphi1:2.*wphi1/(slipNum-1):wphi1;

%单个传感器的指向性函数

AF1 = besselj(1 ,k*ra*sin(theta) )./ (k*ra*sin(theta));   %单个超声波传感器的指向性图

AF1=2*abs(AF1);

%%%----------------单个发射阵元的方向指向性图--------------------------------

figure(1)

polar(theta,AF1,'b');