恒模约束波形设计（Constant Modulus Waveform Design）第2章 MIMO雷达波形设计基础

目录2.1 MIMO雷达系统模型2.1.1 发射-接收信号模型与阵列几何配置2.1.2 信号依赖干扰建模2.1.3 扩展目标与点目标检测模型2.2 波形设计目标函数2.2.1 波束赋形设计准则2.2.2 输出SINR最大化框架2.2.3 CRLB最小化方法2.3 硬件约束与工程实现2.3.1 高功率放大器饱和特性与恒模需求2.3.2 峰均比约束与包络恒定化技术2.3.3 相似性约束的p-范数形式第二部分结构化伪代码算法1恒模波束赋形联合优化算法算法2SINR最大化联合收发设计算法3CRLB最小化线性近端方法算法4PAPR约束联合优化ADMM框架算法5扩展目标检测波形优化2.1 MIMO雷达系统模型2.1.1 发射-接收信号模型与阵列几何配置考虑由Nt​ 个发射天线和Nr​ 个接收天线组成的MIMO雷达系统阵列采用均匀线阵ULA几何配置阵元间距为半波长。发射信号矩阵记为S[s1​,s2​,…,sNt​​]T∈CNt​×L 其中sn​[sn​(1),sn​(2),…,sn​(L)]T 表示第n 个发射阵元在L 个快拍上的发射序列。在远场条件下目标方向θ 处的发射导向矢量为a(θ)[1,ejπsinθ,…,ejπ(Nt​−1)sinθ]T目标处的入射信号可表示为aT(θ)S 。对于位于方向θt​ 的扩展目标其目标冲激响应TIR记为t∈CQ 相应的TIR矩阵构造为T∑i1Q​t(i)Ji−1​ 其中Ji​ 为(QL−1)×L 移位矩阵其元素满足Ji​(m,n){1,0,​m−nim−ni​接收端信号矩阵表示为Yt​b(θt​)aT(θt​)STT其中b(θ) 为接收导向矢量。接收信号经匹配滤波后输出信号矢量包含目标回波、杂波与噪声分量。2.1.2 信号依赖干扰建模信号依赖干扰Signal-Dependent Clutter的统计特性与发射波形直接相关其协方差矩阵构造为Rc​(S)k1∑Kc​​σk2​(b(θk​)aT(θk​)STkT​)(b(θk​)aT(θk​)STkT​)Hσn2​I其中Kc​ 为杂波单元数σk2​ 表示第k 个杂波单元的反射功率θk​ 为杂波入射方向Tk​ 为对应杂波单元的响应矩阵。杂波与目标回波通过相同的传播路径到达接收机形成与发射波形耦合的干扰环境。在空时自适应处理STAP框架下信号依赖杂波呈现角度-多普勒耦合特性其空时协方差矩阵维度为Nr​Np​×Nr​Np​ Np​ 为脉冲数。2.1.3 扩展目标与点目标检测模型扩展目标检测模型假设目标占据多个距离单元其冲激响应服从复高斯分布t∼CN(0,Σt​) 。检测问题转化为二元假设检验{H0​:yyc​nH1​:yyt​yc​n​其中yc​ 为信号依赖杂波分量n 为加性高斯白噪声。对于点目标检测目标回波简化为yt​αb(θt​)aT(θt​)S α 为目标复幅度。扩展目标检测需联合处理多个散射中心的相干积累其检测统计量构造为ΛyHR−1(S)yH0​≷H1​​​η其中R(S)Rc​(S)σn2​I 为总干扰协方差矩阵。2.2 波形设计目标函数2.2.1 波束赋形设计准则波束赋形Beam Pattern Matching设计准则旨在使合成波束方向图逼近理想方向图d(θ) 。发射协方差矩阵Rs​SSH/L 决定方向图特性其合成波束功率为P(θ)aH(θ)Rs​a(θ)优化目标函数构造为加权最小二乘形式Smin​m1∑M​wm​​aH(θm​)Rs​a(θm​)−d(θm​)​2其中wm​ 为角度权重系数{θm​}m1M​ 为采样角度集合。通过引入辅助变量XRs​ 问题转化为半定规划SDP形式利用凸松弛技术求解。对于恒模约束需将解映射至单位圆上保持幅度恒定∣sn​(l)∣c 。2.2.2 输出SINR最大化框架输出信干噪比SINR最大化是波形设计的核心准则。接收滤波器w 与发射波形S 联合优化输出SINR表达式为SINRwHRc​(S)wσn2​wHwwHRt​(S)w​其中Rt​(S)(b(θt​)aT(θt​)STT)(b(θt​)aT(θt​)STT)H 为目标协方差矩阵。固定S 时最优接收滤波器为MVDR波束形成器wopt​bH(θt​)R−1(S)b(θt​)R−1(S)b(θt​)​联合优化问题表述为分式规划S,wmax​wHR(S)wwHRt​(S)w​s.t.∣sn​(l)∣c该非凸优化问题采用Dinkelbach变换与坐标下降法求解迭代更新S 与w 直至收敛。2.2.3 CRLB最小化方法参数估计的克拉美-罗下界CRLB为波形设计提供理论性能边界。对于目标方向θt​ 与距离参数rt​ 的联合估计Fisher信息矩阵FIM元素计算为[J]ij​2Re{tr(∂ξi​∂Ry​​Ry−1​∂ξj​∂Ry​​Ry−1​)}∂ξi​∂myH​​Ry−1​∂ξj​∂my​​其中ξ[θt​,rt​]T 为待估参数矢量my​ 与Ry​ 分别为接收信号的均值与协方差。CRLB矩阵为FIM的逆CJ−1 。波形优化目标为最小化CRLB的迹或最大特征值Smin​tr(J−1)或Smin​λmax​(J−1)通过线性化近似与顺序凸规划SCP将原问题转化为一系列凸子问题求解。引入近端项保证子问题凸性利用KKT条件导出闭式解避免半定松弛带来的高计算复杂度。2.3 硬件约束与工程实现2.3.1 高功率放大器饱和特性与恒模需求高功率放大器HPA的非线性特性要求输入信号保持恒包络以避免幅度-幅度调制AM/AM与幅度-相位调制AM/PM失真。当信号峰均比PAPR过高时HPA工作于饱和区导致频谱再生与误码率恶化。恒模约束强制所有发射样本位于复平面圆上∣sn​(l)∣c,∀n∈{1,…,Nt​},l∈{1,…,L}该约束使优化问题呈现非凸特性属于NP-hard问题类别。常用求解方法包括循环算法Cyclic Algorithm、坐标下降法CD与流形优化Manifold Optimization。在恒模约束下波形仅保留相位自由度参数空间维度从2Nt​L 降至Nt​L 。2.3.2 峰均比约束与包络恒定化技术峰均比PAPR约束为恒模条件的松弛形式定义为PAPR(sn​)L1​∑l1L​∣sn​(l)∣2maxl​∣sn​(l)∣2​≤η其中η∈[1,Nt​L] 为PAPR阈值。过采样情况下PAPR约束表述为∥(FosH​⊗INt​​)X˙∥i,l2​≤NNt​ηPt​​,∀i,lFos​ 为过采样DFT矩阵Υ 为过采样率。包络恒定化技术包括直接相位编码、迭代限幅滤波与基于ADMM的凸优化方法。在ISAC系统中PAPR约束需与通信多用户干扰MUI联合优化通过交替方向乘子法实现通信与感知性能的折衷。2.3.3 相似性约束的p-范数形式相似性约束保证优化波形与参考波形如线性调频信号保持相位相关性以维持脉冲压缩与模糊函数特性。基于p -范数的相似性约束表述为∥S−S0​∥p​≤ξ其中S0​ 为参考波形矩阵ξ 为相似度阈值。当p2 时约束为Frobenius范数形式∥S−S0​∥F2​≤ξ当p∞ 时转化为逐元素幅度约束。p -范数约束可等价表述为相关函数形式​n1∑Nt​​l1∑L​sn​(l)s0,n∗​(l−τ)e−j2πfd​l−χ0​(τ,fd​)​≤ϵ其中χ0​(τ,fd​) 为参考波形模糊函数ϵ 为容许偏差。相似性约束与恒模约束联合作用时形成非凸可行域需采用投影梯度法或序列二阶锥规划SOCP求解。第二部分结构化伪代码算法1恒模波束赋形联合优化算法plain复制Algorithm Constant_Modulus_Beamforming Input: 理想方向图 d(θ), 角度采样网格 {θ_m}_{m1}^M, 阵元数 N_t, 序列长度 L Output: 恒模波形矩阵 S ∈ ℂ^{N_t×L} Initialize: S^(0) ← exp(j·random_phase(N_t, L)) R_s^(0) ← S^(0)(S^(0))^H / L k ← 0 while not converged do // 协方差矩阵优化凸松弛 R_s^(k1) ← argmin_{R_s ⪰ 0} Σ_{m1}^M w_m |a^H(θ_m)R_s a(θ_m) - d(θ_m)|^2 subject to: diag(R_s) c·1 // 功率约束 rank(R_s) 1 // 松弛后忽略 // 恒模映射Cyclic Algorithm for n 1 to N_t do for l 1 to L do s_n^(k1)(l) ← c · exp(j·arg([R_s^(k1)]_{n,l})) end end // 相似性约束投影可选 if ||S^(k1) - S_0||_p ξ then S^(k1) ← S_0 (S^(k1) - S_0) · ξ / ||S^(k1) - S_0||_p end // 收敛判断 if ||S^(k1) - S^(k)||_F / ||S^(k)||_F ε then converged ← true end k ← k 1 end while return S^(k) end Algorithm算法2SINR最大化联合收发设计plain复制Algorithm Joint_SINR_Maximization Input: 目标方向 θ_t, 杂波方向集合 {θ_k}_{k1}^{K_c}, 噪声功率 σ_n^2 Output: 最优波形 S* 与接收滤波器 w* Initialize: S^(0) ← LFM_reference(N_t, L) w^(0) ← b(θ_t) // 初始导向矢量匹配 k ← 0 while k K_max do // 固定S优化接收滤波器闭式解 R(S^(k)) ← Σ_{k1}^{K_c} σ_k^2 [b(θ_k)a^T(θ_k)S^(k)T_k^T][·]^H σ_n^2 I w^(k1) ← R^{-1}(S^(k)) b(θ_t) / ||R^{-1}(S^(k)) b(θ_t)|| // 固定w优化发射波形坐标下降 for n 1 to N_t do for l 1 to L do // 提取与s_n(l)相关的SINR分量 num ← |w^H [b(θ_t)a^T(θ_t)S T^T]_{:,l}|^2 den ← w^H [Σ_k σ_k^2 [b(θ_k)a^T(θ_k)S T_k^T][·]^H σ_n^2 I] w // 相位更新梯度下降 ∇_φ ← ∂(num/den)/∂φ_n(l) φ_n^(k1)(l) ← φ_n^(k)(l) μ · ∇_φ // 恒模投影 s_n^(k1)(l) ← c · exp(j·φ_n^(k1)(l)) end end // Dinkelbach分数规划迭代 λ^(k) ← SINR(w^(k1), S^(k)) if |λ^(k) - λ^(k-1)| δ then break end k ← k 1 end while S* ← S^(k) w* ← w^(k) return {S*, w*} end Algorithm算法3CRLB最小化线性近端方法plain复制Algorithm CRLB_Minimization_LPM Input: 初始波形 S^(0), 参数向量 ξ [θ, r], FIM计算函数 J(S) Output: CRLB最小化波形 S* Initialize: S^(0) ← random_complex(N_t, L) with |s_n(l)| c λ ← 1 // 近端参数 k ← 0 while not converged do // 在当前点线性化FIM J_lin(S; S^(k)) ← J(S^(k)) ⟨∇J(S^(k)), S - S^(k)⟩ // 构造近端子问题 S^(k1) ← argmin_S tr([J_lin(S; S^(k))]^{-1}) (λ/2)||S - S^(k)||_F^2 subject to: |s_n(l)| c, ∀n,l // 恒模约束 ||S - S_0||_p ≤ ξ // 相似性约束 // 利用KKT条件推导闭式解 for each n,l do // 相位更新规则 s_n^(k1)(l) ← c · exp(j·arg( [ (λS^(k) - ∇_S tr(J^{-1})|_{S^(k)}) ]_{n,l} )) end // 检查约束违反 if ||S^(k1) - S_0||_p ξ then S^(k1) ← Project_onto_p_norm_ball(S^(k1), S_0, ξ) end // 收敛判定子线性收敛率O(1/k) if ||S^(k1) - S^(k)||_F ε · (1 ||S^(k)||_F) then converged ← true end k ← k 1 end while return S^(k) end Algorithm算法4PAPR约束联合优化ADMM框架plain复制Algorithm PAPR_Constrained_Joint_Design Input: 通信信道 H, 参考波形 S_0, PAPR阈值 η, 权重 ρ ∈ [0,1] Output: 双功能波形 X ∈ ℂ^{N_t×L} Initialize: X^(0) ← S_0 Z^(0) ← 0 // 辅助变量 U^(0) ← 0 // 对偶变量 ν ← 0.1 // 惩罚参数 k ← 0 while k K_max do // X-update: 最小化雷达-通信加权目标 X^(k1) ← argmin_X ρ·MUI(X; H) (1-ρ)·||X - S_0||_F^2 (ν/2)||X - Z^(k) U^(k)||_F^2 subject to: (1/L)||X||_F^2 P_t // 总功率约束 // Z-update: 处理PAPR约束逐元素投影 for i 1 to N_t·L do if |x_i^(k1)|^2 η·P_t/(N_t·L) then z_i^(k1) ← √(η·P_t/(N_t·L)) · exp(j·arg(x_i^(k1))) else z_i^(k1) ← x_i^(k1) end end // U-update: 对偶变量更新 U^(k1) ← U^(k) X^(k1) - Z^(k1) // 残差计算 r_prim ← ||X^(k1) - Z^(k1)||_F r_dual ← ||ν(Z^(k1) - Z^(k))||_F if r_prim ε_pri and r_dual ε_dual then break end k ← k 1 end while return X^(k) end Algorithm算法5扩展目标检测波形优化plain复制Algorithm Extended_Target_Detection Input: 目标TIR统计特性 Σ_t, 杂波协方差 R_c, 检测门限 η Output: 检测概率最大化波形 S Initialize: S^(0) ← Orthogonal_Waveform(N_t, L) P_D^(0) ← 0 k ← 0 while not converged do // 计算当前波形下的检测概率 R_t^(k) ← E_t{b(θ_t)a^T(θ_t)S^(k)T^T T^* (S^(k))^H a^*(θ_t)b^H(θ_t)} Λ^(k) ← Generalized_Eigenvalue(R_t^(k), R_c σ_n^2 I) P_D^(k) ← Q_1(√(2Λ^(k)), √(-2ln(P_FA))) // Marcum Q函数 // 基于梯度上升更新波形 ∇_S P_D ← ∂P_D/∂S |_{S^(k)} // 流形优化保持恒模 for n 1 to N_t do for l 1 to L do // 切空间投影 grad_n,l ← ∇_S P_D|_{n,l} - Re{s_n^(k)*(l) · ∇_S P_D|_{n,l}} · s_n^(k)(l) // 更新步长回溯线搜索 α ← Armijo_Line_Search(P_D, S^(k), grad) // 流形上的指数映射 s_n^(k1)(l) ← c · exp(j·(arg(s_n^(k)(l)) α·grad_n,l/c)) end end // 验证相似性约束 if ||S^(k1) - S_0||_∞ ε_sim then S^(k1) ← argmin_S ||S - S^(k1)||_F s.t. ||S - S_0||_∞ ≤ ε_sim end if |P_D^(k) - P_D^(k-1)| δ then converged ← true end k ← k 1 end while return S^(k) end Algorithm参考文献索引本文技术内容参考以下国外学术文献体系基于恒模约束的MIMO雷达波形设计方法、低峰均比联合收发滤波器设计、信息论准则下的SINR约束波形优化、基于CRLB最小化的高效线性近端方法、双功能雷达通信系统的恒模波形设计、PAPR约束下的ISAC波形ADMM优化、多频谱约束低PAPR空时自适应处理、以及扩展目标检测的博弈论波形设计等。所有算法伪代码遵循IEEE算法环境规范采用结构化程序设计与数学符号混排形式。