vsfz.m

单个矢量传感器的波束仿真

clc
clear all    
close all                          %二维矢量水听器，平面波条件，p(t)=x(t) vx(t)=cos(thetas)*x(t) vy(t)=sin(thetas)*x(t)%x(t)
thetas=0:0.1:2*pi;                 %为声压波形, thetas为声传播的水平方位角%上式中,某一时刻t的速度vx(t),vy(t)只与方位角
vx=cos(thetas);                    %thetas有关,可画出关系图如下%
subplot(1,2,1);
polar(thetas,abs(vx));
vy=sin(thetas);
subplot(1,2,2);
polar(thetas,abs(vy));

pvx=(1+cos(thetas)).*cos(thetas);  %由公式(p(t)+vx(t))*vx(t)=(1+cos(thetas)).*cos(thetas)*x(t)^2，某一时刻t的组合量(p(t)+vx(t))*vx(t)只与方位角thetas有关
pvy=(1+sin(thetas)).*sin(thetas);  %由公式(p(t)+vy(t))*vy(t)=(1+sin(thetas)).*sin(thetas)*x(t)^2，某一时刻t的组合量(p(t)+vx(t))*vx(t)只与方位角thetas有关
figure
subplot(1,2,1);
polar(thetas,abs(pvx));
subplot(1,2,2);
polar(thetas,abs(pvy));

pvxf=(1+cos(thetas)).^2;            %由公式(p(t)+vx(t))^2=(1+cos(thetas))^2*x(t)^2,某一时刻t的组合量(p(t)+vx(t))^2只与方位角thetas有关
pvyf=(1+sin(thetas)).^2;            %由公式(p(t)+vy(t))^2=(1+sin(thetas))^2*x(t)^2,某一时刻t的组合量(p(t)+vy(t))^2只与方位角thetas有关
figure
subplot(1,2,1);
polar(thetas,abs(pvxf));
subplot(1,2,2);
polar(thetas,abs(pvyf));

fai=pi/4;                   %下面实现波束在水平面内的电子旋转,设旋转角为fai
pvx=(1+cos(thetas-fai)).*cos(thetas-fai);  
pvy=(1+sin(thetas-fai)).*sin(thetas-fai); 
figure
subplot(1,2,1);
polar(thetas,abs(pvx));
subplot(1,2,2);
polar(thetas,abs(pvy));

vx=cos(thetas-fai);            % 通过改变fai值，可以在全方位进行角度搜索，当fai=thetas时,即搜索到目标方位,vx有极大值,同时vy有极小值. 
vy=sin(thetas-fai);
figure
subplot(1,2,1);
polar(thetas,abs(vx));
subplot(1,2,2);
polar(thetas,abs(vy));