clc
clear
syms AF1;
%**************************************************************************
% AF阵列因子
%**************************************************************************
f = 40000; 中心频率为40%的超声传感器kHz
omega = 2*pi*f; %ω角速度
rho_air = 1.205; 摄氏度时的空气密度
c0 = 344; 当摄氏度时，超声波在空气中的传播速度
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');