matlab仿真天线阵代码

matlab仿真天线阵代码

ID:29620867

大小:58.51 KB

页数:10页

时间:2018-12-21

matlab仿真天线阵代码_第1页
matlab仿真天线阵代码_第2页
matlab仿真天线阵代码_第3页
matlab仿真天线阵代码_第4页
matlab仿真天线阵代码_第5页
资源描述:

《matlab仿真天线阵代码》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、天线阵代码http://www.pudn.com/downloads164/sourcecode/math/detail750575.html一、clcclearallf=3e9;N1=4;N2=8;N3=12;a=pi/2;%馈电相位差i=1;%天线电流值lambda=(3e8)/f;%lambda=c/f波长d=lambda/2;beta=2.*pi/lambda;W=-2*pi:0.001:2*pi;y1=sin((N1.*W./2))./(N1.*(sin(W./2)));%归一化阵因子y1=abs(y1);r1=max(y1);y2=sin((N2.*W./2))./(N

2、2.*(sin(W./2)));%归一化阵因子y2=abs(y2);r2=max(y2);y3=sin((N3.*W./2))./(N3.*(sin(W./2)));%归一化阵因子y3=abs(y3);r3=max(y3);%归一化阵因子绘图程序,figure(1)subplot(311);plot(W,y1);gridon;%绘出N=4等幅等矩阵列的归一化阵因子xlabel('f=3GHz,N=4,d=1/2波长,a=π/2')subplot(312);plot(W,y2);gridon;%绘出N=8等幅等矩阵列的归一化阵因子xlabel('f=3GHz,N=8,d=1/2波长,

3、a=π/2')subplot(313);plot(W,y3);gridon;%绘出N=12等幅等矩阵列的归一化阵因子xlabel('f=3GHz,N=12,d=1/2波长,a=π/2')%---------------------%只有参数N改变的天线方向图t=0:0.01:2*pi;W=a+(beta.*d.*cos(t));z1=(N1/2).*(W);z2=(1/2).*(W);W1=sin(z1)./(N1.*sin(z2));%非归一化的阵因子K1K1=abs(W1);%----------------------W=a+(beta.*d.*cos(t));z3=(N2/

4、2).*(W);z4=(1/2).*(W);W2=sin(z3)./(N2.*sin(z4));%非归一化的阵因子K2K2=abs(W2);%-------------------------W=a+(beta.*d.*cos(t));z5=(N3/2).*(W);z6=(1/2).*(W);W3=sin(z5)./(N3.*sin(z6));%非归一化的阵因子K3K3=abs(W3);%--------------------绘图函数figure(2)subplot(131);polar(t,K1);xlabel('f=3GHz,N=4,d=1/2波长,a=π/2');subpl

5、ot(132);polar(t,K2);xlabel('f=3GHz,N=8,d=1/2波长,a=π/2');subplot(133);polar(t,K3);xlabel('f=3GHz,N=12,d=1/2波长,a=π/2');%----------------------%只有阵列单元方向角a改变的天线方向图N4=10;d1=lambda/4;a1=0;a2=pi/2;a3=pi/2+pi/10;A=a1+(beta.*d1.*cos(t));x1=(N4/2).*(A);x2=(1/2).*(A);A1=sin(x1)./(N4.*sin(x2));%非归一化的阵因子K4K

6、4=abs(A1);%---------------------------B=a2+(beta.*d1.*cos(t));y_1=(N4/2).*(B);y_2=(1/2).*(B);B1=sin(y_1)./(N4.*sin(y_2));%非归一化的阵因子K5K5=abs(B1);%----------------------------C=a3+(beta.*d1.*cos(t));v1=(N4/2).*(C);v2=(1/2).*(C);C1=sin(v1)./(N4.*sin(v2));%非归一化的阵因子K6K6=abs(C1);%--------------------

7、------绘图函数figure(3)subplot(131);polar(t,K4);xlabel('f=3GHz,N=10,d=1/4波长,a=0');subplot(132);polar(t,K5);xlabel('f=3GHz,N=10,d=1/4波长,a=π/2');subplot(133);polar(t,K6);xlabel('f=3GHz,N=10,d=1/4波长,a=π/2+π/10');%------------------------------

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。