Simulink PowerGUI FFT数据导出到Origin的完整指南-二趣网

1. 从Simulink到Origin：为什么FFT数据输出难倒了那么多人

做电力电子仿真的人基本都遇到过这个场景：Simulink里搭了个三相整流或逆变模型，仿真跑完，打开PowerGUI的FFT分析工具，看着那张漂亮的幅频特性曲线，心想“终于拿到结果了”。可当你想把这张图放进论文或者项目报告里时，问题来了——PowerGUI FFT分析窗口里既不能直接导出数据，也不能复制坐标轴数值。你只能要么截图，要么用肉眼一条条读数据，然后再到Origin里手动点出几十个点的曲线。

这事我折腾过很久。早期我的做法是：把Scope里的波形用To Workspace模块导出到工作区，然后自己写FFT代码重新分析。但这么做有一个绕不开的坑——你MATLAB手写FFT和PowerGUI FFT分析的结果往往对不上。差别可能来自窗函数、采样点数、基波频率设置、谐波次数范围等一堆细节，哪怕差一个点，你论文里的THD数值就和仿真对不上了，审稿人要是较真，很容易被质疑。

直到后来我找到了一个通用办法：直接从PowerGUI的FFT分析窗口内部下手，把它的底层数据完整抽出来，再按Origin的格式重组。这个方法不需要额外装任何工具箱，只需要一段MATLAB脚本，就能把PowerGUI FFT分析里显示的那条曲线、表格里的幅值/相位/频率数据，完完整整地搬到Origin里，画出来的图和Simulink里显示的一模一样。下面把整个思路和源码逐段讲清楚。

2. 核心思路拆解：PowerGUI FFT分析窗口到底藏了什么数据

2.1 PowerGUI FFT工具的工作原理

PowerGUI（Powergui）是Simulink Simscape Electrical（以前叫SimPowerSystems）里的一个交互式工具模块。你双击模型里的Powergui模块，就能看到FFT Analysis选项卡。它的工作流程是：从模型里选定一个信号（比如三相电流中的某一相），设定好起始时间和周期数，然后调用MATLAB的fft算法，把时域波形变成频域分量。

关键点在于：这个工具本质上就是 MATLAB 内部函数的一个图形化封装。也就是说，它的底层数据全在MATLAB工作区里，只不过图形界面把结果“画”给你看了而已。你没法从界面里复制，但可以通过脚本从工作区里把计算底层的数据元信息捞出来。

具体的实现机制其实是这样的：当你在FFT Analysis界面设置好参数并点击Display之后，PowerGUI会把选定信号的时间序列、采样时间、FFT分析点数等一系列信息，更新到当前工作区的powergui对象以及相关结构体里。MATLAB的内部函数power_fft（注意这个函数名）就是真正干活的函数。

2.2 为什么不能直接按“导出数据”按钮

在较新版本的Simulink里，PowerGUI界面里确实多了一个“Export”相关的选项，但实测下来，它导出到的目标要么是工作区的一个结构体，要么是生成一张figure图。这个导出的结构体里虽然包含了频谱数据，但格式并不友好——比如三相电流的时候，它会把三路信号叠在一个矩阵里，但没有清楚的列标签，你导入Origin以后还得花半天去识别哪列是A相、哪列是频率、哪列是幅值。

而且最麻烦的一点是：导出的频谱数据并不一定和你界面上看到的曲线完全一致。原因在于窗函数和归一化设置，这个我们后面细说。

所以我的思路很简单：绕开PowerGUI的那个“Export”按钮，直接调用power_fft函数，把FFT分析的原始输入信号和参数拿过来自己算一遍，然后把结果重组成一个清晰可用的表格。这样既保证数据和PowerGUI显示的一致性，又能按你想要的格式输出。

2.3 保证“一模一样”的三个前提条件

要实现导出的数据和PowerGUI FFT界面显示完全一致，有三个条件必须同时满足：

同一个输入信号：你脚本里做FFT的信号，必须和你在PowerGUI里选的那个信号完全相同。这个信号通常来自模型的To Workspace导出。
相同的分析参数：起始时间（Start time）、周期数（Number of cycles）、最大频率（Max frequency）、基波频率（Base frequency）这些参数必须和你在PowerGUI界面设置的完全一样。
相同的FFT实现方式：理论上MATLAB的fft函数是统一的，但PowerGUI的power_fft在内部做了一些特殊处理——比如abc三相变换到dq坐标系下的FFT、信号周期延拓、直流分量处理等。如果直接用裸fft，结果会有微小的差异。所以要“一模一样”，最好的办法是调用power_fft函数本身，或者复刻它的计算流程。

如果你只是想导出和自己界面上看到的图完全一样的结果，用power_fft的返回值是最省心的方案。

3. 实操前置准备：Simulink模型与工作区数据导出

3.1 在Simulink里设置好信号导出

在动手写脚本之前，先确保你Simulink模型里已经加好了To Workspace模块。这一步很多人会忽视，导致后面脚本没数据可用。

具体操作建议：

在你需要做FFT分析的信号线（比如三相电流i_a）上引出分支，接到一个To Workspace模块。
To Workspace模块的Variable Name设置成i_a，输出数据类型建议选Timeseries而不是Array。原因在于Timeseries保留了时间向量，这样你后续做周期对齐、起始时间定位就方便得多；Array只有数值矩阵，适合你已经明确知道采样时间的情况。
仿真Stop Time要留足裕量，至少要覆盖你FFT分析需要的时长。比如你要分析2个工频周期（0.04秒），仿真时长建议至少0.1秒，留出前面的启动暂态。

提示：如果你更习惯用Scope模块，也可以直接在Scope里勾选“Logging”属性，将数据记录到工作区。但个人经验是To Workspace更可控，特别是信号名字可以自定义，不容易搞混。

3.2 获取FFT分析所需要的模型参数

在脚本里，你要动态获取以下几个参数：

采样时间Ts：通常就是仿真步长（定步长仿真）或PowerGUI里显示的采样时间。如果你用的是变步长仿真，FFT前需要先把信号重采样到均匀时间间隔，否则power_fft会报错或结果不准。
基波频率f_base：这个一般就是50Hz或60Hz，你模型里用的什么频率就填什么。
起始时间t_start：PowerGUI里默认是0秒。如果你要看稳态后的波形，可以设成比如0.04秒。
周期数n_cycles：通常是1、2、3等整数。假定你设了2个周期，那么FFT分析的数据长度就是2 / f_base所对应的时间跨度。

这些参数在脚本里可以用变量的方式灵活配置，模拟你在PowerGUI界面里的操作。

3.3 关于稳态判定的一个小建议

PowerGUI FFT分析的准确性高度依赖你选的“分析窗口”落在稳态区域。如果在整流器启动过程中做FFT，谐波含量会被启动暂态严重污染。建议你在仿真模型里加一个计时器或者直接用时间条件判断，等系统进入稳态后再导出To Workspace的数据。

比如三相不控整流桥带阻感负载，启动暂态一般持续2~3个工频周期。你可以把t_start设为0.06秒（约3个周期），这时候基波幅值已经基本稳定了。用MATLAB跑FFT前先画一下时域波形肉眼确认，别盲目自信。

4. MATLAB源码解析：从工作区提取FFT数据并格式化输出

4.1 主脚本结构

先给出一段可以直接运行的脚本，这段脚本是我目前最常用、经过多个版本Simulink验证的版本。拿过去改一下变量名和路径就能用。

% ===================================================== % PowerGUI FFT 数据导出脚本 (通用版) % 适用条件：Simulink模型中使用PowerGUI进行FFT分析 % 依赖：电力电子仿真模型已通过To Workspace导出信号 % ===================================================== clear; clc; close all; %% 1. 用户配置区（改成你自己的参数） % 信号名称（对应Simulink中To Workspace导出的变量） signal_name = 'i_a'; % 要分析的单相信号 % 如果分析三相，可以分别设 'i_a','i_b','i_c' % FFT参数（必须与PowerGUI界面设置一致！） f_base = 50; % 基波频率(Hz) t_start = 0.06; % 起始时间(s) n_cycles = 2; % 分析周期数 max_freq = 1000; % 最大分析频率(Hz) THD_base = f_base; % THD计算基准频率 %% 2. 加载数据 % 从工作区读取Simulink导出的Timeseries数据 if exist(signal_name, 'var') ~= 1 error('找不到信号变量 %s，请先在Simulink中运行仿真并导出数据。', signal_name); end data_ts = eval(signal_name); t = data_ts.Time; x = data_ts.Data; x = x(:); % 强制列向量 % 剔除NaN或Inf（变步长仿真偶尔会出现） valid_idx = isfinite(t) & isfinite(x); t = t(valid_idx); x = x(valid_idx); %% 3. 重采样到均匀时间网格 % 如果模型是变步长，这里做线性插值；如果已经是定步长，跳过也行 Ts = 1 / (max_freq * 20); % 采样周期：至少10倍于最大分析频率 % 这里用20倍是为了让FFT结果更平滑，实际也可以用10倍 t_uniform = (t_start : Ts : t_start + n_cycles/f_base - Ts)'; x_uniform = interp1(t, x, t_uniform, 'linear'); %% 4. 调用power_fft函数（核心） % 注意：power_fft是Simscape Electrical的内部函数，名字在不同版本可能略有差异 % 建议先用 which power_fft 确认一下你的MATLAB里是否存在 try % 计算FFT % 输出说明: % f : 频率向量 (Hz) % mags : 幅值向量 (与PowerGUI显示的幅值定义一致) % ph : 相位向量 (度) [f, mags, ph] = power_fft(x_uniform, Ts, f_base, ... t_start, n_cycles/max_freq, ... max_freq, [], 'freq'); % 说明：不同版本power_fft输出顺序可能有区别， % 建议先调用帮助文档 doc power_fft 仔细核对 catch ME warning('power_fft调用失败，改用MATLAB自带fft函数计算(结果可能与PowerGUI略有出入)'); disp(ME.message); % 后备方案：直接用fft（需要自己做窗函数和幅值修正） N = length(x_uniform); Y = fft(x_uniform); P2 = abs(Y / N); P1 = P2(1:floor(N/2)+1); P1(2:end-1) = 2 * P1(2:end-1); f = (0:floor(N/2)) * (1 / (N * Ts)); mags = P1; ph = angle(Y(1:floor(N/2)+1)) * 180 / pi; end %% 5. 控制需要输出到Origin的频点范围 % 只保留从基波到max_freq的范围，同时把直流分量(0Hz)剔除 idx_keep = (f >= f_base - 0.5) & (f <= max_freq + 0.5) & (f > 0); % 如果希望保留0Hz直流分量，可以把条件改成 f <= max_freq f_out = f(idx_keep); mags_out = mags(idx_keep); ph_out = ph(idx_keep); % 计算THD（总谐波畸变率） % 注意：这个THD和PowerGUI界面显示的THD可能略有不同， % 因为PowerGUI通常只算到某一次谐波，需要你手动核对。 fund_idx = find(abs(f_out - f_base) < 1e-3, 1); if ~isempty(fund_idx) fund_mag = mags_out(fund_idx); harm_mags = mags_out; harm_mags(fund_idx) = []; % 去掉基波 THD = sqrt(sum(harm_mags.^2)) / fund_mag * 100; else THD = NaN; end fprintf('基波幅值: %.4f, THD: %.2f%%\n', fund_mag, THD); %% 6. 导出为CSV（Origin可以直接导入） output_table = table(f_out, mags_out, ph_out, ... 'VariableNames', {'Frequency_Hz', 'Magnitude', 'Phase_Deg'}); writetable(output_table, 'fft_export_data.csv'); disp('数据已导出到 fft_export_data.csv'); %% 7. (可选)画一个简易频谱图验证 figure; stem(f_out, mags_out, 'MarkerSize', 4); xlabel('频率 (Hz)'); ylabel('幅值'); title('FFT频谱 (建议与PowerGUI对比)'); grid on;

4.2 关键代码段的逻辑解释

第一处关键：重采样频率选择

Ts = 1 / (max_freq * 20)这行我故意用了20倍而不是理论上的2倍（奈奎斯特采样定理只需要2倍以上即可）。原因很简单：PowerGUI内部对FFT的处理用了线性插值和平滑处理，如果你用刚好2倍采样率，频谱图上会出现轻微的“栅栏效应”导致幅值偏差，和PowerGUI显示的结果对不上。用20倍采样密度，FFT结果会和PowerGUI的结果在误差范围内几乎完全一致。

要注意的是，重采样拉高了采样率之后，FFT点数会增加很多，计算会稍微变慢，但对今天的电脑来说完全不是问题。

第二处关键：power_fft函数的调用方式

power_fft这个函数是Simscape Electrical模块集的内部函数，官方文档里没有完整的说明，但它在各个MATLAB版本里基本都存在。它的输入参数比裸fft多，用处就在于它自己处理了周期延拓、幅值定义、窗函数选择这些细节，保证输出和PowerGUI分析窗口里看到的一致。

这里特别提醒一句：不同MATLAB版本对power_fft的输入参数顺序可能有一点调整。在我测试过的R2020b和R2023a版本里，上述传参方式都能正常工作，但如果你用的是更老或更新的版本，建议先执行：

doc power_fft

或者直接在命令窗口输入：

which power_fft

来确认函数存在以及查看它的帮助信息。

第三处关键：CSV导出格式

表格里输出了三列：频率、幅值、相位。Origin导入CSV时，默认就会创建三个列，直接可以用。如果你后面想在Origin里做对数坐标、双纵轴或者画柱状频谱图，这三列数据已经完全够用了。

4.3 如果你用的是三相系统

三相系统里最常遇到的场景是：你需要在PowerGUI里分别看A相、B相、C相的频谱。脚本的扩展非常直接：把i_a替换成i_b、i_c，分别导出三个CSV，然后在Origin里合并到一个图里。还有一个更简便的做法——在脚本外层写一个循环：

signal_list = {'i_a', 'i_b', 'i_c'}; for k = 1:length(signal_list) signal_name = signal_list{k}; % 后面的代码段一样 end

你还可以顺便把三相FFT的基波幅值和THD打印成一个摘要表，方便写结论。

5. Origin端操作要点：让图表和Simulink里完全一致

5.1 如何从CSV导入数据并快速绘图

Origin导入CSV非常简单：

打开Origin，选择File > Import > CSV...，选中刚才生成的fft_export_data.csv。
在弹出的导入向导里直接点“Finish”就行。默认设置下，第一列成为X轴（Frequency_Hz），第二列和第三列成为Y列（Magnitude与Phase_Deg）。
做频谱图时，通常只需要Magnitude那一列。选中前两列，点击Plot > Line或者Plot > Symbol > Stem，就能得到频谱曲线。

注意：如果你在Simulink里看到的是“幅值（标幺值或实际值）”的离散谱线，那你Origin里也选Stem图或者Stem+Line叠加图，视觉上才能对得上。普通的连续Line图会让人误以为谱线之间连续变化，实际上FFT分析应该是离散的谱线。

5.2 调整坐标轴范围与标签

让Origin里的图和PowerGUI里的图“一模一样”，关键在坐标轴范围。

X轴：在Origin里设置X轴范围从0到你的max_freq，和PowerGUI里设置的最大频率一致。
Y轴：PowerGUI里默认显示幅值，单位取决于你的信号单位（电流就是A，电压就是V）。在Origin里也要把Y轴单位标成一样的，不然会看起来“不一样”。
网格和字体：这个看个人习惯，但有个小细节——PowerGUI的频谱图是线宽较粗的离散谱线（类似stem），如果你在Origin里用默认的细线，打印出来会觉得“不太像”。建议在Origin里把线宽调到1.5或2，符号用方块或圆圈。

5.3 如何把导出的数据做成“论文级”的图

发论文的话，通常不会直接用PowerGUI默认的样子，而是要重新排版。我的习惯是：

用Origin画一个双Y轴图：左边Y轴放电流幅值（或电压幅值），右边Y轴放THD值（如果用柱状图表示各次谐波占比的话）。X轴放谐波次数（1次、2次、3次……直到n次），而不是频率Hz。这样看图的人一眼就能看出各次谐波的分布规律。

要把频率（Hz）换算成谐波次数很简单：用Origin里的Set Column Values功能，添加一列公式Frequency_Hz / 50就得到谐波次数。

5.4 批量处理多组数据的小技巧

如果你要比较不同工况下的频谱（比如轻载、额定负载、过载），每一组都跑一遍Simulink仿真并生成一个CSV，然后把所有CSV都在Origin里画成一个图。这时候有一个小技巧：

用Origin的Import Multiple Files功能，一次选中多个CSV文件，然后逐列添加，再设置成多Y轴图。这样每一条谱线都能区分开，论文里就能按“不同负载条件下的频谱特性”为标题出一张对比图。

6. 常见问题与排查技巧实录

6.1`power_fft`函数不存在或报错

这是遇到最多的问题。原因一般是你的MATLAB安装时没有完整安装Simscape Electrical工具箱，或者版本太老。排查步骤：

命令窗口输入ver，看看是否存在Simscape Electrical这项。
如果存在，试试which power_fft -all，有时这个函数名可能变成了power_fft_fft或者带版本后缀的名字。你可以在MATLAB安装目录下搜索power_fft开头的文件。

如果确实找不到，就用脚本里的后备方案（直接调fft）。但你要做好心理准备：结果和PowerGUI可能有差异，特别是幅值精度上。

6.2 FFT结果和PowerGUI界面上显示的幅值不一致

我排查过的多数情况都是因为重采样率不够高或者窗函数没对齐。PowerGUI内部在计算FFT之前会做一次“周期截断”，也就是把时域信号截取恰好整数个周期。如果你重采样后的时间跨度不是严格等于n_cycles/f_base的整数倍，FFT出来的谱线就会发生频率泄漏，导致基波幅值偏小。

解决方法是严格保证t_uniform的长度N满足：

N = round(n_cycles / (f_base * Ts)); % 确保是整数 t_uniform = t_start + (0:N-1)' * Ts;

这个细节非常重要。我在脚本里写的是t_uniform = (t_start : Ts : t_start + n_cycles/f_base - Ts)'，在浮点数运算下末点可能有一点误差，建议改成上面这种通过整数N来构造的方式。

另外，还有个容易忽略的点：PowerGUI显示的是“单边幅值谱”，也就是正频率分量乘以2后的结果。power_fft函数默认已经做了单边处理，但后备方案的fft代码里也手动做了P1(2:end-1) = 2 * P1(2:end-1)，这个别漏了。

6.3 Origin导入CSV后频率列出现科学计数法

这个不是错误，但看起来不舒服。在Origin里选中频率列，右键Properties，把Display Format改成Decimal（小数位数设成2位或自定义）。这样频率坐标轴上就不会显示4.00E+01这种形式了。

6.4 模型用了Carsim或其他外部联合仿真平台

做Carsim + Simulink联合仿真的朋友经常会遇到一个问题：Simulink里的信号带了外部平台的时间戳或者包含了外部平台的噪声，导致FFT结果异常。

解决办法是在做FFT之前，先把信号用低通滤波或移动平均处理掉高频数值噪声。但注意，滤波本身也会改变信号频谱，所以这个操作要谨慎，且滤波后要和PowerGUI的结果交叉验证一下。

6.5 THD数值和PowerGUI差零点几个百分点

这个我花了很长时间才定位到原因。PowerGUI的THD计算是基于基波分量和2次及以上的各次谐波，但它默认只考虑到指定最大频率对应的最高次谐波。比如你最大频率设了1000Hz，基波50Hz，那么PowerGUI只算到20次谐波。而如果你后台脚本里把所有谐波分量（哪怕超过1000Hz的）都算进去了，THD就会偏大。

修正方法：在写THD计算代码时，过滤掉频率超过max_freq的谱线再做平方和。上面给的脚本里，idx_keep已经限制了f <= max_freq + 0.5，所以只要保持这个条件不变，理论上THD应该和PowerGUI一致。如果还差一丢丢，多半是相位对消或者数值截断导致的，这时候建议以PowerGUI为准，脚本主要用于批量数据提取。

7. 避免踩坑：关于幅值定义和标准化

7.1 实际值、标幺值与dB三种显示模式

PowerGUI的FFT窗口提供了三种幅值表示方式：

实际值（Magnitude）：直接是信号的幅值，比如220V电压的基波幅值就是311V（峰值）。
标幺值（Per Unit）：以你设定的“基波幅值”为1个单位，其他各次分量相对基波的比值。
分贝值（dB）：20*log10(实际值/基准值)。

这些需要和你论文里用的坐标系匹配。如果是描述电流谐波含量，建议用实际值（A）；如果要对比不同工况下的各次谐波占比，用标幺值更直观。

脚本里如果你需要输出标幺值，只需要在导出前做一次归一化：

mags_pu = mags_out / fund_mag;

7.2 相位数据的处理

PowerGUI里除了幅值还会显示相位，但相位角度通常不太稳定（特别是接近0幅值的谱线，相位会随机跳变）。论文里一般不放相位，除非你做的是阻抗分析或者闭环控制设计需要相位信息。

导出到Origin时，相位列保留即可，但在画图时不要随意连接相位散点，容易形成误导。如果要用相位信息，建议用极坐标图或者Bode图样式。

8. 自动批量仿真：把脚本嵌入到工作流里

上面讲的是单次导出的流程。实际项目中，我更推荐把这套脚本嵌入到整个仿真自动化流程里，实现“跑完仿真，自动导出FFT数据，自动生成Origin图”的流水线。

具体可以这样做：

在Simulink模型里用set_param循环修改工况参数（比如负载电阻值、转速给定值）。
每次修改后，用sim命令让模型运行一次。
仿真结束后，调用上面那段脚本，导出当前工况的FFT数据到独立的CSV文件。
所有工况跑完后，写一个小脚本或者用Origin的Batch Import功能，一次性生成所有频谱图。

举个例子，如果你在研究三相整流器在不同负载下的谐波情况，你可以设置负载从10%到100%每次增加10%，一共10种工况。脚本自动跑完仿真、自动导出10个CSV，然后再统一画图，效率远超手动操作。

这里给一个简化的自动化循环框架：

loads = 0.1:0.1:1.0; for k = 1:length(loads) % 修改Simulink模型里的负载参数 load_val = loads(k); % 假设模型里有一个变量R_load表示负载电阻 R_load = 100 / load_val; % 根据具体模型调整 assignin('base', 'R_load', R_load); % 运行仿真 simOut = sim('three_phase_rectifier.slx'); % 导出数据（注意每次先清空上一轮的变量） clear i_a i_b i_c; i_a = simOut.get('i_a'); i_b = simOut.get('i_b'); i_c = simOut.get('i_c'); % 调用FFT导出函数 run_fft_export(signal_name='i_a', f_base=50, ... t_start=0.06, n_cycles=2, max_freq=1000, ... output_filename=sprintf('fft_load_%02d.csv', k)); end

想省事的话，也可以把FFT导出部分封装成一个function文件run_fft_export.m，每次调用只需要传参数，返回一个表格，统一存CSV。

9. 我踩过的坑：关于版本兼容和数据精度

9.1 Simulink版本升级后PowerGUI模块名变化

如果你用过几年前的MATLAB，应该记得，以前这些电力电子模块叫SimPowerSystems，后来改成了Simscape Electrical。PowerGUI模块在新的工具箱下可能不再显示为“Powergui”，而是“powergui”或其他名字，但核心功能没变。如果你的模型打开后PowerGUI不见了，可以去库浏览器里搜powergui重新拖一个进来，不影响已有模型运行。

9.2 采样时间一定要低于奈奎斯特频率

很多初学者在定步长仿真时把采样步长设得不够小，直接导致FFT的高频部分完全乱掉。比如你仿真逆变器输出，开关频率是10kHz，那你至少要保证采样步长低于50微秒（即奈奎斯特频率大于10kHz）才能看到10kHz处的谐波。如果采样步长是100微秒，那10kHz的谐波会直接混叠成你看不懂的频率。

这个比FFT要算多少个点更重要。我在做并网逆变器的谐波分析时，一般把采样时间设置为1微秒（1e-6），这样可以保证看到50kHz以内的频谱，足够覆盖常见的开关谐波。

9.3 Origin绘图时不要直接连接离散点

这个问题我在帮几个同学看论文图的时候经常发现：他们在Origin里把频谱数据画成了一条连续的曲线，看起来和PowerGUI的离散谱线完全不像。正确做法是：

选中数据列，右键 ->Plot As改成Stem（或者Column）。
如果想让图更好看，用Stem + Line混合样式，数据点位置显示竖线和符号，输出效果最接近PowerGUI。

9.4 如何应对数据量巨大、Origin卡顿

有些仿真一个工况跑完，FFT数据导出后可能有成千上万根谱线。导入Origin后操作会明显卡顿。建议在脚本里做一次“谱线稀疏化”处理，只要保留幅值大于基波幅值某百分比的谱线（比如0.1%），其他忽略。这样可以把几千条谱线压缩到几十条，画图更快，也不影响解读。

threshold = max(mags_out) * 0.001; % 0.1%的基波幅值作为阈值 idx_save = mags_out >= threshold; f_out = f_out(idx_save); mags_out = mags_out(idx_save); ph_out = ph_out(idx_save);

这样操作后导出的CSV可能只有几十行，但关键谐波信息全都在。

10. 把脚本整合成一个小工具

用了几次之后，我把这个脚本优化成了一个可以直接在命令窗口调用的函数。在命令行输入exportFFTtoOrigin('i_a', 50, 0.06, 2, 2000)这种格式，就能直接导出数据，非常方便。

一个简洁版函数定义如下：

function exportFFTtoOrigin(signal_name, f_base, t_start, n_cycles, max_freq) % 导出PowerGUI FFT数据到Origin兼容的CSV文件 % 示例: exportFFTtoOrigin('i_a', 50, 0.06, 2, 2000) % 省略了上面详述的数据处理流程，核心逻辑一致 % 输出文件名自动包含信号名和时间戳，防止覆盖 output_file = sprintf('FFT_%s_%s.csv', signal_name, datestr(now, 'yyyymmdd_HHMMSS')); writetable(output_table, output_file); fprintf('数据已导出到 %s\n', output_file); end

如果你想在Origin里一键调用MATLAB脚本，还可以在Origin的Command Window里使用matlab -batch或者用Origin的LabTalk的system命令调用MATLAB。但说实话，最方便的方式还是先在MATLAB里把CSV都生成好，再到Origin里导入，两个软件同步操作反而容易乱。

11. 扩展思路：从“导数据”到“自动出图”

如果你已经能做到FFT数据自动导出，下一步可以尝试用MATLAB直接出图，而不经过Origin。虽然最终论文图一般还是用Origin做，但在前期数据验证阶段，用MATLAB快速画一张频谱图可以省掉很多来回导入导出的时间。

方法很简单：power_fft的结果拿到手之后，用stem(f, mags)画谱线，设置好坐标轴标注和网格，直接就能和PowerGUI对照。这个验证步骤我一般每次跑完仿真都会做一遍，确保脚本和PowerGUI设置没有偏差。

有了脚本基础，你还可以进一步扩展：

把三相的FFT结果放在同一个图里对比。
输出一个“谐波次数-幅值”表，方便直接粘贴到Word文档里。
自动计算指定次谐波的占比，比如3次、5次、7次谐波。

这些扩展都很简单，核心还是那三列数据：频率、幅值、相位。

12. 为什么“一模一样”导出对工程报告这么重要

最后聊聊这个方法的实际价值。

在工程实践中，仿真与论文报告的数据保持一致，是基本要求。如果PowerGUI里看到的THD是5.21%，而你论文里写的THD是5.18%，虽然数值接近，但严格来说就已经是“对不上”了。对严谨的工程报告来说，这种数据不一致很致命，可能会导致返工或者被质疑。

有了这套脚本，每次仿真结束后，导出的数据和PowerGUI显示的数据是同一套计算逻辑产生的，自然就是完全一致的。你再也不需要手动抄录或补画几条线来“凑数”，也不用担心截图里字太小或像素不够。

另一个价值是效率。手动从PowerGUI导出一次数据然后作图，熟练的人至少也要5~10分钟。用脚本自动化之后，从Simulink仿真结束到Origin图出来，整个过程可以控制在几秒内。遇到需要批量对比十几个工况的时候，省下的时间非常可观。

我自己用这个方法处理过三相整流器、逆变器、有源滤波器等多类模型的谐波分析，目前在MATLAB R2020b和R2023a上都验证过，没有发现兼容性问题。如果你用的版本确实比较特殊，建议先在一台机器上测试一遍，确认power_fft函数可调用后再推广到整个项目流程。

企业官网建设流程全解析