自适应波束形成与Matlab程序代码注解VIP

1.均匀线性阵列样式(1)matlab程序ClcClear allClose allimag=sqrt(-1)；Element _ num=32阵列元素%为8d _ lam da=1/2；%阵列元素间距d与波长lamda的关系theta=linspace(-pi/2，pi/2，200)；theta 0=0；% oll波方向w=exp(imag * 2 * pi * d _ lamd * sin(thet a 0)*03360 element _ num-1)；for j=1: length(theta(theta)a=exp(imag * 2 * pi * d _ lamd a* sin(theta(j)*03360 element _ num-1)；p(j)=w * a；Endpatternag=ABS(p)；patternnmagnorm=patternnmag/max(max(patternnmag)；patter ndb=20 * log 10(patter nmag)；patterndbnorm=20 * log 10(patternnmagnorm)；Figure(1)Plot (thet a * 180/pi，patternnmag)；栅格开；Xlabel(theta/radian)Ylabel(amplitude/dB)title(num 2 str(element _ num)阵列元素均匀线性样式，波的方向为num2str(theta0*180/pi)。)；霍尔德温；Figure(2)plot(theta，pattern bnorm，r)；栅格开；Xlabel(theta/radian)Ylabel(amplitude/dB)title(num 2 str(element _ num)阵列元素均匀线性样式，波的方向为num2str(theta0*180/pi)。)；axis(-1 . 5 1 . 5-50)；(2)模拟结果A.波的方向为0非规格化规范化B.波的方向是45非规格化规范化C.随着阵列元素数量的增加，梁宽度变窄，分辨率提高，仿真图表如下：非规格化规范化非规格化规范化梁宽度与DOA方向和阵列元素数的关系(1)matlab程序ClcClear allClose allimag=sqrt(-1)；Element _ num1=16Element _ num2=128Element _ num3=1024兰达=0.1；d=0.5 * lamda；theta=0:0.53336990；for j=1: length(theta(theta)fai(j)=theta(j)* pi/180-asin(sin(theta(j)* pi/180)-lamda/(element _ num 1)psi(j)=theta(j)* pi/180-asin(sin(theta(j)* pi/180)-lamda/(element _ num 2)beta(j)=theta(j)* pi/180-asin(sin(theta(j)* pi/180)-lamda/(element _ num 3)EndFigureplot(theta、pier、r、theta、psi、b、theta、beta、g)；栅格开；xlabel(theta)；Ylabel(径向宽度)；Title(梁宽度与最大波方向和阵列元素数之间的关系)；Legend(N=16，N=128，N=1024)；(2)模拟结果结果3.阵列元素间距时，栅格瓣出现，空间变暗(1)模拟结果非规格化规范化4.与时域过滤类似，天线模式是具有最佳权重的傅立叶变换(1)matlab程序ClcClear allClose allimag=sqrt(-1)；Element _ num=32source _ num=1；D _ lambda=0.5theta=linspace(-pi/2，pi/2，200)；theta 0=0；w=exp(imag * 2 * pi * d _ lambda * sin(thet a 0)*03360 element _ num-1)；for j=1: length(theta(theta)a=exp(imag * 2 * pi * d _ lamb da * sin(theta(j)*03360 element _ num-1)；p(j)=w * a；Endpatternag=ABS(p)；patternnmagnorm=patternnmag/max(max(patternnmag)；patter ndb=20 * log 10(patter nmag)；patterndbnorm=20 * log 10(patternnmagnorm)；Figure(1)Subplot(1，2，1)；plot(theta，pattern bnorm)；栅格开；xlabel(theta/radian)；ylabel(amplitude/dB)；axis(-2 . 0 2 . 0-50)；Subplot(1，2，2)；Pfft=fftshift(fft(w，256)；pfftmag=ABS(pfft)；pfftmagnorm=pfftmag/max(max(pfftmag)；PfftdB=20 * log 10(pfftmagnorm)；pfftdbnorm=20 * log 10(pfftmagnorm)；Plot (linspace (-pi/2、pi/2、256)、pfftdbnorm)；栅格开；xlabel(theta/radian)；ylabel(FFT _ amplitude/dB)；axis(-2 . 0 2 . 0-50)；(2)模拟结果5.最大信噪比参考模式和功率谱(1) MATLAB程序ClcClear allClose allimag=sqrt(-1)；element _ num=8；阵列元素%为8D _ lambda=0.5%间隔是半波长theta=-90:0.53336990；%扫描范围theta 0=0；% oll波位置theta 1=20；%碰撞方向L=512%取样点For i=1:Lmp0=10 * randn(1)；am P1=200 * randn(1)；ampn=1；S (:I)=a mp0 * exp(imag * 2 * pi * 0.5 * sin(thet a 0 * pi/180)*03360 elementJ (:I)=amp 1 * exp(imag * 2 * pi * 0.5 * sin(thet a 1 * pi/180)*03360 elementN (:I)=ampn * exp (randn (element _ num，1) imag * rannn (element _ num，1)EndRS=1/L * s * s；%信号自相关矩阵rnj=1/L *(j * j n * n)；%干扰噪声的自相关矩阵V，D=eig(Rs，Rnj)；%(Rs，Rnj)的广义特征值和特征向量D，I=sort(diag(D)；%唯一向量排序Wopt=V(:I(8)；%最佳权重向量for j=1: length(theta(theta)a=exp(imag * 2 * pi * d _ lamb da * sin(thet a(j)* pi/180)*03360 element _ numf(j)=wop t * a；p(j)=a* RS * a* rnj * a；Endf=20 * log 10(ABS(f)/max(ABS(f)；p=20 * log 10(ABS(p)/max(ABS(p)；Subplot(1，2，1)plot(theta，f)；栅格开；霍尔德温；plot(theta 0，-80:0，)；plot(theta 1，-80:0，)；xlabel(theta/0)；ylabel(F in dB)；Title(max-SNR地图)；axis(-90 90-80)；霍尔德温；Subplot(1，2，2)；plot(theta、p、r)；栅格开；xlabel(theta/0)；Ylabel(电源in dB)；Title(max-SNR功率谱)；栅格开；axis(-90 90-80)；(2)模拟结果卸下Asc旁瓣- MSE说明(1) MATLAB程序ClcClose allClear allimag=sqrt(-1)；M=32%第二天线数D _ lambda=0.5%阵列元素间距theta 0=-30；% oll波方向theta 1=60；%碰撞方向L=512%采样单位数S=zeros(1，512)；%预先划分区域For ii=1:Lmp0=1 * randn(1)；%信号的振幅是随机生成的，并确保信号之间没有关联性am P1=200 * randn(1)；ampn=1；Jam (:ii)=amp 1 * exp(imag * 2 * pi * d _ lambda * sin(thet a 1 * pi/180)*000%干扰噪音s(ii)=a mp0 * exp(imag * 2 * pi * d _ lamb da * sin(thet a 0 * pi/180)amp 1 * exp()%接收信号(信号干扰噪声)s0(ii)=a mp0 * exp(imag * 2 * pi * d _ lamb da * sin(thet a 0 * pi/180)；Endrx=1/L * jam * jam；r _ xd=1/L * jam * s；wopt=pinv(Rx)* r _ xd；delta=s0-(s-Wopt * jam)；Delta 1=ABS (delta .2)-(delta)(2)；theta=linspace(-pi/2，pi/2，200)；For jj=1:length(theta)a=exp(imag * 2 * pi * d _ lamb da * sin(theta(jj)*0336m-1)；f(jj)=wop t * a；Endf=20 * log 10(ABS(f)/max(ABS(f)；Figure(1)Plot(theta*180/pi，F)；栅格开；霍尔德温；plot(theta 0，-5033600，)；plot(theta 1，-5033600，)；xlabel(theta/)；ylabel(F/dB)；Title(基于MSE标准的模式)；axis(-90 90-50)；(2)模拟结果7.线性约束最小方差(LCMV)准则(1) MATLAB程序ClcClear allClose allimag=sqrt(-1)；element _ num=8；阵列元素%数D _ lambda=0.5%阵列元素间距和波长的关系theta=-90:0.53336990；%搜索范围theta 0=0；% 3个信号源的波方向theta 1=30；theta 2=60；L=512%采样单位数For i=1:Lmp0=10 * randn(1)；amp 1=100 * randen(1)；amp 2=10 * randn(1)；Ampn=10X(:I)=a mp0 * Exp(imag * 2 * pi * d _ lambda * sin(thet a 0