UXO
C13 / 电磁场与磁性目标 · 2.0 · 参考阅读 6 分钟,推导与练习另计

静磁极化张量:把材料、形状和姿态连起来

严格区分输入 H 与 B,推导主轴响应和实验坐标变换。

展开本章阅读路线(8 节)

一个标量为什么不够用?

球体沿每个方向响应相同,可以用一个比例系数描述。但椭球沿长轴与短轴不同,外场换方向后,磁矩甚至不再与外场平行。张量就是同时记录“各方向输入如何产生各方向输出”的线性映射。

对均匀各向同性线性椭球,把退磁关系代入材料关系:

M=χ(H0NM),(I+χN)M=χH0.\mathbf M=\chi(\mathbf H_0-N\mathbf M),\quad (I+\chi N)\mathbf M=\chi\mathbf H_0.

乘体积 V 后定义本课程的静磁极化张量 P:

m=PH0,P=V(I+χN)1χI.\mathbf m=P\mathbf H_0,\qquad P=V(I+\chi N)^{-1}\chi I.

P 的单位为 m³,因为 m 的单位 A·m²,H 的单位 A/m。有些文献用 B 作为输入,此时对应张量是 P/μ₀,不能在两种定义之间只改字母不改系数。

在主轴上读出物理含义

主轴坐标下 Pi=Vχ/(1+χNi)P_i=V\chi/(1+\chi N_i)。V 描述尺寸尺度,χ 描述材料线性响应,Nᵢ描述形状反馈。P 是对称矩阵,因此其特征向量给出正交主方向;球体三个特征值相同,没有唯一可辨识主轴。

若 R 把目标坐标分量变到实验坐标,先将外场变回目标系,再施加目标响应,最后把磁矩变到实验系:

mlab=RPbodyRTH0,lab,Bs=G(xx0)RPbodyRTH0.\mathbf m_{lab}=R P_{body}R^T\mathbf H_{0,lab},\quad \mathbf B_s=G(\mathbf x-\mathbf x_0)R P_{body}R^T\mathbf H_0.

这条顺序解释了为什么不是简单的 RP。接收方向还要左乘 nT\mathbf n^T,得到单轴通道读数。

一个二维旋转例子

设 P 的两个平面主值为 2、1(单位相同),旋转 45° 后该平面块为 [1.50.50.51.5]\begin{bmatrix}1.5&0.5\\0.5&1.5\end{bmatrix}。沿实验 x 方向激励时,会同时产生 x、y 方向磁矩。非对角项并不意味着材料具有奇怪的新机制,它可能仅仅来自坐标旋转。

从自洽磁化到真正的张量表达

在均匀椭球、线性材料且无剩磁时,M=χ(H0NM)M=\chi(H_0-NM),所以

P=V(I+χN)1χ,m=PH0.\mathbf P=V(\mathbf I+\boldsymbol\chi\mathbf N)^{-1}\boldsymbol\chi,\qquad\mathbf m=\mathbf P\mathbf H_0.

当 χ 可逆时也可写 P=V(χ1+N)1P=V(\chi^{-1}+N)^{-1}。这个形式显示:若材料磁化率和退磁张量对称,则 P 对称。不能为了方便随意交换两个非对易矩阵的位置。各向同性 χ 退化为标量,才得到熟悉的三个主轴分式。

从六个参数到一条观测行

对称张量取参数向量 p=(Pxx,Pyy,Pzz,Pxy,Pxz,Pyz)Tp=(P_{xx},P_{yy},P_{zz},P_{xy},P_{xz},P_{yz})^T。定义接收传播行向量 gT=eTGg^T=e^TG 和局部激励 h,则标量数据为

d=(gxhxgyhygzhzgxhy+gyhxgxhz+gzhxgyhz+gzhy)p.d=\begin{pmatrix}g_xh_x&g_yh_y&g_zh_z&g_xh_y+g_yh_x&g_xh_z+g_zh_x&g_yh_z+g_zh_y\end{pmatrix}\mathbf p.

每一个实际发射—接收组合提供设计矩阵的一行。非对角项没有统一的额外2;只有 g=h 等特殊情况才简化。若激励永远只沿 z 方向,则 Pxx、Pyy、Pxy 不进入这一观测,重复采样不能恢复它们。

旋转后的不变量与变化量

世界系张量为 Pw=RPbRTP_w=RP_bR^T。迹、行列式和特征值不随坐标旋转变化;各矩阵元素会变化。这些不变量对目标表征很有用,但同时丢掉了方向信息。后续应保留张量估计的不确定性,而不是仅把三条排序曲线当作无误差“指纹”。

深入练习:单测点获得六条数据,就一定能反演六个张量元素吗?

不一定。六条观测行可能线性相关。必须检查设计矩阵的秩和噪声加权后的奇异值;条件数极大时,形式上的满秩也可能没有实用分辨率。

研究拓展

对各向异性材料,起点变为 M=χ(H0NM)\mathbf M=\boldsymbol\chi(\mathbf H_0-N\mathbf M),应保留矩阵乘法顺序。对于导电目标的频域 MPT,实部和虚部一般不保证能由同一个 R 同时对角化;把静磁椭球的共轴结构推广到所有目标必须另加依据。

练习与解析

1. 球体旋转后 P 会改变吗?

不会。P=pI,因此 RPRᵀ=pI。数据无法从这一理想球体响应识别姿态。

2. 为什么单一方向激励通常不足以恢复全部 P?

一次激励只探测 P 对该向量的作用。即便完整知道输出磁矩,也没有测试另外的独立输入方向,需要更多独立激励与足够观测几何。

把推导变成可观察的变化

旋转共轴双指数模型,查看实验坐标中的对角与非对角响应。

打开实验:张量旋转:为什么出现交叉分量 ↗

参考文献与继续学习

外部链接需联网,部分论文全文需要机构访问。页面内容为原创教学解释;文献用于核对理论和进一步阅读。资源链接核对日期:2026-09-12。

把这一章变成自己的理解

笔记保存在当前浏览器本机;建议定期导出备份。
SEARCH / 全文搜索