单脉冲测角-雷达信号处理-半阵法测角-matlab仿真

% 半阵侧向法，和差波束形成 % Author:huasir 2023.11.1 @Beijing N = 8; %阵元数 theta = (-90:0.1:90); %观测角度范围 theta0 =0; %阵列波束指向的方向 theta0 = theta0*pi/180; theta = theta*pi/180; d_lembda = 1/2; %阵元间距比波长 a = exp(j*2*pi*d_lembda*(0:N-1)'*sin(theta0)); %指向theta0处的波束导向矢量 a_phis = exp(j*2*pi*d_lembda*(0:N-1)'*sin(theta)); %观测角度下的波束导向矢量 wSum = a'; %和波束的加权权向量 wDiff = [ones(1,N/2).*(-1),ones(1,N/2)].*(a'); %差波束的加权权向量 ySum = abs(wSum*a_phis); yDiff = abs(wDiff*a_phis); %win = taylorwin(N); % 泰勒窗 %win = chebwin(N,30); %切比雪夫窗，抑制旁瓣 % ww = w.*win; %加窗操作 % y1 = abs(ww'*a); %加窗后的方向图 %% 绘制半阵法和差波束图 figure; plot(theta*180/pi,20*log10(ySum/max(ySum)),'linewidth',1); hold on; plot(theta*180/pi,20*log10(yDiff/max(yDiff)),'--','linewidth',1); xlabel('方位角/°');ylabel('归一化功率和波束方向图/dB'); legend('\Delta 和波束','\Sigma 差波束'); axis tight;% axistight 使得图形框图靠近数据 grid on; %添加栅格线 ylim([-50, 0]); % 为了限制y值范围，使得图像显示的更加合理 title(sprintf('和差波束图，阵元数：%d，波束方向：%.0f°',N,theta0)); %% 绘制单脉冲比曲线 ySum = (wSum*a_phis); %和波束 yDiff = (wDiff*a_phis); %差波束 angleDiscr =imag(yDiff./ySum); %鉴角曲线，取其虚部 thetaDis = (-6:0.1:6); %鉴角曲线取-6°~6° angleDiscr = angleDiscr(841:961); %取-6°-6° figure; plot((-6:0.1:6),angleDiscr,'linewidth',1); axis tight;% axistight 使得图形框图靠近数据 grid on; %添加栅格线 title(sprintf('鉴角曲线，阵元数：%d，波束方向：%.0f°',N,theta0)); %% 取某一个角度，进行验证 angle = 2 %某一个角度 realAnglD = angleDiscr((angle+6)/0.1+1) %该角度对应的单脉冲比 computeAngle = realAnglD/(pi*N*1/2*d_lembda); %由单脉冲比计算出来的角度 computeAngle = computeAngle*180/pi %角度转变为弧度 angleError = computeAngle - angle %计算误差