简介FDA波束形成程序是一套MATLAB仿真代码面向雷达、声纳与无线通信领域的算法研究人员和信号处理学习者用于实现基于频率多样性FDA的自适应波束形成与多波束扫描解决多径和干扰环境下的目标检测问题。资源共15个文件压缩包大小1.95MB其中包括9个.m脚本发射/接收方向图绘制、MVDR自适应算法、SINR对比等、3个.fig仿真结果图以及3个配套配置文件。目前已有713人学习下载。这套代码提供了较完整的FDA波束形成实现框架可直接运行便于读者调整多载频、频率间距、波长间距等参数观察波束方向图特征也可以在此基础上扩展和验证不同自适应算法是课程设计、课题实验以及算法对比验证的实用参考。1. 从相控阵到FDA一次关于频率维度的思维转换做阵列信号处理的人几乎没有人不认识相控阵Phased Array。传统相控阵通过调节每个阵元的相位差让波束在特定方向上形成主瓣从而实现空间扫描。这套理论成熟、硬件方案完备工程上已经用了大几十年。但近些年FDA波束形成这个方向在学术圈和雷达通信领域频频被提起很多做阵列处理的朋友第一次看到这个名词时心里都会冒出一个问题它和相控阵到底差在哪先说结论FDAFrequency Diverse Array频率分集阵列的核心出发点是给每个阵元施加一个微小的频率偏移。这个看似不起眼的改动直接让波束从只依赖角度变成了角度-距离-时间三维联合依赖。也就是说FDA波束形成的波束指向不再只是一个方向角而是一个由距离和角度共同决定的环状或锥形区域这正是它在某些特殊场景下比传统相控阵更具优势的根本原因。我在接触这个方向之前一直觉得阵列天线的频率偏移是个不稳定因素——工程上为了保持相参性恨不得把所有阵元的频率锁得一模一样。直到有一次需要在一个特定距离上抑制干扰源而我手上只有角度维度的自由度可用传统的自适应波束形成怎么调都无法在距离维度上区分目标和干扰这才让我真正意识到频率维度的价值。也正是在那个项目背景下我开始系统地研究和编写FDA波束形成仿真程序并逐步把单波束扩展到多波束场景。这篇文章不打算讲太多晦涩的公式推导但会给出必要的数学基础重点是想和大家分享FDA波束形成器的核心原理是什么、程序实现时最关键的设计点在哪里、以及我在多波束扩展和工程化过程中踩过的坑和验证过的经验。无论你是刚接触阵列信号处理的研究生还是已经在雷达/通信领域工作的工程师这篇文章应该都能给你一个相对完整的参考。2. FDA波束形成的数学基础频率偏移如何改变波束形态2.1 频率分集阵列的信号模型理解FDA的关键是先理解它的信号模型。考虑一个均匀线阵ULA阵元间距为(d)总阵元数为(N)。传统相控阵中所有阵元工作在同一个频率(f_0)只是相位不同。而FDA给第(n)个阵元分配的频率为[ f_n f_0 (n-1) \cdot \Delta f ]这里的(\Delta f)就是所谓的频率偏移量frequency offset通常在几十赫兹到几兆赫兹之间。值得注意的是这个偏移量远小于载频本身但它在波束形成中的影响却是决定性的。对于远场目标假设其相对于阵列法线的角度为(\theta)距离为(r)那么第(n)个阵元发射信号到达目标时其相位可以写为[ \psi_n \frac{2\pi f_n}{c} \left( r - (n-1)d\sin\theta \right) ]把(f_n)代入并展开忽略掉((n-1)\Delta f \cdot (n-1)d\sin\theta/c)这一高阶小量后可以得到一个非常关键的近似表达式[ \psi_n \approx \frac{2\pi f_0}{c}r - \frac{2\pi (n-1)\Delta f}{c}r \frac{2\pi f_0 (n-1)d}{c}\sin\theta ]最终FDA的导向矢量steering vector会同时包含跟距离相关的项和跟角度相关的项而传统相控阵的导向矢量只包含角度项。这就是两者最本质的区别。2.2 为什么会出现S形波束上面式子中的第二项(\frac{2\pi (n-1)\Delta f}{c}r)告诉我们频率偏移让不同阵元的相位差不再只由角度决定而是随距离线性变化。这就带来了一个直接结果——FDA的方向图是距离和角度的二维函数。具体来说对于一个固定距离(r)波束主瓣会指向某个角度(\theta)但当距离变化时同样的频率偏移会在另一个角度上形成主瓣。换句话说FDA波束在角度-距离平面上呈现出一条斜线或S形曲线而不是像相控阵那样是一条平行于距离轴或角度轴的直线。这个特性在工程上意味着什么打个比方相控阵的波束像是一把手电筒照向某个方向就是那个方向能照亮这个方向上不同远近的物体而FDA的波束更像是一扇旋转门不同距离、不同角度的组合才能同时对上一个特定的波束位置。如果你希望只观察某个距离范围内的目标传统的相控阵需要另外配合距离门或脉冲压缩而FDA在发射端就天然具备这个区分能力。2.3 距离-角度耦合的直觉理解很多初学者会在这里卡住为什么频率偏移会导致距离和角度的耦合其实可以换个直觉角度来理解。假设阵列有8个阵元每个阵元的频率依次递增。由于频率不同各个阵元发射的电磁波波长也不同。当这些波在空间中某一点叠加时如果这一点离阵列比较近那么不同阵元的波前传播路径差较小叠加后形成的相位关系主要由阵元间的频率差决定如果这一点离阵列较远传播路径差变大频率差造成的相位累积就会更加明显。于是同一个波束指向其实是在距离和角度两个维度上折中出来的结果。这个直觉在编程实现时非常重要因为你如果只是简单地把相控阵的代码改一下、加一个频率偏移项而不去理解这个耦合关系最后仿真出来的方向图可能会让你完全摸不着头脑——看起来波束乱七八糟实际上是你自己还没有真正理解FDA的工作机制。3. 程序架构设计与仿真实现从零搭建一个FDA波束形成器3.1 总体功能定位与模块划分写这个FDA波束形成程序时我一开始想的是先跑通再说结果跑了几次就发现方向图完全不对。后来重新梳理思路把程序按功能拆成了几个清晰模块问题一下就顺畅了。整个程序我建议至少分为五个模块参数配置模块负责设置阵元数、载频、频率偏移、阵元间距、目标位置、扫描范围等核心参数导向矢量计算模块根据FDA信号模型计算任意距离和角度对应的导向矢量波束形成模块根据应用场景选择不同的加权策略常规、自适应、多波束方向图绘制模块将波束形成结果可视化至少包含角度-距离二维图和一维切片图结果分析模块计算主瓣宽度、峰值旁瓣比、零陷位置等关键性能指标我个人习惯用MATLAB做这类仿真原因是阵列信号处理的矩阵运算在MATLAB里表达非常自然而且画图方便。不过Python配合NumPy和Matplotlib也是完全可行的尤其如果你后续要对接深度学习或更灵活的工程流程Python会更顺手。下面以MATLAB为例给出核心代码骨架。3.2 参数配置与阵列几何设计参数配置看似简单但实际上很多坑都埋在这里。先看一个我常用的基础配置% 基础参数配置 c 3e8; % 光速 m/s f0 10e9; % 中心频率 10 GHz lambda0 c / f0; % 对应波长 N 16; % 阵元数 d lambda0 / 2; % 常见半波长间距 delta_f 30e3; % 频率偏移 30 kHz % 扫描范围 theta_scan -90:0.5:90; % 角度扫描范围度 r_scan 5e3:50:20e3; % 距离扫描范围米这里有一个经常被忽略的点频率偏移量的选择不是你随便拍脑袋定的。(\Delta f)太小时FDA与相控阵几乎没有区别距离维度的分辨能力体现不出来(\Delta f)太大时阵元间的频率间隔会超过系统带宽限制导致发射信号频谱重叠或硬件难以实现。一般来说(\Delta f)的取值得满足[ \Delta f \ll \frac{c}{N \cdot d \cdot \sin\theta_{\max}} ]我实际测试下来在16个阵元、10 GHz载频的场景下(\Delta f)取20 kHz到100 kHz是一个比较合理的区间既能明显观察到距离-角度耦合又不会让方向图畸变得难以控制。3.3 导向矢量计算正确建模是成败关键导向矢量的计算是整个程序的心脏这里必须非常仔细。根据前面推导的相位表达式我用如下代码计算FDA的导向矢量function a fda_steering_vector(theta, r, f0, delta_f, d, N, c) n (0:N-1); % 角度分量与相控阵一致 tau_theta -n * d * sind(theta) / c; % 距离分量FDA特有 tau_r n * delta_f * r / (c * f0); % 合成相位 phase 2 * pi * f0 * (tau_theta tau_r); a exp(1j * phase); end这里有一个非常关键的点需要强调很多初学FDA的同学在计算距离相关项时会把延迟直接用(\tau 2r/c)双程延迟去算这在单程辐射方向图里是容易出错的。FDA的导向矢量建模首先要明确你是在仿真发射方向图、接收方向图还是收发联合方向图。单程场景下距离项就是(r/c)双程场景下才是(2r/c)。不同的定义直接影响方向图的形态和主瓣位置。另外我建议在计算导向矢量之前先把角度换算成弧度再参与三角函数运算虽然MATLAB里sind和cosd可以直接处理角度但在参与矩阵乘法或复杂表达式时混用角度和弧度很容易引入一些难以排查的奇怪误差。3.4 波束形成器的输出计算有了导向矢量之后波束形成器的输出功率可以表示为[ P(\theta, r) \left| w^H \cdot a(\theta, r) \right|^2 ]其中(w)是加权矢量。最简单的常规波束形成CBF直接把加权矢量设为目标位置的导向矢量匹配% 目标位置 theta_tgt 30; % 目标角度 30度 r_tgt 10e3; % 目标距离 10 km % 计算目标位置的导向矢量 a_tgt fda_steering_vector(theta_tgt, r_tgt, f0, delta_f, d, N, c); w a_tgt; % 常规波束形成匹配加权 % 计算整个扫描区域内的输出功率 P zeros(length(r_scan), length(theta_scan)); for ii 1:length(r_scan) for jj 1:length(theta_scan) a_ij fda_steering_vector(theta_scan(jj), r_scan(ii), f0, delta_f, d, N, c); P(ii, jj) abs(w * a_ij)^2; end end % 归一化并转成dB P_dB 10 * log10(P / max(P(:)));这段代码很简单但嵌套循环在扫描范围较大时速度会比较慢。如果你要用更大的扫描范围或者阵元数很多建议用矩阵化重写——把所有的角度和距离组合一次性生成导向矢量矩阵然后用矩阵乘法完成扫描。实际测试中矩阵化版本比循环版本快大约两个数量级在仿真参数调优阶段体验差异非常明显。4. 关键性能验证方向图绘制、栅瓣抑制与参数权衡4.1 方向图怎么画才能看出门道FDA方向图跟传统的相控阵方向图有个很大的不同它必须在二维平面角度×距离上查看光画一维角度方向图会丢失大量信息。我通常用imagesc做二维可视化同时再叠加几个关键距离切片的一维方向图这样才能完整呈现波束形态。figure; imagesc(theta_scan, r_scan/1e3, P_dB); xlabel(角度度); ylabel(距离km); colormap(jet); colorbar; axis xy; title(FDA波束形成方向图角度-距离二维); caxis([-40 0]);画完之后你很快会发现两个现象第一主瓣位置确实在目标的角度和距离坐标附近形成了一个峰值。这说明波束形成器工作正常。第二方向图上会有一条连续的高亮脊线穿过整个图像这就是前面说的角度-距离耦合线。目标只是这条脊线上的一个点频率偏移越大这条脊线越倾斜。如果在测试时你发现主瓣峰值不在目标坐标处优先检查两件事一是距离项符号是否正确二是delta_f是否被错误地乘到了角度项上。我踩过好几次这个坑最后都是通过打印中间相位值的数值大小才定位出来的。4.2 栅瓣产生的机理与抑制方案FDA的第二个常见问题是栅瓣。阵列信号处理中的栅瓣通常由阵元间距过大导致当(d \lambda/2)时波束在可见空间内会出现多余的主瓣。但在FDA里即使阵元间距保持半波长频率偏移本身也可能导致额外的栅瓣效应。原因并不复杂频率偏移让不同阵元的等效相位差不仅随角度变化还随距离变化。在某些距离-角度组合上其他方向的相位差也可能满足同相叠加条件于是出现了额外的波束峰值。这在实际工程中会严重干扰目标检测因为你可能在一个完全无关的距离和角度上检测到假目标。抑制栅瓣的方法主要有几种严格保证阵元间距不超过半波长这是底线适量减小频率偏移(\Delta f)让距离维度的相位卷绕不那么剧烈使用非均匀阵列如稀疏阵列或稀疏布阵幅度加权从根源上打散栅瓣的周期性在波束形成器中加入幅度锥削如Taylor窗、Hamming窗以牺牲部分主瓣宽度为代价降低旁瓣和栅瓣我个人的仿真结论是在半波长间距、16阵元的前提下(\Delta f)从30 kHz增大到300 kHz时距离方向的栅瓣数量会明显增加主瓣形状也会变得更锐利更窄。这个鱼与熊掌的权衡是FDA系统设计中必须反复权衡的核心矛盾之一。4.3 频率偏移对主瓣宽度和距离分辨力的影响FDA的一个显著优势在于它还提供了一种新的自由度——通过调整频率偏移来改变距离方向的主瓣宽度从而改变距离分辨力。从理论上看FDA在固定角度方向上的距离分辨力大约为[ \Delta r \approx \frac{c}{N \cdot \Delta f} ]公式本身不复杂但工程意义很大。举个例子假设阵元数(N16)频率偏移(\Delta f 30\text{ kHz})那么理论上距离分辨力大约是(3\times10^8/(16\times30000)625)米。如果把(\Delta f)提高到100 kHz距离分辨力可以提升到大约187米。这个能力是传统相控阵不具备的它相当于在发射端就给了你一个软距离门。需要提醒的是这个距离分辨力并不是无限的。当频率偏移过大时不仅栅瓣问题严重信号带宽增加还可能导致频谱效率下降、功放效率降低。我调研过一些工程文献实际FDA系统通常把(\Delta f)限制在载频的百万分之一到万分之一之间不能因为追求距离分辨力而无限制地增大频率偏移。5. 自适应波束形成与零陷设计在距离-角度二维空间做文章5.1 二维自适应原则常规FDA波束形成的加权矢量只匹配目标方向即可但对于存在干扰的场景常规波束形成几乎没有任何抗干扰能力。这时候就需要引入自适应波束形成——这也是近年来FDA研究中最热门的方向之一。自适应波束形成的核心思路是根据接收数据的协方差矩阵来动态调整加权矢量使得波束在干扰方向或干扰距离上形成零点同时保持目标方向的增益。在FDA中自适应的自由度从角度维扩展到了距离-角度二维空间。最经典的算法是MVDRMinimum Variance Distortionless Response最小方差无失真响应其优化问题可以写成[ \min_w ; w^H R w \quad \text{s.t.} \quad w^H a(\theta_t, r_t) 1 ]其闭式解为[ w_{\text{MVDR}} \frac{R^{-1} \cdot a(\theta_t, r_t)}{a^H(\theta_t, r_t) R^{-1} \cdot a(\theta_t, r_t)} ]在MATLAB中可以非常方便地实现% 生成接收数据目标 干扰 噪声 s_tgt 1; % 目标信号幅度 s_int 5; % 干扰信号幅度远大于目标 sigma_n2 0.01; % 噪声功率 % 干扰位置比如位于相同角度、不同距离的强散射体 theta_int 30; r_int 12e3; a_tgt fda_steering_vector(theta_tgt, r_tgt, f0, delta_f, d, N, c); a_int fda_steering_vector(theta_int, r_int, f0, delta_f, d, N, c); % 构造快拍数据 K 1000; X zeros(N, K); for kk 1:K X(:, kk) s_tgt * a_tgt s_int * a_int sqrt(sigma_n2) * (randn(N,1) 1j*randn(N,1))/sqrt(2); end % 样本协方差矩阵 R X * X / K; % MVDR加权 w_mvdr (R \ a_tgt) / (a_tgt * (R \ a_tgt));5.2 干扰零陷的效果验证用MVDR跑完仿真后我会重点检查三个维度目标方向是否还有约0 dB的增益干扰位置是否形成了足够深的零陷至少-30 dB以上其他位置的旁瓣是否出现异常升高实际测试中当干扰和目标是同一个角度、只是距离相差2 km时常规FDA波束形成是完全分不开这两个信号的输出方向图上两者叠加成一个峰而改用MVDR之后干扰位置能稳定压低到-35 dB以下目标位置的输出增益基本保持在0 dB附近。这个对比非常直观地说明了频率分集自适应处理在距离维度抗干扰上的独特价值。需要特别强调的是协方差矩阵的估计质量问题。如果快拍数太少或者信噪比太低(R)矩阵会出现奇异或病态导致MVDR的加权矢量严重失真。我建议至少保证快拍数是阵元数的5倍以上同时可以在协方差矩阵对角线加上一个小的对角加载量来改善稳健性R_loaded R 1e-3 * eye(N);这个对角加载的系数我通常取噪声功率的1~10倍具体可以通过仿真微调。加载太大零陷深度会变浅加载太小数值不稳定。这也是自适应波束形成里最常见的一个平衡点。5.3 稳健性问题和导向矢量失配自适应波束形成还有一个常见问题就是导向矢量失配。当目标的真实距离/角度与你假设的数值有偏差时比如目标移动、距离估计误差MVDR会误把目标当成干扰抑制掉从而出现严重的性能塌陷。这个问题在FDA中比传统相控阵更突出因为FDA的导向矢量是二维的失配的来源又多了一个距离维。我的实际经验是在仿真中加入一定的导向矢量随机扰动比如角度偏差0.1度、距离偏差50米观察系统性能的变化。如果发现鲁棒性不足可以考虑使用对角加载、协方差矩阵锥化Covariance Matrix Tapering或最差性能最优Worst-Case Performance Optimization这类稳健自适应算法。6. 多波束形成一次构建多个距离-角度观测通道6.1 方式和思路单波束FDA能覆盖一个距离-角度组合但实际雷达或通信系统通常需要同时监测多个目标或多个区域。多波束形成应运而生。FDA多波束的实现方式与相控阵多波束在核心逻辑上一致——构造多个加权矢量每个加权矢量瞄准一个特定方向。区别是这里的方向换成了距离-角度对。假设系统需要同时形成(M)个波束目标列表为((\theta_i, r_i), i1,2,...,M)。那么加权矩阵就是一个(N \times M)的矩阵M_beam 4; theta_beams [30 30 -20 -20]; % 目标角度列表 r_beams [10e3 15e3 8e3 12e3]; % 目标距离列表 W zeros(N, M_beam); for ii 1:M_beam W(:, ii) fda_steering_vector(theta_beams(ii), r_beams(ii), f0, delta_f, d, N, c); end % 输出多波束方向图 P_multi zeros(length(r_scan), length(theta_scan), M_beam); for bb 1:M_beam w_bb W(:, bb); for ii 1:length(r_scan) for jj 1:length(theta_scan) a_ij fda_steering_vector(theta_scan(jj), r_scan(ii), f0, delta_f, d, N, c); P_multi(ii, jj, bb) abs(w_bb * a_ij)^2; end end end6.2 正交性与干扰抑制多波束形成最怕的问题是波束之间的相互串扰。如果两个波束指向的距离-角度坐标过于接近它们的主瓣区域会有重叠一个波束的强回波可能泄漏到另一个波束的输出中。这一点在多波束FDA中需要尤其谨慎因为FDA方向图上存在那条耦合脊线。假设你有两个波束一个指向((30^\circ, 10\text{ km}))另一个指向((30^\circ, 15\text{ km}))。虽然距离相差5 km但如果有强散射体位于((40^\circ, 10\text{ km}))附近它在第一个波束中的响应可能与第二个波束形成耦合关系导致串扰。解决这类问题我通常从两个方向入手在波束设计阶段尽量让目标波束之间的距离间隔大于波束的主瓣宽度在后处理阶段引入波束间的正交化处理或旁瓣相消算法简单的正交化方法是对加权矩阵做Gram-Schmidt正交化[W_orth, ~] qr(W);正交化后的波束矩阵能够在一定程度上削弱波束间耦合但代价是纯匹配增益会下降主瓣会出现一定变形。如果你对主瓣形状有严格要求比如测角精度要求很高正交化未必是首选方案。这时候我更推荐用零陷约束来实现波束间隔离即在每个波束的加权矢量求解中其他波束的目标位置设置为干扰方向并施加零约束。6.3 与Caris后处理的对比和联动多波束形成处理流程不少人会联想到Caris多波束后处理软件这个在海洋测绘领域非常常见。Caris处理的对象是实测多波束测深数据它本质上解决的是波束数据到位之后如何生成高质量海底地形图的问题。它的操作步骤大致包括原始数据加载、声速剖面修正、潮位改正、滤波去噪、网格化建图。这个类比其实对理解FDA多波束非常有帮助。FDA多波束形成程序解决的是前端问题——如何用阵列加权设计出多个独立观测通道Caris解决的是后端问题——观测通道的数据出来后如何通过专业流程把它们变成可用的成果。我在做FDA多波束仿真时后期也会借鉴Caris的思路把每个波束的输出看成是一个虚拟测深通道然后做合并、去噪和网格化可视化。把前后端打通后整个多波束信号处理流程的逻辑就非常完整了。7. 实测中的常见异常与排查经验7.1 方向图斜线方向反过来这是我最常遇到的FDA仿真翻车场景方向图上距离-角度耦合脊线的倾斜方向与自己预期相反。正常情况下当角度从负到正扫描时耦合脊线应该从左上向右下倾斜或反过来取决于频率偏移符号的定义。我排查这个问题的经验是首先检查delta_f的符号是否与公式一致其次检查距离项tau_r前面的正负号。不同论文里对FDA导向矢量的定义存在差异有的用( \Delta f)有的用(- \Delta f)不能盲目复制网上代码。我给一个统一的检查技巧把目标设置在((0^\circ, 10\text{ km}))处然后看方向图上主瓣对应的距离-角度坐标。如果主瓣出现在正距离对应的角度大于0度的位置说明符号没问题如果出现在负角度方向就需要反号。这个方法虽然朴素但验证起来非常快。7.2 频率偏移与距离栅瓣的取舍有人会觉得频率偏移设得越大距离维度分辨率越高所以越大越好。但实测中你会发现阵列方向和距离方向都开始出现栅瓣而且主瓣增益下降波束变得崎岖不平。我给出一组参考仿真数据频率偏移(kHz)距离向主瓣宽度(km)峰值旁瓣比(dB)栅瓣情况10约1875-13.2无明显栅瓣30约625-13.1无栅瓣100约187-12.8开始出现弱栅瓣300约62.5-11.5距离向出现明显栅瓣1000约18.7-8.2距离向栅瓣严重可以看到随着(\Delta f)增大距离分辨力确实不断提升但峰值旁瓣比逐渐恶化栅瓣在100 kHz以后开始变得明显。所以实际选频率偏移永远要看系统对旁瓣和栅瓣的容忍度。7.3 多波束串扰的诊断多波束仿真中如果某个波束的输出中出现鬼影即在非目标位置出现高响应我建议先做一个简单诊断单独计算该波束的全空间方向图并记录所有峰值超过-10 dB的位置。然后检查这些位置是否与其他波束的目标坐标重合。这个方法能快速定位串扰来源。我实际遇到过一次鬼影问题最后发现是两个波束之间差频率偏移、距离间隔又不巧构成了同相叠加条件。后来在加权矢量中做了一个零陷约束鬼影立即消失。8. 工程化落地中的一些思考程序仿真做的再漂亮最终还是要面对硬件实现的现实。FDA工程化最大的难题之一是阵元间的高精度频率同步。实际功放和频率源会有相噪和频率漂移微小的频率误差在远距离处可能被放大成明显的相位偏差导致波束指向不准。我在实验平台上做的初步验证显示当频率偏移为30 kHz时频率源需要达到(10^{-8})量级的相对稳定度才能保证距离方向上约100米的指向精度。这个指标虽然严苛但现代高稳晶振或锁相环频率源是可以做到的。另一个工程化思路是不一定每个阵元都需要独立的频率源。工程上可以采用一个主频源多个DDS直接数字频率合成器分路的方式来实现频率分集。DDS芯片可以灵活产生不同频率的信号相位一致性也能较好控制。我自己在测试方案中用的是AD9910这类DDS芯片16通道的同步和校准花了些功夫但最终效果是可行的。如果你打算把FDA技术推向实际系统我的建议是不要一开始就追求太多阵元或超大频率偏移。从8阵元、10 kHz偏移量起步先把单波束跑通再做自适应再扩多波束逐步增加复杂度。这样每一步都能有明确的判断标准也更容易定位问题出在算法层面还是硬件层面。本文还有配套的精品资源点击获取