线阵自然指向性函数图Matlab仿真

Water_Sounds

已于 2022-05-05 13:34:20 修改

阅读量2.1k

点赞数 5

分类专栏：学习笔记文章标签： matlab 线性代数矩阵

于 2022-05-05 11:59:05 首次发布

本文链接：https://blog.csdn.net/weixin_45902229/article/details/124585620

版权

学习笔记专栏收录该内容

12 篇文章 7 订阅

订阅专栏

该博客探讨了不同阵列元素数量（5, 8, 12, 15, 20）下，半波间距和大于半波长间距的声纳线阵自然指向性的仿真。通过 MATLAB 代码展示了随着阵元数量变化，主瓣和栅瓣的形成及旁瓣数量的变化，强调了阵列配置对声纳性能的影响。

摘要由CSDN通过智能技术生成

%线阵的自然指向性仿真
clc;clear;close all;
% 阵数目
N = 10;
N1 = 5;
N2 = 8;
N3 = 12;
N4 = 15;
N5 = 20;
% 工作频率
f = 1000;
% 声速
c0 = 1500;
% 波长
lamda = c0/f;

d = lamda;    %布阵间距

n = 1000-1;
% xita = (-pi:2*pi/n:pi);
xita = (-pi:2*pi/n:pi);
fai = pi * d * sin(xita) / lamda;

R_xita = abs(sin(N*fai) ./ (N*sin(fai)));
R_xita1 = abs(sin(N1*fai) ./ (N1*sin(fai)));
R_xita2 = abs(sin(N2*fai) ./ (N2*sin(fai)));
R_xita3 = abs(sin(N3*fai) ./ (N3*sin(fai)));
R_xita4 = abs(sin(N4*fai) ./ (N4*sin(fai)));
R_xita5 = abs(sin(N5*fai) ./ (N5*sin(fai)));

x = sin(xita);
figure;
subplot(2,3,3);plot(x,R_xita);grid minor;
title(['N=10,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,1);plot(x,R_xita1);grid minor;
title(['N=5,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,2);plot(x,R_xita2);grid minor;
title(['N=8,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,4);plot(x,R_xita3);grid minor;
title(['N=12,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,5);plot(x,R_xita4);grid minor;
title(['N=15,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,6);plot(x,R_xita5);grid minor;
title(['N=20,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');

figure;
subplot(2,3,3);polar(xita,R_xita);grid minor;
title(['N=10,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,1);polar(xita,R_xita1);
title(['N=5,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,2);polar(xita,R_xita2);
title(['N=8,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,4);polar(xita,R_xita3);
title(['N=12,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');
subplot(2,3,5);polar(xita,R_xita4);
title(['N=15,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');     
subplot(2,3,6);polar(xita,R_xita5);
title(['N=20,lamda/d=',num2str(lamda/d)]);xlabel('sin(xita)');ylabel('|R(xita)|');

半波间距布阵时，只有主瓣，看不见栅瓣：
在这里插入图片描述
大于半波长间距布阵时，会出现栅瓣：

两个极大值之间有N-1个零点，有N-2个旁瓣