含苞欲放的意思-

2021年1月21日发(作者：江注)
天线阵代码

http://pudn/downloads164/sourcecode/math/

一、

clc

clear all

f=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))./(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);

第
1
页

归一化阵因子


归一化阵因子


归一化阵因子


%
归一化阵因子绘图程序
,

figure(1)

subplot(311);plot(W,y1) grid on; %
绘出
N=4
等幅等矩阵列的
归一化阵因子


xlabel('f=3GHz,N=4,d=1/2
波长,a=π/2')

subplot(312);plot(W,y2) grid on; %
绘出
N=8
等幅等矩阵列的
归一化阵因子


xlabel('f=3GHz,N=8,d=1/2
波长,a=π/2')

subplot(313);plot(W,y3) grid on; %
绘出
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)); %
非归一化的阵因子
K1

K1=abs(W1);

W=a+(beta.*d.*cos(t));

z3=(N2/2).*(W);

z4=(1/2).*(W);

W2=sin(z3)./(N2.*sin(z4)); %
非归一化的阵因子
K2

第
2
页

K2=abs(W2);

W=a+(beta.*d.*cos(t));

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);xl abel('f=3GHz,N=4,d=1/2
波长,a=π/2');

sub plot(132);polar(t,K2);xlabel('f=3GHz,N=8,d=1/2
波长,a=π/2');

subplot(133);polar(t,K3);xlab el('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)); %
非归一化的阵因子
K4

K4=abs(A1);

B=a2+(beta.*d1.*cos(t));

y_1=(N4/2).*(B);

y_2=(1/2).*(B);

第
3
页

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)

subplot(131);polar(t,K4);xl abel('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 ');

%
只有阵列单元间隔
d
改变的天线方向图


N5=20;d2=lambda/4;d3=lambda/2;d4=0.7*lamb da;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));

第
4
页

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)

subplot(131);polar(t,K7);xl abel('f=3GHz,N=20,d=1/4
波
长,a=π/2');
subplot(132);polar(t,K8);xlabel('f=3GHz,N=20,d=1 /2
波
长,a=π/2');

subplot(133);polar( t,K9);xlabel('f=3GHz,N=20,d=0.7
波
长,a=π/2');

%---------------------------------------3 D-
天线方向图


n_tehta = 130; %--------------------
采样视角点的仰角


n_phi = 130; %--------------------
采样点的方向角


[te hta,phi]=meshgrid(eps:pi./(n_tehta-1):pi,... %meshgrid
函
数为矩形区域的设定范围是
epf
第
5
页

0:2*pi./(n_phi-1):2*pi)

t3=tehta;

%-------------
只有参数
N
改变的天线方向
3D
图


M=a+(beta.*d.*cos(t3)); %----N1=4;N2=8;N3=12;

z_1=(N1/2).*(M);

z_2=(1/2).*(M);

M1=sin(z_1)./(N1.*sin(z_2)); %
非归一化的阵因子
K1

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);

M=a+(beta.*d.*cos(t3));

z_3=(N2/2).*(M);

z_4=(1/2).*(M);

M2=sin(z_3)./(N2.*sin(z_4)); %
非归一化的阵因子
K2

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);

M=a+(beta.*d.*cos(t3));

第
6
页

z_5=(N3/2).*(M);

z_6=(1/2).*(M);

M3=sin(z_5)./(N3.*sin(z_6)); %
非归一化的阵因子
K3

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
，采用伪彩色表示曲面的高度


camlight right

light

colorbar

axis image

rotate3D on

title('f=3GHz,N=4,d=1/2
波长,a=π/2');

figure(6)

surf(X2,Y2,Z2);

camlight right

light

colorbar

第
7
页

axis image

rotate3D on

title('f=3GHz,N=8,d=1/2
波长,a=π/2');

figure(7)

surf(X3,Y3,Z3)

camlight right

light

colorbar

axis image

rotate3D on

title('f=3GHz,N=12,d=1/2
波长,a=π/2');
< br>%--------------------
只有阵列单元方向角
a
改变的天 线方向图


N4=10;d1=lambda/4;a1=0;a2=pi/2;a3=pi/2+pi/10;

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)); %
非归一化的阵因子
K4

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);

第
8
页

含苞欲放的意思-

含苞欲放的意思-

含苞欲放的意思-

含苞欲放的意思-

含苞欲放的意思-

含苞欲放的意思-

含苞欲放的意思-

含苞欲放的意思-