【正文】
。2)由于堯化門(mén)2主變跳停，致使3號(hào)機(jī)跳停，金陵分公司供電能力嚴(yán)重不足。由于公司負(fù)荷主要由高低壓電機(jī)組成，低壓電機(jī)一般在電壓小于70%額定電壓時(shí)并持續(xù)幾十毫秒即停運(yùn)。二、事故發(fā)生原因1)南京地區(qū)電網(wǎng)220kv系統(tǒng)為環(huán)網(wǎng)運(yùn)行，龍王山變電站220kV系統(tǒng)發(fā)生故障（接地并短路）使公司兩段110kV系統(tǒng)電壓同時(shí)波動(dòng)，且持續(xù)時(shí)間較長(zhǎng)，約4秒鐘，電壓最低降至40%額定電壓以下。4)烷基苯水廠總降低壓配電室Ⅰ段、Ⅱ段進(jìn)線開(kāi)關(guān)失壓跳閘，母聯(lián)未自投，全廠所有低壓設(shè)備停電。3)烷三車(chē)間的16變高壓Ⅰ段、Ⅱ段進(jìn)線開(kāi)關(guān)未受影響。12變低壓400V2進(jìn)線開(kāi)關(guān)跳閘，母聯(lián)未自投，低壓負(fù)載全部失電。這些裝置一半有電。金陵石化熱電廠電氣一次主接線圖煉油Ⅱ線化肥Ⅰ線煉油Ⅰ線堯金Ⅱ線堯金Ⅰ線煉油Ⅰ線化肥Ⅱ線煉油Ⅱ線高備變主變主變主變高廠變高廠變高廠變廠低變廠低變低備變廠低變廠低變廠低變低備變廠低變輸煤變輸煤變輸煤變除塵變除塵變除塵變除塵變除塵變化水變化水變化水變廠前變輸煤變輸煤變圖1 熱電運(yùn)行部電力系統(tǒng)圖煉油及煤化工區(qū)域6（10）kV變電所供電系統(tǒng)影響情況1)煉油及煤化工區(qū)域53個(gè)6(10)kV分變電所中45個(gè)備自投成功，恢復(fù)供電。5:13：汽機(jī)專(zhuān)業(yè)打跳Ⅰ、Ⅱ機(jī)，Ⅰ、Ⅱ機(jī)主汽門(mén)關(guān)閉。Ⅴ、Ⅵ爐UPS轉(zhuǎn)電池供電，Ⅲ發(fā)電機(jī)公用系統(tǒng)UPS、Ⅲ發(fā)電機(jī)UPS所有指示全無(wú)；至主控檢查時(shí)，發(fā)現(xiàn)Ⅲ機(jī)保護(hù)裝置黑屏。5:09：汽機(jī)專(zhuān)業(yè)就地檢查Ⅰ、Ⅱ機(jī)在運(yùn)行狀態(tài)，4給水泵運(yùn)行，5給水泵跳閘，出口門(mén)未關(guān)閉。隨后南京電網(wǎng)故障切除，但同時(shí)地區(qū)電網(wǎng)堯化門(mén)變電站2主變跳停。5：06分左右，南京地區(qū)電網(wǎng)龍王山變電站220kV系統(tǒng)發(fā)生故障。盡管如此，但是在實(shí)驗(yàn)中依然容易把for循環(huán)跟C語(yǔ)言中的for語(yǔ)句弄混，最后經(jīng)過(guò)不懈努力下，終于弄明白了兩者之間的差別，使我能更好的運(yùn)用這些指令語(yǔ)句。i=i+1。 elsey=1。 elseif n==139。 elseif n==039。if n39。請(qǐng)輸入一個(gè)正數(shù)n=39。!的函數(shù)程序，取n分別為80、10驗(yàn)證其正確性(輸入n為負(fù)數(shù)時(shí)輸出出錯(cuò)信息)。 [65] 39。39。39。[100][40]39。39。39。39。39。[85][76][95]39。39。39。end A A =39。E39。D39。C39。B39。A39。A+39。A(2,:)={85,76,95,100,40,65}。f39。e39。d39。c39。b39。a39。end A=cell(3,6)。c(k)={69+k}。clear,clc for k=1:10a(k)={89+k}。要求：（1）用switch語(yǔ)句實(shí)現(xiàn)。 end，要求輸出成績(jī)等級(jí)A+、A、B、C、D、E。if m0x1=(b+sqrt(m))/(2*a)x2=(bsqrt(m))/(2*a)elseif m==0x=(b)/(2*a)else39。為恒不等式39。m=b^24*a*c。請(qǐng)輸入一個(gè)數(shù)c=39。)。b=input(39。請(qǐng)輸入一個(gè)數(shù)a=39。并輸入幾組典型值加以檢驗(yàn)。二、實(shí)驗(yàn)設(shè)備及條件計(jì)算機(jī)一臺(tái)（）。實(shí)驗(yàn)二 MATLAB程序設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康?。x2y2=1繪制平面曲線，并分析參數(shù)a對(duì)其形狀的影響。2在－54，繪制函數(shù)的曲面圖和等高線z=cosxcosyex2+y24其中x的21個(gè)值均勻分布在[5,5]范圍，y的31個(gè)值均勻分布在[0,10]，要求使用subplot（2，1，1）和subplot（2，1，2）將產(chǎn)生的曲面圖和登高圖畫(huà)在同一個(gè)窗口上。1In(x+1+x2),x0239。e y=237。x163。3，已知236。五、實(shí)驗(yàn)內(nèi)容及步驟：1，設(shè)y=[+3sinx]cosx，在x＝0～2π區(qū)間取101點(diǎn)，繪制函數(shù)曲線。四、預(yù)習(xí)要求：1．復(fù)習(xí)二維與三維圖形的繪圖函數(shù)。為了保證圖象的顯示效果，一般還應(yīng)使用colormap函數(shù)設(shè)置圖象色圖。6，圖像處理1）imread和imwrite函數(shù) 這兩個(gè)函數(shù)分別用于將圖象文件讀入matlab工作空間，以及將圖象數(shù)據(jù)和色圖數(shù)據(jù)一起寫(xiě)入一定格式的圖象文件。一般情況下，x，y，z是維數(shù)相同的矩陣。當(dāng)x，y是同維矩陣時(shí)，則以x，y對(duì)應(yīng)列元素為橫、縱坐標(biāo)分別繪制曲線，曲線條數(shù)等于矩陣的列數(shù)。二、實(shí)驗(yàn)原理：1，二維數(shù)據(jù)曲線圖（1）繪制單根二維曲線 plot(x,y)。=b+x(xc)xa=b==仿真結(jié)果LMS自適應(yīng)濾波器： 仿真圖：遇到的問(wèn)題問(wèn)題的解決辦法：將Spectrum Scope中的buffer input打勾就行了。=(x+x)xamp。仿真結(jié)果，左邊為原始信號(hào)，右邊為積分后的信號(hào)。關(guān)于學(xué)習(xí)的心得就寫(xiě)到這里了，接下來(lái)是一些自己做的仿真。Simulink中的模型都是分級(jí)的，因此可以通過(guò)自上而下或者自下而上的方式建立模型。（4）檢查塊之間的鏈接，是否每一個(gè)塊的輸出端口與它所連接的模塊輸入端口有相同的信號(hào)寬度。（2）模型系統(tǒng)中的各個(gè)層級(jí)將被平展開(kāi)來(lái)，每一個(gè)子系統(tǒng)將被相應(yīng)的模塊所代替。而且里面的SFunction也是比較好的一個(gè)設(shè)計(jì)，用SFunction可以自定義一些組件，使自己的仿真看上去更加清晰與精簡(jiǎn)。不是很麻煩的。Tables（函數(shù)和平臺(tái)模塊）Math（數(shù)學(xué)模塊）Nonlinear（非線性模塊）Signalsamp。神經(jīng)網(wǎng)絡(luò)濾波器等等，鑒于代碼量太大，就不再敖述。Red ComponentGreen ComponentBlue Component將源圖像的RGB分量轉(zhuǎn)換到Y(jié)CbCr上后又轉(zhuǎn)化為RGB得到的三個(gè)分量差異的圖。而且看了一下matlab自帶的rgb2ycbcr和ycbcr2rgb并且用了這兩個(gè)函數(shù)測(cè)試后圖像是一樣的。)。title(39。subplot(1,3,3)。綠色分量差異39。imshow(GD)。)。title(39。subplot(1,3,1)。BD=B(:,:,3)blue(:,:,3)。RD=R(:,:,1)red(:,:,1)。RGB轉(zhuǎn)換為YCrCb后又轉(zhuǎn)換為RGB的圖像39。imshow(Comp2)。Comp2(:,:,3)=blue(:,:,3)。Comp2(:,:,1)=red(:,:,1)。green=*(Cb128)*(Cr128)。Cr=***B+128。Y=*R+*G+*B。position image39。subplot(1,2,2)。Comp(:,:,3)=B(:,:,3)。%合成% Comp(:,:,1)=R(:,:,1)。Blue Component39。figure(2)。B(:,:,2)=zeros(r,c)。%計(jì)算藍(lán)色分量并顯示分解圖% B(:,:,3)=Source(:,:,3)。Green Component39。subplot(1,3,2)。G=uint8(G)。G(:,:,1)=zeros(r,c)。)。imshow(R)title(39。figure(2)。R=uint8(R)。R(:,:,2)=zeros(r,c)。%顯示圖像 [r c d]=size(Source)。original image39。subplot(1,2,1)。)。Source=imread(39。下面是寫(xiě)的程序代碼以及測(cè)試得到的結(jié)果圖像。JPEG解壓時(shí)需要將YCbCr空間的圖像又轉(zhuǎn)化回來(lái)到RGB上。原理：BMP圖像壓縮為JPEG的第一步是將RGB色彩空間通過(guò)這個(gè)公式映射到Y(jié)CbCr空間上。end y=a。lx=length(x)。i=i1。while a0 y(i)=mod(a,2)。a=x。%顯示標(biāo)題 end附帶的兩個(gè)函數(shù)：function y=numdec2bin(x,n)。層后的結(jié)果39。去除第39。%顯示n個(gè)圖像，每層是相應(yīng)的剔除該層的結(jié)果imshow(y)。%將剩下的7個(gè)元素轉(zhuǎn)換為10進(jìn)制數(shù)放入y矩陣中end endy=uint8(y)。%先將每個(gè)像素轉(zhuǎn)換為8位2進(jìn)制序列B(n)=[]。%圖像長(zhǎng)寬之積 for n=1:8。%長(zhǎng) N=size(I,2)。%讀入源圖像I=double(I)。39。要做的就是把每bit的信息提取出來(lái)。,18)。text((x(m)+x(m))/2,(y(m)+y(m))/2,f,39。f=num2str(A(m,m))。legend(e)。)。%連線 end endtitle(39。%將權(quán)值顯示在兩點(diǎn)連線中間 end hold on。Fontsize39。039。for i=1:DEF a=find(A(i,:)0)%將A矩陣每行大于0的數(shù)的在該行的地址找出來(lái)放在a中for j=1:length(a)c=num2str(A(i,j))。r+39。y=100*rand(1,DEF)。A=floor(A)。%設(shè)定一個(gè)東東 方便改變隨機(jī)點(diǎn)的個(gè)數(shù) A=rand(DEF,DEF)。然后畫(huà)出這個(gè)矩陣就行了。以下是我自己寫(xiě)的一些程序的代碼1：網(wǎng)絡(luò)隨機(jī)拓?fù)鋱D目的是要生成一個(gè)度數(shù)隨機(jī)、權(quán)值隨機(jī)的拓?fù)鋱D。P*A = L*U.。1 1] A = [L,U]=lu(A)L =0U =0 我們老師當(dāng)時(shí)的解釋是這是由于我們所用matlab是破解版的緣故，破解版的計(jì)算精度 沒(méi)正版的高所以導(dǎo)致運(yùn)算結(jié)果的錯(cuò)誤，后來(lái)我看了一些LU這個(gè)函數(shù)的help文檔 psychologically lower triangular matrix( product of lower triangular and permutation matrices)in L。1/epusino 1]。01/epusino]。因?yàn)?/e為一個(gè)極大數(shù) 則1—1/e可以看成—1/e；設(shè)e=0.***000000001%10的負(fù)24次方從而L*U=則最后一位1沒(méi)有了 變成了0其中L=U=Matlab 程序?yàn)閑pusino=0.***000000001。235。235。 經(jīng)過(guò)三角分解為L(zhǎng)=101235。 U=234。234。ee1233。1249。10233。其實(shí)這個(gè)事件最好的解決辦法是在3DMAX中，不過(guò)因?yàn)闆](méi)能找到MATLAB與3DMAX是怎么接口的，所以作罷了。然后做出來(lái)的。最后問(wèn)題解決的方式是查著了一些圖書(shū)館的學(xué)報(bào)。當(dāng)然我也不是很懂，只是理解了部分?，F(xiàn)有的許多方法使用了基于RBF的插值技術(shù)，將一般人臉網(wǎng)格變化到特定人臉的形狀。但是，這種參數(shù)模型僅僅局限于那些構(gòu)造參數(shù)已知的情況，并且對(duì)特定人臉參數(shù)的調(diào)整非常困難。確定該方法來(lái)實(shí)現(xiàn)三維人臉變形是可行的，但要構(gòu)造好的基函數(shù)，以及解這些高維方程的解仍是難點(diǎn)。第二種也就是基本的老辦法RBF。分為兩個(gè)過(guò)程內(nèi)部插值和邊緣插值。然后進(jìn)行插值。但是速度挺好。主要有兩種，第一種是face—ls算法。需要通過(guò)插值使其變得光滑，當(dāng)然問(wèn)題的復(fù)雜之處在于，插值會(huì)改變?cè)镜木W(wǎng)格結(jié)構(gòu)。問(wèn)題可以概要為：人臉網(wǎng)格插值。最后的處理辦法是把這些數(shù)據(jù)c os，sin這些東西變換一下后，使其大致規(guī)則，然后再用polar畫(huà)出極坐標(biāo)下的圖形。然后就想到把極坐標(biāo)的數(shù)據(jù)轉(zhuǎn)換的直角坐標(biāo)系下，pol2cart可以實(shí)現(xiàn)這個(gè)想法，但是隨后而來(lái)，也就是最后導(dǎo)致整個(gè)問(wèn)題失敗的關(guān)鍵也在這里。首先是一個(gè)物理實(shí)驗(yàn)的問(wèn)題；在做一個(gè)關(guān)于光的偏振的實(shí)驗(yàn)的時(shí)候，最后的結(jié)果要在一個(gè)極化坐標(biāo)下顯示出來(lái)；因?yàn)閿?shù)據(jù)是離散的，所以顯示出來(lái)的圖像是一個(gè)折來(lái)折去的一個(gè)東東； 然后很自然的想法是對(duì)這個(gè)曲線進(jìn)行插值處理。但是 手動(dòng)用鉛筆去畫(huà)確實(shí)很麻煩，所以用matlab確實(shí)可以省去很大的麻煩。但是matlab畢竟只是一個(gè)工具，原理和一些基本的編程素質(zhì)還是必須有的，否則matlab最多也只能是驗(yàn)證一些別人的東西而已，根本幫不上什么忙的。其實(shí)最初開(kāi)始學(xué)習(xí)matlab的時(shí)候感覺(jué)這個(gè)東西和C沒(méi)什么兩樣，但是后來(lái)具體到一些 東西，比如信號(hào)處理和數(shù)學(xué)建模上以后才感覺(jué)到使用matlab編寫(xiě)程序去驗(yàn)證結(jié)果比C要節(jié) 省很多時(shí)間，而且matlab寫(xiě)東西基本都是按照自己的思路平鋪直敘很少去考慮什么函數(shù)的 嵌套調(diào)用或者指針等等很頭疼的東西。Matlab的確不愧成為是草稿紙上的語(yǔ)言。：基于Matlab/Simulink[M].北京：機(jī)械工業(yè)出版社，2005.[21]薛定語(yǔ)，[M]．北京：機(jī)械工業(yè)出版社，2005.[19]張森，[M].北京：機(jī)械工業(yè)出版社，2004.[15] 張葛祥，[M]．北京：國(guó)防工業(yè)出版社，2004.[16]在此我向胡老師表示我最誠(chéng)摯的感謝！其次，在大學(xué)的四年里，學(xué)校的各位老師在學(xué)習(xí)上和生活上對(duì)我的關(guān)心、支持和厚愛(ài)，都令我銘記于心，在此也表示我由衷的感謝！各位老師嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度，踏踏實(shí)實(shí)做人的品質(zhì)，和藹可親的作風(fēng)以及無(wú)私奉獻(xiàn)的精神給我留下了深刻的印象，使我終生受益。致 謝 本設(shè)計(jì)在老師的精心指導(dǎo)下圓滿完成，從開(kāi)題報(bào)告到資料查詢(xún)?cè)俚秸撐牡淖珜?xiě)，XX老師都耐心地教導(dǎo)，并為我點(diǎn)出其中的錯(cuò)誤。通過(guò)這次論文的撰寫(xiě)，我對(duì)電力系統(tǒng)故障以及Matlab軟件有了更深的認(rèn)識(shí)。高效的數(shù)據(jù)仿真分析，特別是信號(hào)處理和直觀的圖形顯示功能，且Matlab/simulink環(huán)境下的SimPowerSystems模型庫(kù)及Simulink強(qiáng)大的二次開(kāi)發(fā)功能和豐富的工具箱，能快速而準(zhǔn)確的對(duì)電路及更復(fù)雜的電力系統(tǒng)進(jìn)行仿真、計(jì)算。總結(jié) 本文以三相短路為例，介紹了電力系統(tǒng)故障和Matlab/Simulink的基本特點(diǎn)，探索了電力系統(tǒng)故障中最常見(jiàn)的短路計(jì)算一些常用的計(jì)算方法，和Matlab應(yīng)用的基本方法和步驟，在Matlab軟件中電力系統(tǒng)仿真如何應(yīng)用SimPowerSystems模塊庫(kù)構(gòu)建電力系統(tǒng)故障的仿真模型并對(duì)其仿真結(jié)果進(jìn)行分析，得出一下結(jié)果：應(yīng)用Matlab/simulink進(jìn)行仿真分析的結(jié)果和理論計(jì)算的結(jié)果相差不大，Matlab仿真工具是一種很實(shí)用的工具。在用Matlab進(jìn)行電力系統(tǒng)故障仿真得出的短路點(diǎn)電流的標(biāo)幺值為9，短路點(diǎn)電壓為0；，電壓因?yàn)槎搪范鵀榱?。電力系統(tǒng)在線路上發(fā)生三相短路故障時(shí)，由于短路故障發(fā)生器處于斷開(kāi)狀態(tài)所以短路點(diǎn)電流是由零開(kāi)始迅