原创力文档- 1、本文档共1页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
clear
M=64;c=340;%阵元总数
theta=-2:0.05:2;%入射仰角
N=length(theta);
phi=-2:0.05:2;%入射方向角
R=1;%阵列直径
omiga=4000;%入射角频率
dphi=2*pi/M;%相邻阵元间隔角度
for k=1:N
for m=1:N
for n=0:M-1
t(n+1)=omiga*(R*cos(dphi*n)*(sin(theta(k))*cos(phi(m))-1/2)+R*sin(dphi*n)*(sin(theta(k)*sin(phi(m))-1/2)))/c;%第i号阵元相对于参考阵元的时延
A(n+1)=exp(-1i*t(n+1));
end
B(k,m)=1/M*(sum(A));
end
end
[xx,yy]=meshgrid(theta,phi);
surf(xx,yy,abs(B));
title(4000Hz指向性);
xlabel(入射仰角 x/m);
ylabel(入射方向角 y/m);
figure
contour(B,8);
title(4000Hzμè????);
xlabel(入射仰角 x/m);
ylabel(入射方向角 y/m);
文档评论（0）