包含波束形成代码和波束形成的原理解析,主要针对水声以及信号处理相关的设计人员有很大的帮助。
代码片段和文件信息
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 2278581 2018-11-08 23:04 多波束原理.pdf
文件 1514 2018-11-08 22:57 代码.m
----------- --------- ---------- ----- ----
2280095 2
close all;
clear;
clc;
N = 10;%number of sensors
M = 2;
c = 1500;
f = 1000;
d = c/f/2;%distance between sensors
B = 62.5;
snr = 5;
temp = 0:N-1;
theta = [60;45];
% theta = [10;60;50];
tempr = repmat(temp21);
theta = repmat(theta1N);
fs = 8192;
t = 1/fs:1/fs:0.1;
v = exp(j*2*pi*f/c*d.*cos(theta*pi/180).*tempr);
s = zeros(Mlength(temp));
s1 = sqrt(2*10^(snr/10))*exp(j*2*pi*f*t);
s2 = sqrt(2*10^(snr/10))*exp(j*2*pi*f*t);
% s3 = sqrt(2*10^(snr/10))*exp(j*2*pi*f*t);
% s = [s1;s2;s3];
%noise = rand()
s = [s1+awgn(s15‘measured‘);s2+awgn(s25‘measured‘)];
%
% for i = 1:length(v)
% x(i:) = x(i:)+rand(1length(s))+i*rand(1length(s));
% end
x = v‘*s;
r = x*x‘;
lmin = 0;
lmax = 90;
output=zeros(1(lmax-lmin+1));
output1 = output;
jl = 1;
for thets = lmin:lmax
vs = exp(j*2*pi*f/c*d*cos(thets*pi/180)*temp);
output(jl) = vs*r*vs‘;
jl = jl+1;
end
figure;
plot(lmin:lmaxabs(output))
[eigv eig1] = eig(r);
[rol col]=size(eigv);
[eigord point]= sort(diag(eig1));
%[valbpoint] = find(diag(eig1)==eigord);
u = eigv(:point);
% eigord = sort(diag(eig1));
% for l1 = 1:N-M
% [valbpoint] = find(eig1==eigord(l1));
% u(:l1) = eigv(:point);
% end
jl = 1;
P1=zeros(1length(lmin:lmax));
for thets = lmin:lmax
vs = exp(j*2*pi*f/c*d*cos(thets*pi/180)*temp);
output1(jl) = sum((vs*u(:1:(N-M))).^2);
% for n=1:N-M
% P1(jl)= P1(jl)+(vs*u(:n)).^2;
% end
jl = jl+1;
end
figure;
plot(lmin:lmax1./abs(output1))
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 2278581 2018-11-08 23:04 多波束原理.pdf
文件 1514 2018-11-08 22:57 代码.m
----------- --------- ---------- ----- ----
2280095 2
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件举报,一经查实,本站将立刻删除。
评论列表(条)