简介这是一份面向无线通信与集成感知ISAC研究方向的学术论文PDF资源作者围绕液态天线辅助的ISAC系统展开适合高校科研人员、算法工程师以及关注6G通信物理层设计的读者学习。论文指出传统固定位置天线在满足雷达感知波束增益与发射功率约束时难以兼顾下行通信速率因此提出在基站和通信用户处引入可动态调整位置的液态天线并对液态天线位置与发射波束成形进行联合优化。由于原问题具有强非凸性和变量高度耦合作者设计了交替优化算法将联合优化拆分为三个子问题利用凸优化、泰勒展开等数学工具逐一处理最终获得可求解的凸问题仿真结果表明该系统相比基准方案显著提升了ISAC综合性能。资源压缩包仅包含1个PDF文件大小约623KB内容涵盖论文摘要、系统模型、算法推导以及仿真验证便于读者复现和借鉴其中的优化思路。目前已有68人学习下载适合作为ISAC、液态天线、凸优化与波束成形交叉领域的精读资料也可为5G/6G网络中的频谱共享和高性能MIMO设计提供参考。1. 液态天线辅助集成感知与通信系统为什么共孔径方案绕不开液体金属5G-A 和 6G 的讨论里集成感知与通信ISAC已经从论文概念走向了原型验证但真正落地的瓶颈往往不在基带算法而在天线口径。感知任务需要宽波束覆盖来搜索目标通信任务需要高增益窄波束来保证链路预算传统相控阵在单元数受限时很难兼顾这两种模式。液态天线辅助集成感知与通信系统的核心思路是把辐射体本身变成可重构的流体通过改变液态金属柱的高度或形状动态调整工作频率、波束宽度和极化方式。这篇文章面向要做原型验证的工程师从材料选型、等效建模、可复现代码到封装避坑讲清楚一条能直接上手的技术路径。2. 液态天线在ISAC中的物理基础它到底动了哪块电磁边界2.1 感知与通信对天线的需求冲突在哪个频段最尖锐先说一个反直觉的结论ISAC 的频率越低液态天线的收益反而越大。Sub-6GHz 频段下天线尺寸大相控阵要做到 8 通道以上阵面面积和成本都很夸张感知波束和通信波束想共用同一组单元只能靠时分切换这牺牲了感知的持续跟踪能力。到了毫米波频段阵列集成度高波束赋形灵活液态天线的优势更多体现在极化可重构和频段捷变上。从电磁边界看天线的谐振频率主要由辐射体的电长度决定。液态金属柱的高度变化直接改变有效电长度等效于在同一个物理端口上实现了不同的频段覆盖。对于 ISAC 而言这意味着通信频段例如 3.5GHz和感知频段例如 4.9GHz可以通过调整液态柱高度切换而不是靠两套独立的射频链路。这是液态天线辅助集成感知与通信系统的第一层价值频谱复用。第二层价值在波束形状。感知模式偏好较宽的半功率波束宽度以便在低信噪比下捕获目标回波通信模式则需要窄波束来抑制干扰、提升信噪比。液态天线通过改变辐射体的横向尺寸例如液态金属圆盘直径可以在同一个工作频率上改变方向图波束宽度而这一能力在传统贴片天线或者金属腔体天线上几乎不可能在线实现。2.2 液态金属材料选型EGaIn、Galinstan 与电解质溶液的参数对比材料是一切可复现工作的前提。我在实际测试中发现并不是所有能导电的液体都适合做射频天线。下表列出三种常见方案的核心参数这是做仿真和实测前必须定下来的基础。材料类型电导率 (S/m)粘度 (cP)熔点/凝固点驱动方式典型应用场景共晶镓铟 (EGaIn)3.4 × 10⁶~2.015.5°C过冷可到 -19°C气压、电润湿可重构单极子、缝隙天线镓铟锡合金 (Galinstan)3.3 × 10⁶~2.4-19°C气压、介电泳需要低温环境的外场测试离子液体/盐水5 ~ 2010可定制电泵、机械泵低频段、损耗容忍度高的宽带天线实际选型时电导率直接关联辐射效率。EGaIn 的电导率比盐水高出近六个数量级在铜5.8×10⁷ S/m和铝3.5×10⁷ S/m之间适合做谐振类辐射体。盐水虽然也能辐射但损耗太大在 ISAC 系统里用来做频率选择表面或者反射面较多不适合做主辐射单元。Galinstan 的优势是不含铟价格略低凝固点更低适合在北方冬季外场测试。但它对铜和铝的腐蚀性更强封装材料要避开这两种金属。EGaIn 在接触铜时会发生镓致脆化GIMC这是做液态天线复现时第一个要踩的坑后文避坑章节会详述。2.3 从系统级视角看液态天线如何参与 ISAC 共孔径架构ISAC 的共孔径有两种主流架构时分复用和频分复用。时分复用最简单感知和通信分时占用同一组天线但感知刷新率受限于切换速度液态天线的响应时间通常 10-100ms 级在这种架构下完全够用。频分复用则要求天线同时支持两个频段这正好对应液态天线的频段捷变能力。例如一个典型的 ISAC 帧结构里下行通信时隙用 3.6GHz 窄波束做数传帧间空闲时隙切到 4.8GHz 宽波束发射线性调频信号做感知。液态天线的液态柱高度从 20mm 收缩到 15mm谐振点从 3.6GHz 移到 4.8GHz中间不需要额外的匹配网络。这里的关键是液态柱的高度与谐振频率的关系并不是线性而是近似反比。实际建模时需要用到后文的等效电路法。3. 可复现的建模链路从液态柱参数到天线方向图的计算3.1 感应电动势法估算液态单极子天线输入阻抗要做可复现必须先有能算的数学模型。液态天线辅助集成感知与通信系统中最常见的辐射体是液态金属单极子它有解析解适合作为参照基准。这里用感应电动势法EMF写一个 Python 脚本计算不同液态柱高度下的输入阻抗。import numpy as np import matplotlib.pyplot as plt # 感应电动势法计算单极子天线输入阻抗 # 参考文献Balanis, Antenna Theory, Ch. 8 def monopole_z_in(freq_hz, length_m, radius_m): 计算垂直导电单极子在地平面上的输入阻抗 :param freq_hz: 工作频率单位 Hz :param length_m: 液态金属柱高度单位 m :param radius_m: 液态金属柱半径单位 m :return: 输入阻抗 复数 c 2.99792458e8 k 2 * np.pi * freq_hz / c eta 120 * np.pi # 自由空间波阻抗 ~377 欧姆 # 对正弦电流分布做数值积分 N 400 # 分段数步长越小越接近连续电流分布 z np.linspace(0, length_m, N) dz z[1] - z[0] # 计算辐射电阻 Rr # 使用 Pocklington 积分方程的近似形式 # 这里采用经典的感应电动势法 # 先假设电流 I(z) I0 * sin(k * (L - z)) # 再求远场辐射电阻 current np.sin(k * (length_m - z)) # 远场电场积分系数 # E_theta 正比于 sin(theta)通过坡印廷积分得到 Rr # 数值积分公式(参考 Balanis 第8章) integral 0.0 M 1000 theta np.linspace(0.001, np.pi - 0.001, M) dtheta theta[1] - theta[0] for th in theta: # 方向图函数 F(theta) Fth np.abs((np.cos(k * length_m * np.cos(th)) - np.cos(k * length_m)) / np.sin(th)) integral Fth**2 * np.sin(th) * dtheta Rr (eta / (2 * np.pi)) * integral # 单极子仅上部空间辐射系数修正 # 输入电抗 Xm 近似(基于Esin的近似公式) # 这里用tan函数校正当趋近谐振点时精度会下降 if abs(np.mod(k * length_m, np.pi)) 0.05: Xm 0.0 # 正好谐振时电抗趋于0 else: # 近似计算 Xm -120 * (np.log(length_m / radius_m) - 1) / np.tan(k * length_m) return complex(Rr, Xm) # 扫描液态柱高度固定频率为 3.5GHz freq 3.5e9 radius 0.002 # 2mm 半径 heights_mm np.linspace(8, 25, 50) z_series [] for h_mm in heights_mm: z monopole_z_in(freq, h_mm / 1000, radius) z_series.append(z) z_series np.array(z_series) # 找出最接近50欧姆的高度 impedance_mag np.abs(z_series - 50) best_idx np.argmin(impedance_mag) best_h heights_mm[best_idx] print(f在 3.5GHz 下最接近 50Ω 匹配的液态柱高度: {best_h:.1f} mm) print(f对应阻抗: 实部 {z_series[best_idx].real:.1f} Ω, 虚部 {z_series[best_idx].imag:.1f} Ω)这段代码用感应电动势法将液态单极子的长度映射到输入阻抗。best_h是首次匹配的参考高度实际工程中因为地平面寄生参数和封装的影响会存在 5%-10% 的偏差但这已经足够确定初始物理尺寸。注意Xm的计算用了对数近似在液态柱高度接近半波长时精度最好短于四分之一波长时误差变大后续仿真要用全波软件修正。3.2 从电流分布推导方向图液态柱高度对波束宽度的影响知道阻抗还不够ISAC 系统需要评估波束覆盖范围。单极子天线的方向图只取决于电流分布而电流分布又取决于电长度。下面这段代码画出不同液态柱高度对应的归一化方向图。import numpy as np import matplotlib.pyplot as plt def monopole_pattern(length_m, theta, freq_hz): 计算理想地平面上单极子的归一化辐射方向图 :param length_m: 液态金属柱高度 :param theta: 仰角0为天顶方向90度为水平方向弧度 :param freq_hz: 工作频率 :return: 归一化场强幅度 c 2.99792458e8 k 2 * np.pi * freq_hz / c L length_m # 单极子方向图公式 F(theta) [cos(kL cos(theta)) - cos(kL)] / sin(theta) # 分子分母在 theta0 或 pi 时同时趋于0需要特殊处理 with np.errstate(divideignore, invalidignore): F np.abs((np.cos(k * L * np.cos(theta)) - np.cos(k * L)) / np.sin(theta)) # 处理 theta 接近 0 或 pi 时的不定型此时方向图趋于0 F[np.isnan(F)] 0 F[F 1e-6] 0 # 消除数值噪声 # 归一化 if np.max(F) 0: F F / np.max(F) return F theta np.linspace(0, np.pi, 181) # 上半空间 0~180 度 freq 3.5e9 heights_to_plot [0.5 * (299.79 / 3500), 0.25 * (299.79 / 3500)] # 半波长和四分之一波长对应的高度 labels [0.5 波长, 0.25 波长] plt.figure(figsize(8, 5)) for H, lab in zip(heights_to_plot, labels): F_val monopole_pattern(H, theta, freq) plt.plot(np.degrees(theta), 20 * np.log10(F_val 1e-6), labellab) plt.xlabel(仰角 (度)) plt.ylabel(归一化方向图幅度 (dB)) plt.title(液态单极子高度对 ISAC 波束覆盖的影响) plt.ylim(-30, 0) plt.xlim(0, 180) plt.grid(True, alpha0.3) plt.legend() plt.show()运行这段脚本可以清楚看到液态柱为半波长时方向图在水平方向90度增益最大这是因为半波单极子等效为端射阵变到四分之一波长时波束明显展宽水平增益下降但覆盖范围更广。在 ISAC 系统里这就是感知搜素模式和通信精确指向模式之间的方向图切换机制。3.3 ISAC 激励下液态天线的联合波束赋形权重计算实际 ISAC 系统不只有单根天线。液态天线阵列需要和传统相控阵一样计算波束赋形权重只不过液态天线额外提供了频段和极化自由度。下面给出一个感知-通信联合权重生成的示例假设阵元间距为半波长目标在 30 度方向。import numpy as np def isac_beamformer(freq_hz, target_deg, num_elements8, d_ratio0.5): 生成 ISAC 联合波束赋形权重 :param freq_hz: 工作频率 :param target_deg: 目标方向角相对法线 :param num_elements: 阵元数 :param d_ratio: 阵元间距与波长比例 :return: 复数权重向量 c 2.99792458e8 lam c / freq_hz d d_ratio * lam # 相控阵导向向量 theta_rad np.radians(target_deg) steering np.exp(-1j * 2 * np.pi * d * np.arange(num_elements) * np.sin(theta_rad) / lam) # 液态天线带来的幅度扰动液态柱高度不一致引起的相位偏移 # 这里加入一个可通过液态驱动的修正项以 0.15*pi 为默认 liquid_phase np.array([0.0, 0.1, 0.2, 0.15, -0.05, -0.1, 0.0, 0.25]) * np.pi liquid_compensation np.exp(-1j * liquid_phase) weights steering * liquid_compensation # 归一化使总功率为1 weights weights / np.linalg.norm(weights) return weights # 通信目标角度 comm_weights isac_beamformer(3.5e9, target_deg30) # 感知目标角度更宽的覆盖范围例如 -20 度 sense_weights isac_beamformer(4.8e9, target_deg-20) print(通信模式权重 (3.5GHz):, np.round(comm_weights, 3)) print(感知模式权重 (4.8GHz):, np.round(sense_weights, 3))这段代码的意义在于它模拟了液态天线阵列中每个单元的独立调谐。由于封装工艺的误差液态柱高度不可能完全一致liquid_phase参数就是用来建模这种不一致性的。在实际系统中可以通过给每个液态单元加装位置传感器实时闭环补偿阵列相位。这是液态天线辅助集成感知与通信系统在波束赋形层面区别于传统相控阵的关键特征传统阵列相位是由移相器决定的而液态阵列可以通过物理变形来修正阵元之间的互耦影响。4. 液态天线辅助ISAC复现避坑从仿真到封装的5个教训4.1 镓致脆化导致金属结构断裂现象液态金属填充微通道后钎焊点或 PCB 走线在数小时内变脆稍微弯折就断裂。原因EGaIn 和 Galinstan 中的镓元素会沿铜、铝的晶界渗透生成镓铜合金力学强度急剧下降这在微电子封装领域是公认的失效机理。解决所有直接接触液态金属的金属部分改用钼、钨或钛制作这三个金属对镓的耐受性较好。PCB 上如果必须走铜线要用氮化硅或氧化铝涂层做物理隔离。4.2 液态金属氧化结皮导致阻抗漂移现象天线 S11 参数在静置数小时后发生变化谐振频率偏移超过 10%辐射效率下降。原因EGaIn 暴露在空气中会迅速形成一层厚度为纳米级的氧化皮主要成分是 Ga₂O₃这层氧化皮改变了电磁边界条件等效于在液态金属表面增加了一层损耗介质。解决注入液态金属时在手套箱内完成或者在封装腔内通入氮气保护。使用注射器加 PE 管方式填充确保液面不与空气接触。如果条件不允许唯一可行的办法是在液态金属表面覆盖一层硅油实测将阻抗漂移时间从小时级延长到数周级。4.3 电磁仿真不收敛且结果依赖网格尺寸现象在 HFSS/CST 中绘制液态金属天线剖分网格从 0.1mm 变到 0.05mm谐振频率结果变化超过 15%仿真始终不收敛。原因液态金属的介电常数不是简单实数它作为导体有趋肤深度问题。在毫米波频段趋肤深度很小但液态金属表面形状受表面张力影响存在微米级的倒角和曲面这导致网格剖分时在尖角处产生奇异点。解决在仿真中不要直接建模液体金属的自由曲面而是把它固定为圆柱体并且把圆柱体与微通道壁面的接触边设成半径为 0.02mm 的圆角。这一项改动能让收敛速度提升数倍结果也更接近实测。4.4 重力导致液态柱高度不稳定现象天线在水平放置时谐振频率正常垂直放置或倒置时谐振频率漂移且伴随机械晃动。原因液态金属密度在 6g/cm³ 左右重力引起的静压变化会导致液态柱形状改变特别是气压驱动方式下气液界面在失重或反置时会变形。解决在微通道内部填充多孔介质比如聚氨酯海绵利用毛细力锁定液态金属的位置。实测在 0.5mm 孔径的海绵中EGaIn 可以在任意姿态下保持形状这是目前最可靠的无运动部件方案。4.5 接触电阻不稳定现象天线效率在多次调谐后下降网络分析仪测量发现射频馈点处的阻抗时大时小。原因液态金属与金属触点之间的接触面积因为液体流动而变化接触电阻直接从毫欧级跳到欧姆级这部分能量损耗在回波损耗上几乎不可接受。解决改用螺纹机械触点让探针尖端嵌入液态金属内部。这里可以借鉴电化学测试中的三电极设计把铂丝加工成弹簧针形式插入液态柱中心位置。靠液体表面张力和机械力共同保持接触实测接触电阻稳定在 0.5Ω 以内。5. 进阶验证微波暗室里的相位校准与感知回波测试当仿真和封装都稳定之后液态天线辅助集成感知与通信系统要真正证明价值必须上暗室做三件事近场方向图扫描、S参数在线监测、以及感知-通信联合回波验证。前两项是天线常规测试这里重点讲感知回波验证。一个可行的验证方法是搭一套桌面级 ISAC 系统使用矢量网络分析仪和射频开关切换两副天线。液态天线作为收发共用天线配合一个角反射器。在感知模式下用 VNA 的时域门功能测量回波峰值在通信模式下切换信号源速率。关键指标是切换时间对感知中断的影响。切换液态天线的物理响应时间可以用 30ms 作为基准。FMCW 波形设计时chirp 时间需要小于液态天线稳定时间的一半才能避免回波信号展宽。这个参数跟液态金属粘度、微通道直径都有关系。我的实际经验是在验证液态天线的频段捷变能力时千万不要直接把 VNA 的电校准当成系统校准。液态天线切换时射频线的机械弯折和接头磨损会引入额外的幅度相位误差。最好每次切换后都做一次直通校准或者在测试链路里加入实时的参考通道归一化。这算不上什么高深理论但真能省下一天时间。希望帮到你。本文还有配套的精品资源点击获取