【正文】 。從文開始到現(xiàn)在論文的完成，戶老師在這一過程中在專業(yè)視野，尤其是使學(xué)生對電力系統(tǒng)穩(wěn)定這方面有了系統(tǒng)的了解。到現(xiàn)在完成畢業(yè)設(shè)計后回頭看，這一階段的努力，應(yīng)該是在老師點撥下，獨立完成一個項目的過程，是一個把書本知識轉(zhuǎn)化為應(yīng)用技術(shù)的過程。fclose(fid)。],stability39。,1,size(stability,2)),39。fprintf(fid,[repmat(39。穩(wěn)態(tài)有功\t穩(wěn)態(tài)無功\t穩(wěn)態(tài)電壓\t電壓相角\n39。)。\n39。%f\t39。)。39。fprintf(fid,39。],b139。,1,size(b1,2)),39。fprintf(fid,[repmat(39。\nB\n39。)。\n39。%f\t39。\n電納\n39。)。\n39。%f\t39。)。fprintf(fid,39。 stability(:,3)=voltage_stability(:,1)。 stability(:,1)=p_stability(:,1)。)。\njudge_resutl:1:secure。)。%d\t%f\t%f \t%f\t%d \t%f\t%d \t%f\t%d\n39。)。 fprintf(fid,39。result(:,9)=judge_result(:,3)。result(:,5)=judge_result(:,1)。result(:,6)=analysis_data_v(:,2)。result(:,3)=。result(:,1)=1:5。disp(judge_result)。],q_revise39。,1,size(q_revise,2)),39。fprintf(fid,[repmat(39。\n支路1_2開斷\t支路1_3開斷\t支路2_3開斷\n39。)。fprintf(fid,39。],p_revise39。,1,size(p_revise,2)),39。fprintf(fid,[repmat(39。\n支路1_2開斷\t支路1_3開斷\t支路2_3開斷\n39。)。 fprintf(fid,39。], q_inject39。, 1, size(q_inject,2)), 39。fprintf(fid,[repmat(39。\n支路1_2\t支路1_3\t支路2_3\t支路2_4\t支路3_5\n39。)。%轉(zhuǎn)置并適時(寫完一列后)添
加換行fprintf(fid,39。], p_inject39。, 1, size(p_inject,2)), 39。fprintf(fid,[repmat(39。支路1_2\t支路1_3\t支路2_3\t支路2_4\t支路3_5\n39。)。%轉(zhuǎn)置并適時(寫完一列后)
添加換行fprintf(fid,39。], branch_q39。, 1, size(branch_q,2)), 39。fprintf(fid,[repmat(39。支路1_2\t支路1_3\t支路2_3\t支路2_4\t支路3_5\n39。)。%轉(zhuǎn)置并適時(寫完一列后)
添加換行fprintf(fid,39。], branch_p39。, 1, size(branch_p,2)), 39。fprintf(fid,[repmat(39。支路1_2\t支路1_3\t支路2_3\t支路2_4\t支路3_5\n39。)。 %t模式按照文本而非二進制模式讀寫 fprintf(fid,39。wt39。d:\39。)。 else fprintf(fid,39。secure\n39。amp。)。 end %% 對4號發(fā)電機節(jié)點無功校驗 fprintf(fid,39。danger\n39。)。 analysis_branch_q(2*m,5)5 fprintf(fid,39。 if analysis_branch_q(2*m,5) amp。qg5\t39。)。 else fprintf(fid,39。secure\n39。amp。,m,node_f
(m,1),node_e(m,1))。 endend%% 對5號節(jié)點發(fā)電機功率的校驗for m=1:3 fprintf(fid,39。 result=0。danger\n39。 result=1。secure\n39。amp。,n)。 for n=1:5 fprintf(fid,39。\n第 %d 條支路:%d_%d開斷時節(jié)點電壓校驗 \n39。judge_result=zeros(5,3)。\n靜態(tài)安全分析結(jié)果\n39。)。,39。endanalysis_branch_panalysis_branch_q%% 靜態(tài)安全分析,同時比較結(jié)果,存入文檔fid=fopen(39。analysis_branch_q=zeros(6,5)。 %圖形存放地址% draw_voltage_break(analysis_data_v(:,m),address)。endanalysis_data_v %得到的開斷后的節(jié)點電壓幅值% %% 畫3次線路開斷后的電壓分布% for m=1:3% frame2(v)。%將數(shù)據(jù)整成兩個5*3的數(shù)據(jù)矩陣其存取速度是否會快于讀
取一個5*6的矩陣analysis_data_ang=zeros(5,3)。支路傳輸功率,發(fā)
電機的發(fā)出功率。[p_inject,q_inject]=simulink_dynamic(p,q,branch_p,branch_q,b0,b1,b2,y0,node_f,node_e)。 %圖形存放地址% draw_voltage(v,address)。q=q_stability。ang=ang_stability。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()。 p=p0。 q(n2,1)=q(n2,1)+q_inject(2,m)。 p_revise(:,m)=p。 p(n1,1)=p(n1,1)+p_inject(1,m)。for m=1:3 n1=node_f(m)。 endend攀枝花學(xué)院本科畢業(yè)設(shè)計（論文） 附錄E:修改節(jié)點注入功率子程序 附錄E：修改節(jié)點注入功率子程序%輸入穩(wěn)態(tài)時,節(jié)點注入功率,%輸出斷線后,節(jié)點注入功率function[p_revise,q_revise]=pq_revise(p,q,p_inject,q_inject,node_f,node_e) p0=p。1,c]*[b2(n1,n1),b2(n1,n2)。 c=1+2*y0(n1,1)/b0(n1,n2)。1,1]*[b1(n1,n2),b1(n1,n2)。 n2=node_e(m,1)。q_inject=zeros(2,5)。 %qjibranch_q(2,m)=v(n2,1)^2*(y0(n2,1)+b0(n2,n1))v(n2,1)*v(n1,1)*(g(n2,n1)*sin(
ang_mid)b0(n2,n1)*cos(ang_mid))。%此處由于網(wǎng)絡(luò)拓?fù)鋽?shù)據(jù)的特殊，沒有電納，故此處公式為簡化后的公式 %qijbranch_q(1,m)=v(n1,1)^2*(y0(n1,1)+b0(n1,n2))v(n1,1)*v(n2,1)*(g(n1,n2)*sin
(ang_mid)b0(n1,n2)*cos(ang_mid))。 ang_mid=ang(n1,1)ang(n2,1)。為方便模擬短線注入功率的求解將pij和qij設(shè)置為2*5的格式 n1=node_f(m,1)。2）默認(rèn)i為z矩陣中的第一列元素，為低壓側(cè)數(shù)據(jù)，即默認(rèn)branch_p中第一列總是存放：從低壓側(cè)到高壓側(cè)的支路功率%以上的默認(rèn)是對數(shù)據(jù)的結(jié)構(gòu)體化，對數(shù)據(jù)的結(jié)構(gòu)體一樣固定其格式，便于對程序的優(yōu)化，從中受益已不止一次。iteration=iteration+1。v_correction=b2*q_correction./v2。%跳出對電壓相角的計算 endv2=v(1:3,1)。, iteration)。39。 if normF tol fprintf(39。 q_correction ]。q_correction=q(1:3,1)q_input。 end q_input(m,1)=voltage(m,1)*q_input(m,1)。q_input=zeros(3,1)。ang_correction=b1*p_correction./v./v。p_correction=p(1:4,1)p_input。 end p_input(m)=p_input(m)*voltage(m)。 p_input=zeros(4,1)。tol=1e6。[voltage,ang]=PQ(voltage,ang,p,q,g,b0,b1,b2)。node_f=Z(:,1)。 %對地導(dǎo)納的共軛矩陣node_f=zeros(5,1)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%計算支路潮流所需的參數(shù)y0 =Z(:,4)。p=real(s)。 %對電壓列向量初始化for m=1:5 s(m)=var(m,1)var(m,2)。ang =zeros(5,1)。%/1平衡節(jié)點，2PQ節(jié)點，3pv節(jié)點voltage=zeros(5,1)。g=conductance。b2=inv(b2)。 %4*4b1=inv(b1)。 %電導(dǎo)英文:conductancesusceptance=imag(admittance)。 %是否添加對地導(dǎo)納參數(shù)%計算高壓側(cè)自導(dǎo)納admittance(p,p)=admittance(p,p)+1./Z(m,3)+Z(m,4)。 admittance(q,p)=admittance(p,q)。 %p=1,1,2,4,5 q=Z(m,2)。admittance=zeros(5)。 攀枝花學(xué)院本科畢業(yè)設(shè)計（論文） 參考文獻參 考 文 獻 ［1］KUNDUR. P． Power System Stability and Control [M]. 北京：中國電力出版社，2002.［2］吳際舜．電力系統(tǒng)靜態(tài)安全分析［M］．上海交通大學(xué)出版社，1985.［3］王錫凡, 方萬良, 杜正春．現(xiàn)代電力系統(tǒng)分析［M］．科學(xué)出版社，2003.［4］吳天鵬，謝小介，彭斌. MATLAB電力系統(tǒng)設(shè)計與分析［M］. 北京:國防工業(yè)出版社,2004.［5］王正林，劉明. 精通MATLAB7［M］.電子工業(yè)出版社，2006.［6］何仰贊．電力系統(tǒng)分析［M］．武漢:武漢華中科技大學(xué)出版社，1985.［7］于永源．電力系統(tǒng)分析［M］．北京：中國電力出版社，2007.［8］陳珩．電力系統(tǒng)穩(wěn)態(tài)分析(第三版)［M］. 北京：中國電力出版社，2007.［9］李光琦. 電力系統(tǒng)暫態(tài)分析［M］．北京:中國電力出版社，1995.［10］李文沅. 電力系統(tǒng)風(fēng)險評估模型方法和應(yīng)用［M］．北京:科學(xué)出版社，2006.［11］李文沅. 電力系統(tǒng)經(jīng)濟安全運行模型與方法［M］．重慶:重慶大學(xué)出版社，1989.［12］周榮光．電力系統(tǒng)故障分析［M］．北京:北京清華大學(xué)出版社，1988.［13］中華人民共和國國家經(jīng)濟貿(mào)易委員會. 電力系統(tǒng)安全穩(wěn)定導(dǎo)則DL7552001.［14］許高達．電力系統(tǒng)之計算機應(yīng)用［M］．新興圖書公司，1979.［15］朱仁峰. 精通MATLAB7［M］.清華大學(xué)出版社，2006.［16］Math Works. 產(chǎn)品和服務(wù)［EB/OL］. ，2014513.［17］劉瀏. 百度百科［EB/OL］．:// ,2014513.43攀枝花學(xué)院本科畢業(yè)設(shè)計（論文） 附錄A:基礎(chǔ)數(shù)據(jù)生成子程序附錄A:基礎(chǔ)數(shù)據(jù)生成子程序程序_基礎(chǔ)數(shù)據(jù)生成子程序function[p,q,voltage,an