MATLAB/Simulink電力系統與仿真，第七章輸電線路故障行波的仿真經驗

本人在學習simulink時參考此書，按照書中教程和參數搭建潮流計算模型，但是書中并未詳細給出是以的設定參數，對元件應用描述還沒有很細，仿真過程未較長的描述，是以在模組化過程中碰到了一些錯誤，可能是本人自身原因所造成，此書還是非常值得大家參考的。以下列舉一些我的模組化經驗：

simulink圖

使用的是Matlab2016b版本

首先發電機，用“Three-phase source”模型，三個發電機中E1相位設為0，E2相位設為30，E3相位設為60，其餘參數一樣，EM參數如圖：

測量子產品，用“Three-phase V-I measurement”，測量端口接示波器

輸電線路仿真子產品采用“Distributed Parameters Line”分布參數子產品，Line1、2、3、4的線路長度分别為100km,100km,150km,260km，其餘參數一樣，Line1參數如圖：

通過“To File”子產品轉變為M檔案格式，參數如圖：

設定仿真時間起始為0s,終止時間0.1s，采用步長為Ode23tb算法，設定故障為A相接地短路，故障時間為：[0.035 0.1]

示波器結果如圖：

U：

I：

輸電線路故障行波的提取：在matlab腳本中輸入如下代碼

%設定仿真模型在0.035s時發生故障
%故障分量取為從故障後的0.035s~0.039s減去故障前的0.015s~0.019s

clc
clear
load xingbo.mat;
m=n';
ua=m(3501:3900,2)-m(1501:1900,2);
ia=m(3501:3900,5)-m(1501:1900,5);
ub=m(3501:3900,3)-m(1501:1900,3);
ib=m(3501:3900,6)-m(1501:1900,6);
uc=m(3501:3900,4)-m(1501:1900,4);
ic=m(3501:3900,7)-m(1501:1900,7);
Q=1/3*[2 -1 -1 
    0 sqrt(3) -sqrt(3) 
    1 1 1];
um1=Q(1,:)*[ua ub uc]';
im1=Q(1,:)*[ia ib ic]';   %進行Clarke變換得到電壓電流的模量
Lm1=0.8984e-3;
Cm1=12.94e-9;
Zcm1=sqrt(Lm1/Cm1);       %求波阻抗
uf=(um1+im1*Zcm1);
ur=(um1-im1*Zcm1);        %求出正反向行波
uf1=uf';
ur1=ur';
t1=0:10:3990;
t=t1';
plot(t,uf1,'r',t,ur1,'b--');
xlabel('t/us');ylabel('u/V');
legend('正向行波','反向行波','location','northwest');  %Legend位置在左上角（西北方）
           

正向行波和反向行波結果如圖：

方向行波的知識書中有簡要描述，公式也有列出。

有任何問題歡迎留言或私信我。

這是我營運的公衆号，裡面書的PDF和此模型檔案以及各種幹貨後續更新：

聲明：本文隻為了學習使用，如有侵權，請聯系我撤回！