阵列信号低副瓣处理matlab程序

共形阵列低副瓣处理matlab程序，里面包括圆阵和半圆阵。

clc
clear all
close all
%%%%%% 初始参数设置 %%%%%%
f = 3e9; % 中心频率 1300MHz
c = 3e8; % 光速
lambda = c/f; % 波长
k = 2*pi/lambda;

sensor_num=20;% 阵元数
central_ang = [0:sensor_num-1]‘*pi/（sensor_num-1）; % 阵元对应的中心角
%radius = lambda*sensor_num/4/pi; % 圆阵半径，保证阵元间距为半波长
radius=1.5*lambda;
x = radius*cos（central_ang）;
y = radius*sin（central_ang）;
position=[xy];
theta_n=central_ang*180/pi;
%%%%%%%%%%%%%%%%期望信号导向矢量%%%%%%%%%%%%%%%
theta_d=0;
v_dirction_d = [sind（theta_d）;cosd（theta_d）]; % 方向矢量%%即入射信号的单位方向矢量%
U_d= exp（1i*k*position*v_dirction_d）;%%%方向导向矢量%

%%%%%% 虚拟干扰方向导向矢量 %%%%%%%%%%%%%%%%%%%
theta_i =-90:3:90;%仰角干扰
v_dirction_i = [sind（theta_i）;cosd（theta_i）]; % 干扰方向矢量
U_i= exp（1i*k*position*v_dirction_i）; % 干扰方向导向矢量
%%加入阵元方向图

%%%%%%%%%%%%%%%%扫描信号导向矢量%%%%%%%%%%%%%%%
theta_scan=-90:1:90;  % 一维方位角扫描范围（与x轴正向夹角）
v_dirction = [sind（theta_scan）;cosd（theta_scan）]; % 方向矢量
U_scan=exp（1i*k*position*v_dirction）;


%%%%%%%%% 期望波形图 %%%%%%%%%
SSL = 40; % 期望副瓣电平;
P_opt = 10^（-SSL/2

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-11-15 09:19  UCA\
     文件        3449  2013-10-21 09:07  UCA\S_UCA.asv
     文件        2685  2013-10-22 15:17  UCA\S_UCA.m
     文件        3567  2013-10-23 20:25  UCA\S_UCAf.m
     文件        2769  2013-10-21 21:22  UCA\UCA.asv
     文件        2776  2013-10-23 19:43  UCA\UCA.m
     文件        3506  2013-10-22 14:49  UCA\UCAf.asv
     文件        3525  2013-10-23 20:18  UCA\UCAf.m
     文件        4673  2013-10-10 10:49  UCA\UCA_d.asv
     文件        4860  2013-10-11 09:57  UCA\UCA_d.m
     文件        4076  2013-10-15 08:57  UCA\UCA_zhuban.asv
     文件        4077  2013-10-15 09:00  UCA\UCA_zhuban.m
     文件        1514  2013-10-13 16:02  UCA\UCA_zishiying.m