shengqiangfafzh1.m

矢量传感器的声强法定向程序。利用矢量传感器输出生涯和阵速的相关性

clear all
close all
clc
m=20;n=4096;
% Generate a uniform distribution of random numbers on a specified interval [a,b]. 
a1 = -pi/2; b1 = pi/2;
alfa1 = a1 + (b1-a1) * rand(1,m);
% Generate a uniform distribution of random numbers on a specified interval [a,b]. 
a2 = 0; b2 = 2*pi;
theta1 = a2 + (b2-a2) * rand(1,m);
%generates an m-by-n matrix of white Gaussian noise. 
ro0=1000;c=1500;

pni=zeros(m,n);
for i=1:m
    pni(i,:)=wgn(1,n,0);
end
pn=sum(pni,1)/m;

vxi=zeros(m,n);
for i=1:m
    vxi(i,:)=pni(i,:)*cos(alfa1(i))*cos(theta1(i))/(ro0*c);
end
vnx=sum(vxi,1)/m;

vyi=zeros(m,n);
for i=1:m
    vyi(i,:)=pni(i,:)*cos(alfa1(i))*sin(theta1(i))/(ro0*c);
end
vny=sum(vxi,1)/m;

vzi=zeros(m,n);
for i=1:m
    vzi(i,:)=pni(i,:)*sin(alfa1(i))/(ro0*c);
end
vnz=sum(vxi,1)/m;
% Generate cw numbers . 
f=500;fs=2000;
t=0:1/fs:(n-1)/fs;
ps=sin(2*pi*f*t);
plot(ps);
alfa=5*pi/180;theta=50*pi/180;
vsx=ps*cos(alfa)*cos(theta)/(ro0*c);
vsy=ps*cos(alfa)*sin(theta)/(ro0*c);
vsz=ps*sin(theta)/(ro0*c);
% Generate diejia signal
p=ps+pn;
vx=vsx+vnx;
vy=vsy+vny;
vz=vsz+vnz;
Ix=mean(p.*vx);
Iy=mean(p.*vy);
thetaq=atan(Iy/Ix)*180/pi
Ix1=mean(ps.*vsx);
Iy1=mean(ps.*vsy);
thetaq1=atan(Iy1/Ix1)*180/pi