基于 MATLAB 的电力系统动态分析研究【IEEE9、IEEE68系节点】

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

2.1 IEEE9节点

2.2 IEEE68节点

🎉3 参考文献

🌈4 Matlab代码、数据、文章

💥1 概述

本文介绍了为电力系统动态分析开发的基于 MATLAB 的程序。可以获得时域仿真、系统线性化、模态分析、参与因子分析和可视化、控制器的优化放置、反馈信号选择、频率响应分析和控制设计。除了解决电力系统问题外，该软件包还提供模型在时域和状态空间中的符号和矢量化表示。该软件包充分利用了 MATLAB 强大的求解器的优势，用于求解非刚性和刚性问题。显式和隐式技术都用于求解微分代数方程（DAE）。假设同步电机配备了励磁器、涡轮和稳定器。负载可以建模为电压相关负载和独立负载。本文使用的测试系统是IEEE 9节点和68节点系统，以及德克萨斯州的2007节点合成电源系统。不同类型的干扰应用于系统，包括发电机侧和网络侧干扰。

📚2 运行结果

2.1 IEEE9节点

2.2 IEEE68节点

部分代码：

%%%----------------------------% LL-1--------------------------------------
T1_LL1 = realp('T1_LL1',1); % T1 coefficient (name and initial value)
T2_LL1 = realp('T2_LL1',0.1); % T2 coefficient
T1_LL1.Minimum = 0.1; T1_LL1.Maximum = 1; % Set min-max values for T1
T2_LL1.Minimum = 0.01; T2_LL1.Maximum = 0.1; % Set min-max values for T2
T1_LL1.Free = true;
T2_LL1.Free = true;
LL1=tf([T1_LL1 1],[T2_LL1 1]);
%%%----------------------------% LL-2--------------------------------------
T3_LL2 = realp('T3_LL2',1); % T3 coefficient
T4_LL2 = realp('T4_LL2',0.1); % T4 coefficient
T3_LL2.Minimum = 0.1; T3_LL2.Maximum = 1; % Set min-max values for T3
T4_LL2.Minimum = 0.01; T4_LL2.Maximum = 0.1; % Set min-max values for T4
T3_LL2.Free = true;
T4_LL2.Free = true;
LL2=tf([T3_LL2 1],[T4_LL2 1]);
%%%--------------------------Washout---------------------------------------
Tw = realp('Tw',1); % Tw coefficient
Tw.Minimum = 1; Tw.Maximum = 10;
Tw.Free = true;
WO=tf([Tw 0],[Tw 1]);
%%%--------------------------Gain------------------------------------------
Ck = realp('Ck',1);
Ck.Minimum = 1; Ck.Maximum = 50;
Kg=tf(Ck);
%%%==========================System tuning=================================
CL0 = feedback(LL1*LL2*Ck*G,1, -1); % Closed-loop TF (with PSS)
CL0.InputName = 'ws';
CL0.OutputName = 'w';
Req1 = TuningGoal.Poles(0,0.2,Inf); % [min decay, min damping ratio, max freq].
options = systuneOptions('Display','iter');
[CL,fSoft] = systune(CL0,Req1, options);
sys2=CL.Blocks;
OptimizedParam=[sys2.T1_LL1 sys2.T2_LL1 sys2.T3_LL2 sys2.T4_LL2 sys2.Ck]; OptimizedParam=double(OptimizedParam);
figure (1)
viewGoal(Req1,CL); xlim([-50 1]); hold on; % Plot the results with the desired goal to check if it is satisfactory
%=================================END======================================