模拟平面波入射和相关法测向过程

clear all;clc;close all;
%% 初始化
d=1;%阵间距
N=10;%传感器数量
l=d*(N-1);%阵长
c=343;%声速
fs=2000;%采样率
t=0:1/fs:2-1/fs;
s=zeros(length(t),N);%每通道一列
%% 模拟任意方向来的平面波
f=100;%信号频率
w=f2pi;%角频率
lamda=c/f;%波长(保证两倍阵间距小于波长
k=w/c;%波数
theta=56;%1~180 正前方90度
for i=1:N
phi=(i-1)cos(theta/180pi)dk;
s(:,i)=exp(1j*(w*t-phi));% -kx
end

%% 检测平面波方向
xiang=zeros(1,180);
g0=(N-1)d/cfs;%最大偏移量
g0=ceil(g0);%向上取整
dth=0.5;%角度分辨率
n=1;
for m=1:dth:180
theta=m/180pi;
% su=ones(1,32);
% ss=su
s;%求和
ss1=zeros(length(t)-2*g0,1);

dd=dcos(theta);
tao=dd/c;
gg=tao
fs;
for i=1:N
g=round((i-1)*gg);
ss1=ss1+s(g+g0+1:g+length(t)-g0,i);
end
xiang(n)=max(xcorr(s(g0+1:length(t)-g0,1),ss1));
n=n+1;

end
plot(1:dth:180,xiang)
[~,m]=max(xiang);
m=m*dth;
fprintf(‘平面波来自%.1f’,m);

  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值