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

本文详述了一种基于MATLAB的电力系统动态分析工具,涵盖了时域仿真、系统线性化、控制器设计等技术,以IEEE9和68节点系统为例,展示了如何进行系统调优和干扰处理。MATLAB代码提供了关键环节的实现细节。
摘要由CSDN通过智能技术生成

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

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

⛳️座右铭:行百里者,半于九十。

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

目录

💥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======================================
 

🎉3 参考文献

文章中一些内容引自网络,会注明出处或引用为参考文献,难免有未尽之处,如有不妥,请随时联系删除。

🌈4 Matlab代码、数据、文章

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值