array antenna

(Local) In[]:=
λ=1;(*[m]*)
(Local) In[]:=
k=2Pi/λ;
(Local) In[]:=
N_elem=5;(*thenumberoftheelements*)

broadside array antenna

(Local) In[]:=
D_bs[d_,δ_,r_,θ_,φ_]=Block[{ψ},​​ψ=k*d*Cos[θ]+δ;​​Cos[Cos[θ]*Pi/2]/Sin[θ]*Sin[N_elem*ψ/2]/(N_elem*Sin[ψ/2])​​]
(Local) Out[]=
1
5
Cos
1
2
πCos[θ]Csc[θ]Csc
1
2
(δ+2dπCos[θ])Sin
5
2
(δ+2dπCos[θ])
(Local) In[]:=
Block[{d=0.4λ,θ0=Pi/3},​​Plot[D_bs[d,-k*d*Cos[θ0],1,θ,0],{θ,0,Pi},AxesLabel{Style["θ",16],Style["D_bs[θ]",16]}]​​]
(Local) Out[]=
(Local) In[]:=
Block[{d=0.4λ,θ0=Pi/3,D_bs_cartesian},​​D_bs_cartesian[x_,y_,z_]=TransformedField["Spherical"->"Cartesian",D_bs[d,-k*d*Cos[θ0],r,θ,φ],{r,θ,φ}->{x,y,z}];​​SliceDensityPlot3D[D_bs_cartesian[x,y,z],"CenterSphere",{x,-1,1},{y,-1,1},{z,-1,1},AxesLabel{Style["x",16],Style["y",16],Style["z",16]}]​​]
(Local) Out[]=

end-fire array antenna

(Local) In[]:=
D_ef[d_,δ_,r_,θ_,φ_]=Block[{ψ},​​ψ=k*d*Sin[θ]Cos[φ]+δ;​​Cos[Cos[θ]*Pi/2]/Sin[θ]*Sin[N_elem*ψ/2]/(N_elem*Sin[ψ/2])​​]
(Local) Out[]=
1
5
Cos
1
2
πCos[θ]Csc[θ]Csc
1
2
(δ+2dπCos[φ]Sin[θ])Sin
5
2
(δ+2dπCos[φ]Sin[θ])
(Local) In[]:=
Block[{d=0.4λ,θ0=Pi/2,φ0=0},​​Plot[D_ef[d,-k*d*Sin[θ0]Cos[φ0],1,Pi/2,φ],{φ,0,2Pi},AxesLabel{Style["φ",16],Style["D_ef[π/2,φ]",16]}]​​]
(Local) Out[]=
(Local) In[]:=
Block[{d=0.4λ,θ0=Pi/2,φ0=0,D_ef_cartesian},​​D_ef_cartesian[x_,y_,z_]=TransformedField["Spherical"->"Cartesian",D_ef[d,-k*d*Sin[θ0]Cos[φ0],r,θ,φ],{r,θ,φ}->{x,y,z}];​​SliceDensityPlot3D[D_ef_cartesian[x,y,z],"CenterSphere",{x,-1,1},{y,-1,1},{z,-1,1},AxesLabel{Style["x",16],Style["y",16],Style["z",16]}]​​]
(Local) Out[]=

Playground