用途
该技能用于估算电磁场在介质组织体模内部的比吸收率（SAR），面向射频与天线工程中的合规评估、组织吸收与功率沉积分析。适用场景包括：计算天线在生物组织附近或内部的 SAR；评估 FCC/ICNIRP 射频暴露合规性；建模鸟笼线圈用于 MRI SAR 分析；计算点 SAR 或质量平均 SAR（1g/10g）；设计植入式天线并估算 SAR。
不适用场景：设计天线本身（应使用 matlab-design-antenna）；RF 传播/覆盖分析（应使用 matlab-analyze-rf-propagation）；RCS 或散射分析（应使用 matlab-analyze-rcs）。
SAR 公式
SAR = σ|E|² / (2ρ)，单位 W/kg。其中 σ 为组织电导率（S/m），|E| 为组织内电场幅值（V/m），ρ 为组织质量密度（kg/m³）。电导率与损耗正切的关系为：σ = ω ε₀ ε_r tanδ。
三种方法
方法一：birdcage + Phantom（完整组织属性）
birdcage 天线是唯一支持体介质体的对象，通过 Phantom 属性将体模直接纳入 MoM 求解器。
Phantom 数据格式：
phantom = struct( ...
Points=vertices, ... % N x 3 顶点坐标（米）
Tetrahedra=elements, ... % M x 4 四面体单元连接
EpsilonR=epsR, ... % 相对介电常数（标量）
LossTangent=tanD); % 介质损耗正切（标量，无上限）
随附体模：humanheadcoarse.mat（584 顶点，2818 四面体，较快）；humanheadfine.mat（更高分辨率，更精确但更慢）。两者提供变量 P（顶点）与 T（四面体），需应用 scaleFactor = 0.003 得到物理尺寸。
工作流示例（128 MHz，3T MRI）：设置组织属性（灰质：EpsilonR=77，Sigma=0.51 S/m，Rho=1040 kg/m³），由 sigma 换算 tanD；加载体模并缩放；用 birdcage 创建天线并 show；用 createObservationGrid 生成观测网格；调用 EHfields 获取电场；按公式计算点 SAR。
方法二：conformalArray + shape.Custom3D（任意天线）
使用 conformalArray 将任意天线与介质头模并列放置。shape.Custom3D 作为无源介质散射体参与 MoM 问题，不需要馈电，可直接放入 conformalArray.Element，无需包装为 customAntenna。
工作流示例（2.4 GHz）：加载并缩放头模（缩放到 180mm 宽度）；用 triangulation 与 freeBoundary 提取表面网格；创建 headShape = shape.Custom3D(headTri)，设置 Dielectric = "TMM10"；设计偶极子；构建 conformalArray，将 {d, headShape} 作为 Element，设置 ElementPosition 与 Reference；计算电场并按公式计算 SAR（使用 TMM10 属性：epsR=9.8，tanD=0.0022，rho=2270）。
观测网格生成：通过 createObservationGrid 对三维网格过滤，仅保留头体积内的点；使用 delaunayTriangulation 与 pointLocation 判断点是否在体内。注意 EHfields 期望 3×M 矩阵（列为点），因此传入 obsPoints'。
功率归一化：默认馈电电压为 1V。用 impedance 获取输入阻抗（birdcage 多端口时取 Z_in(1)），P_accepted = 0.5*real(Z_in)/abs(Z_in)²，再缩放到 1W 接受功率。
质量平均 SAR（10g/1g）：取平均质量 0.010 kg（ICNIRP）或 0.001 kg（FCC），由 avgMass/tissueRho 开立方得到立方体边长，对每个点在其半边长范围内取均值。
功率平衡校验：方法一用 efficiency 得到 P_absorbed_eff = P_accepted*(1-eta)；方法二用 SAR 积分得到 P_absorbed_SAR = sum(SAR_point)tissueRhodeltaV。相对误差 <20% 为 PASS，20–50% 为 MARGINAL，>50% 需细化网格。
可视化：用 scatter3 展示观测点与 SAR_1W，添加 colormap、colorbar 与轴标签，view(30,20)。
介质目录与限制
DielectricCatalog 可通过 add(catalog, "BrainTissue900MHz", 52, 0.36, 900e6) 添加材料。但求解器实例化材料时会调用 dielectric("name")，强制 LossTangent ≤ 0.03，超过该上限的自定义材料无法用于分析。dielectric 类本身也强制 LossTangent ≤ 0.03，真实组织（如 900 MHz 脑组织 tanδ≈0.36）超过该上限；只有 birdcage 的 Phantom 通过自定义 struct 绕过此限制，直接 EHfields 方法则通过在 SAR 公式中应用组织电导率来绕开该上限。
可用目录材料：TMM10（εr=9.8，tanδ=0.0022，最高介电常数）、TMM6（6.3，0.0023）、FR4（4.8，0.026，最高损耗正切）、Teflon（2.1，0.0002，极低损耗）。SAR 演示建议使用 TMM10。
组织属性参考
脑（128 MHz）：εr=77，σ=0.51 S/m，tanδ=0.93，ρ=1040 kg/m³；脑（900 MHz）：52，0.94，0.36，1040；脑（2.4 GHz）：48，1.8，0.28，1040；肌肉（900 MHz）：55，0.94，0.34，1040；皮肤（900 MHz）：41，0.87，0.42，1100。
法规 SAR 限值
FCC（美国）：1.6 W/kg，1g，头部/身体；ICNIRP（欧盟）：2.0 W/kg，10g，头部/躯干；ICNIRP（欧盟）：4.0 W/kg，10g，四肢；IEC 60601-2-33（MRI）：3.2 W/kg，10g，头部。
方法三：植入式天线（直接 EHfields）
对于嵌入组织内的天线（植入式医疗设备），conformalArray 不可用，因为求解器拒绝相交几何。解决方法是使用独立天线的 EHfields，再在后处理中套用 SAR 公式。
限制：EM 求解器按自由空间辐射计算场，未捕捉组织对天线阻抗的加载；组织电导率仅在 SAR 公式中应用。这会低估失谐效应，但能正确演示 SAR 工作流。
conformalArray 失败原因：其在网格划分时调用 checkIntersection，若天线几何与 shape.Custom3D 组织体相交，会抛出 "Intersection detected in specified geometry."。这是硬性约束，天线必须位于介质体之外。
方法选择
使用 birdcage + Phantom：建模 MRI 线圈 SAR；需要真实组织属性（高 LossTangent）；天线为鸟笼线圈。
使用 conformalArray + shape.Custom3D：建模任意天线类型的手机/设备 SAR；天线位于组织体之外（非嵌入）；可接受 LossTangent 上限（0.03）。
使用直接 EHfields（方法三）：天线植入组织内；任意天线类型（pcbStack、目录、阵列）；组织电导率在后处理中应用（真实值，无上限）。
频率解析
支持 MHz、GHz、Hz，无单位时默认 Hz。常见 SAR 频率：128 MHz（3T MRI）、900 MHz（GSM）、1.8 GHz（LTE）、2.4 GHz（Wi-Fi）。按参考表使用频率相关的组织属性。
编码规范
4 空格缩进，变量使用 lowerCamelCase，Name-Value 参数使用 UpperCamelCase；字符串使用双引号；Antenna Toolbox 绘图（show、pattern、impedance、rfplot）不加标题；手动 plot() 与 scatter3() 图形加标题；用 fprintf 输出格式化数值；所有图形显示在独立 figure 中；输出包含单位（米、欧姆、W/kg、dB）。
使用准则
始终用功率平衡检查验证 SAR 结果；始终归一化到已知输入功率（通常 1W 接受功率）；使用 EHfields 时传入转置后的点（3×M，不是 M×3）；birdcage 取 Z_in(1)（impedance 返回多端口向量）；网格间距影响精度，快速演示用 20mm，发表质量用 5mm；除非用户明确需要 birdcage/MRI SAR，默认使用 conformalArray 方法；用户请求组织属性时提示 LossTangent 上限；说明 TMM10 的局限（演示工作流但低估真实组织吸收）；不要过度解释电磁理论。
版权 2026 The MathWorks, Inc.