【正文】 K=line(k,7)。 if (K==0)amp。(J~=0) % 普通線路: K=0 YtYm(k,3) = Yt。 YtYm(k,4) = Ym。 YtYm(k,5) = Ym。 end if (K==0)amp。(J==0) % 對地支路: K=0,J=0,R=X=0 YtYm(k,4) = Ym。 end if K0 % 變壓器線路: Zt和Ym為折算到i側(cè)的值,K在j側(cè) YtYm(k,3) = Yt/K。 YtYm(k,4) = Ym+Yt*(K1)/K。 YtYm(k,5) = Yt*(1K)/K/K。 end if K0 % 變壓器線路: Zt和Ym為折算到K側(cè)的值,K在i側(cè) YtYm(k,3) = Yt*K。 YtYm(k,4) = Ym+Yt*(1+K)。 YtYm(k,5) = Yt*(K^2+K)。 endend 計(jì)算并返回節(jié)點(diǎn)數(shù)據(jù)結(jié)果function bus_res = bus_res(bus)[nb,mb]=size(bus)。bus_res = zeros(nb,4)。 % bus_res矩陣儲存著節(jié)點(diǎn)計(jì)算結(jié)果bus_res(:,1:2) = bus(:,1:2)。bus_res(:,3) = bus(:,3) *180 / pi。 % 相角采用角度制bus_res(:,4) = bus(:,4) + (sqrt(1))*bus(:,5)。 % 注入功率end 計(jì)算并返回線路潮流function S_res = S_res(bus,line,YtYm)[nl,ml]=size(line)。S_res = zeros(nl,5)。 % S_res矩陣儲存著線路潮流計(jì)算結(jié)果S_res(:,1:2) = line(:,1:2)。 % 前兩列為節(jié)點(diǎn)編號for k=1:nl I = S_res(k,1)。 J = S_res(k,2)。 if (J~=0)amp。(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(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))。 S_res(k,5)=S_res(k,3) + S_res(k,4)。 % 利用公式計(jì)算非接地支路的潮流 else if(J==0 ) S_res(k,3)=bus(I,2)^2*conj(YtYm(k,4))。 S_res(k,5)=S_res(k,3)+S_res(k,4)。 else S_res(k,4)=bus(J,2)^2*conj(YtYm(k,5))。 S_res(k,5)=S_res(k,3)+S_res(k,4)。 % 利用公式計(jì)算接地支路的潮流 end endendend 作業(yè)y=pifun(10000)y = clear all tic。 y=pifun(100000)。 toc。Elapsed time is seconds.function y= pifun(n)y=0。for i=1:n y=y+(1)^(i1)*1/(2*i1)。end y=y*4。end a=0:pi/12:2*pi。 x=5*cos(a)。 y=3*sin(a)。 plot(x,y) clearformat longmax1=100。 eps1=。eps2=。 x0=[ ]39。 x=x0。for i=1:max1, A=Fd(x)。 b=F(x)。 dx=A\b。 %矩陣左除，即b除以A x=xdx。 sprintf(39。第%d次迭代結(jié)果：\n39。,i) [x dx F(x)] %在屏幕上輸出每次的x(i),dx(i),F(x(i)) if (max(abs(dx))eps1)amp。(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))。 2*x(1), 162*(x(2)+),cos(x(3))。 0,x(3)*exp(x(2)*x(3)),x(2)*exp(x(2)*x(3))+20]。endfunction y=F(x) y(1)=3*x(1)cos(x(2)*x(3))。 y(2)=x(1)^281*(x(2)+)^2+sin(x(3))+。 y(3)=exp(x(2)*x(3))+20*x(3)+(10*pi3)/3。 y=[y(1)。y(2)。y(3)]。 end format longmax1=100。eps=10^10。x0=[ ]39。x=x0。for i=1:max1,A=Fd(x)。b=F(x)。dx=A\b。 %矩陣左除，即b除以Ax=xdx。sprintf(39。第%d次迭代結(jié)果：\n39。,i)[x dx F(x)] %在屏幕上輸出每次的x(i),dx(i),F(x(i))if (max(abs(dx))eps)amp。(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。 y(2)=3*x(1)*x(2)+x(1)*x(3)11*x(1)。 y(3)=2*x(2)*x(3)+40*x(1)。 y=[y(1)。y(2)。y(3)]。 endfunction J=Fd(x) J=[1,10*x(2),14*x(3)。 3*x(2)+x(3)11, 3*x(1),x(1)。 40,2*x(3),2*x(2)]。end