【正文】 end 。 3*x(2)+x(3)11, 3*x(1),x(1)。y(3)]。 y=[y(1)。 y(2)=3*x(1)*x(2)+x(1)*x(3)11*x(1)。,i) [x dx F(x)] %在屏幕上輸出每次的 x(i),dx(i),F(x(i)) if (max(abs(dx))eps)amp。 sprintf(39。 dx=A\b。 for i=1:max1, A=Fd(x)。 x0=[ ]39。 end format long max1=100。y(2)。 y(3)=exp(x(2)*x(3))+20*x(3)+(10*pi3)/3。 end 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 function y=F(x) y(1)=3*x(1)cos(x(2)*x(3))。 2*x(1), 162*(x(2)+),cos(x(3))。,i) [x dx F(x)] %在屏幕上輸出每次的 x(i),dx(i),F(x(i)) if (max(abs(dx))eps1)amp。 sprintf(39。 dx=A\b。 for i=1:max1, A=Fd(x)。 x0=[ ]39。 eps1=。 y=3*sin(a)。 end a=0:pi/12:2*pi。 for i=1:n y=y+(1)^(i1)*1/(2*i1)。 toc。 % 利用公式計算接地支路的潮流 end end end end 作業(yè) y=pifun(10000) 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 y = clear all tic。 else S_res(k,4)=bus(J,2)^2*conj(YtYm(k,5))。 % 利用公式計算非接地支路的潮流 else if(J==0 ) S_res(k,3)=bus(I,2)^2*conj(YtYm(k,4))。 S_res(k,4)=bus(J,2)^2*(conj(YtYm(k,3))+conj(YtYm(k,5)))bus(I,2)*bus(J,2)*(conj(cos(bus(I,3))+j*sin(bus(I,3))))*(cos(bus(J,3))+j*sin(bus(J,3)))*conj(YtYm(k,3))。 if (J~=0)amp。 % 前兩列為節(jié)點編號 for k=1:nl I = S_res(k,1)。 S_res = zeros(nl,5)。 % 相角采用角度制 bus_res(:,4) = bus(:,4) + (sqrt(1))*bus(:,5)。 % bus_res矩陣儲存著節(jié)點計算結(jié)果 bus_res(:,1:2) = bus(:,1:2)。 end end 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 計算并返回節(jié)點數(shù)據(jù)結(jié)果 function bus_res = bus_res(bus) [nb,mb]=size(bus)。 YtYm(k,4) = Ym+Yt*(1+K)。 YtYm(k,5) = Yt*(1K)/K/K。 end if K0 % 變壓器線路 : Zt和 Ym為折算到 i側(cè)的值 ,K在 j側(cè) YtYm(k,3) = Yt/K。 end if (K==0)amp。 YtYm(k,4) = 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)。 % 對 YtYm矩陣賦初值 0 YtYm(:,1:2) = line(:,1:2)。 break end end end end % 恢復 line的編號 作用為計算線路的等效 Yt和 Ym，以計算線路潮流 function YtYm = YtYm(line) [nl,ml]=size(line)。 end end end % 利用 bus矩陣的首列 編號重新對 bus矩陣排序并存入bus_temp矩陣中 bus = bus_temp。 for i = 1 :nb for j = 1 : nb if bus(j,1) == k bus_temp(k,:) = bus(j,:)。 bus_temp = zeros(nb,mb)。 end % 計算平衡節(jié)點的無功及有功注入 作用為對節(jié)點和線路數(shù)據(jù)恢復編號 function [bus,line] = ReNum(bus,line,nodenum) [nb,mb]=size(bus)。 end end % 利用公式計算 PV節(jié)點的無功注入 for j =1:n bus(n,4)=bus(n,4)+bus(n,2)*bus(j,2)*(real(Y(n,j))*cos(bus(n,3)bus(j,3))+imag(Y(n,j))*sin(bus(n,3)bus(j,3)))。 % 生成雅克比矩陣 作用為計算每個節(jié)點的功率注入 function bus = PQ(bus,Y,nPQ,nPV) n = nPQ+nPV+1。 end end % 分別計算 L矩陣的對角及非對角元素 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 end end end J = [H N。 end % 分別計算 K矩陣的對角及非對角元素 if j nPQ+1 if i~=j L(i,j)=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)))。 end end % 分別計算 N矩陣的對角及非對角元素 if i nPQ+1 if i~=j K(i,j)=bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)bus(j,3))+imag(Y(i,j))*sin(bus(i,3)bus(j,3)))。 end % 分別計算 H矩陣的對角及非對角元素 if j nPQ+1 if i~=j N(i,j)=bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)bus(j,3))+imag(Y(i,j))*sin(bus(i,3)bus(j,3)))。 end end % 初始化并計算 Qi和 Pi for i = 1:nb1 for j = 1:nb1 if i~=j H(i,j)=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)))。 for i = 1:nb1 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 for j = 1:nb Pi(i,1)=Pi(i,1)+bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)bus(j,3))+imag(Y(i,j))*sin(bus(i,3)bus(j,3)))。 K L]，并初始化 Qi = zeros(nb1,1)。 L = zeros(nPQ,nPQ)。 N = zeros(nb1,nPQ)。 end end end % 利用循環(huán)計算求取 dP和 dQ 作用為計算雅克比矩陣 function J = Jac(bus,Y,nPQ) [nb,mb]=size(bus)。 % 對 dP和 dQ賦初值 PV節(jié)點不需計算 dQ 平衡節(jié)點不參與計算 for i = 1:n1 for j = 1:n dP(i,1) = dP(i,1)bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)bus(j,3))+imag(Y(i,j))*sin(bus(i,3)bus(j,3)))。 % 總節(jié)點個數(shù) dP = bus(1:n1,4)。 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。 Y(I,I)=Y(I,I)+Ym。 end if (K==0)amp。 Y(I,J)=Y(I,J)Yt。 Y(I,I)=Y(I,I)+Yt+Ym。 if (K==0)amp。 Ym=line(k,5)+j*line(k,6)。 Zt=line(k,3)+j*line(k,