freepeople性欧美熟妇, 色戒完整版无删减158分钟hd, 无码精品国产vα在线观看DVD, 丰满少妇伦精品无码专区在线观看,艾栗栗与纹身男宾馆3p50分钟,国产AV片在线观看,黑人与美女高潮,18岁女RAPPERDISSSUBS,国产手机在机看影片

正文內(nèi)容

畢業(yè)論文-高壓輸電網(wǎng)潮流的計算機算法程序設(shè)計(參考版)

2025-01-15 19:47本頁面
  

【正文】 在主程序中打開上面的 m文件,運行結(jié)果如下: 節(jié)點導納矩陣 Y +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() 第 1 次迭代結(jié)果 第 1 次迭代的雅比矩陣 J 貴州大學課程設(shè)計 第 頁 29 。 3 4 65 100 0。 1 5 3 110 。 5 2]。 3 1。 貴州大學課程設(shè)計 第 頁 28 第六章 運行結(jié)果分析 正常運行結(jié)果分析 MATLAB 運行結(jié)果 用 MATLAB把 m文件,文件內(nèi)容如下: % (bus)(volt) (ang) (p) (q) (bus type) bus=[ 1 3。, lPQ(I,1),lPQ(I,2),real(lPQ(I,3)),imag(lPQ(I,3)),real(lPQ(I,4)),imag(lPQ(I,4)),real(lPQ(I,5)),imag(lPQ(I,5)))。 for I=1:nl, fprintf(myf, 39。節(jié)點 I 節(jié)點 J 線路功率 S(I,J) 線路功率 S(J,I) 線路損耗 dS(I,J)\n39。)。 end fprintf(myf, 39。 %2d %10f %10f %10f+j%10f \n39。)。 fprintf(myf, 39。 節(jié)點計算 結(jié)果: \n39。)。 fprintf(myf, 39。 global lPQ。 global nl。 for I=1:nl for J=1:2 for k=1:nb if lPQ(I,J)==nodenum(k,1) lPQ(I,J)=nodenum(k,2)。 tem(I,1)=I。 for I=1:nb for J=1:nb if nodenum(J,2)==I break。 global nodenum。 global line。 global nb。 lPQ(k,3:5)=[Sij Sji delS]。 Sji=Uj*conj(Iji)。 Iji=(UjUi)*Yt*K+Uj*Yt*K*(K1)。 end if K0 K=K。 end if K0 Iij=(UiUj)*Yt/K+Ui*(Ym+Yt*(K1)/K)。(J==0) Iij=Ui*Ym。 Iji=Uj*(Yt+Ym)Ui*Yt。 end if (K==0)amp。 Ui=bus(I,2)*(cos(bus(I,3))+j*sin(bus(I,3)))。 end Ym=line(k,5)+j*line(k,6)。 貴州大學課程設(shè)計 第 頁 26 if J~=0 Zt=line(k,3)+j*line(k,4)。 lPQ(k,1)=I。 for k=1:nl I=line(k,1)。 global bus。 end 計算線路功率及損耗: global nl。 sum=sum+bus(I,2)*bus(J,2)*B。 for J=1:nb ang=bus(I,3)bus(J,3)。 end bus(I,4)=sum。 A=real(Y(I,J))*cos(ang)+imag(Y(I,J))*sin(ang)。 for I=nPQ+1:nb if bus(I,6)==3 sum=0。 global line。 global nPQ。)。(max(abs(dPQ(Y,bus)))eps2) break end end if i==max1 error(39。, I,u(I,1))。 fprintf(myf, 39。, I,angl(I,1))。 fprintf(myf, 39。 u=x(nb:nb+nPQ1,1)。,i)。 fprintf(myf,39。 end for I=nb:nb+nPQ1 dx(I,1)=dx(I,1)*x(I,1)。\n39。, dx(I,1))。 for I=1:nb+nPQ1 fprintf(myf, 39。 第 %d 次迭代的節(jié)點相角和電壓的偏差dx\n39。)。 fprintf(myf, 39。%10e 39。,i)。 end fprintf(myf,39。\n39。, Jac(I,J))。 for I=1:nb+nPQ1 for J=1:nb+nPQ1 fprintf(myf, 39。第 %d次迭代的雅比矩陣 J\n39。,i)。 fprintf(myf,39。 del=dPQ(Y,bus)。u]。 u(:,1)=bus(1:nPQ,2)。)。 end fprintf(myf, 39。%10f+j*(%10f) 39。)。 fprintf(myf, 39。 貴州大學課程設(shè)計 第 頁 24 eps1=。 global myf。 global line。 global nPQ。 Y(J,I)=Y(I,J)。 Y(J,J)=Y(J,J)+K*K*Yt。 Y(J,I)=Y(I,J)。 Y(J,J)=Y(J,J)+Yt/K/K。(J==0) Y(I,I)=Y(I,I)+Ym。 Y(J,I)=Y(I,J)。 Y(J,J)=Y(J,J)+Yt+Ym。 if (K==0)amp。 end% Ym=line(k,5)+j*line(k,6)。 Zt=line(k,3)+j*line(k,4)。 for k=1:nl I=line(k,1)。 global myf。 global line。 global nb。 for I=1:nl for J=1:2 for k=1:nb if line(I,J)==nodenum(k,2) line(I,J)=nodenum(k,1)。 nodenum=[newbus bus(:,1)]。SW]。 end end bus=[PQ。 else nPQ = nPQ + 1。 elseif type == 2, nPV = nPV +1。 if type == 3, nSW = nSW + 1。 nPQ = 0。 nSW = 0。 global line。 global nPQ。 節(jié)點重新編號: function change global nb。w39。)。 貴州大學課程設(shè)計 第 頁 22 end outfile=strcat(outfile,39。 break。 for I=1:length(dfile) if dfile(I)==39。 outfile=39。 end [nl,ml]=size(line)。) else lfile =length(dfile)。 if pathname == 0 error(39。Select Data File39。*.m39。 global myf。 global line。 output。 flow。 NR。 change。 global myf。 global nodenum。 global line。 global nl。 global nSW nPQ nPV。 程序框圖如下: 輸 入 原 始 數(shù) 據(jù)形 成 節(jié) 點 導 納 矩 陣設(shè) 電 壓 初 值 e( 0 )、 f( 0 )設(shè) 迭 代 次 數(shù) k = 0計 算 誤 差 向 量? P( k )、 ? Q( k )、 ? U2( k )收 斂 否求 雅 可 比 矩 陣 元 素解 修 正 方 程 ,求 解 ? e( k )、 ? f( k )修 正 節(jié) 點 電 壓e( k + 1 )= e( k ) ? e( k )f( k + 1 )= f( k ) ? f( k )K = K + 1K K m a x不 收 斂 停 機求 P U 節(jié) 點 無 功 功 率 ,求 平 衡 節(jié) 點 功 率求 支 路 功 率 分 布 和 損 耗停 機是否 貴州大學
點擊復制文檔內(nèi)容
環(huán)評公示相關(guān)推薦
文庫吧 www.dybbs8.com
備案圖鄂ICP備17016276號-1