【正文】 n2=node_e(m)。 p(n1,1)=p(n1,1)+p_inject(1,m)。 p(n2,1)=p(n2,1)+p_inject(2,m)。 p_revise(:,m)=p。 q(n1,1)=q(n1,1)+q_inject(1,m)。 q(n2,1)=q(n2,1)+q_inject(2,m)。 q_revise(:,m)=q。 p=p0。 q=q0。endend攀枝花學(xué)院本科畢業(yè)設(shè)計（論文） 附錄F: 靜態(tài)安全分析主程序 附錄F: 靜態(tài)安全分析主程序clearclc%% 求穩(wěn)態(tài)解,并用變量存儲,作為基礎(chǔ)數(shù)據(jù)庫[p_stability,q_stability,voltage_stability,ang_stability,g,b0,b1,b2,y0,node_f,node_e,name_cell
]=basic_data()。v=voltage_stability。ang=ang_stability。p=p_stability。q=q_stability。% %% 畫穩(wěn)態(tài)時的電壓圖% frame %圖形框架% address=name_cell{1,1}。 %圖形存放地址% draw_voltage(v,address)。 %畫圖同時存放圖片%%% 求支路數(shù)據(jù)[branch_p,branch_q]=branch_s(v,ang,g,b0,y0,node_f,node_e)。[p_inject,q_inject]=simulink_dynamic(p,q,branch_p,branch_q,b0,b1,b2,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é)點電壓幅值,相位。支路傳輸功率,發(fā)
電機(jī)的發(fā)出功率。%第一組數(shù)據(jù):電壓相角數(shù)據(jù)analysis_data_v=zeros(5,3)。%將數(shù)據(jù)整成兩個5*3的數(shù)據(jù)矩陣其存取速度是否會快于讀
取一個5*6的矩陣analysis_data_ang=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)。endanalysis_data_v %得到的開斷后的節(jié)點電壓幅值% %% 畫3次線路開斷后的電壓分布% for m=1:3% frame2(v)。% address=name_cell{1,1+m}。 %圖形存放地址% draw_voltage_break(analysis_data_v(:,m),address)。 %畫圖同時存放圖片% end% analysis_data_ang %得到的開斷后的節(jié)點電壓幅值%數(shù)據(jù)整合.analysis_branch_p=zeros(6,5)。analysis_branch_q=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)。endanalysis_branch_panalysis_branch_q%% 靜態(tài)安全分析,同時比較結(jié)果,存入文檔fid=fopen(39。D:\39。,39。wt39。)。fprintf(fid,39。\n靜態(tài)安全分析結(jié)果\n39。)。judge_result=zeros(5,3)。for m=1:3 fprintf(fid,39。\n第 %d 條支路:%d_%d開斷時節(jié)點電壓校驗 \n39。,m,node_f(m,1),node_e
(m,1))。 for n=1:5 fprintf(fid,39。%d 節(jié)點: \t39。,n)。 if analysis_data_v(n,m)= amp。amp。 analysis_data_v(n,m)= %% 電壓校驗及結(jié)果
輸出 fprintf(fid,39。secure\n39。)。 result=1。 else fprintf(fid,39。danger\n39。)。 result=0。 end judge_result(n,m)=result。 endend%% 對5號節(jié)點發(fā)電機(jī)功率的校驗for m=1:3 fprintf(fid,39。\n第 %d 條支路:%d_%d開斷時發(fā)電機(jī)功率校驗 \npg5\t39。,m,node_f
(m,1),node_e(m,1))。 %% 發(fā)電機(jī)有功判斷分析 if analysis_branch_p(2*m,5)1 amp。amp。 analysis_branch_p(2*m,5)8 fprintf(fid,39。secure\n39。)。 else fprintf(fid,39。danger\n39。)。 end %% 對5號節(jié)點發(fā)電機(jī)無功功的校驗 fprintf(fid,39。qg5\t39。)。 if analysis_branch_q(2*m,5) amp。amp。 analysis_branch_q(2*m,5)5 fprintf(fid,39。secure\n39。)。 else fprintf(fid,39。danger\n39。)。 end %% 對4號發(fā)電機(jī)節(jié)點無功校驗 fprintf(fid,39。pg4\t39。)。 if analysis_branch_q(2*m,4)3 amp。amp。 analysis_branch_q(2*m,4)3 fprintf(fid,39。secure\n39。)。 else fprintf(fid,39。danger\n39。)。 endend%% 輸出數(shù)據(jù),以文本形式輸出.fid = fopen(39。d:\39。, 39。wt39。)。 %t模式按照文本而非二進(jìn)制模式讀寫 fprintf(fid,39。支路傳輸有功率功率\n\n39。)。fprintf(fid,39。支路1_2\t支路1_3\t支路2_3\t支路2_4\t支路3_5\n39。)。fprintf(fid,[repmat(39。%f\t39。, 1, size(branch_p,2)), 39。\n39。], branch_p39。)。%轉(zhuǎn)置并適時(寫完一列后)
添加換行fprintf(fid,39。\n\n支路傳輸無功率功率\n39。)。fprintf(fid,39。支路1_2\t支路1_3\t支路2_3\t支路2_4\t支路3_5\n39。)。fprintf(fid,[repmat(39。%f\t39。, 1, size(branch_q,2)), 39。\n39。], branch_q39。)。%轉(zhuǎn)置并適時(寫完一列后)
添加換行fprintf(fid,39。\n被開斷支路兩節(jié)點注入有功率\n39。)。fprintf(fid,39。支路1_2\t支路1_3\t支路2_3\t支路2_4\t支路3_5\n39。)。fprintf(fid,[repmat(39。%f\t39。, 1, size(p_inject,2)), 39。\n39。], p_inject39。)。%轉(zhuǎn)置并適時(寫完一列后)添
加換行fprintf(fid,39。\n\n被開斷支路兩節(jié)點注入無功率\n39。)。fprintf(fid,39。\n支路1_2\t支路1_3\t支路2_3\t支路2_4\t支路3_5\n39。)。fprintf(fid,[repmat(39。%f\t39。, 1, size(q_inject,2)), 39。\n39。], q_inject39。)。 fprintf(fid,39。\n\n修改后的節(jié)點注入有功率\n39。)。fprintf(fid,39。\n支路1_2開斷\t支路1_3開斷\t支路2_3開斷\n39。)。fprintf(fid,[repmat(39。%f\t39。,1,size(p_revise,2)),39。\n39。],p_revise39。)。fprintf(fid,39。\n\n修改后的節(jié)點注入無功率\n39。)。fprintf(fid,39。\n支路1_2開斷\t支路1_3開斷\t支路2_3開斷\n39。)。fprintf(fid,[repmat(39。%f\t39。,1,size(q_revise,2)),39。\n39。],q_revise39。)。disp(judge_result)。result=zeros(5,8)。result(:,1)=1:5。result(:,2)=。result(:,3)=。result(:,4)=analysis_data_v(:,1)。result(:,6)=analysis_data_v(:,2)。result(:,8)=analysis_data_v(:,3)。result(:,5)=judge_result(:,1)。result(:,7)=judge_result(:,2)。result(:,9)=judge_result(:,3)。 disp(result)。 fprintf(fid,39。\n\n節(jié)點\t電壓上限\t電壓下限\t1_2線路開斷電壓數(shù)據(jù)\t判斷結(jié)果\t1_3線
路開斷\t結(jié)果\t2_3線路開斷\t判斷結(jié)果\n39。)。fprintf(fid,39。%d\t%f\t%f \t%f\t%d \t%f\t%d \t%f\t%d\n39。,result39。)。fprintf(fid,39。\njudge_resutl:1:secure。\t0:danger39。)。%% 穩(wěn)態(tài)值 stability=zeros(5,4)。 stability(:,1)=p_stability(:,1)。 stability(:,2)=q_stability(:,1)。 stability(:,3)=voltage_stability(:,1)。 stability(:,4)=ang_stability(:,1)。fprintf(fid,39。\n電導(dǎo)\n39。)。fprintf(fid,[repmat(39。%f\t39。,1,size(g,2)),39。\n39。],g39。)。fprintf(fid,39。\n電納\n39。)fprintf(fid,[repmat(39。%f\t39。,1,size(b0,2)),39。\n39。],b039。)。fprintf(fid,39。\nB\n39。)。fprintf(fid,[repmat(39。%f\t39。,1,size(b1,2)),39。\n39。],b139。)。fprintf(fid,39。\nB39。39。\n39。)。fprintf(fid,[repmat(39。%f\t39。,1,size(b2,2)),39。\n39。],b239。)。fprintf(fid,39。穩(wěn)態(tài)有功\t穩(wěn)態(tài)無功\t穩(wěn)態(tài)電壓\t電壓相角\n39。)。fprintf(fid,[repmat(39。%f\t39。,1,size(stability,2)),39。\n39。],stability39。)。fclose(fid)。攀枝花學(xué)院本科畢業(yè)設(shè)計（論文） 致謝致 謝 斷斷續(xù)續(xù)的畢業(yè)設(shè)計中，從一開始對題目的膚淺認(rèn)識而有所畏懼。到現(xiàn)在完成畢業(yè)設(shè)計后回頭看，這一階段的努力，應(yīng)該是在老師點撥下，獨立完成一個項目的過程，是一個把書本知識轉(zhuǎn)化為應(yīng)用技術(shù)的過程。首先感謝指導(dǎo)老師，戶秀瓊老師的盡心點撥。從文開始到現(xiàn)在論文的完成，戶老師在這一過程中在專業(yè)視野，尤其是使學(xué)生對電力系統(tǒng)穩(wěn)定這方面有了系統(tǒng)的了解。同時感謝數(shù)學(xué)與計算機(jī)學(xué)院的林宗兵老師，林老師在MATLAB編程上的兩次指點，讓學(xué)生對MATLAB編程中關(guān)于子函數(shù)調(diào)用和數(shù)據(jù)類型不再困惑。