基于牛頓拉夫遜法潮流計(jì)算的matlab實(shí)驗(yàn)報(bào)告定稿-wenkub.com
2025-06-24 20:13 本頁(yè)面
【正文】 40,2*x(3),2*x(2)]。y(2)。(max(abs(F(x)))eps)breakendendans =第1次迭代結(jié)果：ans = +002 * ans =第2次迭代結(jié)果：ans = ans =第3次迭代結(jié)果：ans = ans =第4次迭代結(jié)果：ans = ans =第5次迭代結(jié)果：ans = ans =第6次迭代結(jié)果：ans = 0 0 0 x = 1 5 4function y=F(x) y(1)=x(1)5*x(2)^2+7*x(3)^2+12。 %矩陣左除，即b除以Ax=xdx。x=x0。y(3)]。 y(2)=x(1)^281*(x(2)+)^2+sin(x(3))+。(max(abs(F(x)))eps2) break endend ans =第1次迭代結(jié)果：ans = ans =第2次迭代結(jié)果：ans = ans =第3次迭代結(jié)果：ans = ans =第4次迭代結(jié)果：ans = ans =第5次迭代結(jié)果：ans = ans =第6次迭代結(jié)果：ans = ans =第7次迭代結(jié)果：ans = ans =第8次迭代結(jié)果：ans = 0 0 x = function J=Fd(x) J=[3,x(3)*sin(x(2)*x(3)),x(2)*sin(x(2)*x(3))。 %矩陣左除，即b除以A x=xdx。 x=x0。 plot(x,y) clearformat longmax1=100。end y=y*4。 y=pifun(100000)。 S_res(k,5)=S_res(k,3)+S_res(k,4)。(I~=0) S_res(k,3)=bus(I,2)^2*(conj(YtYm(k,3))+conj(YtYm(k,4)))bus(I,2)*bus(J,2)*(cos(bus(I,3))+j*sin(bus(I,3)))*(conj(cos(bus(J,3))+j*sin(bus(J,3))))*conj(YtYm(k,3))。 % S_res矩陣儲(chǔ)存著線路潮流計(jì)算結(jié)果S_res(:,1:2) = line(:,1:2)。bus_res(:,3) = bus(:,3) *180 / pi。 YtYm(k,5) = Yt*(K^2+K)。 YtYm(k,4) = Ym+Yt*(K1)/K。 YtYm(k,5) = Ym。 K=line(k,7)。 J=line(k,2)。YtYm = zeros(nl,5)。 k = k + 1。[nl,ml]=size(line)。 % n為總節(jié)點(diǎn)數(shù)for i = nPQ+1:n1 for j = 1:n bus(i,5)=bus(i,5)+bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)bus(j,3))imag(Y(i,j))*cos(bus(i,3)bus(j,3)))。 else L(i,j)=Qi(i,1)+imag(Y(i,j))*((bus(i,2))^2)。 else N(i,j)=Pi(i,1)real(Y(i,j))*((bus(i,2))^2)。 Qi(i,1)=Qi(i,1)+bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)bus(j,3))imag(Y(i,j))*cos(bus(i,3)bus(j,3)))。 % 將雅克比矩陣分塊，即：J = [H N。H = zeros(nb1,nb1)。dQ = bus(1:nPQ,5)。Y(I,J)=Y(I,J)+K*Yt。Y(I,J)=Y(I,J)Yt/K。(J==0) % 對(duì)地支路: K=0,J=0,R=X=0。Y(J,J)=Y(J,J)+Yt+Ym。K=line(k,7)。 %讀入線路參數(shù)J=line(k,2)。breakendendendend 作用為計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣function Y = y(bus,line)[nb,mb]=size(bus)。newbus=[1:nb]39。endend。elseif type == 2,nPV = nPV +1。 % number of PV bus counternPQ = 0。 % 迭代結(jié)束后繼續(xù)在“”寫(xiě)入節(jié)點(diǎn)計(jì)算結(jié)果和線路計(jì)算結(jié)果 程序結(jié)束 作用為對(duì)節(jié)點(diǎn)重排序，并修改相應(yīng)的線路數(shù)據(jù)function [bus,line,nPQ,nPV,nodenum] = Num(bus,line)[nb,mb]=size(bus)。fprintf(myf,39。,real(S_res(i,3)),imag(S_res(i,3)))。% 39。for i = 1:nlfprintf(myf,39。,real(bus_res(i,4)),imag(bus_res(i,4)))。% 39。fprintf(myf,39。)。a39。 % 計(jì)算節(jié)點(diǎn)數(shù)據(jù)結(jié)果的子程序S_res = S_res(bus,line,YtYm)。(max(abs(dAng))EPS)break end % 判斷是否滿足精度誤差，如滿足則跳出，否則返回繼續(xù)迭代endbus = PQ(bus,Y,nPQ,nPV)。 % 計(jì)算相角修正量dU = UD*(dAngU(nb:nb+nPQ1,1))。for i = 1:nPQUD(i,i) = bus(i,2)。 % 在當(dāng)前目錄下生成“”文件，寫(xiě)入節(jié)點(diǎn)導(dǎo)納矩陣 format long EPS = 。fprintf(myf,39。39。end [nb,mb]=size(bus)。if pathname == 0 error(39。*.m39。 20 19 0 0 ] 。 35 22 0 0 。 31 6 0 0 。 5 0 0 0 0 0。 28 26 0 0。 23 22 0 0。 24 16 0 0。 16 15 0 0。 11 10 0 0。 11 6 0 0。 14 4 0 0。 25 2 0 0。2) 負(fù)荷數(shù)據(jù)見(jiàn)表D1表D1 負(fù)荷 數(shù)據(jù)節(jié)點(diǎn)號(hào)有功（MW）無(wú)功（Mvar）節(jié)點(diǎn)號(hào)有功（MW）無(wú)功（Mvar）BUS 3BUS 4BUS 7BUS 8BUS 12BUS 15BUS 16BUS 18BUS 20BUS 21BUS 23BUS 24BUS 25BUS 26BUS
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號(hào)-1