MATLAB仿真天线阵代码
MATLAB仿真天线阵代码
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 、 y1=abs(y1);r1=max(y1);y2=sin((N2 、 *W 、 /2)) 、/(N2 、 y2=abs(y2);r2=max(y2);y3=sin((N3 、 *W 、 /2)) 、/(N3 、y3=abs(y3);r3=max(y3);%归一化阵因子绘图程序figure(1) subplot(311);plot(W,y1) ; grid on; % 阵因子xlabel('f=3GHz,N=4,d=1/2 波长,a= n /2') subplot(312);plot(W,y2) ; grid on; % 阵因子xlabel('f=3GHz,N=8,d=1/2 波长,a= n /2')subplot(313);plot(W,y3) ; grid on; % 化阵因子 xlabel('f=3GHz,N=12,d=1/2 波长,a= n /2')% --------------- %只有参数N 改变的天线方向图 t=0:0 、0 1 :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);% ---------------天线阵代码波长*(sin(W 、/2))); % *(sin(W 、/2))); % *(sin(W 、/2))); % 归一化阵因子 归一化阵因子 归一化阵因子 绘出N=4等幅等矩阵列的归一化 绘出N=8等幅等矩阵列的归一化 绘出N=12等幅等矩阵列的归一W=a+(beta 、 *d 、*cos(t)); z3=(N2/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 subplot(132);polar(t,K2);xlabel('f=3GHz,N=8,d=1/2 subplot(133);polar(t,K3);xlabel('f=3GHz,N=12,d=1/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);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);xlabel('f=3GHz,N=10,d=1/4 subplot(132);polar(t,K5);xlabel('f=3GHz,N=10,d=1/4 subplot(133);polar(t,K6);xlabel('f=3GHz,N=10,d=1/4 长,a= n /2+ n /10');% ------------------------------------------------波长,a= n /2'); 波长,a= n /2'); 波长,a= n 波长,a=0'); 波长,a= n /2');波%只有阵列单元间隔 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)); % 非归一化的阵因子 K7K7=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)subplot(131);polar(t,K7);xlabel('f=3GHz,N=20,d=1/4 波长,a= n /2'); subplot(132);polar(t,K8);xlabel('f=3GHz,N=20,d=1/2波长,a= n /2'); subplot(133);polar(t,K9);xlabel('f=3GHz,N=20,d=0% ------------------------------------------------ % ---------------------------- 3D- 天线方向图n_tehta = 130; % ------------- 采样视角点的仰角n_phi = 130; % --------------- 采样点的方向角[tehta,phi]=meshgrid(eps:pi 、/(n_tehta-1):pi, 、、、 %meshgrid为矩形区域的设定范围就是epfvtehtav n 0<phi<2 n0: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)); % 非归一化的阵因子 K1K_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);7 波长,a= n /2');函数M2=sin(z_3) 、/(N2 、*sin(z_4)); % 非归一化的阵因子K2K_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));z_5=(N3/2) 、*(M);z_6=(1/2) 、*(M);M3=sin(z_5) 、/(N3 、*sin(z_6)); % 非归一化的阵因子K3K_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 lightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=4,d=1/2 波长,a= n /2');figure(6) surf(X2,Y2,Z2);camlight rightlightcolorbaraxis imagerotate3D on title('f=3GHz,N=8,d=1/2 波长,a= n /2');figure(7) surf(X3,Y3,Z3)camlight rightlightcolorbaraxis imagerotate3D on title('f=3GHz,N=12,d=1/2 波长,a= n /2');% ---------------------------------% -------------- 只有阵列单元方向角 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); %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)); % 非归一化的阵因子K5 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)); % 非归一化的阵因子K6 K_6=abs(C_1);radio_6 =K_6;X6=radio_6 、*sin(tehta) 、*cos(phi);Y6=radio_6 、*sin(tehta) 、*sin(phi);Z6=radio_6 、*cos(tehta);% -------------------------figure(8)surf(X4,Y4,Z4); % 三维绘图函数surf, 采用伪彩色表示曲面的高度camlight right lightcolorbaraxis imagerotate3D on title('f=3GHz,N=10,d=1/4 波长,a=0');figure(9) surf(X5,Y5,Z5);camlight rightlight colorbar axis image rotate3D on title('f=3GHz,N=10,d=1/4 figure(10)surf(X6,Y6,Z6) camlight rightlight colorbar axis image rotate3D ontitle('f=3GHz,N=10,d=1/4 波长,a= n /2+ n /10'); %% ------------ 只有阵列单元波长,a= n /2');间隔d改变的天线方向3D图N5=20;d2=lambda/4;d3=lambda/2;d4=0 、7*lambda;a4=pi/2;D_3d=a4+(beta、*d2、*cos(t3));p_1=(N5/2) 、*(D_3d);p_2=(1/2) 、*(D_3d);D_1=sin(p_1) 、/(N5、*sin(p_2)); % 非归一化的阵因子K7K_7=abs(D_1);radio_7 =K_7;X7=radio_7 、*sin(tehta) 、*cos(phi);Y7=radio_7 、*sin(tehta) 、*sin(phi);Z7=radio_7 、*cos(tehta);% ----------------------------E_3d=a4+(beta、*d3、*cos(t3));q_1=(N5/2) 、*(E_3d);q_2=(1/2) 、*(E_3d);E_1=sin(q_1) 、/(N5 、*sin(q_2)); % 非归一化的阵因子K8K_8=abs(E_1);radio_8 =K_8;X8=radio_8 、*sin(tehta) 、*cos(phi);Y8=radio_8 、*sin(tehta) 、*sin(phi);Z8=radio_8 、*cos(tehta);% ------------------------------F_3d=a4+(beta 、*d4、*cos(t3));r_1_3d=(N5/2) 、*(F_3d);r_2_3d=(1/2) 、*(F_3d);F_1=sin(r_1_3d) 、/(N5 、*sin(r_2_3d)); % 非归一化的阵因子K9K_9=abs(F_1);radio_9 =K_9;X9=radio_9 、*sin(tehta) 、*cos(phi);Y9=radio_9 、*sin(tehta) 、*sin(phi);Z9=radio_9 、*cos(tehta);% -------------------------figure(11)surf(X7,Y7,Z7); % 三维绘图函数surf, 采用伪彩色表示曲面的高度camlight right lightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=1/4 波长,a= n /2'); figure(12)surf(X8,Y8,Z8);camlight rightlightcolorbaraxis imagerotate3D on title('f=3GHz,N=20,d=1/2 波长,a= n /2');figure(13)surf(X9,Y9,Z9)camlight rightlightcolorbaraxis imagerotate3D on title('f=3GHz,N=20,d=0 、7 波长,a= n /2');二、% ------------ 均匀直线阵列天线的应用之一: 边射阵clcclear allf=3e10; %30GH z,厘米波i=1; % 天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2、*pi/lambda;N=15;t=0:0 、01:2*pi;d1=lambda/4; %没有栅瓣效应的边射阵，即间隔d＜波长W1=beta、*d1 、*cos(t); % 定义kdcos( 方向角) z1=(N/2) 、*W1; z2=(1/2) 、*W1;F1=sin(z1) 、/(N、*sin(z2));K1=abs(F1);d2=lambda*1、5; %有栅瓣现象的边射阵,即间隔d＞波长W2=beta、*d2、*cos(t); % 定义kdcos( 方向角)z3=(N/2) 、*W2;z4=(1/2) 、*W2;F2=sin(z3) 、/(N 、*sin(z4));K2=abs(F2);figure(1)subplot(121);polar(t,K1);title('subplot(122);polar(t,K2);title(' 波长');% ------------ 均匀直线阵列天线的应用之二 : 普通端射阵clc clear all f=3e10; %30GHz, 厘米波 i=1; % 天线电流值lambda=(3e8)/f; %lambda=c/f 波长 beta=2 、 *pi/lambda;N=15;t=0:0 、01:2*pi;d1=lambda/4; % 没有栅瓣效应的普通端射阵 , 即间隔 d<1/2 波长W1=beta 、 *d1 、 *cos(t); % 定义 kdcos( 方向角 ) z1=((N/2) 、*W1)+N/2*beta*d1;z2=((1/2) 、 *W1)+1/2*beta*d1;F1=sin(z1) 、/(N 、*sin(z2));K1=abs(F1); d2=lambda*0、7; %有栅瓣现象的普通端射阵 ,即间隔 d>1/2 波长 W2=beta 、 *d2 、 *cos(t); % 定义 kdcos( 方向角 ) z3=((N/2) 、*W2)+N/2*beta*d2;z4=((1/2) 、 *W2)+1/2*beta*d2;F2=sin(z3) 、/(N 、*sin(z4));K2=abs(F2); figure(2) subplot(121);polar(t,K1);title(' 普通端射阵 f=30GHz,N=15,d=1/4 波长 ');subplot(122);polar(t,K2);title(' 普通端射阵 ( 有栅瓣) f=30GHz,N=15,d=0、7倍波长 ');四% ------------ 均匀直线阵列天线的应用之三 : 强方向性端射阵clc clear all f=3e10; %30GHz, 厘米波 i=1; % 天线电流值lambda=(3e8)/f; %lambda=c/f 波长 beta=2 、 *pi/lambda;N=15;t=0:0 、01:2*pi;d1=lambda/4; % 没有栅瓣效应的强方向性端射阵 , 即间隔 d<(1/2 波 长)*(1-1/N)W1=beta 、 *d1 、 *cos(t); % 定义 kdcos( 方向角 ) z1=((N/2) 、*W1)+N/2*(beta*d1+pi/N);z2=((1/2) 、*W1)+1/2*(beta*d1+pi/N);F1=sin(pi/2/N) 、*sin(z1) 、/(sin(z2));K1=abs(F1);d2=lambda*0、5; %有栅瓣现象的强方向性端射阵，即间隔d>1/2波长*(1-1/N) W2=beta 、 *d2、 *cos(t); % 定义 kdcos( 方向角) z3=((N/2) 、 边射阵 f=30GHz,N=15,d=1/4 波长 '); 边射阵(有栅瓣) f=30GHz,N=15,d=1 、5倍*W2)+N/2*(beta*d2+pi/N);z4=((1/2) 、*W2)+1/2*(beta*d2+pi/N);F2=sin(pi/2/N) 、*sin(z3) 、/(sin(z4));K2=abs(F2);figure(3)subplot(121);polar(t,K1);title(' 强方向性端射阵f=30GHz,N=15,d=1/4 波长');subplot(122);polar(t,K2);title(' 强方向性端射阵(有栅瓣)f=30GHz,N=15,d=0、5倍波长');。
MATLAB中V-Blast的ZF和MMSE检测算法仿真代码
MATLAB中V-Blast的ZF和MMSE检测算法仿真代码发送天线m1，接受天线m2，高斯白噪声，ZF时的仿真程序%发射天线数t某,接收天线数r某,发射矩阵长度L(帧长)t某=m1;r 某=m2;L=10000;Modulation='BPSK';EbN0=[0:5:20];B=30000;T=1/24300;%建立EbN0与SNR之间的换算关系SNR=EbN0-10某log10(T某B);%信源AA=randint(t某某L,1);%经过BPSK调制的V-Blat发射矩阵某某=zero(t某,L);fork=1:t某某(k,:)=(-1).^(A(k:t某:end)+1);end%信道传输============================================================%快衰落高斯信道HH=qrt(1/2)某(randn(r某,t某,L));%均值为0方差为1的高斯白噪声nn=qrt(1/2)某(randn(r某,L));%未叠加噪声的接收信号RR=zero(r某,L);fork=1:LR(:,k)=qrt(1/t某)某H(:,:,k)某某(:,k);end%检测%ZF========================================================= =========dip('berz');berz=[];%在不同的信噪比下计算ZF接收机误比特率berzform=SNRm%每个子信道的平均信噪比为nr的接受信号R_noiednr=10^(m/10);R_noied=awgn(R,m,1);某=[];a=zero(t某某L,1);%逐时隙对接收符号矢量进行检测，合并得到一帧发射矩阵某的估计某fort=1:Lr=R_noied(:,t);%迫零矩阵GG=pinv(H(:,:,t));y=G某r;某temp=(y>=0)-(y<0)+0;某=[某,某temp];end%从某求A的估计afork=1:t某a(k:t某:end)=(某(k:t某:end)+1)/2;end%比较A和a计算错值率temp_ber[errbit,temp_ber]=biterr(A,a);berz=[berz,temp_ber];end emilogy(EbN0,berz,'o-r'),gridon某label('Eb/N0(dB)');ylabel('误比特率');title('不同天线配置ZF系统仿真传输特性')te某t(9,0.15,'\\fontize{14}\\color{red}\\fontname{隶书}2某2')holdon发送天线m1，接受天线m2，高斯白噪声，MMSE时的仿真程序%发射天线数t某,接收天线数r某,发射矩阵长度L(帧长)t某=m1;r 某=m2;L=10000;Modulation='BPSK';EbN0=[0:5:20];B=30000;T=1/24300;%建立EbN0与SNR之间的换算关系SNR=EbN0-10某log10(T某B);%信源AA=randint(t某某L,1);%经过BPSK调制的V-Blat发射矩阵某某=zero(t某,L);fork=1:t某某(k,:)=(-1).^(A(k:t某:end)+1);end%信道传输============================================================%快衰落高斯信道HH=qrt(1/2)某(randn(r某,t某,L));%均值为0方差为1的高斯白噪声nn=qrt(1/2)某(randn(r某,L));%未叠加噪声的接收信号RR=zero(r某,L);fork=1:LR(:,k)=qrt(1/t某)某H(:,:,k)某某(:,k);end%检测%MMSE======================================================= ========dip('berz');berz=[];%在不同的信噪比下计算ZF接收机误比特率berzform=SNRm%每个子信道的平均信噪比为nr的接受信号R_noiednr=10^(m/10);R_noied=awgn(R,m,3);某=[];a=zero(t某某L,1);%逐时隙对接收符号矢量进行检测，合并得到一帧发射矩阵某的估计某fort=1:Lr=R_noied(:,t);HH=H(:,:,t);某temp=zero(t某,1);w=inv(HH'某HH+(1/nr)某eye(t某))某HH';y=w某r;某temp=(y>=0)-(y<0)+0;某=[某,某temp];end%从某求A的估计afork=1:t某a(k:t某:end)=(某(k:t某:end)+1)/2;end%比较A和a计算错值率temp_ber[errbit,temp_ber]=biterr(A,a);berz=[berz,temp_ber];end emilogy(EbN0,berz,'o-b'),gridon某label('Eb/N0(dB)');ylabel('误比特率');title('不同天线配置MMSE系统仿真传输性能')te某t(8,0.08,'\\fontize{14}\\color{blue}\\fontname{隶书}1某2')holdon。
(两个辐射单元)偶极子天线阵辐射场图matlab仿真工程电磁场小论文
偶极子天线辐射场图——MATLAB动态仿真【摘要】天线遍布于生活中的每一个角落，为了更好地学习天线，本文对直线天线的简单模型——半波偶极子进行分析。
应用MATLAB这个学习软件，对偶极子天线进行了动态仿真，通过结果分析，很好地符合书本中的实验结论，对抽象的天线理论很好地结合到了实际理解当中。
【关键字】偶极子天线元辐射场MATLAB动态仿真偶极子（dipole）定义：指相距很近的符号相反的一对电荷或“磁荷”。
在电磁学的概念里，有两种偶极子：电偶极子和磁偶极子。
电偶极子是两个分隔一段距离，电量相等，正负相反的电荷。
应用有偶极子天线。
磁偶极子是一圈封闭循环的电流，例如一个有常定电流运行的线圈，称为载流回路。
偶极子的性质可以用它的偶极矩描述。
电偶极矩由负电荷指向正电荷，大小等于正电荷量乘以正负电荷之间的距离。
磁偶极矩的方向，根据右手法则，是大拇指从载流回路的平面指出的方向，而其它拇指则指向电流运行方向，磁偶极矩的大小等于电流乘以线圈面积。
而将两个辐射单元（天线元或者阵元），也就是偶极子，按照一定方式排列的列阵天线，如果排列在直线上，称线阵天线（图一），如果排列在一个平面上，则称为面阵天线。
而这里媒质是线性的，根据线性系统的叠加定理，列阵天线的辐射场就是这两个天线元辐射场的矢量和。
并且适当地各天线元激励电流的大小和相位，就可以得到所需的辐射特性。
从而也很好地讨论由相似天线元组成的线阵天线的方向性。
偶极子天线用来发射和接收固定频率的信号。
虽然在平时的测量中都使用宽带天线，但在场地衰减和天线系数的测量中都需要使用偶极子天线。
SCHWARZBECK 偶极子天线的频率范围由30MHz~4GHz。
其中的VHAP和UHAP是一套精确偶极子天线，特别适用于场地衰减和天线系数的测量。
同时该天线为日本VCCI等标准机构指定的电波暗室和开阔场场地衰减测量等的唯一专用天线。
该天线为众多实验室所采用，作为实验室的天线标准。
垂直天线实际上是一种偶极子天线。
圆形相控阵天线matlab程序
圆形相控阵天线是一种常见的天线布局形式，它能够实现空间波束赋形和指向，适用于雷达、通信和无线电导航等领域。
在本文中，我将从圆形相控阵天线的原理、设计和matlab程序实现三个方面展开详细的讨论。
1. 圆形相控阵天线的原理圆形相控阵天线是由多个单元天线组成的阵列，每个单元天线均具有相位调控能力。
通过调节各个单元天线的相位，可以实现对不同方向的波束进行形成和指向。
在雷达和通信系统中，这种特性可以实现目标探测和数据传输的定向性。
2. 圆形相控阵天线的设计在设计圆形相控阵天线时，需要考虑阵元间的距离、阵元的数量、工作频率和辐射特性等因素。
通过合理的设计，可以实现对特定方向的高增益波束形成，并且抑制其他方向的干扰信号。
还需要考虑相控阵天线的布局方式、控制电路和信号处理算法等方面。
3. 圆形相控阵天线的matlab程序实现借助matlab等仿真工具，可以方便地进行圆形相控阵天线的仿真分析和性能评估。
在编写matlab程序时，需要考虑阵列的结构、波束赋形算法和相位控制策略。
通过调用matlab中的信号处理和天线阵列工具箱，可以快速实现圆形相控阵天线的仿真和性能优化。
总结圆形相控阵天线作为一种重要的天线结构，在雷达、通信和导航等领域具有广泛的应用前景。
通过深入理解其原理和设计方法，以及借助matlab等工具进行仿真分析，可以更好地实现对目标的探测和数据传输。
在未来的研究和应用中，圆形相控阵天线将发挥更加重要的作用。
个人观点在圆形相控阵天线的研究和应用中，我认为需要更加注重对其在多个领域的交叉应用，例如在无人机导航和卫星通信中的应用。
借助人工智能和大数据分析等技术手段，可以进一步提升圆形相控阵天线的性能和应用范围。
通过本文的阐述，我希望读者能够对圆形相控阵天线有一个更加深入和全面的了解，从而促进其在工程技术和科学研究中的应用和发展。
感谢阅读！圆形相控阵天线的原理及其工作原理圆形相控阵天线是由许多具有相位调制能力的单元天线组成的阵列。
手把手教你天线设计——用MATLAB仿真天线方向图
手把手教你天线设计——用MATLAB仿真天线方向图吴正琳天线是一种变换器，它把传输线上传播的导行波，变换成在无界媒介（通常是自由空间）中传播的电磁波，或者进行相反的变换。
在无线电设备中用来发射或接收电磁波的部件。
无线电通信、广播、电视、雷达、导航、电子对抗、遥感、射电天文等工程系统，凡是利用电磁波来传递信息的，都依靠天线来进行工作。
此外，在用电磁波传送能量方面，非信号的能量辐射也需要天线。
一般天线都具有可逆性，即同一副天线既可用作发射天线，也可用作接收天线。
同一天线作为发射或接收的基本特性参数是相同的。
这就是天线的互易定理。
天线的基本单元就是单元天线。
1、单元天线对称振子是一种经典的、迄今为止使用最广泛的天线，单个半波对称振子可简单地单独立地使用或用作为抛物面天线的馈源，也可采用多个半波对称振子组成天线阵。
两臂长度相等的振子叫做对称振子。
每臂长度为四分之一波长、全长为二分之一波长的振子，称半波对称振子。
对称振子是一种经典的、迄今为止使用最广泛的天线，单个半波对称振子可简单地单独立地使用或用作为抛物面天线的馈源，也可采用多个半波对称振子组成天线阵。
两臂长度相等的振子叫做对称振子。
每臂长度为四分之一波长、全长为二分之一波长的振子，称半波对称振子。
1.1用MATLAB画半波振子天线方向图主要是说明一下以下几点：1、在Matlab中的极坐标画图的方法：polar(theta,rho,LineSpec)；theta：极坐标坐标系0-2*pirho：满足极坐标的方程LineSpec：画出线的颜色2、在方向图的过程中如果rho不用abs(f)，在polar中只能画出正值。
也就是说这时的方向图只剩下一半。
3、半波振子天线方向图归一化方程：Matlab程序：clear alllam=1000;%波长k=2*pi./lam;L=lam/4;%天线臂长theta=0:pi/100:2*pi;f1=1./(1-cos(k*L));f2=(cos(k*L*cos(theta))-cos(k*L))./sin(theta);rho=f1*f2;polar(theta,abs(rho),'b');%极坐标系画图2、线性阵列天线2.1方向图乘积定理阵中第i 个天线单元在远区产生的电场强度为：2(,)ij i i i i ie E K If r πλθϕ-=式中，i K 为第i 个天线单元辐射场强的比例常数，i r 为第i 个天线单元至观察点的距离，(,)i f θϕ为第i 个天线单元的方向图函数，i I 为第i 个天线单元的激励电流，可以表示成为：Bji i i I a e φ-∆=式中，i a 为幅度加权系数，B φ∆为等间距线阵中，相邻单元之间的馈电相位差，亦称阵内相移值。
MATLAB仿真天线阵代码共8页
天线阵代码http://pudn/downloads164/sourcecode/math/detail750575.html一、clcclear allf=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);%归一化阵因子绘图程序,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)); %非归一化的阵因子K1K1=abs(W1);W=a+(beta.*d.*cos(t));z3=(N2/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'); subplot(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)); %非归一化的阵因子K4K4=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);%--------------------------绘图函数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');%只有阵列单元间隔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)); %非归一化的阵因子K7K7=abs(D1);E=a4+(beta.*d3.*cos(t));q1=(N5/2).*(E);q2=(1/2).*(E);E1=sin(q1)./(N5.*sin(q2)); %非归一化的阵因子K8K8=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)); %非归一化的阵因子K9K9=abs(F1);%-----------------------绘图函数figure(4)subplot(131);polar(t,K7);xlabel('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');%---------------------------------------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;N2=8;N3=12;z_1=(N1/2).*(M);z_2=(1/2).*(M);M1=sin(z_1)./(N1.*sin(z_2)); %非归一化的阵因子K1K_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)); %非归一化的阵因子K2K_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));z_5=(N3/2).*(M);z_6=(1/2).*(M);M3=sin(z_5)./(N3.*sin(z_6)); %非归一化的阵因子K3K_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 rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=4,d=1/2波长,a=π/2');figure(6)surf(X2,Y2,Z2);camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=8,d=1/2波长,a=π/2');figure(7)surf(X3,Y3,Z3)camlight rightlightcolorbaraxis imagerotate3D ontitle('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_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)); %非归一化的阵因子K4K_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)); %非归一化的阵因子K5K_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)); %非归一化的阵因子K6K_6=abs(C_1);radio_6 =K_6;X6=radio_6.*sin(tehta).*cos(phi);Y6=radio_6.*sin(tehta).*sin(phi);Z6=radio_6.*cos(tehta);figure(8)surf(X4,Y4,Z4); %三维绘图函数surf，采用伪彩色表示曲面的高度camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=0');figure(9)surf(X5,Y5,Z5);camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=π/2');figure(10)surf(X6,Y6,Z6)camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=π/2+π/10');%-------------------只有阵列单元间隔d改变的天线方向3D图N5=20;d2=lambda/4;d3=lambda/2;d4=0.7*lambda;a4=pi/2;D_3d=a4+(beta.*d2.*cos(t3));p_1=(N5/2).*(D_3d);p_2=(1/2).*(D_3d);D_1=sin(p_1)./(N5.*sin(p_2)); %非归一化的阵因子K7K_7=abs(D_1);radio_7 =K_7;X7=radio_7.*sin(tehta).*cos(phi);Y7=radio_7.*sin(tehta).*sin(phi);Z7=radio_7.*cos(tehta);E_3d=a4+(beta.*d3.*cos(t3));q_1=(N5/2).*(E_3d);q_2=(1/2).*(E_3d);E_1=sin(q_1)./(N5.*sin(q_2)); %非归一化的阵因子K8K_8=abs(E_1);radio_8 =K_8;X8=radio_8.*sin(tehta).*cos(phi);Y8=radio_8.*sin(tehta).*sin(phi);Z8=radio_8.*cos(tehta);F_3d=a4+(beta.*d4.*cos(t3));r_1_3d=(N5/2).*(F_3d);r_2_3d=(1/2).*(F_3d);F_1=sin(r_1_3d)./(N5.*sin(r_2_3d)); %非归一化的阵因子K9K_9=abs(F_1);radio_9 =K_9;X9=radio_9.*sin(tehta).*cos(phi);Y9=radio_9.*sin(tehta).*sin(phi);Z9=radio_9.*cos(tehta);figure(11)surf(X7,Y7,Z7); %三维绘图函数surf，采用伪彩色表示曲面的高度camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=1/4波长,a=π/2');figure(12)surf(X8,Y8,Z8);camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=1/2波长,a=π/2');figure(13)surf(X9,Y9,Z9)camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=0.7波长,a=π/2');二、%-----------------均匀直线阵列天线的应用之一：边射阵clcclear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的边射阵，即间隔d<波长W1=beta.*d1.*cos(t); %定义kdcos(方向角)z1=(N/2).*W1;z2=(1/2).*W1;F1=sin(z1)./(N.*sin(z2));K1=abs(F1);d2=lambda*1.5; %有栅瓣现象的边射阵，即间隔d>波长W2=beta.*d2.*cos(t); %定义kdcos(方向角)z3=(N/2).*W2;z4=(1/2).*W2;F2=sin(z3)./(N.*sin(z4));K2=abs(F2);figure(1)subplot(121);polar(t,K1);title('边射阵 f=30GHz,N=15,d=1/4波长'); subplot(122);polar(t,K2);title('边射阵（有栅瓣） f=30GHz,N=15,d=1.5倍波长');三、%-----------------均匀直线阵列天线的应用之二：普通端射阵clcclear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的普通端射阵，即间隔d<1/2波长W1=beta.*d1.*cos(t); %定义kdcos(方向角)z1=((N/2).*W1)+N/2*beta*d1;z2=((1/2).*W1)+1/2*beta*d1;F1=sin(z1)./(N.*sin(z2));K1=abs(F1);d2=lambda*0.7; %有栅瓣现象的普通端射阵，即间隔d>1/2波长W2=beta.*d2.*cos(t); %定义kdcos(方向角)z3=((N/2).*W2)+N/2*beta*d2;z4=((1/2).*W2)+1/2*beta*d2;F2=sin(z3)./(N.*sin(z4));K2=abs(F2);figure(2)subplot(121);polar(t,K1);title('普通端射阵 f=30GHz,N=15,d=1/4波长'); subplot(122);polar(t,K2);title('普通端射阵（有栅瓣） f=30GHz,N=15,d=0.7倍波长');四%-----------------均匀直线阵列天线的应用之三：强方向性端射阵clcclear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的强方向性端射阵，即间隔d<(1/2波长)*(1-1/N)W1=beta.*d1.*cos(t); %定义kdcos(方向角)z1=((N/2).*W1)+N/2*(beta*d1+pi/N);z2=((1/2).*W1)+1/2*(beta*d1+pi/N);F1=sin(pi/2/N).*sin(z1)./(sin(z2));K1=abs(F1);d2=lambda*0.5; %有栅瓣现象的强方向性端射阵，即间隔d>1/2波长*(1-1/N) W2=beta.*d2.*cos(t); %定义kdcos(方向角)z3=((N/2).*W2)+N/2*(beta*d2+pi/N);z4=((1/2).*W2)+1/2*(beta*d2+pi/N);F2=sin(pi/2/N).*sin(z3)./(sin(z4));K2=abs(F2);figure(3)subplot(121);polar(t,K1);title('强方向性端射阵 f=30GHz,N=15,d=1/4波长');subplot(122);polar(t,K2);title('强方向性端射阵（有栅瓣）f=30GHz,N=15,d=0.5倍波长');。
matlab仿真天线阵代码
matlab仿真天线阵代码天线阵代码tail750575.html一、clcclear allf=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);%归一化阵因子绘图程序,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)); %非归一化的阵因子K1K1=abs(W1);%---------------------- W=a+(beta.*d.*cos(t));z3=(N2/2).*(W);z4=(1/2).*(W);W2=sin(z3)./(N2.*sin(z4)); %非归一化的阵因子K2 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);xlabel('f=3GHz,N=4,d=1/2波长,a=π/2');subplot(132);polar(t,K2);xlabel('f=3GHz,N=8,d=1/2波长,a=π/2');,a=π/2'); subplot(133);polar(t,K3);xlabel('f=3GHz,N=12,d=1/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);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);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');%----------------------------------------------------------------------%只有阵列单元间隔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)); %非归一化的阵因子K7K7=abs(D1);%------------------------------ E=a4+(beta.*d3.*cos(t));q1=(N5/2).*(E);q2=(1/2).*(E);E1=sin(q1)./(N5.*sin(q2)); %非归一化的阵因子K8K8=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)); %非归一化的阵因子K9K9=abs(F1);%-----------------------绘图函数figure(4)subplot(131);polar(t,K7);xlabel('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');%--------------------------------------------------------------------------%---------------------------------------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;N2=8;N3=12;z_1=(N1/2).*(M);z_2=(1/2).*(M);M1=sin(z_1)./(N1.*sin(z_2)); %非归一化的阵因子K1K_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));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 rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=4,d=1/2波长,a=π/2'); figure(6)surf(X2,Y2,Z2);camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=8,d=1/2波长,a=π/2'); fi gure(7)surf(X3,Y3,Z3)camlight rightlightcolorbaraxis imagerotate3D ontitle('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_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)); %非归一化的阵因子K4K_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)); %非归一化的阵因子K5 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)); %非归一化的阵因子K6K_6=abs(C_1);radio_6 =K_6;X6=radio_6.*sin(tehta).*cos(phi); Y6=radio_6.*sin(tehta).*sin(phi);Z6=radio_6.*cos(tehta);%----------------------------------- figure(8)surf(X4,Y4,Z4); %三维绘图函数surf，采用伪彩色表示曲面的高度 camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=0');figure(9)surf(X5,Y5,Z5);camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=π/2');figure(10)surf(X6,Y6,Z6)camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=π/2+π/10');%----------------------------------------------------------------只有阵列单元间隔d改变的天线方向3D图 % N5=20;d2=lambda/4;d3=lambda/2;d4=0.7*lambda;a4=pi/2;D_3d=a4+(beta.*d2.*cos(t3)); p_1=(N5/2).*(D_3d);p_2=(1/2).*(D_3d);D_1=sin(p_1)./(N5.*sin(p_2)); %非归一化的阵因子K7K_7=abs(D_1);radio_7 =K_7;X7=radio_7.*sin(tehta).*cos(phi); Y7=radio_7.*sin(tehta).*sin(phi);Z7=radio_7.*cos(tehta);%--------------------------------------- E_3d=a4+(beta.*d3.*cos(t3)); q_1=(N5/2).*(E_3d);q_2=(1/2).*(E_3d);E_1=sin(q_1)./(N5.*sin(q_2)); %非归一化的阵因子K8K_8=abs(E_1);radio_8 =K_8;X8=radio_8.*sin(tehta).*cos(phi); Y8=radio_8.*sin(tehta).*sin(phi);Z8=radio_8.*cos(tehta);%------------------------------------------F_3d=a4+(beta.*d4.*cos(t3)); r_1_3d=(N5/2).*(F_3d);r_2_3d=(1/2).*(F_3d);F_1=sin(r_1_3d)./(N5.*sin(r_2_3d)); %非归一化的阵因子K9 K_9=abs(F_1);radio_9 =K_9;X9=radio_9.*sin(tehta).*cos(phi);Y9=radio_9.*sin(tehta).*sin(phi);Z9=radio_9.*cos(tehta); %-----------------------------------figure(11)surf(X7,Y7,Z7); %三维绘图函数surf，采用伪彩色表示曲面的高度camlight rightlightcolorbaraxis imagerotate3D on,a=π/2'); title('f=3GHz,N=20,d=1/4波长figure(12)surf(X8,Y8,Z8); camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=1/2波长,a=π/2'); figure(13)surf(X9,Y9,Z9)camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=0.7波长,a=π/2');二、%-----------------均匀直线阵列天线的应用之一:边射阵 clc clear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的边射阵，即间隔d<波长 W1=beta.*d1.*cos(t); %定义kdcos(方向角) z1=(N/2).*W1;z2=(1/2).*W1;F1=sin(z1)./(N.*sin(z2)); K1=abs(F1);d2=lambda*1.5; %有栅瓣现象的边射阵，即间隔d>波长 W2=beta.*d2.*cos(t); %定义kdcos(方向角) z3=(N/2).*W2;z4=(1/2).*W2;F2=sin(z3)./(N.*sin(z4)); K2=abs(F2);figure(1)subplot(121);polar(t,K1);title('边射阵 f=30GHz,N=15,d=1/4波长');subplot(122);polar(t,K2);title('边射阵(有栅瓣) f=30GHz,N=15,d=1.5倍波长');三、%-----------------均匀直线阵列天线的应用之二:普通端射阵 clcclear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的普通端射阵，即间隔d<1/2波长W1=beta.*d1.*cos(t); %定义kdcos(方向角)z1=((N/2).*W1)+N/2*beta*d1; z2=((1/2).*W1)+1/2*beta*d1;F1=sin(z1)./(N.*sin(z2)); K1=abs(F1);d2=lambda*0.7; %有栅瓣现象的普通端射阵，即间隔d>1/2波长W2=beta.*d2.*cos(t); %定义kdcos(方向角)z3=((N/2).*W2)+N/2*beta*d2; z4=((1/2).*W2)+1/2*beta*d2;F2=sin(z3)./(N.*sin(z4)); K2=abs(F2);figure(2)subplot(121);polar(t,K1);title('普通端射阵 f=30GHz,N=15,d=1/4波长');subplot(122);polar(t,K2);title('普通端射阵(有栅瓣)f=30GHz,N=15,d=0.7倍波长');四%-----------------均匀直线阵列天线的应用之三:强方向性端射阵 clcclear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的强方向性端射阵，即间隔d<(1/2波长)*(1-1/N)W1=beta.*d1.*cos(t); %定义kdcos(方向角)z1=((N/2).*W1)+N/2*(beta*d1+pi/N); z2=((1/2).*W1)+1/2*(beta*d1+pi/N); F1=sin(pi/2/N).*sin(z1)./(sin(z2)); K1=abs(F1);d2=lambda*0.5; %有栅瓣现象的强方向性端射阵，即间隔d>1/2波长*(1-1/N) 方向角) W2=beta.*d2.*cos(t); %定义kdcos(z3=((N/2).*W2)+N/2*(beta*d2+pi/N); z4=((1/2).*W2)+1/2*(beta*d2+pi/N); F2=sin(pi/2/N).*sin(z3)./(sin(z4)); K2=abs(F2);figure(3)subplot(121);polar(t,K1);title('强方向性端射阵 f=30GHz,N=15,d=1/4波长');subplot(122);polar(t,K2);title('强方向性端射阵(有栅瓣)f=30GHz,N=15,d=0.5倍波长');。
matlab无线通信仿真代码
matlab无线通信仿真代码
本文分享了一份基于MATLAB平台的无线通信仿真代码。
该代码包含了常用的调制解调算法、信道模型、编码解码等模块，通过调用这些模块可以实现多种无线通信系统的仿真。
具体来说，该代码可以用于仿真以下几种无线通信系统：
1. 传统的模拟调制系统，包括AM、FM、PM等模式的调制解调。
2. 数字调制系统，包括PSK、QAM、FSK等模式的调制解调。
3. 多径衰落信道模型，包括AWGN信道、瑞利衰落信道、莱斯衰落信道等。
4. 信道编码解码，包括卷积码、Turbo码、LDPC码等。
除此之外，该代码还包含了多个演示程序，可以直观地展示无线通信系统的仿真结果。
这些演示程序涉及到信道容量、误码率、信噪比等指标，可以帮助用户深入理解无线通信系统的性能。
总之，该MATLAB无线通信仿真代码简单易用，功能齐全，可用于学术研究、教学实验、工程应用等多个领域。
- 1 -。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
天线阵代码一、clcclear allf=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);%归一化阵因子绘图程序,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)); %非归一化的阵因子K1K1=abs(W1);%----------------------W=a+(beta.*d.*cos(t));z3=(N2/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'); subplot(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)); %非归一化的阵因子K4K4=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);%--------------------------绘图函数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');%----------------------------------------------------------------------%只有阵列单元间隔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)); %非归一化的阵因子K7K7=abs(D1);%------------------------------E=a4+(beta.*d3.*cos(t));q1=(N5/2).*(E);q2=(1/2).*(E);E1=sin(q1)./(N5.*sin(q2)); %非归一化的阵因子K8K8=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)); %非归一化的阵因子K9K9=abs(F1);%-----------------------绘图函数figure(4)subplot(131);polar(t,K7);xlabel('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');%--------------------------------------------------------------------------%---------------------------------------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;N2=8;N3=12;z_1=(N1/2).*(M);z_2=(1/2).*(M);M1=sin(z_1)./(N1.*sin(z_2)); %非归一化的阵因子K1K_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)); %非归一化的阵因子K2K_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));z_5=(N3/2).*(M);z_6=(1/2).*(M);M3=sin(z_5)./(N3.*sin(z_6)); %非归一化的阵因子K3K_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 rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=4,d=1/2波长,a=π/2');figure(6)surf(X2,Y2,Z2);camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=8,d=1/2波长,a=π/2');figure(7)surf(X3,Y3,Z3)camlight rightlightcolorbaraxis imagerotate3D ontitle('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_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)); %非归一化的阵因子K4K_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)); %非归一化的阵因子K5K_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)); %非归一化的阵因子K6K_6=abs(C_1);radio_6 =K_6;X6=radio_6.*sin(tehta).*cos(phi);Y6=radio_6.*sin(tehta).*sin(phi);Z6=radio_6.*cos(tehta);%-----------------------------------figure(8)surf(X4,Y4,Z4); %三维绘图函数surf，采用伪彩色表示曲面的高度camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=0');figure(9)surf(X5,Y5,Z5);camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=π/2');figure(10)surf(X6,Y6,Z6)camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=10,d=1/4波长,a=π/2+π/10');%---------------------------------------------%-------------------只有阵列单元间隔d改变的天线方向3D图N5=20;d2=lambda/4;d3=lambda/2;d4=0.7*lambda;a4=pi/2;D_3d=a4+(beta.*d2.*cos(t3));p_1=(N5/2).*(D_3d);p_2=(1/2).*(D_3d);D_1=sin(p_1)./(N5.*sin(p_2)); %非归一化的阵因子K7K_7=abs(D_1);radio_7 =K_7;X7=radio_7.*sin(tehta).*cos(phi);Y7=radio_7.*sin(tehta).*sin(phi);Z7=radio_7.*cos(tehta);%---------------------------------------E_3d=a4+(beta.*d3.*cos(t3));q_1=(N5/2).*(E_3d);q_2=(1/2).*(E_3d);E_1=sin(q_1)./(N5.*sin(q_2)); %非归一化的阵因子K8K_8=abs(E_1);radio_8 =K_8;X8=radio_8.*sin(tehta).*cos(phi);Y8=radio_8.*sin(tehta).*sin(phi);Z8=radio_8.*cos(tehta);%------------------------------------------F_3d=a4+(beta.*d4.*cos(t3));r_1_3d=(N5/2).*(F_3d);r_2_3d=(1/2).*(F_3d);F_1=sin(r_1_3d)./(N5.*sin(r_2_3d)); %非归一化的阵因子K9K_9=abs(F_1);radio_9 =K_9;X9=radio_9.*sin(tehta).*cos(phi);Y9=radio_9.*sin(tehta).*sin(phi);Z9=radio_9.*cos(tehta);%-----------------------------------figure(11)surf(X7,Y7,Z7); %三维绘图函数surf，采用伪彩色表示曲面的高度camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=1/4波长,a=π/2');figure(12)surf(X8,Y8,Z8);camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=1/2波长,a=π/2');figure(13)surf(X9,Y9,Z9)camlight rightlightcolorbaraxis imagerotate3D ontitle('f=3GHz,N=20,d=0.7波长,a=π/2');二、%-----------------均匀直线阵列天线的应用之一：边射阵clcclear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的边射阵，即间隔d<波长W1=beta.*d1.*cos(t); %定义kdcos(方向角)z1=(N/2).*W1;z2=(1/2).*W1;F1=sin(z1)./(N.*sin(z2));K1=abs(F1);d2=lambda*1.5; %有栅瓣现象的边射阵，即间隔d>波长W2=beta.*d2.*cos(t); %定义kdcos(方向角)z3=(N/2).*W2;z4=(1/2).*W2;F2=sin(z3)./(N.*sin(z4));K2=abs(F2);figure(1)subplot(121);polar(t,K1);title('边射阵 f=30GHz,N=15,d=1/4波长'); subplot(122);polar(t,K2);title('边射阵（有栅瓣） f=30GHz,N=15,d=1.5倍波长');三、%-----------------均匀直线阵列天线的应用之二：普通端射阵clcclear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的普通端射阵，即间隔d<1/2波长W1=beta.*d1.*cos(t); %定义kdcos(方向角)z1=((N/2).*W1)+N/2*beta*d1;z2=((1/2).*W1)+1/2*beta*d1;F1=sin(z1)./(N.*sin(z2));K1=abs(F1);d2=lambda*0.7; %有栅瓣现象的普通端射阵，即间隔d>1/2波长W2=beta.*d2.*cos(t); %定义kdcos(方向角)z3=((N/2).*W2)+N/2*beta*d2;z4=((1/2).*W2)+1/2*beta*d2;F2=sin(z3)./(N.*sin(z4));K2=abs(F2);figure(2)subplot(121);polar(t,K1);title('普通端射阵 f=30GHz,N=15,d=1/4波长'); subplot(122);polar(t,K2);title('普通端射阵（有栅瓣） f=30GHz,N=15,d=0.7倍波长');四%-----------------均匀直线阵列天线的应用之三：强方向性端射阵clcclear allf=3e10; %30GHz,厘米波i=1; %天线电流值lambda=(3e8)/f; %lambda=c/f 波长beta=2.*pi/lambda;N=15;t=0:0.01:2*pi;d1=lambda/4; %没有栅瓣效应的强方向性端射阵，即间隔d<(1/2波长)*(1-1/N)W1=beta.*d1.*cos(t); %定义kdcos(方向角)z1=((N/2).*W1)+N/2*(beta*d1+pi/N);z2=((1/2).*W1)+1/2*(beta*d1+pi/N);F1=sin(pi/2/N).*sin(z1)./(sin(z2));K1=abs(F1);d2=lambda*0.5; %有栅瓣现象的强方向性端射阵，即间隔d>1/2波长*(1-1/N) W2=beta.*d2.*cos(t); %定义kdcos(方向角)z3=((N/2).*W2)+N/2*(beta*d2+pi/N);z4=((1/2).*W2)+1/2*(beta*d2+pi/N);F2=sin(pi/2/N).*sin(z3)./(sin(z4));K2=abs(F2);figure(3)subplot(121);polar(t,K1);title('强方向性端射阵 f=30GHz,N=15,d=1/4波长');subplot(122);polar(t,K2);title('强方向性端射阵（有栅瓣）f=30GHz,N=15,d=0.5倍波长');。