MATLAB仿真天线阵代码17页Word下载.docx
- 文档编号:14592339
- 上传时间:2022-10-23
- 格式:DOCX
- 页数:12
- 大小:16.89KB
MATLAB仿真天线阵代码17页Word下载.docx
《MATLAB仿真天线阵代码17页Word下载.docx》由会员分享,可在线阅读,更多相关《MATLAB仿真天线阵代码17页Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
2*pi;
y1=sin((N1.*W./2))./(N1.*(sin(W./2)));
%归一化阵因子
y1=abs(y1);
r1=max(y1);
y2=sin((N2.*W./2))./(N2.*(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);
%绘出N=8等幅等矩阵列的归一化阵因子
f=3GHz,N=8,d=1/2波长,a=π/2'
subplot(313);
plot(W,y3);
%绘出N=12等幅等矩阵列的归一化阵因子
f=3GHz,N=12,d=1/2波长,a=π/2'
%只有参数N改变的天线方向图
t=0:
0.01:
W=a+(beta.*d.*cos(t));
z1=(N1/2).*(W);
z2=(1/2).*(W);
W1=sin(z1)./(N1.*sin(z2));
%非归一化的阵因子K1
K1=abs(W1);
z3=(N2/2).*(W);
z4=(1/2).*(W);
W2=sin(z3)./(N2.*sin(z4));
%非归一化的阵因子K2
K2=abs(W2);
z5=(N3/2).*(W);
z6=(1/2).*(W);
W3=sin(z5)./(N3.*sin(z6));
%非归一化的阵因子K3
K3=abs(W3);
%--------------------绘图函数
figure
(2)
subplot(131);
polar(t,K1);
xlabel('
);
subplot(132);
polar(t,K2);
subplot(133);
polar(t,K3);
%只有阵列单元方向角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));
%非归一化的阵因子K4
K4=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));
%非归一化的阵因子K5
K5=abs(B1);
C=a3+(beta.*d1.*cos(t));
v1=(N4/2).*(C);
v2=(1/2).*(C);
C1=sin(v1)./(N4.*sin(v2));
%非归一化的阵因子K6
K6=abs(C1);
%--------------------------绘图函数
figure(3)
polar(t,K4);
f=3GHz,N=10,d=1/4波长,a=0'
polar(t,K5);
f=3GHz,N=10,d=1/4波长,a=π/2'
polar(t,K6);
f=3GHz,N=10,d=1/4波长,a=π/2+π/10'
%只有阵列单元间隔d改变的天线方向图
N5=20;
d2=lambda/4;
d3=lambda/2;
d4=0.7*lambda;
a4=pi/2;
D=a4+(beta.*d2.*cos(t));
p1=(N5/2).*(D);
p2=(1/2).*(D);
D1=sin(p1)./(N5.*sin(p2));
%非归一化的阵因子K7
K7=abs(D1);
E=a4+(beta.*d3.*cos(t));
q1=(N5/2).*(E);
q2=(1/2).*(E);
E1=sin(q1)./(N5.*sin(q2));
%非归一化的阵因子K8
K8=abs(E1);
F=a4+(beta.*d4.*cos(t));
r_1=(N5/2).*(F);
r_2=(1/2).*(F);
F1=sin(r_1)./(N5.*sin(r_2));
%非归一化的阵因子K9
K9=abs(F1);
%-----------------------绘图函数
figure(4)
polar(t,K7);
f=3GHz,N=20,d=1/4波长,a=π/2'
polar(t,K8);
f=3GHz,N=20,d=1/2波长,a=π/2'
polar(t,K9);
f=3GHz,N=20,d=0.7波长,a=π/2'
%---------------------------------------3D-天线方向图
n_tehta=130;
%--------------------采样视角点的仰角
n_phi=130;
%--------------------采样点的方向角
[tehta,phi]=meshgrid(eps:
pi./(n_tehta-1):
pi,...%meshgrid函数为矩形区域的设定范围是epf<
tehta<
π0<
phi<
2π
0:
2*pi./(n_phi-1):
2*pi);
t3=tehta;
%-------------只有参数N改变的天线方向3D图
M=a+(beta.*d.*cos(t3));
%----N1=4;
z_1=(N1/2).*(M);
z_2=(1/2).*(M);
M1=sin(z_1)./(N1.*sin(z_2));
K_1=abs(M1);
radio_1=K_1;
X1=radio_1.*sin(tehta).*cos(phi);
Y1=radio_1.*sin(tehta).*sin(phi);
Z1=radio_1.*cos(tehta);
z_3=(N2/2).*(M);
z_4=(1/2).*(M);
M2=sin(z_3)./(N2.*sin(z_4));
K_2=abs(M2);
radio_2=K_2;
X2=radio_2.*sin(tehta).*cos(phi);
Y2=radio_2.*sin(tehta).*sin(phi);
Z2=radio_2.*cos(tehta);
z_5=(N3/2).*(M);
z_6=(1/2).*(M);
M3=sin(z_5)./(N3.*sin(z_6));
K_3=abs(M3);
radio_3=K_3;
X3=radio_3.*sin(tehta).*cos(phi);
Y3=radio_3.*sin(tehta).*sin(phi);
Z3=radio_3.*cos(tehta);
%------------------------------------3D绘图函数
figure(5)
surf(X1,Y1,Z1);
%三维绘图函数surf,采用伪彩色表示曲面的高度
camlightright
light
colorbar
axisimage
rotate3Don
title('
figure(6)
surf(X2,Y2,Z2);
figure(7)
surf(X3,Y3,Z3)
%--------------------只有阵列单元方向角a改变的天线方向图
A_3d=a1+(beta.*d1.*cos(t3));
x_1=(N4/2).*(A_3d);
x_2=(1/2).*(A_3d);
A_1=sin(x_1)./(N4.*sin(x_2));
K_4=abs(A_1);
radio_4=K_4;
X4=radio_4.*sin(tehta).*cos(phi);
Y4=radio_4.*sin(tehta).*sin(phi);
Z4=radio_4.*cos(tehta);
B_3d=a2+(beta.*d1.*cos(t3));
y_1_3d=(N4/2).*(B_3d);
y_2_3d=(1/2).*(B_3d);
B_1=sin(y_1_3d)./(N4.*sin(y_2_3d));
K_5=abs(B_1);
radio_5=K_5;
X5=radio_5.*sin(tehta).*cos(phi);
Y5=radio_5.*sin(tehta).*sin(phi);
Z5=radio_5.*cos(tehta);
C_3d=a3+(beta.*d1.*cos(t3));
v_1=(N4/2).*(C_3d);
v_2=(1/2).*(C_3d);
C_1=sin(v_1)./(N4.*sin(v_2));
K_6=abs(C_1);
ra
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 仿真 天线阵 代码 17