简介：面向MIMO（多入多出）雷达与通信波形设计研究者和工程师的Matlab源码包，由李建团队编写，围绕主动感知系统波形设计的前七章内容展开，涵盖波形生成、性能评估、参数扫描等核心环节。压缩包共8个文件，以.mat数据文件与.m主程序脚本为主，体积约16.6MB；其中chapter5_stopband相关文件专门演示了阻止带控制场景，便于理解滤波器频段约束对雷达探测距离、角度分辨率、速度分辨率及抗干扰能力的影响。该资源已有418人学习下载，适合课程设计、课题预研或工程复现时参考。读者可结合源码系统梳理MIMO雷达中的空间分集、时间分集、频率分集与编码技术，熟悉主动波形设计中的多目标区分、低截获概率、带宽利用等实际问题；还可以通过修改扫描变量、对比不同.mat数据结果，掌握从理论到仿真的完整实验方法。这些代码既展示了李建团队在MIMO雷达波形设计上的具体实现思路，也为后续在军事监控、航空导航、交通管理等领域开发或优化雷达系统提供了可扩展的基础工程代码，具有较强的研究参考价值。
1. MIMO 雷达主动波形设计：解压这份代码包先别急着运行 main.m
解压《Matlab Code.rar》之后，目录里躺着 main.m、chapter5_stopband、twobands.mat，以及 scan1.mat、scan2.mat、scan_continuous.mat 这一串结果文件。如果直接双击 main.m，多半会碰到矩阵维度不一致或者找不到自定义函数；这不是坏包，而是 MIMO 雷达主动波形设计里很常见的多阶段产出结构。李建团队的这套源码把发射波形设计从正交基线一直推到 stopband 连续扫描，核心回答了一个问题：在给定发射通道数和采样点数的前提下，如何同时压低指定频段的能量，又让各通道波形保持足够低的互相关旁瓣。适合正在做认知雷达、雷达通信一体化，或者想从正交波形走向频谱成形设计的人。新手可以按文件依赖关系还原实验，熟手可以对比 rho/lambda 两个权重参数的取舍。
2. 发射波形自由度与正交性约束：信号模型怎么落到代码
2.1 先把 MIMO 雷达的发射端矩阵推清楚
MIMO 雷达与相控阵的差别在于发射阵元不必发同一个波形。假设发射阵元数为 M，每个脉冲采样 N 点，那么一个脉冲内的发射波形可以写成 M×N 复矩阵 W。第 m 行是第 m 个阵元发出的采样序列，第 n 列是全部阵元在同一时刻的空间快照。接收端对每个发射波形做匹配滤波，等价于在时域把回波与各发射序列相关。要恢复出每一路发射分量的独立信息，波形矩阵需要满足理想正交条件：
W * W' / N = I
这里 W' 是共轭转置，I 是 M×M 单位阵。从空间角度看，正交波形能形成各向同性的发射协方差矩阵 R = W*W'/N，也就是方向图接近均匀的发射功率分布；从目标检测看，正交性保证了多通道匹配滤波后干扰项不会互相叠加。可惜这是理想极限，实际波形还要受到恒模放大器约束和频谱占用限制，完全正交波形往往无法同时满足。
因此工程上的做法是放弃完全正交，改成“低自相关旁瓣 + 低互相关旁瓣 + 指定频谱形状”的联合优化，把波形设计变成一个高维非线性规划问题。前面提到的 rho 和 lambda 两个扫描参数，分别对应旁瓣项与频谱项在代价函数里的权重，这也是为什么 main.m 反复加载 scan_rho.mat、scan_lambda.mat。
2.2 用 Matlab 快速生成一组正交基线
要理解这些约束，第一步先在 Matlab 里生成一个最理想的正交集。常见做法是用离散傅里叶变换（DFT）基构造波形矩阵，每一行是一个不同频率的复正弦，采样点数为 N 时各行天然正交。代码如下：
% gen_ortho_baseline.m M = 4; % 发射阵元数 N = 256; % 每个波形的采样点数 t = (0:N-1).' / N; % 归一化时间轴 base = exp(1j * 2 * pi * (0:M-1).' * t.'); % M×N 复正弦矩阵 base = base ./ sqrt(sum(abs(base).^2, 2)); % 对每行做能量归一化 corr_mat = (base * base') / N; % 互相关矩阵 disp('互相关矩阵幅度：'); disp(abs(corr_mat));这段代码先构造一个 4×256 的复指数矩阵，第 m 行频率是 m/N，能量归一化后每行总功率为 1。corr_mat 用来检查通道间相关：理想情况下对角元是 1，非对角元接近 0。DFT 基的好处是构造简单、数值稳定，不用迭代；坏处是每一行都是单频信号，频谱只有一根谱线，带宽利用率极低，完全不能当实际雷达波形用。它只适合作为后续迭代算法的初始点，或者做正交性判定的参照物。
如果把 DFT 基换成随机相位波形，比如exp(1j * randn(M,N))，互相关会变成类似噪声的小值，但最差旁瓣仍然不够低。实际雷达更常用恒模相位编码序列（如 P4 码、Gold 序列）做初始值，因为它们本身已经具备较好的自相关旁瓣。
2.3 用一张表量化波形性能
正交性和频谱成形效果不能只靠肉眼，通常要盯住几个指标。这里整理成一张速查表：
| 指标 | 定义 | Matlab 计算 | 对雷达的影响 |
|---|---|---|---|
| 自相关旁瓣 PSL | 主瓣外最大旁瓣与主瓣峰值之比 | max(abs(xcorr(w(1,:), 'coeff'))) | 弱目标被强目标旁瓣淹没 |
| 互相关峰值 | 不同通道波形相关结果的峰值 | max(abs(xcorr(w(1,:), w(2,:)))) | 目标跨通道串扰 |
| 阻带能量比 | 阻带频段内功率与总功率之比 | sum(abs(fft(w, 2^nextpow2(N))).^2, 对应下标) | 对带外设备的干扰 |
设计目标一般包含三项：自相关集成旁瓣电平尽量小，通道间正交性尽量高，阻带能量尽量低于某个门限。三者互相矛盾，所以没有“唯一正确波形”，只有针对当前雷达工作场景的 Pareto 最优波形。这就是为什么代码包里 scan_rho.mat 和 scan_lambda.mat 是分开保存的：扫描权重而不是追求单一解。
2.4 为什么主动波形设计不能停在“正交”
如果目标是单纯的正交，DFT 基就结束了，整套代码也就没有理由出现 stopband 和 twobands 文件。主动波形设计的关键在于“主动”两个字：系统知道频谱环境里有固定频段受干扰，便主动在波形上切出阻带，而不是等接收端用滤波器处理。接收滤波只能抑制接收路径上的干扰，发射端若持续在干扰频带辐射能量，不仅浪费功率，还可能干扰同址或邻近设备。把频谱约束前置到发射波形设计里，是 MIMO 雷达获得频谱共存能力的核心手段。
从算法实现上看，这意味着优化问题的可行域被强行切掉一块，正交波形集合不再闭合。每一步迭代都要在保持恒模、压低自相关旁瓣的同时，把频谱幅度在指定频段内压低。下一章讲的 stopband 与 twobands，就是这一约束在 Matlab 里的具体表达。
3. 阻带抑制与多频段连续扫描：stopband 和 twobands 的算法逻辑
3.1 把 stopband 写成代价函数里的约束
从文件名 chapter5_stopband 看，这对应书稿的第五章，主题是阻带约束下的波形设计。雷达频谱上存在受保护频段，例如民航、通信基站和卫星下行链路所在频段，任何本地辐射都不能携带显著能量穿过这些区间。但要直接在优化变量上施加严格为零的约束，波形自由度会被压得过小，迭代难以收敛。常见做法是把阻带能量作为软约束加入代价函数，并给一个足够大的权重 lambda：
J(W) = (1 - rho) * 正交性损失 + rho * 自相关旁瓣损失 + lambda * 阻带能量
其中 rho 控制自相关/互相关与恒模之间的平衡，lambda 控制阻带抑制度。lambda 取 0 时是不加频谱约束的正交波形优化；lambda 取太大时，波形频谱会被削出很深的凹槽，但自相关旁瓣和互相关指标会急剧恶化。优秀的设计是 lambda 与 rho 都落在某一个平台上，指标变化平缓，这正是扫描这两个参数的实际意义。
3.2 交替投影是最快的验证基线
这类问题的可复现基线通常用交替投影：在“时域恒模/低旁瓣”和“频域阻带约束”两个集合之间来回投影。我自己做仿真时，会先写一个最小函数，不追求最优性能，只验证 stopband 约束是否真的能压进去：
function w = design_stopband_min(M, N, stop_idx, lambda, nIter) % M: 发射阵元数; N: 采样点数; stop_idx: 阻带频点下标 % lambda: 阻带幅度衰减因子; nIter: 迭代次数 w = exp(1j * 2 * pi * rand(M, N)); % 随机相位初始 S = fft(w, [], 2); % 每行做时域到频域 for k = 1:nIter % 投影 1：频域阻带压低 S(:, stop_idx) = lambda * S(:, stop_idx); % 回时域 w = ifft(S, [], 2); % 投影 2：恒模映射 w = exp(1j * angle(w)); % 重新进频域 S = fft(w, [], 2); end end代码里 lambda 是 0 到 1 之间的衰减因子，表示阻带里的复频谱幅度被压缩到原来的 lambda 倍，工程上常取 0.01~0.1。投影 2 的恒模映射保证了每个采样点幅度等于 1，避免波形经过功率放大器后失真。正式源码一般会在恒模映射后加一个时域旁瓣优化步骤，但原理相同：频域压低、时域恢复恒模、再频域压低，绕圈收敛。这是理解 chapter5_stopband 的最小路径。
这里有一个需要注意的细节：迭代时如果每次都完整 FFT/IFFT，长波形的计算量不小。M=16、N=1024 时一次迭代大约是 M×N×logN 的复杂度，300 次迭代在普通台式机上要几十秒。所以 scan_continuous.mat 这类文件往往是一次性离线算好的结果，主程序加载它做后处理和可视化，而不是每次打开 main.m 都重新扫描。
3.3 twobands 与连续扫描的实验布局
twobands.mat 从名字看保存的是“两个连续阻带”的边界，典型情形是雷达工作频段两端各有一个保护带，或者中间夹着两个通信频点。实际载入后可能是一个 2×2 矩阵：第一行是第一个阻带的下、上界，第二行是第二个阻带的下、上界。扫频时把边界换算成 FFT 下标，例如 N 点 FFT 下频率分辨率是 fs/N，下标等于 round(f*N/fs)。
scan_continuous.mat 从命名看是“持续扫描”的结果。我做类似实验时，会按阻带中心频率从低到高移动，每个中心频率记录一个波形文件，最后堆成三维矩阵。扫描的目的不是找某一个波形，而是观察阻带位置对旁瓣的影响。参数调节时可以参照下表：
| 参数 | 推荐范围 | 作用 | 注意 |
|---|---|---|---|
| lambda | 0.01~0.1 | 阻带抑制深度 | 太小凹槽浅，太大旁瓣抬升 |
| rho | 0.1~0.9 | 旁瓣项权重 | 大则正交性好，压不住阻带 |
| nIter | 200~500 | 迭代次数 | 超过 500 后收益递减 |
% scan_continuous_example.m fc_list = 0.10 : 0.01 : 0.45; % 归一化中心频率 scan_result = zeros(M, N, length(fc_list)); for ii = 1:length(fc_list) [lo, hi] = freq2idx(fc_list(ii), 0.05, N); % 带宽 5% wc = design_stopband_min(M, N, lo:hi, 0.02, 200); scan_result(:, :, ii) = wc; end save('scan_continuous.mat', 'scan_result', 'fc_list');这里 freq2idx 是自己写的下界到 FFT 下标换算函数。为什么连续扫而不是只测 twobands 两个位置？因为靠近频谱边缘的阻带容易压，靠近中心频带的阻带往往需要更大 lambda，否则凹槽深度不够。连续扫描数据可以画出“阻带深度-旁瓣电平”随中心频率的变化曲线，直接暴露哪些频段是安全的、哪些频段会牺牲太多正交性。wescan.mat 和 scan2.mat 则大概率是扫描过程中保存的中间波形，用于对比不同迭代次数的收敛状态；由于没有配套说明文档，稳妥的做法是先用whos -file查变量名和维度，再决定重放方式。
4. main.m 与 scan_*.mat 的组合关系：从解压到出图的复现路径
4.1 文件角色与数据流
资源包里的文件不是随机散落的。按命名可以分成三组：入口脚本、设置文件、扫描结果。下表是一个可以快速入手的映射关系：
| 文件/文件夹 | 推断角色 | 使用时机 |
|---|---|---|
| main.m | 主流程，加载 .mat 并调用设计函数 | 第一步尝试运行 |
| chapter5_stopband | 第五章配套脚本/函数，处理阻带约束 | 从 main.m 中被调用 |
| twobands.mat | 两个连续阻带的上下界 | 设计前的频谱掩码 |
| scan1.mat / scan2.mat | 两个不同初始化或权重下的波形结果 | 对比实验 |
| wescan.mat | 波形扫描结果保存集 | 后处理画图 |
| scan_rho.mat / scan_lambda.mat | rho 与 lambda 参数扫描记录 | 画折中曲线 |
| scan_continuous.mat | 阻带中心连续扫描结果矩阵 | 观察阻带移动影响 |
这些推断以文件名为准。没有 docs 说明时，依赖关系要通过matlab -batch逐步探测。我的习惯是先运行whos -file scan_continuous.mat，确认每个变量名、维度和数据类型，再决定分析脚本怎么写。
4.2 修正路径之后重跑 main.m
直接双击运行最可能遇到的问题的是路径不对：调用load('twobands.mat')时当前工作目录不在源码目录。用如下命令统一设置路径并运行：
cd /path/to/MatlabCode matlab -batch "addpath(pwd); main"在 Matlab 内则推荐：
cd(fileparts(mfilename('fullpath'))) % 切到源码根目录 addpath('chapter5_stopband') % 把函数目录加进来 main然后观察命令行是否提示缺少函数。MIMO 雷达波形设计一般只依赖信号处理工具箱，极少调用深度学习或通信工具箱，因此在 R2023b 这类常见版本上可以直接跑；如果报错来自freqz、xcorr，检查 Signal Processing Toolbox 是否安装。main.m 运行完成后通常会在工作区留下几个 scan 矩阵，但这些矩阵不一定自动保存，源码包里已有的 scan1.mat、scan2.mat 应该是作者当时跑完留下的。
4.3 读 .mat 遇到版本问题怎么办
Matlab 从 R2012 开始默认保存格式是 HDF5 变体，旧版打开新版 .mat 会报 “Unable to read MAT-file”。更隐蔽的是load在有大矩阵时把全部数据读进内存，scan_continuous.mat 如果达到几百 MB，普通电脑会卡死。这时改用matfile做局部读取：
m = matfile('scan_continuous.mat'); [rows, cols, pages] = size(m, 'scan_result'); % 只读第 50 页对应的波形切片 slice = m.scan_result(1:4, :, 50);matfile不会立刻加载整个文件，而是只把需要的切片读到工作区，处理连续扫描三维矩阵非常有用。若 .mat 来自 Linux 下旧版 Matlab，需要检查 endian 字段，但通常load会自动处理，只有极少数跨平台大文件会出问题。遇到Unsupported MAT-file version 7.3时，用对方电脑重新save(...,'-v7')再拷贝即可。
4.4 从波形结果画第一张频谱图
把算法生成的波形与 twobands 边界放在同一张图上，是验证是否压住停止带的必要一步。常见做法如下：
% plot_stopband_check.m load('twobands.mat', 'twobands'); % 假设保存为 twobands load('scan1.mat', 'wav'); % 假设保存为 wav N = size(wav, 1); % 波形长度 f = (0:N-1) / N; % 归一化频率 spec = abs(fft(wav)); % 单通道频谱 plot(f, mag2db(spec)); hold on; for k = 1:size(twobands, 1) xline(twobands(k, 1), 'r--', 'stopband'); xline(twobands(k, 2), 'r--'); end xlabel('归一化频率 (cycles/sample)'); ylabel('幅度 (dB)');这段代码用xline标出阻带边界，快速看凹槽是否落在两个边界之间。凹槽深度至少要比通带低 20 dB 才基本达标，否则要把 lambda 增加两个数量级重跑。如果主程序输出的是 scan2.mat 而不是 scan1.mat，把第二行文件路径替换即可，整体绘图逻辑不变。
5. rho/lambda 折中曲线的验证技巧与几个反直觉结论
5.1 用网格扫描画出旁瓣-阻带抑制前沿
rho 控制旁瓣项权重，lambda 控制阻带项权重，单点调参很难判断算法是否到了极限。建议把两个参数各取 10 个 log 间隔值，跑 100 组短迭代，记录每组结果的阻带抑制度与互相关旁瓣峰值，然后画散点图：
rho_list = logspace(-2, 0, 10); lam_list = logspace(-3, -1, 10); metric = zeros(10, 10, 2); for i = 1:10 for j = 1:10 w = design_stopband_min(4, 128, stop_idx, lam_list(j), 50); metric(i,j,1) = stopband_power(w, stop_idx); metric(i,j,2) = max_abs_xcorr(w); end end scatter(metric(:,:,1), metric(:,:,2), 20, 'filled');这里的 design_stopband_min 可以用第三章的最小实现，stopband_power 自己写两行 FFT 能量统计。把两个指标画在二维平面上，能定义出一条 Pareto 前沿：落在右下方的点比左上方的好。直接在 log 轴上搜索，比盯着单个 lambda 盲调效率高得多。
5.2 恒模约束和频谱凹槽并不互斥
很多第一次接触的人会认为恒模波形频谱形状被锁死，无法在某个频段压出深凹槽。实际上恒模只限制幅度恒定，相位序列的自由度仍然很高；经过几百次交替投影，频谱凹槽能压到 -30 dB 以下，恒模偏差控制在浮点误差内。关键不是恒模本身，而是旁瓣项的权重 rho 不能太大，否则迭代过程会被旁瓣梯度占满，没余量学频谱形状。
5.3 越“白”的波形往往越不耐看
用功率谱看设计结果，若频谱过度平坦，阻带凹槽附近会出现明显的“振铃”式旁瓣抬升。把 lambda 从 0.01 加到 0.1，阻带能量确实下降，但通带边缘会出现额外 3~5 dB 的波纹，进而抬高自相关旁瓣。因此收敛后要用阻带边缘一个窄过渡带的宽度做验收，而不是只看凹槽深度。用sum(abs(diff(spec(stop_idx))).^2)量化过渡带爬升，这一指标也应写进扫描结果。
5.4 验证方向图比验证波形本身更可靠
设计完波形不要只盯着波形矩阵，应立刻计算发射协方差 R 并画方向图plot(abs(steering'*R*steering))。MIMO 雷达的探测性能最终由角度功率分布决定，阻带优化很可能改善了频谱但让方向图在某个角度出现凹陷。两件事同时验收，才算一次完整的主动波形设计闭环；若发现方向图异常，回到扫描结果里找对应 rho/lambda 组合，比重新设计更快。
本文还有配套的精品资源，点击获取
