【正文】 * 1 4 2 j* 0 5 3 j* 0 ]。 %初始化導納矩陣for m=1:5 p=Z(m,1)。 %q=2,3,3,2,3 admittance(p,q)=admittance(p,q)1./(Z(m,3)*Z(m,5))。 %互導納相等 %計算自導納admittance(q,q)=admittance(q,q)+1./(Z(m,3)*Z(m,5)^2)+Z(m,4)。 %是否添加對地導納參數(shù)%計算低壓側(cè)自導納endconductance=real(admittance)。 %電納英文:susceptanceb1=susceptance(1:4,1:4)。%PQ分解法中對應(yīng)p中的電納矩陣//僅用作求逆b2=susceptance(1:3,1:3)。%PQ分解法中對應(yīng)q的電納矩陣b0=susceptance。var=[%發(fā)電機功率 負荷功率 電壓初值 節(jié)點類型 0 * 1 2 0 2j 1 2 0 * 1 2 5 0 3 0 0 1]。 %創(chuàng)建電壓列向量 v=zeros(4,1)。 %創(chuàng)建相角列向量voltage=var(:,3)。ends=s39。q=imag(s)。 %取對地導納矩陣y0 =imag(y0)。node_e=zeros(5,1)。node_e=Z(:,2)。Return攀枝花學院本科畢業(yè)設(shè)計（論文） 附錄B:潮流計算子程序附錄B：潮流計算子程序%,輸入voltage,ang,p,q,g,b0,b1,b2,輸出voltage,angfunction [voltage,ang]=PQ(voltage,ang,p,q,g,b0,b1,b2)%volatge,ang,為迭代計算的初值；iteration=1。 %%收斂容許誤差while(iteration20)v=voltage(1:4,1)。for m=1:4 for n=1:5 p_input(m)=voltage(n)*(g(m,n)*cos(ang(m)ang(n))+b0(m,n)*sin(ang(m)ang(n)))+p_input(m)。endp_correction=zeros(4,1)。ang_correction=zeros(4,1)。ang(1:4,1)=ang(1:4,1)ang_correction。for m=1:3 for n=1:5 q_input(m,1)=voltage(n,1)*(g(m,n)*sin(ang(m)ang(n))b0(m,n)*cos(ang(m)ang(n)))+q_input(m)。endq_correction=zeros(3,1)。%% check for convergence F = [ p_correction。 normF = norm(F, inf)。\nNewton39。s method power flow converged in %d iterations.\n39。 return。v_correction=zeros(3,1)。voltage(1:3,1)=voltage(1:3,1)v_correction。end攀枝花學院本科畢業(yè)設(shè)計（論文） 附錄C:支路功率計算子程序附錄C：支路傳輸功率計算子程序%功能求支路流過的功率%注明：1）pij和pji在數(shù)值上是不同的。function [branch_p,branch_q]=branch_s(v,ang,g,b0,y0,node_f,node_e)for m=1:5 %pij i=node_f, j=node_f %為程序的可讀性，設(shè)置如下3個中間變量。 n2=node_e(m,1)。 %pijbranch_p(1,m)= v(n1,1)^2*g(n1,n2)v(n1,1)*v(n2,1)*(g(n1,n2)*sin(ang_mid)+b0
(n1,n2)*sin(ang_mid))。 %pjibranch_p(2,m)=v(n2,1)^2*g(n2,n1)v(n2,1)*v(n1,1)*(g(n1,n2)*cos(ang_mid)+b0
(n2,n1)*sin(ang_mid))。endend攀枝花學院本科畢業(yè)設(shè)計（論文） 附錄D:支路開斷時首末節(jié)點注入功率增量計算子函數(shù) 附錄D：支路開斷時首末節(jié)點注入功率增量計算子函數(shù)function [p_inject,q_inject]=simulink_dynamic(p,q,branch_p,branch_q,b0,b1,b2,y0,node_f,node_e) p_inject=zeros(2,5)。for m=1:3 n1=node_f(m,1)。 p_inject(:,m)=inv(eye(2)[1,1。b1(n2,n1),b1(n2,n2)]*b1(n1,n2))*branch_p(:,m)。 q_inject(:,m)=inv(eye(2)[c,1。b2(n2,n1),b2(n2,n2)*b0(n1,n2)])*branch_q(:,m)。 q0=q。 n2=node_e(m)。 p(n2,1)=p(n2,1)+p_inject(2,m)。 q(n1,1)=q(n1,1)+q_inject(1,m)。 q_revise(:,m)=q。 q=q0。v=voltage_stability。p=p_stability。% %% 畫穩(wěn)態(tài)時的電壓圖% frame %圖形框架% address=name_cell{1,1}。 %畫圖同時存放圖片%%% 求支路數(shù)據(jù)[branch_p,branch_q]=branch_s(v,ang,g,b0,y0,node_f,node_e)。p_injectq_inject[p_revise,q_revise]=pq_revise(p,q,p_inject,q_inject,node_f,node_e)p_reviseq_revise%% 以下為靜態(tài)安全分析所需的數(shù)據(jù),所需數(shù)據(jù)有節(jié)點電壓幅值,相位。%第一組數(shù)據(jù):電壓相角數(shù)據(jù)analysis_data_v=zeros(5,3)。%% 求功率注入修改后的潮流分布for m=1:3 [analysis_data_v(:,m),analysis_data_ang(:,m)]=pq(v,ang,p_revise(:,m),q_revise
(:,m),g,b0,b1,b2)。% address=name_cell{1,1+m}。 %畫圖同時存放圖片% end% analysis_data_ang %得到的開斷后的節(jié)點電壓幅值%數(shù)據(jù)整合.analysis_branch_p=zeros(6,5)。for m=1:3 [analysis_branch_p(2*m1:2*m,:),analysis_branch_q(2*m1:2*m,:)]=branch_s
(analysis_data_v(:,m),analysis_data_ang(:,m),g,b0,y0,node_f,node_e)。D:\39。wt39。fprintf(fid,39。)。for m=1:3 fprintf(fid,39。,m,node_f(m,1),node_e
(m,1))。%d 節(jié)點: \t39。 if analysis_data_v(n,m)= amp。 analysis_data_v(n,m)= %% 電壓校驗及結(jié)果
輸出 fprintf(fid,39。)。 else fprintf(fid,39。)。 end judge_result(n,m)=result。\n第 %d 條支路:%d_%d開斷時發(fā)電機功率校驗 \npg5\t39。 %% 發(fā)電機有功判斷分析 if analysis_branch_p(2*m,5)1 amp。 analysis_branch_p(2*m,5)8 fprintf(fid,39。)。danger\n39。 end %% 對5號節(jié)點發(fā)電機無功功的校驗 fprintf(fid,39。)。amp。secure\n39。 else fprintf(fid,39。)。pg4\t39。 if analysis_branch_q(2*m,4)3 amp。 analysis_branch_q(2*m,4)3 fprintf(fid,39。)。danger\n39。 endend%% 輸出數(shù)據(jù),以文本形式輸出.fid = fopen(39。, 39。)。支路傳輸有功率功率\n\n39。fprintf(fid,39。)。%f\t39。\n39。)。\n\n支路傳輸無功率功率\n39。fprintf(fid,39。)。%f\t39。\n39。)。\n被開斷支路兩節(jié)點注入有功率\n39。fprintf(fid,39。)。%f\t39。\n39。)。\n\n被開斷支路兩節(jié)點注入無功率\n39。fprintf(fid,39。)。%f\t39。\n39。)。\n\n修改后的節(jié)點注入有功率\n39。fprintf(fid,39。)。%f\t39。\n39。)。\n\n修改后的節(jié)點注入無功率\n39。fprintf(fid,39。)。%f\t39。\n39。)。result=zeros(5,8)。result(:,2)=。result(:,4)=analysis_data_v(:,1)。result(:,8)=analysis_data_v(:,3)。result(:,7)=judge_result(:,2)。 disp(result)。\n\n節(jié)點\t電壓上限\t電壓下限\t1_2線路開斷電壓數(shù)據(jù)\t判斷結(jié)果\t1_3線
路開斷\t結(jié)果\t2_3線路開斷\t判斷結(jié)果\n39。fprintf(fid,39。,result39。fprintf(fid,39。\t0:danger39。%% 穩(wěn)態(tài)值 stability=zeros(5,4)。 stability(:,2)=q_stability(:,1)。 stability(:,4)=ang_stability(:,1)。\n電導\n39。fprintf(fid,[repmat(39。,1,size(g,2)),39。],g39。fprintf(fid,39。)fprintf(fid,[repmat(39。,1,size(b0,2)),39。],b039。fprintf(fid,39。)。%f\t39。\n39。)。\nB39。\n39。fprintf(fid,[repmat(39。,1,size(b2,2)),39。],b239。fprintf(fid,39。)。%f\t39。\n39。)。攀枝花學院本科畢業(yè)設(shè)計（論文） 致謝致 謝 斷斷續(xù)續(xù)的畢業(yè)設(shè)計中，從一開始對題目的膚淺認識而有所畏懼。首先感謝指導老師，戶秀瓊老師的盡心點撥。同時感謝數(shù)學與計算機學院的林宗兵老師，林老師在MATLAB編程上的兩次指點，讓學生對MATLAB編程中關(guān)于子函數(shù)調(diào)用和數(shù)據(jù)類型不再困