第一篇：MATLAB計(jì)算24點(diǎn)游戲代碼
clear,closeall clc a=5;b=7;c=10;
d=4;%這里輸入需要計(jì)算的四個(gè)數(shù)字a,b,c,d f=[a b c d];tic;g=perms(f);[m,n]=size(g);h='+-*/';fori=1:24 for k1=1:4 for k2=1:4 for k3=1:4
str11=[num2str(g(i,1)),h(k1),num2str(g(i,2)),h(k2),num2str(g(i,3)),h(k3),num2str(g(i,4))];
str22=['(',num2str(g(i,1)),h(k1),num2str(g(i,2)),')',h(k2),num2str(g(i,3)),h(k3),num2str(g(i,4))];
str33=['(',num2str(g(i,1)),h(k1),num2str(g(i,2)),h(k2),num2str(g(i,3)),')',h(k3),num2str(g(i,4))];
str44=['(',num2str(g(i,1)),h(k1),num2str(g(i,2)),')',h(k2),'(',num2str(g(i,3)),h(k3),num2str(g(i,4)),')',];A=str2num(str11);B=str2num(str22);C=str2num(str33);D=str2num(str44);if A==24||B==24||C==24||D==24 break else end end
if A==24||B==24||C==24||D==24 break else end end
if A==24||B==24||C==24||D==24 break else end end
if A==24||B==24||C==24||D==24 break else end end
if A==24
answer=str11;elseif B==24
answer=str22;elseif C==24
answer=str33;elseif D==24
answer=str44;else
answer='無(wú)解';end
disp(['計(jì)算方法',num2str(answer)])time=toc;
disp(['計(jì)算耗時(shí)',num2str(time),'s'])
第二篇：24點(diǎn)游戲代碼
// +-* /---/ // 0 1 2 3 4 5
#include
int treat(float a,float b,float c,float d);float myF(int flag,float m,float n);void myPrint(int type,int i,int j,int k,float a,float b,float c,float d);
int time,temp=0;void main(){ int i,j,k,t,again,res,flag;float num[4];again=1;while(again==1){
printf(“nPlease Enter 4 nums(1~13):n”);
i=0;
flag=0;
while(flag==0)
{
i++;
// printf(“Input num-%dn”,i);
for(i=0;i<4;i++)
{
scanf(“%f”,&num[i]);
if(num[i]<1 || num[i]>13 || num[i]!=int(num[i]))
flag++;
}
if(flag!=0)
{
printf(“Error input againn”,i);
flag=0;
}
else
flag=1;
}
for(i=0;i<4;i++)
for(j=0;j<4;j++)
if(j!=i)
for(k=0;k<4;k++)
if(k!=j && k!=i)
for(t=0;t<4;t++)
if(t!=i && t!=j && t!=k)
{
res=treat(num[i],num[j],num[k],num[t]);
} if(res==0)
printf(“nNo answern”);else;// printf(“time=%dnn”,time);printf(“n1: Go onn2: Quitn”);scanf(“%d”,&again);} }
int treat(float a,float b,float c,float d){ int i,j,k;float sum1,sum2,sum3;for(i=0;i<4;i++)
for(j=0;j<6;j++)
for(k=0;k<6;k++)
{
if((!(i==3 && b==0))&&(!(j==3 && c==0))&&(!(k==3 && d==0)))
{
sum1=myF(i,a,b);
sum2=myF(j,sum1,c);
sum3=myF(k,sum2,d);
if(fabs(sum3-24)<0.1)
{
temp++;
myPrint(1,i,j,k,a,b,c,d);
// printf(“sum1:myF(%d,%2.0f,%2.0f)sum1=%fn”,i,a,b,sum1);
// printf(“sum2:myF(%d,%2.0f,%2.0f)sum2=%fn”,j,c,d,sum2);
// printf(“1:myF(%d,myF(%d,myF(%d,%2.0f,%2.0f),%2.0f),%2.0f)
sum3=%fnn”,k,j,i,a,b,c,d,sum3);
}
}
if(k==2)
{
sum1=myF(i,a,b);
sum2=myF(j,c,d);
sum3=sum1*sum2;
if(fabs(sum3-24)<0.1)
{
temp++;
myPrint(2,i,j,k,a,b,c,d);
// printf(“sum1:myF(%d,%2.0f,%2.0f)sum1=%fn”,i,a,b,sum1);
// printf(“sum2:myF(%d,%2.0f,%2.0f)sum2=%fn”,j,c,d,sum2);
// printf(“2:myF(%d,myF(%d,%2.0f,%2.0f),myF(%d,%2.0f,%2.0f))
sum3=%fnn”,k,i,a,b,j,c,d,sum3);
}
}
if(k==3)
{
sum1=myF(i,a,b);
sum2=myF(j,c,d);
if(sum2!=0)
{
sum3=sum1/sum2;
if(fabs(sum3-24)<0.1)
{
temp++;
myPrint(3,i,j,k,a,b,c,d);
// printf sum1=%fn“,i,a,b,sum1);
// printf sum2=%fn”,j,c,d,sum2);
// printf(“3:myF(%d,myF(%d,%2.0f,%2.0f),myF(%d,%2.0f,%2.0f))
}
}
}
} if(temp==0)
return 0;else
return 1;}
float myF(int flag,float m,float n){
// time++;if(flag==0)
return(m+n);if(flag==1)
return(m-n);if(flag==2)
(”sum1:myF(%d,%2.0f,%2.0f)(“sum2:myF(%d,%2.0f,%2.0f)sum3=%fnn”,k,i,a,b,j,c,d,sum3);
return(m*n);if(flag==3)
if(n==0)
return 30000;
else
return(m/n);if(flag==4)
return(n-m);if(flag==5)
if(m==0)
return 30000;
else
return(n/m);return 0;}
void myPrint(int type,int i,int j,int k,float a,float b,float c,float d){ char sigle[6];
sigle[0]='+';
sigle[1]='-';
sigle[2]='*';
sigle[3]='/';
sigle[4]='-';
sigle[5]='/';if(type==1){
if(j==4 || j==5)
{
if(k==4 || k==5)
printf(“%2.0f %c(%2.0f %c =24n”,d,sigle[k],c,sigle[j],a,sigle[i],b);
else
printf(“(%2.0f %c(%2.0f %c =24n”,c,sigle[j],a,sigle[i],b,sigle[k],d);
}
else if(k==4 || k==5)
{
printf(“%2.0f %c((%2.0f %c =24n”,d,sigle[k],a,sigle[i],b,sigle[j],c);
}
else
printf(“((%2.0f %c %2.0f)%c =24n”,a,sigle[i],b,sigle[j],c,sigle[k],d);}
(%2.0f %2.0f))%2.0f)%2.0f)%c %2.0f))
%c %2.0f
%c %2.0f)
%c %2.0f if(type==2 || type==3){ // if(k==4 || k==5)// printf(“(%2.0f %c(%2.0f %c %2.0f)=24n”,c,sigle[j],d,sigle[k],a,sigle[i],b);// else
printf(“(%2.0f %c %2.0f)%c =24n”,a,sigle[i],b,sigle[k],c,sigle[j],d);} }
%2.0f)%c
(%2.0f %c %2.0f)
第三篇：Q學(xué)習(xí)走迷宮MATLAB代碼
gamma = 0.4;%%設(shè)置初值 Q = zeros(6, 6);R = [-1-1-1-1 0-1;-1-1-1 0-1 100;-1-1-1 0-1-1;-1 0 0-1 0-1;0-1-1 0-1 100;-1 0-1-1 0 100];now = 2;seq = [];for now = 1:1:6
for i = 0:1:20 %%打亂路徑順序，選取隨機(jī)的路徑 rcolumn = R(now, :);
rcolumn(rcolumn ==-1)= [];num = randperm(size(rcolumn,2));
next = find(R(now,:)== rcolumn(num(1)), 1);seq = [seq next];
Qmax = max(Q(next,:));%%更新Q函數(shù)
Q(now, next)= R(now, next)+ gamma * Qmax gamma = gamma*0.99;%%由隨即策略漸漸變?yōu)殡S機(jī)策略 now = next;% for i = 0:1:5
% Qmax = max(Q(now, :));
% next = find(Q(now, :)== Qmax, 1);% seq = [seq next];% now = next;% end
第四篇：用matlab電力系統(tǒng)潮流計(jì)算
題目：潮流計(jì)算與matlab
教學(xué)單位 電氣信息學(xué)院
姓 名 學(xué) 號(hào)
年 級(jí)
專 業(yè) 電氣工程及其自動(dòng)化
指導(dǎo)教師
職 稱 副教授
摘 要
電力系統(tǒng)穩(wěn)態(tài)分析包括潮流計(jì)算和靜態(tài)安全分析。本文主要運(yùn)用的事潮流計(jì)算，潮流計(jì)算是電力網(wǎng)絡(luò)設(shè)計(jì)與運(yùn)行中最基本的運(yùn)算，對(duì)電力網(wǎng)絡(luò)的各種設(shè)計(jì)方案及各種運(yùn)行方式進(jìn)行潮流計(jì)算，可以得到各種電網(wǎng)各節(jié)點(diǎn)的電壓，并求得網(wǎng)絡(luò)的潮流及網(wǎng)絡(luò)中的各元件的電力損耗，進(jìn)而求得電能損耗。本位就是運(yùn)用潮流計(jì)算具體分析，并有MATLAB仿真。
關(guān)鍵詞： 電力系統(tǒng) 潮流計(jì)算 MATLAB
Abstract Electric power system steady flow calculation and analysis of the static safety analysis.This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power.The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation.Key words: Power system;Flow calculation;MATLAB simulation
目 錄 任務(wù)提出與方案論證....................................................................................................................................2 2 總體設(shè)計(jì)........................................................................................................................................................3 2.1潮流計(jì)算等值電路.............................................................................................................................3 2.2建立電力系統(tǒng)模型.............................................................................................................................3 2.3模型的調(diào)試與運(yùn)行.............................................................................................................................3 3 詳細(xì)設(shè)計(jì)........................................................................................................................................................4 3.1 計(jì)算前提............................................................................................................................................4 3.2手工計(jì)算.............................................................................................................................................7 4設(shè)計(jì)圖及源程序...........................................................................................................................................11 4.1MATLAB仿真.......................................................................................................................................11 4.2潮流計(jì)算源程序...............................................................................................................................11 5 總結(jié).............................................................................................................................................................31 參考文獻(xiàn)..........................................................................................................................................................32 任務(wù)提出與方案論證
潮流計(jì)算是在給定電力系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)、參數(shù)和決定系統(tǒng)運(yùn)行狀態(tài)的邊界條件的情況下確定系統(tǒng)穩(wěn)態(tài)運(yùn)行狀態(tài)的一種基本方法,是電力系統(tǒng)規(guī)劃和運(yùn)營(yíng)中不可缺少的一個(gè)重要組成部分?？梢哉f(shuō),它是電力系統(tǒng)分析中最基本、最重要的計(jì)算,是系統(tǒng)安全、經(jīng)濟(jì)分析和實(shí)時(shí)控制與調(diào)度的基礎(chǔ)。常規(guī)潮流計(jì)算的任務(wù)是根據(jù)給定的運(yùn)行條件和網(wǎng)路結(jié)構(gòu)確定整個(gè)系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。潮流計(jì)算的結(jié)果是電力系統(tǒng)穩(wěn)定計(jì)算和故障分析的基礎(chǔ)。在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式時(shí)，采用離線潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線潮流計(jì)算。是電力系統(tǒng)研究人員長(zhǎng)期研究的一個(gè)課題。它既是對(duì)電力系統(tǒng)規(guī)劃設(shè)計(jì)和運(yùn)行方式的合理性、可靠性及經(jīng)濟(jì)性進(jìn)行定量分析的依據(jù)，又是電力系統(tǒng)靜態(tài)和暫態(tài)穩(wěn)定計(jì)算的基礎(chǔ)。
潮流計(jì)算經(jīng)歷了一個(gè)由手工到應(yīng)用數(shù)字電子計(jì)算機(jī)的發(fā)展過(guò)程，現(xiàn)在的潮流算法都以計(jì)算機(jī)的應(yīng)用為前提用計(jì)算機(jī)進(jìn)行潮流計(jì)算主要步驟在于編制計(jì)算機(jī)程序，這是一項(xiàng)非常復(fù)雜的工作。對(duì)系統(tǒng)進(jìn)行潮流分析，本文利用 MATLAB中的SimpowerSystems工具箱設(shè)計(jì)電力系統(tǒng)，在simulink 環(huán)境下，不僅可以仿真系統(tǒng)的動(dòng)態(tài)過(guò)程，還可以對(duì)系統(tǒng)進(jìn)行穩(wěn)態(tài)潮流分析。
總體設(shè)計(jì)
SimpowerSystems使用Simulink環(huán)境，可以將該系統(tǒng)中的發(fā)電機(jī)、變壓器，線路等模型聯(lián)結(jié)起來(lái)，形成電力系統(tǒng)仿真模擬圖。在加人測(cè)量模塊，并對(duì)各元件的參數(shù)進(jìn)行設(shè)置后，用measurement和sink中的儀器可以觀察各元件的電壓、電流、功率的大小。
2.1潮流計(jì)算等值電路
10MWYN,d114?63MWVA15MWGGGG120MW10kV??p0?15.7kW??ps?73kW?I%?0.5?0?Vs%?10.5YN,d1116MWVA4?63MW“?xd?0.134?x2?0.161?x?0.06?0?cos?N?0.8510kV??p0?11kW??ps?50kW?I%?0.55?0?Vs%?10.5YN,d112?10MWVA35kV32km25MW110kV80km25MW110kV70km110kVYN,d112?20MWVA20MWGGG4?15MW”?xd?0.136?x2?0.16?x?0.073?0?cos?N?0.8??p0?18.6kW??ps?89kW?I%?0.530MW?0?Vs%?10.510kV??p0?15.7kW??ps?73kW?I%?0.5GG?0V%?10.5s?YN,d112?16MWVA63MWVA??p0?44kW??ps?121kW10kV?I%?0.35?0?Vs%?10.5GG35MWYN,Y,d112?10MVA10kVGGG3?12MW1?50MW“?xd?0.128?x2?0.154?x?0.054?0?cos?N?0.852?25MW”?xd?0.128?x2?0.157?x?0.0591?0?cos?N?0.8?x?0.136?x2?0.161?x?0.075?0?cos?N?0.8“d
2.2建立電力系統(tǒng)模型
在Simulink中按照電力系統(tǒng)原型選擇元件進(jìn)行建模。所建立的模型和建立的方法在詳細(xì)設(shè)計(jì)中詳述。
在電力系統(tǒng)模型的建立工程中主要涉及到的是：元器件的選擇及其參數(shù)的設(shè)置；發(fā)電機(jī)選型；變壓器選擇；線路的選擇；負(fù)荷模型的選擇；母線選擇。
2.3模型的調(diào)試與運(yùn)行
建立系統(tǒng)模型，并設(shè)置好參數(shù)以后，就可以在Simulink環(huán)境下進(jìn)行仿真運(yùn)行。運(yùn)行的具體結(jié)果和分析也在詳細(xì)設(shè)計(jì)中詳述。
30km35kV0km31??p0?44kW??ps?121kW?I%?0.35?0?Vs%?10.5??p0?13.2kW??ps?63kW?I0%?0.55?V%?10.5s(1?2)80MW?Vs(2?3)%?6.55MW??Vs(1?3)%?17.53 詳細(xì)設(shè)計(jì)
3.1 計(jì)算前提
首先是發(fā)電機(jī)的參數(shù)計(jì)算，先對(duì)5個(gè)發(fā)電廠簡(jiǎn)化為5臺(tái)發(fā)電機(jī)來(lái)計(jì)算。發(fā)電機(jī)G1：
P1?4?15?60MWQ1?60?tan(arccos0.8)?45MVar發(fā)電機(jī)G2：
P2?4?63?252MWQ2?252?tan(arccos0.85)?156MVar發(fā)電機(jī)G3：
P3?3?12?36MWQ3?36?tan(arccos0.8)?27MVar發(fā)電機(jī)G4：
P4?1?50?50MWQ4?50?tan(arccos0.85)?31MVar發(fā)電機(jī)G5：
P5?2?25?50MWQ5?50?tan(arccos0.8)?37.5MVar
其次是變電站的參數(shù)計(jì)算，我們還是對(duì)7個(gè)變電站簡(jiǎn)化為7臺(tái)變壓器來(lái)計(jì)算。變壓器T1：
2?ps?VN73?11023RT1??10??103?3.450?232SN(16?10)2Vs%?VN10.5?1102XT1??10??10?79.406?SN16?103?S01??p0?j變壓器T2：（雙并聯(lián)）
I0%?SN?(0.0157?j0.0800)MVA 100RT2XT221?ps?VN189?11023???10???103?1.346?2322SN2(20?10)21Vs%?VN110.5?1102???10???10?31.7625? 32SN220?10?S02?2?(?p0?j變壓器T3：（四并聯(lián)）
I0%?SN)?(0.0372?j0.2000)MVA100 1?ps?VN21121?11023RT3???10???103?0.092?2324SN4(63?10)XT31Vs%?VN2110.5?1102???10???10?5.042? 4SN463?103I0%?SN)?(0.1760?j0.8820)MVA100?S03?4?(?p0?j變壓器T4：（雙并聯(lián)）
1RT1?1.7250?21 XT4?XT1?39.7030?2?S04?2?S01?(0.0314?j0.1600)MVART4?變壓器T5：
RT5?4RT3?0.3680?XT5?4XT3?20.168??S05?1?S03?(0.0440?j0.2205)MVA4163?3523???10?0.386? 322(10?10)
變壓器T6：（兩個(gè)三繞組變壓器并聯(lián)）
RT6?1?RT6?2?RT6?31?[Vs(1?2)%?Vs(1?3)%?Vs(2?3)%]?10.7521Vs2%??[Vs(1?2)%?Vs(2?3)%?Vs(1?3)%]??0.25
21Vs3%??[Vs(1?3)%?Vs(2?3)%?Vs(1?2)%]?6.752Vs1%?21Vs1%?VNXT6?1???10?6.584?2SNXT6?2XT6?321Vs2%?VN???10??0.153?2SN21Vs3%?VN???10?4.134? 2SN?S06?2?(?P06?j變壓器T7：（雙并聯(lián)）
I0%?10)?(0.0264?j0.1100)MVA 100 RT7XT721?ps?VN150?3523???10???103?0.306?2322SN2(10?10)21Vs%?VN110.5?352???10???10?6.431?2SN210?103
?S07?2?(?p0?jI0%?SN)?(0.0220?j0.1100)MVA100再次是傳輸線參數(shù)計(jì)算，5條傳輸線的具體計(jì)算如下。
根據(jù)教材查得r0?0.21?/km x0?0.4?/km b0?2.8?10S/km ?6線路L1：
線路L2：
線路L3：（雙回路）
線路L4：
線路L5：（雙回路）RL1?r0?l1?0.21?40?8.4?XL1?x0?l1?0.4?40?16?B?6L1?b0?l1?2.8?10?40?1.12?10?4S ?Q1L1??2BL1V2N??0.6776MVarRL2?r0?l2?0.21?130?27.3?XL2?x0?l2?0.4?130?52?B?6L2?b0?l2?2.8?10?130?3.64?10?4S ?Q1L2??2BL2V2N??2.2022MVarR?12?rl1L30?3?2?0.21?70?7.35?X11L3?2?x0?l3?2?0.4?70?14? BL3?2?b?40?l3?2?2.8?10?6?70?3.92?10S?Q1L3??2B2L3VN??2.3716MVarRL4?r0?l4?0.21?60?12.6?XL4?x0?l4?0.4?60?24?BL4?b0?l4?2.8?10?6?60?1.68?10?4S ?Q12L4??2BL4VN??1.0164MVar
11RL5??r0?l5??0.21?20?2.1?2211XL5??x0?l5??0.4?20?4? 22BL5?2?b0?l5?2?2.8?10?6?20?1.12?10?4S1?QL5??BL3VN2??0.0686MVar23.2手工計(jì)算
FLR1：
P2102?ST1?2(RT1?jXT1)?(3.450?j74.406)?(0.0285?j0.6562)MVA2VN110Sa?10MW??ST1??S01?j?QL1?(10.0442?j0.1142)MVAP2?Q210.04422?0.11422?SL1?(RL1?jXL1)?(8.4?j16)?(0.070?j0.1334)MVAVN21102?ST2P2?Q2402?452?(RT2?jXT2)?(1.346?j31.7625)?(0.4032?j9.5156)MVAVN21102FLR2Sb?SG1?20??ST2?60?j45?20?0.4032?j9.5156?(39.5968?j35.4844)MVASc?Sb?Sa?25?jQL1??SL1?(4.4826?j35.9144)MVA：
?ST3P2?Q22522?1562?(RT3?jXT3)?(0.092?j5.042)?(0.6679?j36.6024)MVA22VN110P?Q4.4931?34.1048(R?jX)?(27.3?j52)?(2.67?j5.0854)MVAL2L222VN1102222Sc'?(4.4931?j34.1048)MVA?SL2?FLRSd?SG2?Sc'?120??ST3??S03?jQL2??SL2?(132.9792?j149.229)MVA3：
?ST4P2?Q262?272?(RT4?jXT4)?(1.725?j39.703)?(0.1091?j2.5101)MVAVN21102P2?Q2133.59552?149.99562?(RL3?jXL3)?(7.35?j14)?(24.51?j46.682)MVA22VN110'Sd?(133.5955?j149.9956)MVA?SL3'Se?SG3?Sd?30?25??ST4??S04?jQL3??SL3?(89.945?j130.0151)MVAFLR4： ?ST5P2?Q2502?312?(RT5?jXT5)?(0.368?j20.168)?(0.1052?j5.7687)MVA22VN110P2?Q292.74872?133.99372?(RL4?jXL4)?(12.6?j24)?(27.654?j52.674)MVA22VN110Se'?(92.7481?j133.9937)MVA?SL4Sf?SG4?Se'?80??ST5??S05?jQL4??SL4?(34.9449?j107.3469)MVAFLR5: 152?ST7?2?(0.306?j6.431)?(0.0562?j1.1812)MVA35Sh?15??ST7??S07?j?QL5?(15.0782?j0.3422)MVA15.07822?0.34222?SL5??(2.1?j4)?(0.3899?j0.743)MVA352Si?Sh??SL5??S06?j?QL5?5?(20.4945?j1.1266)MVA 152?37.52?ST6?3??(0.386?j4.34)?(0.514?j5.7793)MVA35220.65052?0.54512?ST6?2??(0.386?j0.153)?(0.1345?j0.0533)MVA23526.3362?98.73692?ST6?1??(0.386?j6.584)?(3.2905?j56.1256)MVA352Sg?Sf??ST6?1?SG5??ST6?2??ST6?3?Si?35?(25.5114?j194.12)MVA計(jì)算每一個(gè)FLR的功率分布和電壓分布計(jì)算如下： FLR1：
?VT2?PR?QX40?1.346?45?31.7625??12.8970kVVN115 Vb?115??VT2?102.1030kVPR?QX10.0442?8.4?0.1442?16?VL1???0.8489kVVb102.1030Va?Vb??VL1?101.2541kVFLR2：
功率分布：
SL2?ZZ?Z*T3*L2*Sd?T3(Vb?VN)Z?ZL2****?VN?(0.092?j5.042)?(132.9792?j149.229)?1418.6727.392?j57.042T3?(4.8812?j13.8097)MVAST3?ZZ?Z*L2*L2*Sd?T3(Vb?VN)Z?ZL2?VN?(27.3?j52)?(132.9792?j149.229)?1418.6727.392?j57.042T3?(108.687?j122.62)MVA 電壓分布：
Sc1?SL2??SL2?(4.8812?j13.8097)?(2.67?j5.0854)?(7.5512?j8.7243)MVA7.5512?27.3?8.7243?52??2.424kV102.1030Vd?Vb??VL2?102.103?(?2.424)?104.527kV?VL2?功率分布：
FLR3：
SL3?ZZ?Z*T4*L3*Se?T4(VG3?Vd)Z?ZL3****?VN?(1.725?j39.703)?(89.945?j130.0151)?1037.9279.075?j53.73T4?(59.444?j16.846)MVAST4?ZZ?Z*L3*L3*Se?T4(Vb?VN)Z?ZL3?VN?(7.35?j14)?(89.945?j130.0151)?1037.9279.075?j53.73T4?(31.811?j60.1256)MVA 電壓分布：
Se1?SL3??SL3?(59.444?j19.846)?(24.51?j46.682)?(83.954?j26.836)MVA83.954?7.35?26.836?14?9.404kV105.5643Ve?Vd??VL3?96.16kV?VL3?功率分布：
FLR4：
SL4?ZZ?Z*T5*L4*Sf?T5(VG3?Vd)Z?ZL4****?VN=(0.368?j20.168)?(34.9449?j107.3469)?1037.92712.968?j44.168T5?(20.843?j19.689)MVAST4?ZZ?Z*L4*L4*Sf?T5(VG3?Vd)Z?ZL4?VN=(12.6?j24)?(34.9449?j107.3469)?1037.92712.968?j44.168T5?(1.398?j44.389)MVA 電壓分布：
Se1?SL3??SL3?(59.444?j16.846)?(24.51?j46.682)?(83.954?j63.528)MVA83.954?12.6?63.528?24?24.464kV105.5643Ve?Vd??VL3?81.10kV?VL4?FLR5: 這里我們先將f點(diǎn)和發(fā)電機(jī)G5當(dāng)做電源，經(jīng)過(guò)ZT61和ZT63構(gòu)成兩端供電網(wǎng)絡(luò)以g點(diǎn)作為運(yùn)算負(fù)荷進(jìn)行計(jì)算。ST6?ST4(0.386?j4.134)?(20.2656?j70.9293)?(22.0938?37)?35?(3.900?j25.1175)MVA0.772?j10.718(0.386?j6.584)?(20.2656?j70.9293)?(22.0938?37)?35??(16.5061?j91.7905)MVA0.772?j10.718電壓分布：
ST631?ST63??ST63?(16.6421?j97.5698)MVA16.6421?0.386?97.5698?4.134?10.9186kV37Vg?37??VT63?26.0814V?VT63?20.2656?0.386?70.9293?(?0.153)??0.1162kV
26.0814Vi?Vg??VT62?26.1976?VT62?20.4945?2.1?1.1266?4?1.815kV26.1976Vh?Vi??VL5?24.3826?VL5?
4設(shè)計(jì)圖及源程序
4.1MATLAB仿真
相關(guān)的原始數(shù)據(jù)輸入格式如下：
1、B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點(diǎn)編號(hào)。節(jié)點(diǎn)編號(hào)由小到大編寫(xiě)。
2、對(duì)于含有變壓器的支路，第一列為低壓側(cè)節(jié)點(diǎn)編號(hào)，第二列為高壓側(cè)節(jié)點(diǎn)編號(hào)，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理，第三列為支路的串列阻抗參數(shù)，第四列為支路的對(duì)地導(dǎo)納參數(shù)，第五烈為含變壓器支路的變壓器的變比，第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，“0”為不含有變壓器。
3、B2為節(jié)點(diǎn)參數(shù)矩陣，其中第一列為節(jié)點(diǎn)注入發(fā)電功率參數(shù)；第二列為節(jié)點(diǎn)負(fù)荷功率參數(shù)；第三列為節(jié)點(diǎn)電壓參數(shù)；第六列為節(jié)點(diǎn)類型參數(shù)，其中“1”為平衡節(jié)點(diǎn)，“2”為PQ節(jié)點(diǎn)，“3”為PV節(jié)點(diǎn)參數(shù)。
4、X為節(jié)點(diǎn)號(hào)和對(duì)地參數(shù)矩陣。其中第一列為節(jié)點(diǎn)編號(hào)，第二列為節(jié)點(diǎn)對(duì)地參數(shù)。
4.2潮流計(jì)算源程序
%本程序的功能是用牛頓——拉夫遜法進(jìn)行11節(jié)點(diǎn)潮流計(jì)算 clear;n=11;%input('請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=');nl=11;%input('請(qǐng)輸入支路數(shù)：nl=');isb=1;%input('請(qǐng)輸入平衡母線節(jié)點(diǎn)號(hào)：isb=');pr=0.00001;%input('請(qǐng)輸入誤差精度：pr=');B1=[1
0.03512+0.08306i
0.13455i
0;
0.0068+0.18375i
0
1.02381
1;
0.05620+0.13289i
0.05382i
0;
0.00811+0.24549i
0
1.02381
1;
0.05620+0.13289i
0.05382i
0;
0.04215+0.09967i
0.04037i
0;
0.0068+0.18375i
0
1.02381
1;
0.02810+0.06645i
0.10764i
0;
0.05620+0.13289i
0.05382i
0;0.00811+0.24549i
0
1;
0.03512+0.08306i
0.13455i
0] B2=[0
0
1.1
1.1
0
1;
0
0
0
0
2;
0
0.343+0.21256i
0
0
2;
0
0
0
0
2;
0
0.204+0.12638i
0
0
2;
0
0
0
0
2;
0
0.306+0.18962i
0
0
2;
0
0
0
0
2;
0.5
0
1.1
1.1
0
3;
0
0.343+0.21256i
0
0
2;
0
0
0
0
2];% B1矩陣：
1、支路首端號(hào)；
2、末端號(hào)；
3、支路阻抗；
4、支路對(duì)地電納 %
5、支路的變比；
6、支路首端處于K側(cè)為1，1側(cè)為0 % B2矩陣：
1、該節(jié)點(diǎn)發(fā)電機(jī)功率；
2、該節(jié)點(diǎn)負(fù)荷功率；
3、節(jié)點(diǎn)電壓初始值 %
4、PV節(jié)點(diǎn)電壓V的給定值；
5、節(jié)點(diǎn)所接的無(wú)功補(bǔ)償設(shè)備的容量 %
6、節(jié)點(diǎn)分類標(biāo)號(hào)：1為平衡節(jié)點(diǎn)（應(yīng)為1號(hào)節(jié)點(diǎn)）；2為PQ節(jié)點(diǎn)； %
3為PV節(jié)點(diǎn)；
%input('請(qǐng)輸入各節(jié)點(diǎn)參數(shù)形成的矩陣： B2=');Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);% % %--------------------for i=1:nl
%支路數(shù)
if B1(i,6)==0
%左節(jié)點(diǎn)處于1側(cè)
p=B1(i,1);q=B1(i,2);
else
%左節(jié)點(diǎn)處于K側(cè)
p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非對(duì)角元
Y(q,p)=Y(p,q);
%非對(duì)角元
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;%對(duì)角元K側(cè)
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
%對(duì)角元1側(cè)
end %求導(dǎo)納矩陣
disp('導(dǎo)納矩陣 Y=');disp(Y)%---------------------------G=real(Y);B=imag(Y);
%分解出導(dǎo)納陣的實(shí)部和虛部
for i=1:n
%給定各節(jié)點(diǎn)初始電壓的實(shí)部和虛部
e(i)=real(B2(i,3));
f(i)=imag(B2(i,3));
V(i)=B2(i,4);
%PV節(jié)點(diǎn)電壓給定模值
end for i=1:n
%給定各節(jié)點(diǎn)注入功率
S(i)=B2(i,1)-B2(i,2);
%i節(jié)點(diǎn)注入功率SG-SL
B(i,i)=B(i,i)+B2(i,5);%i節(jié)點(diǎn)無(wú)功補(bǔ)償量
end %===================== P=real(S);Q=imag(S);
%分解出各節(jié)點(diǎn)注入的有功和無(wú)功功率 ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;%迭代次數(shù)ICT1、a；不滿足收斂要求的節(jié)點(diǎn)數(shù)IT2 while IT2~=0
% N0=2*n 雅可比矩陣的階數(shù)；N=N0+1擴(kuò)展列
IT2=0;a=a+1;
for i=1:n
if i~=isb
%非平衡節(jié)點(diǎn)
C(i)=0;D(i)=0;
for j1=1:n
C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)
D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)
end
P1=C(i)*e(i)+f(i)*D(i);%節(jié)點(diǎn)功率P計(jì)算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)
Q1=C(i)*f(i)-e(i)*D(i);%節(jié)點(diǎn)功率Q計(jì)算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)%求i節(jié)點(diǎn)有功和無(wú)功功率P',Q'的計(jì)算值
V2=e(i)^2+f(i)^2;
%電壓模平方
%========= 以下針對(duì)非PV節(jié)點(diǎn)來(lái)求取功率差及Jacobi矩陣元素 =========
if B2(i,6)~=3
%非PV節(jié)點(diǎn)
DP=P(i)-P1;
%節(jié)點(diǎn)有功功率差
DQ=Q(i)-Q1;
%節(jié)點(diǎn)無(wú)功功率差
%=============== 以上為除平衡節(jié)點(diǎn)外其它節(jié)點(diǎn)的功率計(jì)算 ================= %================= 求取Jacobi矩陣 ===================
for j1=1:n
if j1~=isb&j1~=i
%非平衡節(jié)點(diǎn)&非對(duì)角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i);% dP/de=-dQ/df
X2=B(i,j1)*e(i)-G(i,j1)*f(i);% dP/df=dQ/de
X3=X2;
% X2=dp/df X3=dQ/de
X4=-X1;
% X1=dP/de X4=dQ/df
p=2*i-1;q=2*j1-1;
J(p,q)=X3;J(p,N)=DQ;m=p+1;
% X3=dQ/de J(p,N)=DQ節(jié)點(diǎn)無(wú)功功率差
J(m,q)=X1;J(m,N)=DP;q=q+1;
% X1=dP/de J(m,N)=DP節(jié)點(diǎn)有功功率差
J(p,q)=X4;J(m,q)=X2;
% X4=dQ/df X2=dp/df
elseif j1==i&j1~=isb %非平衡節(jié)點(diǎn)&對(duì)角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dQ/de
X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df
p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%擴(kuò)展列△Q
m=p+1;
J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%擴(kuò)展列△P
J(m,q)=X2;
end
end
else
%=============== 下面是針對(duì)PV節(jié)點(diǎn)來(lái)求取Jacobi矩陣的元素 ===========
DP=P(i)-P1;
% PV節(jié)點(diǎn)有功誤差
DV=V(i)^2-V2;
% PV節(jié)點(diǎn)電壓誤差
for j1=1:n
if j1~=isb&j1~=i
%非平衡節(jié)點(diǎn)&非對(duì)角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i);
% dP/de
X2=B(i,j1)*e(i)-G(i,j1)*f(i);
% dP/df
X5=0;X6=0;
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;% PV節(jié)點(diǎn)電壓誤差
m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;
% PV節(jié)點(diǎn)有功誤差
J(m,q)=X2;
elseif j1==i&j1~=isb %非平衡節(jié)點(diǎn)&對(duì)角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
X5=-2*e(i);
X6=-2*f(i);
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;% PV節(jié)點(diǎn)電壓誤差
m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;
% PV節(jié)點(diǎn)有功誤差
J(m,q)=X2;
end
end
end
end
end %========= 以上為求雅可比矩陣的各個(gè)元素及擴(kuò)展列的功率差或電壓差 =====================
for k=3:N0
% N0=2*n（從第三行開(kāi)始，第一、二行是平衡節(jié)點(diǎn)）
k1=k+1;N1=N;
% N=N0+1 即 N=2*n+1擴(kuò)展列△P、△Q 或 △U
for k2=k1:N1
% 從k+1列的Jacobi元素到擴(kuò)展列的△P、△Q 或 △U
J(k,k2)=J(k,k2)./(J(k,k)+eps);% 用K行K列對(duì)角元素去除K行K列后的非對(duì)角元素進(jìn)行規(guī)格化
end
J(k,k)=1;
% 對(duì)角元規(guī)格化K行K列對(duì)角元素賦1
%==================== 回代運(yùn)算
=======================================
if k~=3
% 不是第三行
k > 3
k4=k-1;
for k3=3:k4
% 用k3行從第三行開(kāi)始到當(dāng)前行的前一行k4行消去
for k2=k1:N1 %
k3行后各行上三角元素
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算（當(dāng)前行k列元素消為0）
end
%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K14 列元素
J(k3,k)=0;%當(dāng)前行第k列元素已消為0
end
if k==N0
%若已到最后一行
break;
end
%================== 前代運(yùn)算
==================================
for k3=k1:N0
% 從k+1行到2*n最后一行
for k2=k1:N1
% 從k+1列到擴(kuò)展列消去k+1行后各行下三角元素
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算
end
%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素
J(k3,k)=0;%當(dāng)前行第k列元素已消為0
end
else
%是第三行k=3
%====================== 第三行k=3的前代運(yùn)算 ========================
for k3=k1:N0
%從第四行到2n行（最后一行）
for k2=k1:N1
%從第四列到2n+1列（即擴(kuò)展列）
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算（當(dāng)前行3列元素消為0）
end
%用當(dāng)前行K2列元素減去當(dāng)前行3列元素乘以第三行K2列元素
J(k3,k)=0;%當(dāng)前行第3列元素已消為0
end
end
end %====上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣=====
for k=3:2:N0-1
L=(k+1)./2;
e(L)=e(L)-J(k,N);
%修改節(jié)點(diǎn)電壓實(shí)部
k1=k+1;
f(L)=f(L)-J(k1,N);
%修改節(jié)點(diǎn)電壓虛部
end
%------修改節(jié)點(diǎn)電壓-----------
for k=3:N0
DET=abs(J(k,N));
if DET>=pr
%電壓偏差量是否滿足要求
IT2=IT2+1;%不滿足要求的節(jié)點(diǎn)數(shù)加1
end
end
ICT2(a)=IT2;
%不滿足要求的節(jié)點(diǎn)數(shù)
ICT1=ICT1+1;
%迭代次數(shù) end %用高斯消去法解”w=-J*V“
disp('迭代次數(shù)：');disp(ICT1);disp('沒(méi)有達(dá)到精度要求的個(gè)數(shù)：');disp(ICT2);for k=1:n
V(k)=sqrt(e(k)^2+f(k)^2);
%計(jì)算各節(jié)點(diǎn)電壓的模值
sida(k)=atan(f(k)./e(k))*180./pi;
%計(jì)算各節(jié)點(diǎn)電壓的角度
E(k)=e(k)+f(k)*j;
%將各節(jié)點(diǎn)電壓用復(fù)數(shù)表示 end %=============== 計(jì)算各輸出量 =========================== disp('各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(E);
%顯示各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E用復(fù)數(shù)表示 disp('----------------------');disp('各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(V);
%顯示各節(jié)點(diǎn)的電壓大小V的模值 disp('----------------------');disp('各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(sida);
%顯示各節(jié)點(diǎn)的電壓相角 for p=1:n
C(p)=0;
for q=1:n
C(p)=C(p)+conj(Y(p,q))*conj(E(q));%計(jì)算各節(jié)點(diǎn)的注入電流的共軛值
end
S(p)=E(p)*C(p);
%計(jì)算各節(jié)點(diǎn)的功率 S = 電壓 X 注入電流的共軛值 end disp('各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(S);
%顯示各節(jié)點(diǎn)的注入功率
disp('----------------------');disp('各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：');for i=1:nl
p=B1(i,1);q=B1(i,2);
if B1(i,6)==0
Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))...-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
Siz(i)=Si(p,q);
else
Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))...-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
Siz(i)=Si(p,q);
end
disp(Si(p,q));
SSi(p,q)=Si(p,q);
ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];
disp(ZF);
disp('----------------------');end disp('各條支路的末端功率Sj為(順序同您輸入B1時(shí)一致)：');
for i=1:nl
p=B1(i,1);q=B1(i,2);
if B1(i,6)==0
Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))...-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
Sjy(i)=Sj(q,p);
else
Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))...-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
Sjy(i)=Sj(q,p);
end
disp(Sj(q,p));
SSj(q,p)=Sj(q,p);
ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];
disp(ZF);
disp('----------------------');end disp('各條支路的功率損耗DS為(順序同您輸入B1時(shí)一致)：');for i=1:nl
p=B1(i,1);q=B1(i,2);
DS(i)=Si(p,q)+Sj(q,p);
disp(DS(i));
DDS(i)=DS(i);
ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];
disp(ZF);
disp('----------------------');end
%本程序的功能是用牛頓——拉夫遜法進(jìn)行10節(jié)點(diǎn)潮流計(jì)算 %本程序的功能是用牛頓——拉夫遜法進(jìn)行潮流計(jì)算 clear;n=10;%input('請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=');nl=10;%input('請(qǐng)輸入支路數(shù)：nl=');isb=1;%input('請(qǐng)輸入平衡母線節(jié)點(diǎn)號(hào)：isb=');pr=0.00001;%input('請(qǐng)輸入誤差精度：pr=');B1=[1
0.03512+0.08306i
0.13455i
0;
0.0068+0.18375i
0
1.02381
1;
0.05620+0.13289i
0.05382i
0;
0.00811+0.24549i
0
1.02381
1;
0.05620+0.13289i
0.05382i
0;
0.04215+0.09967i
0.04037i
0;
0.0068+0.18375i
0
1.02381
1;
0.02810+0.06645i
0.10764i
0;0.00811+0.24549i
0
1;
0.03512+0.08306i
0.13455i
0] B2=[0
0
1.1
1.1
0
1;
0
0
0
0
2;
0
0.343+0.21256i
0
0
2;
0
0
0
0
2;
0
0.204+0.12638i
0
0
2;
0
0
0
0
2;
0
0.306+0.18962i
0
0
2;
0
0
0
0
2;
0.5
0
1.1
1.1
0
3;
0
0.343+0.21256i
0
0
2];% B1矩陣：
1、支路首端號(hào)；
2、末端號(hào)；
3、支路阻抗；
4、支路對(duì)地電納 %
5、支路的變比；
6、支路首端處于K側(cè)為1，1側(cè)為0 % B2矩陣：
1、該節(jié)點(diǎn)發(fā)電機(jī)功率；
2、該節(jié)點(diǎn)負(fù)荷功率；
3、節(jié)點(diǎn)電壓初始值 %
4、PV節(jié)點(diǎn)電壓V的給定值；
5、節(jié)點(diǎn)所接的無(wú)功補(bǔ)償設(shè)備的容量 %
6、節(jié)點(diǎn)分類標(biāo)號(hào)：1為平衡節(jié)點(diǎn)（應(yīng)為1號(hào)節(jié)點(diǎn)）；2為PQ節(jié)點(diǎn)； %
3為PV節(jié)點(diǎn)；
%input('請(qǐng)輸入各節(jié)點(diǎn)參數(shù)形成的矩陣： B2=');Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);% % %--------------------for i=1:nl
%支路數(shù)
if B1(i,6)==0
%左節(jié)點(diǎn)處于1側(cè)
p=B1(i,1);q=B1(i,2);
else
%左節(jié)點(diǎn)處于K側(cè)
p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非對(duì)角元
Y(q,p)=Y(p,q);
%非對(duì)角元
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;%對(duì)角元K側(cè)
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
%對(duì)角元1側(cè)
end %求導(dǎo)納矩陣
disp('導(dǎo)納矩陣 Y=');disp(Y)%---------------------------G=real(Y);B=imag(Y);
%分解出導(dǎo)納陣的實(shí)部和虛部
for i=1:n
%給定各節(jié)點(diǎn)初始電壓的實(shí)部和虛部
e(i)=real(B2(i,3));
f(i)=imag(B2(i,3));
V(i)=B2(i,4);
%PV節(jié)點(diǎn)電壓給定模值
end for i=1:n
%給定各節(jié)點(diǎn)注入功率
S(i)=B2(i,1)-B2(i,2);
%i節(jié)點(diǎn)注入功率SG-SL
B(i,i)=B(i,i)+B2(i,5);%i節(jié)點(diǎn)無(wú)功補(bǔ)償量
end %===================== P=real(S);Q=imag(S);
%分解出各節(jié)點(diǎn)注入的有功和無(wú)功功率
ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;%迭代次數(shù)ICT1、a；不滿足收斂要求的節(jié)點(diǎn)數(shù)IT2 while IT2~=0
% N0=2*n 雅可比矩陣的階數(shù)；N=N0+1擴(kuò)展列
IT2=0;a=a+1;
for i=1:n
if i~=isb
%非平衡節(jié)點(diǎn)
C(i)=0;D(i)=0;
for j1=1:n
C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)
D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)
end
P1=C(i)*e(i)+f(i)*D(i);%節(jié)點(diǎn)功率P計(jì)算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)
Q1=C(i)*f(i)-e(i)*D(i);%節(jié)點(diǎn)功率Q計(jì)算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)%求i節(jié)點(diǎn)有功和無(wú)功功率P',Q'的計(jì)算值
V2=e(i)^2+f(i)^2;
%電壓模平方
%========= 以下針對(duì)非PV節(jié)點(diǎn)來(lái)求取功率差及Jacobi矩陣元素 =========
if B2(i,6)~=3
%非PV節(jié)點(diǎn)
DP=P(i)-P1;
%節(jié)點(diǎn)有功功率差
DQ=Q(i)-Q1;
%節(jié)點(diǎn)無(wú)功功率差
%=============== 以上為除平衡節(jié)點(diǎn)外其它節(jié)點(diǎn)的功率計(jì)算 ================= %================= 求取Jacobi矩陣 ===================
for j1=1:n
if j1~=isb&j1~=i
%非平衡節(jié)點(diǎn)&非對(duì)角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i);% dP/de=-dQ/df
X2=B(i,j1)*e(i)-G(i,j1)*f(i);% dP/df=dQ/de
X3=X2;
% X2=dp/df X3=dQ/de
X4=-X1;
% X1=dP/de X4=dQ/df
p=2*i-1;q=2*j1-1;
J(p,q)=X3;J(p,N)=DQ;m=p+1;
% X3=dQ/de J(p,N)=DQ節(jié)點(diǎn)無(wú)功功率差
J(m,q)=X1;J(m,N)=DP;q=q+1;
% X1=dP/de J(m,N)=DP節(jié)點(diǎn)有功功率差
J(p,q)=X4;J(m,q)=X2;
% X4=dQ/df X2=dp/df
elseif j1==i&j1~=isb %非平衡節(jié)點(diǎn)&對(duì)角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dQ/de
X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df
p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%擴(kuò)展列△Q
m=p+1;
J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%擴(kuò)展列△P
J(m,q)=X2;
end
end
else
%=============== 下面是針對(duì)PV節(jié)點(diǎn)來(lái)求取Jacobi矩陣的元素
===========
DP=P(i)-P1;
% PV節(jié)點(diǎn)有功誤差
DV=V(i)^2-V2;
% PV節(jié)點(diǎn)電壓誤差
for j1=1:n
if j1~=isb&j1~=i
%非平衡節(jié)點(diǎn)&非對(duì)角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i);
% dP/de
X2=B(i,j1)*e(i)-G(i,j1)*f(i);
% dP/df
X5=0;X6=0;
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;% PV節(jié)點(diǎn)電壓誤差
m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;
% PV節(jié)點(diǎn)有功誤差
J(m,q)=X2;
elseif j1==i&j1~=isb %非平衡節(jié)點(diǎn)&對(duì)角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
X5=-2*e(i);
X6=-2*f(i);
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;% PV節(jié)點(diǎn)電壓誤差
m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;
% PV節(jié)點(diǎn)有功誤差
J(m,q)=X2;
end
end
end
end
end %========= 以上為求雅可比矩陣的各個(gè)元素及擴(kuò)展列的功率差或電壓差 =====================
for k=3:N0
% N0=2*n（從第三行開(kāi)始，第一、二行是平衡節(jié)點(diǎn)）
k1=k+1;N1=N;
% N=N0+1 即 N=2*n+1擴(kuò)展列△P、△Q 或 △U
for k2=k1:N1
% 從k+1列的Jacobi元素到擴(kuò)展列的△P、△Q 或 △U
J(k,k2)=J(k,k2)./J(k,k);% 用K行K列對(duì)角元素去除K行K列后的非對(duì)角元素進(jìn)行規(guī)格化
end
J(k,k)=1;
% 對(duì)角元規(guī)格化K行K列對(duì)角元素賦1
%==================== 回代運(yùn)算
=======================================
if k~=3
% 不是第三行
k > 3
k4=k-1;
for k3=3:k4
% 用k3行從第三行開(kāi)始到當(dāng)前行的前一行k4行消
去
for k2=k1:N1 %
k3行后各行上三角元素
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算（當(dāng)前行k列元素消為0）
end
%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素
J(k3,k)=0;%當(dāng)前行第k列元素已消為0
end
if k==N0
%若已到最后一行
break;
end
%================== 前代運(yùn)算
==================================
for k3=k1:N0
% 從k+1行到2*n最后一行
for k2=k1:N1
% 從k+1列到擴(kuò)展列消去k+1行后各行下三角元素
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算
end
%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素
J(k3,k)=0;%當(dāng)前行第k列元素已消為0
end
else
%是第三行k=3
%====================== 第三行k=3的前代運(yùn)算 ========================
for k3=k1:N0
%從第四行到2n行（最后一行）
for k2=k1:N1
%從第四列到2n+1列（即擴(kuò)展列）
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算（當(dāng)前行3列元素消為0）
end
%用當(dāng)前行K2列元素減去當(dāng)前行3列元素乘以第三行K2列元素
J(k3,k)=0;%當(dāng)前行第3列元素已消為0
end
end
end %====上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣=====
for k=3:2:N0-1
L=(k+1)./2;
e(L)=e(L)-J(k,N);
%修改節(jié)點(diǎn)電壓實(shí)部
k1=k+1;
f(L)=f(L)-J(k1,N);
%修改節(jié)點(diǎn)電壓虛部
end
%------修改節(jié)點(diǎn)電壓-----------
for k=3:N0
DET=abs(J(k,N));
if DET>=pr
%電壓偏差量是否滿足要求
IT2=IT2+1;%不滿足要求的節(jié)點(diǎn)數(shù)加1
end
end
ICT2(a)=IT2;
%不滿足要求的節(jié)點(diǎn)數(shù)
ICT1=ICT1+1;
%迭代次數(shù) end %用高斯消去法解”w=-J*V“ disp('迭代次數(shù)：');disp(ICT1);disp('沒(méi)有達(dá)到精度要求的個(gè)數(shù)：');disp(ICT2);for k=1:n
V(k)=sqrt(e(k)^2+f(k)^2);
%計(jì)算各節(jié)點(diǎn)電壓的模值
sida(k)=atan(f(k)./e(k))*180./pi;
%計(jì)算各節(jié)點(diǎn)電壓的角度
E(k)=e(k)+f(k)*j;
%將各節(jié)點(diǎn)電壓用復(fù)數(shù)表示 end %=============== 計(jì)算各輸出量 =========================== disp('各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(E);
%顯示各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E用復(fù)數(shù)表示 disp('----------------------');disp('各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(V);
%顯示各節(jié)點(diǎn)的電壓大小V的模值 disp('----------------------');disp('各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(sida);
%顯示各節(jié)點(diǎn)的電壓相角 for p=1:n
C(p)=0;
for q=1:n
C(p)=C(p)+conj(Y(p,q))*conj(E(q));%計(jì)算各節(jié)點(diǎn)的注入電流的共軛值
end
S(p)=E(p)*C(p);
%計(jì)算各節(jié)點(diǎn)的功率 S = 電壓 X 注入電流的共軛值 end disp('各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(S);
%顯示各節(jié)點(diǎn)的注入功率
disp('----------------------');disp('各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：');for i=1:nl
p=B1(i,1);q=B1(i,2);
if B1(i,6)==0
Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))...-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
Siz(i)=Si(p,q);
else
Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))...-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
Siz(i)=Si(p,q);
end
disp(Si(p,q));
SSi(p,q)=Si(p,q);
ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];
disp(ZF);
disp('----------------------');end disp('各條支路的末端功率Sj為(順序同您輸入B1時(shí)一致)：');for i=1:nl
p=B1(i,1);q=B1(i,2);
if B1(i,6)==0
Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))...-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
Sjy(i)=Sj(q,p);
else
Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))...-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
Sjy(i)=Sj(q,p);
end
disp(Sj(q,p));
SSj(q,p)=Sj(q,p);
ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];
disp(ZF);
disp('----------------------');end disp('各條支路的功率損耗DS為(順序同您輸入B1時(shí)一致)：');for i=1:nl
p=B1(i,1);q=B1(i,2);
DS(i)=Si(p,q)+Sj(q,p);
disp(DS(i));
DDS(i)=DS(i);
ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];
disp(ZF);
disp('----------------------');end
%本程序的功能是用牛頓——拉夫遜法進(jìn)行12節(jié)點(diǎn)潮流計(jì)算 %本程序的功能是用牛頓——拉夫遜法進(jìn)行潮流計(jì)算 clear;n=12;%input('請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=');nl=12;%input('請(qǐng)輸入支路數(shù)：nl=');isb=1;%input('請(qǐng)輸入平衡母線節(jié)點(diǎn)號(hào)：isb=');pr=0.00001;%input('請(qǐng)輸入誤差精度：pr=');B1=[1
0.03512+0.08306i
0.13455i
0;
0.0068+0.18375i
0
1.02381
1;
0.05620+0.13289i
0.05382i
0;
0.00811+0.24549i
0
1.02381
1;
0.05620+0.13289i
0.05382i
0;
0.04215+0.09967i
0.04037i
0;
0.0068+0.18375i
0
1.02381
1;
0.02810+0.06645i
0.10764i
0;
0.05620+0.13289i
0.05382i
0;0.00811+0.24549i
0
1;
0.03512+0.08306i
0.13455i
0;
0.03512+0.08306i
0.13455i
0] B2=[0
0
1.1
1.1
0
1;
0
0
0
0
2;
0
0.343+0.21256i
0
0
2;
0
0
0
0
2;
0
0.204+0.12638i
0
0
2;
0
0
0
0
2;
0
0.306+0.18962i
0
0
2;
0
0
0
0
2;
0.5
0
1.1
1.1
0
3;
0
0.343+0.21256i
0
0
2;
0
0
0
0
2;
0
0
0
0
2];% B1矩陣：
1、支路首端號(hào)；
2、末端號(hào)；
3、支路阻抗；
4、支路對(duì)地電納 %
5、支路的變比；
6、支路首端處于K側(cè)為1，1側(cè)為0 % B2矩陣：
1、該節(jié)點(diǎn)發(fā)電機(jī)功率；
2、該節(jié)點(diǎn)負(fù)荷功率；
3、節(jié)點(diǎn)電壓初始值 %
4、PV節(jié)點(diǎn)電壓V的給定值；
5、節(jié)點(diǎn)所接的無(wú)功補(bǔ)償設(shè)備的容量 %
6、節(jié)點(diǎn)分類標(biāo)號(hào)：1為平衡節(jié)點(diǎn)（應(yīng)為1號(hào)節(jié)點(diǎn)）；2為PQ節(jié)點(diǎn)； %
3為PV節(jié)點(diǎn)；
%input('請(qǐng)輸入各節(jié)點(diǎn)參數(shù)形成的矩陣： B2=');Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);% % %--------------------for i=1:nl
%支路數(shù)
if B1(i,6)==0
%左節(jié)點(diǎn)處于1側(cè)
p=B1(i,1);q=B1(i,2);
else
%左節(jié)點(diǎn)處于K側(cè)
p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非對(duì)角元
Y(q,p)=Y(p,q);
%非對(duì)角元
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;%對(duì)角元K側(cè)
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
%對(duì)角元1側(cè)
end %求導(dǎo)納矩陣
disp('導(dǎo)納矩陣 Y=');disp(Y)%---------------------------G=real(Y);B=imag(Y);
%分解出導(dǎo)納陣的實(shí)部和虛部
for i=1:n
%給定各節(jié)點(diǎn)初始電壓的實(shí)部和虛部
e(i)=real(B2(i,3));
f(i)=imag(B2(i,3));
V(i)=B2(i,4);
%PV節(jié)點(diǎn)電壓給定模值
end for i=1:n
%給定各節(jié)點(diǎn)注入功率
S(i)=B2(i,1)-B2(i,2);
%i節(jié)點(diǎn)注入功率SG-SL
B(i,i)=B(i,i)+B2(i,5);%i節(jié)點(diǎn)無(wú)功補(bǔ)償量
end %===================== P=real(S);Q=imag(S);
%分解出各節(jié)點(diǎn)注入的有功和無(wú)功功率 ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;%迭代次數(shù)ICT1、a；不滿足收斂要求的節(jié)點(diǎn)數(shù)IT2 while IT2~=0
% N0=2*n 雅可比矩陣的階數(shù)；N=N0+1擴(kuò)展列
IT2=0;a=a+1;
for i=1:n
if i~=isb
%非平衡節(jié)點(diǎn)
C(i)=0;D(i)=0;
for j1=1:n
C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)
D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)
end
P1=C(i)*e(i)+f(i)*D(i);%節(jié)點(diǎn)功率P計(jì)算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)
Q1=C(i)*f(i)-e(i)*D(i);%節(jié)點(diǎn)功率Q計(jì)算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)%求i節(jié)點(diǎn)有功和無(wú)功功率P',Q'的計(jì)算值
V2=e(i)^2+f(i)^2;
%電壓模平方
%========= 以下針對(duì)非PV節(jié)點(diǎn)來(lái)求取功率差及Jacobi矩陣元素 =========
if B2(i,6)~=3
%非PV節(jié)點(diǎn)
DP=P(i)-P1;
%節(jié)點(diǎn)有功功率差
DQ=Q(i)-Q1;
%節(jié)點(diǎn)無(wú)功功率差
%=============== 以上為除平衡節(jié)點(diǎn)外其它節(jié)點(diǎn)的功率計(jì)算 ================= %================= 求取Jacobi矩陣 ===================
for j1=1:n
if j1~=isb&j1~=i
%非平衡節(jié)點(diǎn)&非對(duì)角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i);% dP/de=-dQ/df
X2=B(i,j1)*e(i)-G(i,j1)*f(i);% dP/df=dQ/de
X3=X2;
% X2=dp/df X3=dQ/de
X4=-X1;
% X1=dP/de X4=dQ/df
p=2*i-1;q=2*j1-1;
J(p,q)=X3;J(p,N)=DQ;m=p+1;
% X3=dQ/de J(p,N)=DQ節(jié)點(diǎn)無(wú)功功率差
J(m,q)=X1;J(m,N)=DP;q=q+1;
% X1=dP/de J(m,N)=DP節(jié)點(diǎn)有功功率差
J(p,q)=X4;J(m,q)=X2;
% X4=dQ/df X2=dp/df
elseif j1==i&j1~=isb %非平衡節(jié)點(diǎn)&對(duì)角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dQ/de
X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df
p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%擴(kuò)展列△Q
m=p+1;
J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%擴(kuò)展列△P
J(m,q)=X2;
end
end
else
%=============== 下面是針對(duì)PV節(jié)點(diǎn)來(lái)求取Jacobi矩陣的元素 ===========
DP=P(i)-P1;
% PV節(jié)點(diǎn)有功誤差
DV=V(i)^2-V2;
% PV節(jié)點(diǎn)電壓誤差
for j1=1:n
if j1~=isb&j1~=i
%非平衡節(jié)點(diǎn)&非對(duì)角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i);
% dP/de
X2=B(i,j1)*e(i)-G(i,j1)*f(i);
% dP/df
X5=0;X6=0;
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;% PV節(jié)點(diǎn)電壓誤差
m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;
% PV節(jié)點(diǎn)有功誤差
J(m,q)=X2;
elseif j1==i&j1~=isb %非平衡節(jié)點(diǎn)&對(duì)角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
X5=-2*e(i);
X6=-2*f(i);
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;% PV節(jié)點(diǎn)電壓誤差
m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;
% PV節(jié)點(diǎn)有功誤差
J(m,q)=X2;
end
end
end
end
end %========= 以上為求雅可比矩陣的各個(gè)元素及擴(kuò)展列的功率差或電壓差 =====================
for k=3:N0
% N0=2*n（從第三行開(kāi)始，第一、二行是平衡節(jié)點(diǎn)）
k1=k+1;N1=N;
% N=N0+1 即 N=2*n+1擴(kuò)展列△P、△Q 或 △U
for k2=k1:N1
% 從k+1列的Jacobi元素到擴(kuò)展列的△P、△Q
或 △U
J(k,k2)=J(k,k2)./(J(k,k)+eps);% 用K行K列對(duì)角元素去除K行K列后的非對(duì)角元素進(jìn)行規(guī)格化
end
J(k,k)=1;
% 對(duì)角元規(guī)格化K行K列對(duì)角元素賦1
%==================== 回代運(yùn)算
=======================================
if k~=3
% 不是第三行
k > 3
k4=k-1;
for k3=3:k4
% 用k3行從第三行開(kāi)始到當(dāng)前行的前一行k4行消去
for k2=k1:N1 %
k3行后各行上三角元素
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算（當(dāng)前行k列元素消為0）
end
%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素
J(k3,k)=0;%當(dāng)前行第k列元素已消為0
end
if k==N0
%若已到最后一行
break;
end
%================== 前代運(yùn)算
==================================
for k3=k1:N0
% 從k+1行到2*n最后一行
for k2=k1:N1
% 從k+1列到擴(kuò)展列消去k+1行后各行下三角元素
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算
end
%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素
J(k3,k)=0;%當(dāng)前行第k列元素已消為0
end
else
%是第三行k=3
%====================== 第三行k=3的前代運(yùn)算 ========================
for k3=k1:N0
%從第四行到2n行（最后一行）
for k2=k1:N1
%從第四列到2n+1列（即擴(kuò)展列）
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算（當(dāng)前行3列元素消為0）
end
%用當(dāng)前行K2列元素減去當(dāng)前行3列元素乘以第三行K2列元素
J(k3,k)=0;%當(dāng)前行第3列元素已消為0
end
end
end %====上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣=====
for k=3:2:N0-1
L=(k+1)./2;
e(L)=e(L)-J(k,N);
%修改節(jié)點(diǎn)電壓實(shí)部
k1=k+1;
f(L)=f(L)-J(k1,N);
%修改節(jié)點(diǎn)電壓虛部
end
%------修改節(jié)點(diǎn)電壓-----------
for k=3:N0
DET=abs(J(k,N));
if DET>=pr
%電壓偏差量是否滿足要求
IT2=IT2+1;%不滿足要求的節(jié)點(diǎn)數(shù)加1
end
end
ICT2(a)=IT2;
%不滿足要求的節(jié)點(diǎn)數(shù)
ICT1=ICT1+1;
%迭代次數(shù) end %用高斯消去法解”w=-J*V" disp('迭代次數(shù)：');disp(ICT1);disp('沒(méi)有達(dá)到精度要求的個(gè)數(shù)：');disp(ICT2);for k=1:n
V(k)=sqrt(e(k)^2+f(k)^2);
%計(jì)算各節(jié)點(diǎn)電壓的模值
sida(k)=atan(f(k)./e(k))*180./pi;
%計(jì)算各節(jié)點(diǎn)電壓的角度
E(k)=e(k)+f(k)*j;
%將各節(jié)點(diǎn)電壓用復(fù)數(shù)表示 end %=============== 計(jì)算各輸出量 =========================== disp('各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(E);
%顯示各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E用復(fù)數(shù)表示 disp('----------------------');disp('各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(V);
%顯示各節(jié)點(diǎn)的電壓大小V的模值 disp('----------------------');disp('各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(sida);
%顯示各節(jié)點(diǎn)的電壓相角 for p=1:n
C(p)=0;
for q=1:n
C(p)=C(p)+conj(Y(p,q))*conj(E(q));%計(jì)算各節(jié)點(diǎn)的注入電流的共軛值
end
S(p)=E(p)*C(p);
%計(jì)算各節(jié)點(diǎn)的功率 S = 電壓 X 注入電流的共軛值 end disp('各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(S);
%顯示各節(jié)點(diǎn)的注入功率
disp('----------------------');disp('各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：');for i=1:nl
p=B1(i,1);q=B1(i,2);
if B1(i,6)==0
Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))...-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
Siz(i)=Si(p,q);
else
Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))...-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
Siz(i)=Si(p,q);
end
disp(Si(p,q));
SSi(p,q)=Si(p,q);
ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];
disp(ZF);
disp('----------------------');end disp('各條支路的末端功率Sj為(順序同您輸入B1時(shí)一致)：');for i=1:nl
p=B1(i,1);q=B1(i,2);
if B1(i,6)==0
Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))...-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
Sjy(i)=Sj(q,p);
else
Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))...-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
Sjy(i)=Sj(q,p);
end
disp(Sj(q,p));
SSj(q,p)=Sj(q,p);
ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];
disp(ZF);
disp('----------------------');end disp('各條支路的功率損耗DS為(順序同您輸入B1時(shí)一致)：');for i=1:nl
p=B1(i,1);q=B1(i,2);
DS(i)=Si(p,q)+Sj(q,p);
disp(DS(i));
DDS(i)=DS(i);
ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];
disp(ZF);
disp('----------------------');end
如果源程序的運(yùn)行結(jié)果需要作圖可用下面的程序
figure(1);subplot(1,2,1);plot(V);xlabel('節(jié)點(diǎn)號(hào)');ylabel('電壓標(biāo)幺值');grid on;subplot(1,2,2);plot(sida);xlabel('節(jié)點(diǎn)號(hào)');ylabel('電壓角度');grid on;figure(2);subplot(2,2,1);P=real(S);Q=imag(S);bar(P);xlabel('節(jié)點(diǎn)號(hào)');ylabel('節(jié)點(diǎn)注入有功');grid on;subplot(2,2,2);bar(Q);xlabel('節(jié)點(diǎn)號(hào)');ylabel('節(jié)點(diǎn)注入無(wú)功');grid on;subplot(2,2,3);P1=real(Siz);Q1=imag(Siz);bar(P1);xlabel('支路號(hào)');ylabel('支路首端注入有功');grid on;subplot(2,2,4);bar(Q1);xlabel('支路號(hào)');ylabel('支路首端注入無(wú)功');grid on;
總結(jié)
通過(guò)本次課程設(shè)計(jì)讓我有復(fù)習(xí)了一次潮流計(jì)算的相關(guān)知識(shí)，跟家清晰了什么事潮流計(jì)算以及潮流計(jì)算的在電力系統(tǒng)的重要性。電力系統(tǒng)的穩(wěn)定運(yùn)行狀況即是正常運(yùn)行狀況，是指電力系統(tǒng)在穩(wěn)定運(yùn)行條件下電壓、功率的分布，也稱為潮流分布。電力系統(tǒng)分析的潮流計(jì)算是電力系統(tǒng)分析的一個(gè)重要的部分。通過(guò)對(duì)電力系統(tǒng)潮流分布的分析和計(jì)算，可進(jìn)一步對(duì)系統(tǒng)運(yùn)行的安全性，經(jīng)濟(jì)性進(jìn)行分析、評(píng)估，提出改進(jìn)措施。同時(shí)潮流分布也是電力系統(tǒng)規(guī)劃設(shè)計(jì)的一項(xiàng)基礎(chǔ)工作。
整個(gè)計(jì)算過(guò)程的模型建立并不是十分復(fù)雜，但計(jì)算過(guò)程十分繁瑣、計(jì)算量相當(dāng)?shù)拇?，而且由于枝?jié)太多很容易算錯(cuò)。不過(guò)在計(jì)算潮流計(jì)算的過(guò)程中卻對(duì)以往學(xué)過(guò)的電力系統(tǒng)分析的相關(guān)知識(shí)進(jìn)行了一次較為深入的復(fù)習(xí)。而且整個(gè)計(jì)算對(duì)計(jì)算量的要求很大，鍛煉了我們的計(jì)算能力。而且對(duì)細(xì)節(jié)的把握也得到了鍛煉，做題的精細(xì)程度得到了提高。
參考文獻(xiàn)
[1]何仰贊, 溫增銀《電力系統(tǒng)分析》（第三版）［M].華中科技大學(xué)，2002 [2]http://baike.baidu.com/view/627420.[3]王守相，劉玉田 電力系統(tǒng)潮流計(jì)算研究現(xiàn)狀--《山東電力技術(shù)》1996年05期
[4]何仰贊，溫增銀.電力系統(tǒng)分析（上冊(cè)）第三版［M].湖北：華中科技大學(xué)出版社，2002 [5] 劉同娟.MATLAB在電路分析中的應(yīng)用.電氣電子教學(xué)學(xué)報(bào).2002 [6] 西安交通大學(xué)等.電力系統(tǒng)計(jì)算[M].北京：水利電力出版社，1993.12 [7] 李光琦.電力系統(tǒng)暫態(tài)分析[M].北京： 水利電力出版社，2002.5 [8]何仰贊，溫增銀.電力系統(tǒng)分析（下冊(cè)）第三版［M].湖北：華中科技大學(xué)出版社，2002 [9]韋化，李濱，杭乃善，等．大規(guī)模水一火電力系統(tǒng)最優(yōu)潮流的現(xiàn)代內(nèi)點(diǎn)算法實(shí)現(xiàn)[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2003．23(6)：13一l8．
[10]Chen Luo—nan，Suzuki Hideki，Katou Kazuo．Mean fieldtheory for optimal power flow[J]．IEEE Transactions OilPower Systems，1997，12(4)：1481·1486
第五篇：基于MATLAB的電力系統(tǒng)潮流計(jì)算
基于MATLAB的電力系統(tǒng)潮流計(jì)算
%簡(jiǎn)單潮流計(jì)算的小程序，相關(guān)的原始數(shù)據(jù)數(shù)據(jù)數(shù)據(jù)輸入格式如下：
%B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點(diǎn)編號(hào)。節(jié)點(diǎn)編號(hào)由小到大編寫(xiě) %對(duì)于含有變壓器的支路，第一列為低壓側(cè)節(jié)點(diǎn)編號(hào)，第二列為高壓側(cè)節(jié)點(diǎn) %編號(hào)，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理。%第三列為支路的串列阻抗參數(shù)。%第四列為支路的對(duì)地導(dǎo)納參數(shù)。
%第五烈為含變壓器支路的變壓器的變比
%第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，%“0”為不含有變壓器。
%B2為節(jié)點(diǎn)參數(shù)矩陣，其中第一列為節(jié)點(diǎn)注入發(fā)電功率參數(shù)；第二列為節(jié)點(diǎn) %負(fù)荷功率參數(shù)；第三列為節(jié)點(diǎn)電壓參數(shù)；第六列為節(jié)點(diǎn)類型參數(shù)，其中 %“1”為平衡節(jié)點(diǎn)，“2”為PQ節(jié)點(diǎn)，“3”為PV節(jié)點(diǎn)參數(shù)。
%X為節(jié)點(diǎn)號(hào)和對(duì)地參數(shù)矩陣。其中第一列為節(jié)點(diǎn)編號(hào)，第二列為節(jié)點(diǎn)對(duì)地 %參數(shù)。
n=input('請(qǐng)輸入節(jié)點(diǎn)數(shù):n=');n1=input('請(qǐng)輸入支路數(shù):n1=');isb=input('請(qǐng)輸入平衡節(jié)點(diǎn)號(hào):isb=');pr=input('請(qǐng)輸入誤差精度:pr=');B1=input('請(qǐng)輸入支路參數(shù):B1=');B2=input('請(qǐng)輸入節(jié)點(diǎn)參數(shù):B2=');X=input('節(jié)點(diǎn)號(hào)和對(duì)地參數(shù):X=');Y=zeros(n);Times=1;%置迭代次數(shù)為初始值 %創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣 for i=1:n1 if B1(i,6)==0 %不含變壓器的支路 p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);else %含有變壓器的支路 p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3);Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));end end Y OrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1);%將OrgS、DetaS初始化
%創(chuàng)建OrgS，用于存儲(chǔ)初始功率參數(shù) h=0;j=0;for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理 if i~=isb&B2(i,6)==2 h=h+1;for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));end end end for i=1:n %對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0 if i~=isb&B2(i,6)==3 h=h+1;for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));end end end OrgS %創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓 PVU=zeros(n-h-1,1);t=0;for i=1:n if B2(i,6)==3 t=t+1;PVU(t,1)=B2(i,3);end end PVU %創(chuàng)建DetaS，用于存儲(chǔ)有功功率、無(wú)功功率和電壓幅值的不平衡量 h=0;for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理 if i~=isb&B2(i,6)==2 h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);end end t=0;for i=1:n %對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0 if i~=isb&B2(i,6)==3 h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;end end DetaS %創(chuàng)建I，用于存儲(chǔ)節(jié)點(diǎn)電流參數(shù) i=zeros(n-1,1);h=0;for i=1:n if i~=isb h=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));end end I %創(chuàng)建Jacbi(雅可比矩陣)Jacbi=zeros(2*n-2);h=0;k=0;for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理 if B2(i,6)==2 h=h+1;for j=1:n if j~=isb k=k+1;if i==j %對(duì)角元素的處理
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));else %非對(duì)角元素的處理
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);end if k==(n-1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行
k=0;end end end end end k=0;for i=1:n %對(duì)PV節(jié)點(diǎn)的處理 if B2(i,6)==3 h=h+1;for j=1:n if j~=isb k=k+1;if i==j %對(duì)角元素的處理
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));else %非對(duì)角元素的處理
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;end if k==(n-1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行
k=0;end end end end end Jacbi %求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量 DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU %修正節(jié)點(diǎn)電壓 j=0;for i=1:n %對(duì)PQ節(jié)點(diǎn)處理 if B2(i,6)==2 j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);end end for i=1:n %對(duì)PV節(jié)點(diǎn)的處理 if B2(i,6)==3 j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);end end B2 %開(kāi)始循環(huán)********************************************************************** while abs(max(DetaU))>pr OrgS=zeros(2*n-2,1);%！!初始功率參數(shù)在迭代過(guò)程中是不累加的，所以在這里必須將其初始化為零矩陣 h=0;j=0;for i=1:n if i~=isb&B2(i,6)==2 h=h+1;for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));end end end for i=1:n if i~=isb&B2(i,6)==3 h=h+1;for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));end end end OrgS %創(chuàng)建DetaS h=0;for i=1:n if i~=isb&B2(i,6)==2 h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);end end t=0;for i=1:n if i~=isb&B2(i,6)==3 h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;end end DetaS %創(chuàng)建I i=zeros(n-1,1);h=0;for i=1:n if i~=isb h=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));end end I %創(chuàng)建Jacbi Jacbi=zeros(2*n-2);h=0;k=0;for i=1:n if B2(i,6)==2 h=h+1;for j=1:n if j~=isb k=k+1;if i==j Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));else Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);end if k==(n-1)k=0;end end end end end k=0;for i=1:n if B2(i,6)==3 h=h+1;for j=1:n if j~=isb k=k+1;if i==j Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));else Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;end if k==(n-1)k=0;end end end end end Jacbi DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU %修正節(jié)點(diǎn)電壓 j=0;for i=1:n if B2(i,6)==2 j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);end end for i=1:n if B2(i,6)==3 j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);end end B2 Times=Times+1;%迭代次數(shù)加1 end Times 一個(gè)原始數(shù)據(jù)的例子 節(jié)點(diǎn)數(shù) 5 支路數(shù) 5平衡節(jié)點(diǎn)編號(hào) 5 精度pr 0.000001 B1(支路參數(shù)矩陣)[1 2 0.04+0.25i 0.5i 1 0;1 3 0.1+0.35i 0 1 0;2 3 0.08+0.30i 0.5i 1 0;4 2 0.015i 0 1.05 1;5 3 0.03i 0 1.05 1] B2(節(jié)點(diǎn)參數(shù)矩陣)[0-1.6-0.8i 1 0 0 2;0-2-1i 1 0 0 2;0-3.7-1.3i 1 0 0 2;0 5+0i 1.05 1.05 0 3;0 0 1.05 1.05 0 1] X(節(jié)點(diǎn)號(hào)和對(duì)地參數(shù))[1 0;2 0;3 0;4 0;5 0]
電力系統(tǒng)潮流計(jì)算
——9結(jié)點(diǎn)算例-PQ法
原始數(shù)據(jù)錄入data.txt文檔:
標(biāo)號(hào),起始結(jié)點(diǎn),終止結(jié)點(diǎn),支路電阻參數(shù),支路電抗參數(shù),支路對(duì)地導(dǎo)納參數(shù) 1,2,5,0.0,0.063,0.0, 2,5,9,0.019,0.072,0.075, 3,6,9,0.012,0.101,0.105, 4,3,6,0.0,0.059,0.0, 5,6,8,0.039,0.17,0.179, 6,4,8,0.017,0.092,0.079, 7,5,7,0.032,0.161,0.153, 8,4,7,0.01,0.085,0.088, 9,1,4,0.0,0.058,0.0, 潮流程序chaoliu.txt文檔: #include ; float x1[N-1],x2 ;for(i=1;i guass(1,N-1,y1,B,x1);for(i=1;i guass(N-M,M,y2,B,x2);for(i=N-M;i else { kp=0;if(kq==0)val(u,g,b,r,ku,kr,h);else goto top;} } } void val(float u[N],float g[N][N],float b[N][N],float r[N],int ku, int kr,float h[N][N]){ float ps=0,pv1=0,pv2=0;float qs=0,qv1=0,qv2=0;float p[N][N]={0};float q[N][N]={0};float s[N][N];float dp[N][N]={0};float dq[N][N]={0};float ds[N][N];float dSp=0,dSq=0;int i,j;FILE *fp1;printf(”n=====ping heng jie dian gong lv =====n“);getch();for(i=0;i printf(”n=======shu ju bao cun=====n“);fp1=fopen(”jieguo.txt“,”w+“);{ fprintf(fp1,”xian lu cao liu:n“);for(i=0;i
點(diǎn)擊咨詢 