1.一种基于截断修正SL0算法的MIMO雷达参数估计方法，其特征在于：包括如下步骤：步骤一：建立MIMO雷达接收信号向量模型；

步骤二：通过对MIMO雷达感知矩阵的奇异值进行二次截断和一次修正处理来改善感知矩阵的病态性，并利用SVD反变换获得非病态感知矩阵；

步骤三：利用SL0算法对MIMO雷达目标参数进行估计，并在计算初值和梯度投影值的步骤中，将病态感知矩阵的伪逆由上面获得的非病态感知矩阵的伪逆代替；

步骤四：最终输出求解后的目标场景向量α，根据最终输出求解后的目标场景向量α中的非零元素位置来确定各个目标的角度、距离和多普勒频移值；

所述步骤一包括如下步骤：

1a：建立MIMO雷达的接收信号为 式中，

(·)*和(·)T分别表示矩阵或矢量的复共轭转置和转置；P是雷达目标探测区域内所划分的距离单元个数，K是雷达目标探测区域内所划分的角度单元个数，H是雷达目标探测区域内所划分的多普勒单元个数；E是加性噪声；αp,k,h，p＝1,...,P；k＝1,...,K；h＝1,...,H，表示目标的复散射系数，其中；JP是用来描述从不同距离单元返回的信号所用的转移矩阵；aT(θk)和aR(θk)分别是发射阵列和接收阵列的导向向量，其中，θk表示第k，1≤k≤K，个角度单元对应的目标角度； 是补零发射信号矩阵，其中，ωh为第h个多普勒单元对应的目标多普勒频移，h＝1,...,H，Sd(ωh)为对应ωh的多普勒频移信号矩阵，为维度(P-1)×Mt的零矩阵，Mt为发射阵元个数，N为发射信号长度；

1b：将上式简化成向量形式，定义

式中， vec(·)表示矩阵向量化运算，A是感知矩阵；α是目标场景向量；N为发射信号长度；Mr为接收阵元的数量；

1c：得到接收信号的向量模型为y＝Aα+e，式中，e＝vec(E)；

所述步骤二包括如下步骤：

2a：对矩阵A作奇异值分解:A＝UΣV*，其中 为由奇异值σi，i＝1,2,…,Mr(N+P-1)，构成的对角矩阵，U、V分别为奇异值对应的左奇异矩阵、右奇异矩阵，A是感知矩阵；

2b：选取一次截断门限 保留大于等于 的奇异值，剔除小于的奇异值，即

2c：对保留的大于等于 的奇异值再次截断并修正，选择二次截断门限将保留的奇异值分成大于等于 和小于 的两部分，利用Tikhonov正则化准则对大于等于 的奇异值进行修正，将小于 的奇异值修正为与 接近的常数值，即q1和q2为修正参数；

2d：定义U1和V1分别为修正后的奇异值所对应的左奇异矩阵和右奇异矩阵，其中U1＝U(:,1:t1)由矩阵U的前t1列组成的矩阵，V1＝V(1:t1,:)由矩阵V的前t1行组成的矩阵，Σ1是修正后的奇异值σ'm，1≤m≤t1，构成的对角矩阵；利用SVD反变换获得非病态感知矩阵A1＝U1Σ1V1*，矩阵A1的伪逆

2.根据权利要求1所述的基于截断修正SL0算法的MIMO雷达参数估计方法，其特征在于：所述步骤三包括如下步骤：

3a：初始化:

(1)令初值 内循环次数为L、步长为μ，其中L,μ＞0，U1和V1分别为修正后的奇异值所对应的左奇异矩阵和右奇异矩阵；y为向量形式表示的MIMO雷达接收信号；

Σ1是修正后的奇异值σ'm，1≤m≤t1，构成的对角矩阵；

(2)选取一个序列[δ1,δ2,...,δJ]，且δ1＞δ2＞...＞δJ，其中，

3b：算法迭代：

forj＝1,2,...,J

(1)令δ＝δj

(2)利用最速下降法求解Fδ(α)的最小值，并将其投影到可行集上；

初始化：

forl＝1,...,L

(a)设φ＝[α1exp(-α12/2δ2),...,αZexp(-αZ2/2δ2)]T，其中Z＝P·K·H，P、K和H分别是在雷达目标探测区域中所划分的距离单元个数,角度单元个数，以及多普勒单元个数；

(b)

(c)

(3)令

3c：当 或δj≤δJ时，转至步骤四，否则，令j＝j+1，δj＝ρδj-1，返回步骤

3b；其中η为一个正数，用于判断算法是否终止，ρ为衰减因子，0＜ρ＜1。

3.根据权利要求2所述的基于截断修正SL0算法的MIMO雷达参数估计方法，其特征在于：所述η＝0.01，ρ＝0.8。

4.根据权利要求2所述的基于截断修正SL0算法的MIMO雷达参数估计方法，其特征在于：所述内循环次数L＝50，步长μ＝2。

5.根据权利要求2所述的基于截断修正SL0算法的MIMO雷达参数估计方法，其特征在于：所述修正参数q1＝115，q2＝95，C＝25。