简介：电磁仿真CST场分布MATLAB后处理代码定位于CST仿真后的数据可视化工具，主要面向天线设计、微波器件分析等场景的电磁工程师与相关专业学生，用于解决场监视器导出的ASCII数据难直接解读的问题。压缩包内包含1个.m脚本文件，整体仅约955B，体积轻量、便于快速部署，只需将CST场监视器导出文件放入同一目录即可运行。脚本通过textscan完成数据导入与时间步提取，利用surf、contourf、quiver等函数绘制二维/三维幅值场与矢量方向图，支持复数场数据的相位分布呈现，并可进一步计算场均值、峰值或扩展频率特性分析，帮助用户直观对比不同设计条件下的场变化，定位潜在问题并辅助优化。目前已有8675人学习使用，适合希望用MATLAB对CST结果进行高效后处理与定制化分析的工程技术人员。
1. 为什么我劝你别在 CST 里直接截图：场分布后处理的真实工作量
做电磁仿真的人只要碰过 CST 场分布，就迟早会走到写 MATLAB 后处理代码这一步。CST 自带的 2D/3D 场可视化适合人工判断“哪里强、哪里弱”，但当你需要跨频点比较场分布峰值、沿某条路径提取场强曲线、把近场数据换算成功率密度，或者把多组仿真结果放进同一套标准里做量化分析时，截图完全不够用，手里得握着坐标和场值的原始数据。这篇笔记就是围绕“电磁仿真 CST 场分布 + MATLAB 后处理代码”这条主线，讲清楚从 CST 导出场数据开始，到 MATLAB 里完成解析、插值、绘图和验证的全过程。适合已经用 CST 跑完仿真、但被数据导出和后处理卡住的微波与天线工程师，也适合正在搭自家仿真后处理流程的团队。
2. 从 CST 导出场分布数据：ASCII、CSV 与二进制的取舍
CST 的场监视器算完之后，数据并不会自动给你一份干干净净的表格，需要你主动导出。而这第一步选什么格式、怎么设参数，直接决定后面 MATLAB 解析脚本的复杂度和运行速度。我见过太多人在这里随手点了个导出，回头在 MATLAB 里解析时才发现列顺序不对、单位不对、注释格式不对，来回折腾半天。这一章先把导出这步的决策讲透。
2.1 三种导出格式在 MATLAB 侧的代价对比
CST 里能导出的场数据格式不止一种，常见的至少包括 ASCII 文本、CSV、二进制 .bin，以及面向第三方可视化的 VTP。对 MATLAB 后处理来说，它们各有各的脾气，我按实际使用频率排个序。
ASCII 文本是最稳的选择。CST 导出的 .txt 文件通常是按行组织的列数据，每行包含 x、y、z 坐标和若干场分量，MATLAB 的 textscan 或 readmatrix 都能直接处理，不需要额外工具箱，遇到格式问题也能打开文件肉眼排查。代价是文件体积大，网格密集时单个文件动辄上百 MB，读取慢，而且列顺序在导出模板之间可能有差异。
CSV 本质上是 ASCII 的变种，只是分隔符固定为逗号，readmatrix 可以直接吃，脚本更简洁。它的缺点和 ASCII 类似，文件也偏大。二进制 .bin 读取最快、体积最小，但前提是你完全清楚文件内部的数据排布——CST 导出的二进制数据没有统一标准，不同版本和模板的排列可能不同，一旦读错就是整片数据错乱，而且错乱后很难靠肉眼发现。VTP 是为 ParaView 生态准备的，MATLAB 侧读写要写额外的解析代码，除非你已经引入 VTK 类库，否则我不建议在这条路上去找麻烦。
对大多数天线和微波无源器件项目，我的默认方案是：网格总量在百万点以内用 ASCII 或 CSV，百万点以上优先检查能否把导出区域裁剪小一点，而不是直接上二进制。二进制留给确实跑不动大文件、且你愿意花一个下午调试读取函数的场景。
2.2 导出前必设的三个参数：单位、坐标系与监视器频率
CST 场导出的对话框里通常有几个直接影响 MATLAB 侧数据的参数，我每次导出前都会强制自己检查一遍。
第一个是单位。CST 建模默认使用 mm，但场导出的单位选项可能最终以 m 输出，也可能保留 mm，取决于你选择的导出模板和当时激活的单位体系。单位不一致会让所有坐标和场分布整体偏移一到三个数量级，在 MATLAB 里画出的图看起来形状完全正确，就是尺寸奇奇怪怪，这是最隐蔽的坑之一。我建议导出前把单位固定成你建模时用的单位，并在文件名或者数据说明里标注清楚。
第二个是坐标系。CST 支持全局坐标和局部坐标，阵列模型或者带子坐标系的结构里，如果选成局部坐标，导出的场分布位置会和实际模型错位。后处理阶段我统一用全局坐标，局部坐标只在极少数专门分析子结构的场景里才用。
第三个是频率监视器的选择。CST 里场监视器是按频点设置的，导出时你要指明导出哪几个监视器结果，如果模型有多个频点，批量导出比一个个手动点要可靠得多。CST 的宏录制功能可以把你手动导出一次的操作录下来，然后批量回放，我一般先把单个监视器的导出参数调好，再录宏跑完所有频点。
下面是几个设置项的建议值，按我常用的工程配置整理：
| 设置项 | 推荐值 | 原因 |
|---|---|---|
| 导出单位 | 与建模单位一致（常用 mm） | 避免坐标整体缩放 |
| 坐标系 | 全局坐标 | 保证与模型装配位置一致 |
| 数据列顺序 | 坐标在前，实部虚部成对出现 | 匹配 MATLAB 解析脚本 |
| 频率监视器 | 按频点批量导出 | 避免手动操作漏项 |
| 导出区域 | 尽量裁剪到目标区域 | 控制文件体积和读取时长 |
2.3 导出后的文件体检：用 5 分钟看清文件头，省后面 50 分钟
拿到 CST 导出的文件后，我强烈建议先做一次“体检”再写解析脚本。做法很简单：用任意文本编辑器打开文件，看前 20 行内容。真正常见的坑都写在文件头里——注释标记是 # 还是 %、每行几个数、坐标单位有没有写进注释、列顺序到底是不是你预期的那样。
在 MATLAB 里快速看文件头可以这样写：
filename = 'E_field_2.4GHz.txt'; fid = fopen(filename, 'r'); for k = 1:20 line = fgetl(fid); if ~ischar(line) break; end fprintf('第 %d 行: %s\n', k, line); end fclose(fid);这段代码的逻辑是：打开文件后循环读取前 20 行并打印到命令行，看到文件实际内容后再决定解析脚本怎么写。fgetl 每次读一行，读到文件末尾会返回 -1，所以用 ischar 做判断来跳出循环。参数说明：20 这个数字一般够用，因为 CST 的注释头通常不超过 10 行；如果文件头特别长，可以改成 50。
体检时我主要确认三件事。第一，是否有注释行、注释符号是什么，这决定脚本里用不用 CommentStyle 跳过注释。第二，每行数据的列数，如果模板输出 9 列但脚本按 13 列解析，后面全乱。第三，确认坐标列是第几到第几列、场分量是实部虚部还是幅度相位，这两者的处理方式完全不同。我曾经拿到一份模板输出的是幅度和相位，没做体检直接按实部虚部解析，结果画出来的场图破破烂烂，检查了半小时才发现文件头注释里写着数据格式是 Magnitude/Phase。
3. MATLAB 读数与数据结构化：CST 的离散场装进可用变量
CST 导出文件的解析是整个后处理流程的地基。很多新手拿到文件后直接 load 或者 importdata 一把梭，数据是读进来了，但不知道每列是什么含义，后面每一步都在猜。这一章我会给出一个经过反复使用的读取函数，并把数据组织成结构体，让后续的插值和绘图代码不用再关心原始文件的格式问题。
3.1 写一个通用读函数：支持注释跳过与列数校验
先给一个我自己一直在用的 CST 场数据读取函数。它不像 readmatrix 那样一行全搞定，但胜在能把错误提示做清楚，什么时候卡住一眼就能看出来。
function fd = readCstField(filename, nCols) % 读取 CST 导出的 ASCII 场分布文件 % 输入: % filename - CST 导出的文件名，支持 .txt 或 .csv % nCols - 每行数据列数，常见为 9 或 13 % 9 列: x y z Re(Ex) Im(Ex) Re(Ey) Im(Ey) Re(Ez) Im(Ez) % 13列: x y z Re(Ex) Im(Ex) Re(Ey) Im(Ey) Re(Ez) Im(Ez) |E| theta phi % 输出: % fd - 结构体，包含 x/y/z 坐标与复数场分量 Ex/Ey/Ez % % 示例: % fd = readCstField('E_field_2.4GHz.txt', 9); fid = fopen(filename, 'r'); if fid == -1 error('无法打开文件: %s', filename); end % 用 textscan 读取全部数值列，CST 注释通常在行首有 # 或 % raw = textscan(fid, repmat('%f', 1, nCols), ... 'Delimiter', {' ', '\t', ','}, ... 'CommentStyle', {'#', '%'}, ... 'MultipleDelimsAsOne', true); fclose(fid); % 检查读到的数据是否为空 if isempty(raw{1}) error('文件为空或格式无法解析: %s', filename); end % 校验实际读到的列数是否符合预期 actualCols = size(raw{1}, 1); if actualCols ~= nCols error('列数不匹配: 期望 %d 列, 实际读到 %d 列', nCols, actualCols); end % 组装输出结构体 fd.x = raw{1}; fd.y = raw{2}; fd.z = raw{3}; if nCols == 9 || nCols >= 13 fd.Ex = complex(raw{4}, raw{5}); fd.Ey = complex(raw{6}, raw{7}); fd.Ez = complex(raw{8}, raw{9}); else error('目前仅支持 9 列或 13 列格式，请检查导出模板'); end end这段代码的核心是 textscan 的配置。Delimiter 同时指定空格、制表符和逗号三种分隔符，覆盖 CST 不同导出模板的习惯；CommentStyle 指定 # 和 % 开头的行自动跳过，避免注释混入数据；MultipleDelimsAsOne 为 true 可以把连续多个空格当作一个分隔符处理，防止空列影响解析。
参数说明里最需要注意的是 nCols 这一项。9 列数据对应坐标加 E 场三方向实部虚部，这是最常见的模板。13 列通常额外加了总场幅度和两个相位角，读取时我们只取前 9 列，后面多出来的列不影响。如果以后碰到 6 列格式（只导出了某个单一分量），直接把分支加进 if 判断即可。这个函数的代价是每次调用都会读取一遍文件，大文件场景下会有点慢，但作为后处理脚本足够可靠。
3.2 用结构体组织场分量数据：坐标和场值不再互相走散
函数里已经把数据放进了结构体，这里单独说一下为什么不用多个独立变量。电磁场数据的特点是坐标和分量强关联，如果拆成 x、y、z、Ex、Ey、Ez 六个向量放在工作区里，后续传给插值函数时得一个一个当参数传，很容易传错顺序。用结构体一次打包，插值、绘图、切片时都是 fd 开头的字段引用，代码可读性和可维护性都高一个档次。
结构体里的 Ex 是复数类型，这是刻意为之。CST 导出的实际场是复数，包含幅度和相位信息，很多后处理操作比如计算功率密度、合成圆极化场、做相控阵加权，都需要保留相位。如果读取时直接把实部虚部分开存成两个实数向量，后面用到相位时就要重新组装，白白多写代码。
一个建议是读取后立刻给 fd 结构体补上频点信息。频点不在 CST 导出的数据列里，但多频点处理时，它是最容易混淆的元数据。手动标记一行代码，后面批量处理时能减少很多糊涂账。
fd.freq = 2.4e9; % 手动标记，单位 Hz fd.unit = 'mm'; % 标记导出单位，后处理换算时用3.3 多监视器批量读入：一个循环读完所有频点
实际工程项目里通常不会只导出一个频点，尤其是天线调试阶段，一次仿真跑五个频点很常见。批量读取的套路不复杂，核心是用 dir 匹配文件名模式，再循环调用读取函数。
files = dir(fullfile('.', 'E_field_*GHz.txt')); fdArr = cell(length(files), 1); for k = 1:length(files) fdArr{k} = readCstField(fullfile(files(k).folder, files(k).name), 9); % 从文件名里解析频率，比如 'E_field_2.4GHz.txt' 取 2.4 nameParts = strsplit(files(k).name, {'_', 'G'}); fdArr{k}.freq = str2double(nameParts{3}) * 1e9; fprintf('已读取 %s，共 %d 个采样点\n', files(k).name, length(fdArr{k}.x)); end这里用 cell 数组而不是 struct 数组，因为每个频点的采样点数不一样，struct 数组在字段长度不一致时处理起来很别扭。频率解析用 strsplit 按下划线和 G 切分，文件名里 E_field_2.4GHz.txt 切出来第 3 段就是 2.4，转成 double 再乘 1e9 就是赫兹。如果你的文件名命名规则不同，按实际情况调整索引即可。
读取完这一步，数据的组织工作就算全部完成。后续不管是对比多频点峰值，还是做插值切片，都只需要从 fdArr 里按索引取数据就行了。
4. 把离散场画出来：规则网格插值、切片图与沿线取场
CST 导出的场数据是定义在网格节点上的，网格本身不一定规则，尤其是四面体网格或者自适应加密过的模型，节点分布完全没有规律。直接拿这些散点画图，PATCH 或 scatter3 也能出图，但做切片、做沿线提取、做后处理数值计算都很别扭。这一章讲怎么把散点场变成规则网格场，然后做三种最常见也是最有用的可视化：切片图、路径曲线和 dB 动态范围图。
4.1 用 scatteredInterpolant 把散点场插值到均匀网格
在 MATLAB 里把无规则散点插值到规则网格，有 griddata 和 scatteredInterpolant 两条路。griddata 是老牌函数，用起来直接，但每次调用都要重新构建三角剖分。scatteredInterpolant 可以一次性创建插值对象，然后对任意多组查询点求值，效率明显更高，代码也更清晰。我一般固定用 scatteredInterpolant。
% 取原始数据里的 E 场幅度 magEx = abs(fd.Ex); % 创建插值对象，插值方法 linear，外推设为 none（不做外推） Fex = scatteredInterpolant(fd.x, fd.y, fd.z, magEx, 'linear', 'none'); % 生成规则网格坐标，网格密度按需调整 Nx = 80; Ny = 80; Nz = 40; xg = linspace(min(fd.x), max(fd.x), Nx); yg = linspace(min(fd.y), max(fd.y), Ny); zg = linspace(min(fd.z), max(fd.z), Nz); [XX, YY, ZZ] = ndgrid(xg, yg, zg); % 对规则网格逐点求插值后的场值 EXq = Fex(XX, YY, ZZ);这段代码里最值得注意的参数是 scatteredInterpolant 的第三个输入 'none'。它表示查询点落在已知点凸包之外时不进行外推，直接返回 NaN。这一点在电磁场里特别重要，因为场分布只存在于求解域内部，求解域外没有物理意义，如果默认外推，边界外会冒出各种奇怪数值，后面画图时会出现一大片蓝色色块覆盖到空气区域。
网格密度 Nx、Ny、Nz 的取值也直接影响结果质量。网格太疏会把峰值抹平，网格太密则计算量和内存占用都会大幅上升。我通常先按 80×80 起步，看一下峰值位置是否与 CST 原始视图一致，峰值偏差超过 2% 就加大密度。对内存敏感的大模型，可以先只对关注区域做局部加密网格。
4.2 用 slice 画三个方向的切片图：先看整体再找热点
拿到规则网格化的数据后，slice 函数是最直接的切片工具。它可以一次性画出 X、Y、Z 三个方向的任意切片，用来快速评估场分布的整体形态非常高效。
figure; % 在三个位置各切一刀，空数组代表该方向不切 slice(XX, YY, ZZ, EXq, [], [], zg(20)); colorbar; colormap(jet); shading interp; xlabel('x (mm)'); ylabel('y (mm)'); zlabel('z (mm)'); title('|Ex| 场分布切片');slice 的参数格式是 slice(X, Y, Z, V, sx, sy, sz)，其中 sx、sy、sz 可以传标量坐标、向量坐标或空数组。传空数组表示该方向不切。这段代码传了 zg(20)，意思是只在 Z 方向上第 20 个网格点处切一刀水平面，适合看某个高度截面上的场分布。如果你需要看器件中心截面，可以把 zg(20) 换成 0 或者其他具体坐标值。
切片图在工程里的价值是快速定位热点区域。看到整体分布之后再针对热点所在坐标做沿线提取或者局部加密看细节，这种先粗后细的做法比一开始就对整个模型做精细可视化高效得多。
4.3 沿两点路径提取场强曲线：传输线和波端口分析必备
天线设计和滤波器调试里，经常需要看场沿某条线的变化趋势，比如沿微带线中心线看电场幅度、穿过谐振腔中心看磁场衰减。这个操作本质上就是在三维空间里做线段参数化，然后用插值对象逐点求值。
% 定义路径起点和终点，单位与场数据一致 x1 = -20; y1 = 0; z1 = 5; x2 = 20; y2 = 0; z2 = 5; % 沿直线取 101 个点 t = linspace(0, 1, 101); px = x1 + (x2 - x1) * t; py = y1 + (y2 - y1) * t; pz = z1 + (z2 - z1) * t; % 计算路径上各点的插值场强 Eline = Fex(px, py, pz); % 转成 dB 显示，加 eps 避免取对数时出现负无穷 ElineDB = 20 * log10(Eline + eps); figure; plot(t, ElineDB); xlabel('沿线归一化距离'); ylabel('|Ex| (dB)'); grid on;这里用参数 t 在 0 到 1 之间均匀取值，把线段分成 101 个点。101 这个数字是经验值，在绝大多数场景下足够平滑，计算开销也小。如果场强沿路径变化剧烈，可以加大到 201 或 401，但要注意路径点太密对插值对象来说也只是多算几十次，不会有什么压力。
dB 转换时加 eps 是一个实用细节。插值结果里如果恰好出现 0 值，log10(0) 会变成负无穷，整条曲线在 MATLAB 里显示时会断裂。加上 eps，数值上几乎不改变真实结果，但能保证曲线连续。
4.4 dB 刻度与动态范围设置：为什么线性色标往往会骗人
电磁场幅度动态范围经常跨好几个数量级，线性色标下，最大值旁边全是同一个颜色，热点细节根本看不出来。dB 刻度是行业默认选择，但直接 20*log10 之后会出现另一个问题——极小值区域会拉低整个色标范围，让真正关心的那部分场分布对比度变差。
工程里常见的做法是设定一个显示下限，把低于这个下限的 dB 值全部裁剪成下限值，这样色标从下限到 0 dB 就能覆盖真正有意义的范围。
% 设定动态范围 40 dB floorDB = -40; ElineDB_clipped = max(ElineDB, floorDB); % 三维切片显示也可以同样处理 EXqDB = 20 * log10(EXq + eps); EXqDB_clipped = max(EXqDB, floorDB);floorDB 取多少取决于你关心的区域和背景噪声的差距。大多数场监视器结果里，-40 dB 已经能看清主路径上的场分布；如果是高 Q 值腔体，动态范围可能要拉到 -60 dB 以上才能看到弱场区域的泄漏路径。这个值不是死的，建议先用 -40 看整体，再逐步往下调，直到图中的细节足够丰富又不会出现大片噪声斑块。
5. CST 场分布后处理避坑：五个常见问题与排查方法
这一章写的全是我在真实项目里踩过或者帮别人排查过的坑，每一个都对应着具体现象和具体解法。做 CST 场分布后处理，算法本身不难，难的是这些藏在数据背后的“暗信息”。我把最常见的五类问题按频率排个序，每条都按“现象 → 原因 → 解决”的路子来讲。
5.1 坐标整体偏移或缩放：毫米与米混用的典型症状
有次处理某滤波器的场分布，数据读进来后画切片图，怎么看怎么不对劲——器件在 CST 里明明只有几十毫米，画出来却是几十米的尺寸，整体形状倒是没变。反复核对导出设置后发现，文件头注释里写着单位是 m，但我脚本里假设它是 mm，坐标对应的场点位置全都偏移了三个数量级。
原因就是导出时单位没固定。CST 项目单位可能保持 mm，但部分导出模板会用 SI 基础单位把坐标输出成米。解决方法是每次导出前把单位固定成建模单位，并在读函数里加一道单位校验。更保险的做法是，在结构体里显式记录单位，后处理时如果需要换算长度，全部通过这个字段换算，而不是在代码里写死。
if fd.unit == "m" fd.x = fd.x * 1000; % 转成 mm，方便后续处理 fd.y = fd.y * 1000; fd.z = fd.z * 1000; fd.unit = "mm"; end这个换算虽然简单，但能省掉后续所有绘图和几何操作里的单位疑云。
5.2 场分量顺序不对导致图像撕裂或出现异常条纹
第二种高频问题是文件列顺序和脚本预期不一致。CST 导出模板有好几种，有的按 Ex、Ey、Ez 的实部虚部排列，有的会把幅度和相位放在实部虚部后面，还有个别版本的模板输出顺序是 Ez 排在最前面。如果你拿 9 列脚本去读一个 13 列文件，解析出来的 Ex 实际是原来的 Ey，后续所有计算都是错的。
这类问题的现象非常迷惑人。场的量级是对的，峰值位置也对，但偏振方向乱了，画出来的图有强烈的方向性条纹。排查方法就是回到我第 2 章说的文件体检，看文件头注释里对每列的描述，然后对照脚本里的列定义逐条核对。
5.3 金属边界处场强异常大：奇点还是插值污染
金属边缘和拐角处，电磁场往往有几何奇异性，CST 网格加密能逼近真实分布，但导出的节点数据在边界附近数值很大。如果直接用 scatteredInterpolant 对整个区域的散点做插值，奇点附近的场会被“抹”到相邻区域，造成热点扩大一圈的假象。
我碰到过一例：一个微带贴片的边缘场强在 CST 原始视图里是很尖锐的窄峰，插值之后变成了一整片高场区。起初以为是插值密度不够，后来才发现是插值把边界奇点的值扩散到了空气一侧。解决方法是把求解域外的点和边界外扩区域的点全部截掉，只对内部点做插值。判断哪些点该截掉，可以用 CST 导出的区域范围配合模型几何边界做一个粗掩膜，或者直接缩小插值范围到不含金属边界的子区域。
5.4 复数共轭约定不一致导致功率方向异常
时谐场的相位约定在不同软件、不同导出模板里可能有差异。CST 默认使用某个相位约定，如果你在 MATLAB 里计算坡印廷矢量时用了相反的共轭顺序，得到的功率流向和实虚部符号都会出问题，典型表现是高 Q 结构里计算出的净功率为负值，或者损耗功率大于输入功率。
这个问题最隐蔽，因为单个标量看没有异常，只有算功率、算能量流这类二次量时才暴露。处理思路是：先确认你的场数据代表的是相量还是瞬时值，再确认导出模板文档里的相位约定，最后在做功率计算前，用已知的简单模型比如矩形波导主模做一次标定。标定时让波导一端输入已知功率，在 MATLAB 里对截面做坡印廷矢量积分，看结果与 CST 报告的波端口功率是否一致，不一致就把共轭取法反过来再测。
5.5 对称边界下数据只导出一半：合并与镜像的麻烦
用了对称边界条件的模型，CST 求解域只有整个结构的四分之一或一半，导出的场分布自然也只有这部分。很多人做后处理时忘了这回事，直接拿半结构数据去算整个器件的总功率或总损耗，结果差一倍还多。
解决这个问题没什么技巧，就是心里绷根弦：凡是用了对称面的模型，导出后第一件事确认求解域范围，然后决定后处理是否要做镜像扩展。镜像扩展本质上就是把坐标取反再复制一份场值，但要注意镜像面上的场分量方向也要跟着变，电场和磁场一个取反一个不取反，搞错了镜像出来的场分布会对不上。
这五个坑覆盖了我见过的绝大多数后处理事故。它们之间的共同点是：问题都不出在 MATLAB 代码逻辑上，而出在对 CST 导出数据本身的假设上。每拿到一个新模板、新项目，先花五分钟确认单位、列顺序、求解域和相位约定，比写完脚本再排错要快得多。
6. 把这套后处理流程固化成自己的工具箱：约定、验证与复用习惯
后处理代码写一次是脚本，写三次就是工具。当你手上有三四个项目都在做 CST 场分布分析时，与其每次复制粘贴再改参数，不如花半天时间把前面几章的代码整理成一套固定的工具箱。整理的关键不是把代码写得多漂亮，而是把输入、输出和验证的约定固定下来。
6.1 输入与输出约定：一份配置说清楚单位、模板与路径
我习惯在每个仿真项目文件夹里放一个配置文件，用简单的 MATLAB 脚本形式记录后处理要用到的所有元信息。字段不需要多，但每个都必须有。
% 项目专用后处理配置 cfg.unit = "mm"; % 建模与导出单位 cfg.nCols = 9; % 数据列数 cfg.coord = "global"; % 全局或局部坐标 cfg.phaseConv = "CST_default"; % 相位约定备注 cfg.exportDir = "CST_export"; % 导出数据目录读取函数和绘图脚本全部以 cfg 结构体为前置输入，这样换项目时只需要改配置文件，主程序一行都不用动。单位、列数、坐标和相位约定这四项是血的教训换来的，缺哪一个都可能让你排查半天还不一定想到。
6.2 验证后处理正确性的两个低成本手段
每次写完后处理代码，我会做两个快速验证，基本能覆盖绝大多数低级错误。第一个是峰值保持校验，把插值前后的场强最大值做对比，正常情况下偏差应该在 2% 以内，超过这个数往往意味着网格插值太疏或者原始数据里有离群点。第二个是对称性校验，模型有对称面时，对比对称位置两点的场强，偏差大说明坐标解析或镜像处理出了问题。
peakRaw = max(magEx); peakInterp = max(EXq(:)); errPct = abs(peakInterp / peakRaw - 1) * 100; fprintf('插值前后峰值偏差: %.2f%%\n', errPct);这两个验证不涉及复杂的电磁理论，只是纯粹的数据一致性检查，能抓住大多数因为单位、列序、坐标范围导致的系统性错误。
6.3 新项目里的一句“老规矩”
现在每次接手新的 CST 后处理任务，第一件事永远是先跑一遍配置确认和文件体检，再谈任何代码。单位是什么、列数是多少、有没有对称边界、相位约定是什么，这四件事在导出数据的那一刻就决定了，晚一步确认，后面所有工作都可能白做。把你自己的检查项固化成一段 5 分钟的例行脚本，你会发现所谓“玄学”般的后处理翻车，绝大多数都是这类元信息问题。希望这篇笔记能帮你把 CST 场分布和 MATLAB 后处理这条路走顺，少踩几次我当年踩过的坑。
本文还有配套的精品资源，点击获取
