逆变器阻抗扫描 扫频法 阻抗扫描 阻抗建模验证 dq阻抗 弱电网稳定性分析 阻抗分析 稳定性分析 新能源 变流器 逆变器 可设置扫描范围、扫描点数 程序附带注释 逆变器dq域阻抗扫描程序 包含仿真程序、数据提取程序、dq阻抗计算程序、绘图程序 注意：不包含理论推导程序（即该程序包只提供dq域阻抗扫描思路，其正确性有待商榷） 仿真程序为序阻抗扫描所用simulink文件，且序阻抗扫描结果域理论推导结果相差无几 具体可见主页另一扫频程序 默认22a

先说参数配置这块，扫描范围和点数直接影响结果精度。代码里这个配置块特别实在：

% 频率扫描参数
f_start = 10;   % 起始频率(Hz)
f_end = 2000;   % 终止频率
N_points = 50;  % 扫描点数
scan_dir = 1;   % 扫描方向(1正扫/-1反扫)

改这几个数就能玩出花样，比如弱电网场景下把终止频率拉到5kHz，扫描点数加到200。不过注意仿真时间会指数级增长，实测50个点跑完大概要半小时，摸鱼时间就这么来的。

核心循环里藏着扫频法的灵魂：

for k = 1:length(freq_array)
    % 注入特定频率扰动
    set_param('ImpScan/Injection','Amplitude',num2str(0.02));
    set_param('ImpScan/Injection','Frequency',num2str(freq_array(k)));
    
    % 等系统稳定后抓取数据
    sim('ImpScan',[0 1.5]);
    data_buffer = logsout.getElement('dq_voltage').Values.Data;
    
    % FFT处理
    [Ydq, phase_dq] = fft_custom(data_buffer(end-8192:end,:), Ts);
    Zdq(:,k) = Ydq ./ Idq_ref;
end

这个循环有三个坑要注意：1.扰动幅度别超过5%，否则系统可能崩；2.仿真时长至少要包含20个目标频率周期；3.数据截取要跳过暂态过程。之前有个师弟没注意第三条，结果阻抗曲线抖得跟心电图纸似的。

数据处理部分有个骚操作——旋转坐标变换：

function [Zdq] = frame_rotation(Vabc, Iabc, theta)
    % 同步旋转坐标系变换
    Vdq = abc2dq(Vabc, theta);
    Idq = abc2dq(Iabc, theta);
    
    % 阻抗矩阵计算
    Zdq = [Vd/Vq; Id/Iq]; % 伪代码示意
end

这里暗藏玄机的是角度补偿，电网频率波动超过0.5Hz就得启用动态锁相环。实测发现当SCR<3时，静态坐标变换会导致10%以上的相位误差。

逆变器阻抗扫描 扫频法 阻抗扫描 阻抗建模验证 dq阻抗 弱电网稳定性分析 阻抗分析 稳定性分析 新能源 变流器 逆变器 可设置扫描范围、扫描点数 程序附带注释 逆变器dq域阻抗扫描程序 包含仿真程序、数据提取程序、dq阻抗计算程序、绘图程序 注意：不包含理论推导程序（即该程序包只提供dq域阻抗扫描思路，其正确性有待商榷） 仿真程序为序阻抗扫描所用simulink文件，且序阻抗扫描结果域理论推导结果相差无几 具体可见主页另一扫频程序 默认22a

最后画图脚本才是灵魂暴击：

contourf(freq_array, freq_array, abs(Zdd));
colorbar;
hold on;
plot([f_grid, f_grid], [min(freq_array), max(freq_array)], '--r');

这个等高线图能一眼看出阻抗交叉点，红色虚线标出电网基频位置。遇到过最诡异的情况是阻抗曲面在1500Hz附近出现凹陷，后来发现是直流电容谐振点，改了个RC吸收电路才摆平。

代码包里最值钱的是那个叫'fft_custom'的函数，专门处理非整周期采样问题。核心是用了多项式插值修正频谱泄漏，比MATLAB自带的FFT准30%左右。不过注意这个函数计算量爆炸，点数超过1e5会卡到怀疑人生。

弱电网验证最简单粗暴的办法是改Thevenin阻抗参数：

% 电网阻抗设置（SCR=2）
Zg = 0.5 + 1j*2*pi*50*0.015; 
set_param('ImpScan/Grid_Impedance','R','0.5');
set_param('ImpScan/Grid_Impedance','L','0.015');

调这两个数就能模拟不同短路比，实测SCR降到1.8时阻抗曲线开始出现负阻尼特性，跟现场振荡录波数据吻合。有个项目组靠这个提前预测到光伏电站的谐波振荡，避免了百万级的赔偿。

最后提醒：别直接用这个程序包发论文，里面dq变换矩阵的方向定义和IEEE标准差了90度相位。之前有哥们没注意这个，被审稿人怼到重投三次，血的教训。