【18X18圆孔径相控阵天线】自编matlab程序

j921y

嗯，程序没有考虑元因子方向性，认为是由全向天线组成的阵列，由方向图乘法原理可以得有向阵元阵列方向图。
程序我没有加注释，大家看起来可能有点费劲，这主要是防止大家不加吸收直接拿来做论文，呵呵。。
有兴趣可以研究下，程序里面主要是分为3个部分：
1、参数设置(更改Nr可以得到更大阵元数量的圆孔径阵列；d是阵元间最小间距)
2、圆孔径阵面形成(里面是圆阵排布的核心,方法有很多中,个人认为这样编程显得简单)
3、方向图的计算及显示
如果大家要看3D的方向图，需要更改phi的范围，将phi改成数组。显示命令用mesh即可。
使用mesh要将theta、phi由向量变成矩阵。
比如在程序第一部分，theta(范围0~pi/2)是1x900,phi(范围0~2*pi)是1x3600长度的向量，需将theta变成900x3600，phi也变成900x3600的矩阵形式。
最简洁的matlab语句是：
theta = theta';
THETA = theta(:,ones(1,3600));
PHI = phi(ones(900,1),:);
mesh(THETA,PHI,AF);
zlim([-80 0]);
xlabel('\theta :rad');
ylabel('\phi :rad');
zlabel('旁瓣电平：dB');
就可以了显示3D方向图了。(这几行代码是加在画图部分的)
下图是低精度的示意图，不是按照900x3600仿真的


当然你也可以使用以下语句来绘制3D方向图:
首先在第4、29行更改:
phi = 0:0.01:2*pi;
theta = 0:0.01:pi/2;  
最后绘图用:
U = sin(theta)'*cos(phi);
V = sin(theta)'*sin(phi);
mesh(U,V,AF);
zlim([-80 0]);
xlabel('U');
ylabel('V');
zlabel('旁瓣电平：dB');

如果你的计算机CPU不是4核的推荐使用低精度的theta、phi画3D方向图，
900x3600的数据量会计算一定时间，尤其是当Nr取值较大的时候。

j921y

soccer11 发表于 2012-5-10 14:56

                               
登录/注册后可看大图

感谢楼主分享原创！

不客气，有空我再传些上来

j921y

00d44 发表于 2012-5-11 08:58

                               
登录/注册后可看大图

支持下原创，欢迎更多原创性文章交流

OK,下了一些论坛的资料,自己也传点,互相交流~~

haoyangang

厉害嘎

shangshan

虽然现在还没发下载，还是很感谢楼主分享。

地狱猎狼人

太贵了，下不起啊

xjfl64

可惜不能下载

dzf2003

好东西啊

supmet

擦还要技术分的啊，没有啊，怎么弄到技术分？

hguang77

thank you for sharing