Диссертация (1144079), страница 15
Текст из файла (страница 15)
and then wecan get variable power flow, correlatively.n_row = 1000;for n = 1:nvarData(:,n)=Xmin(n)+(Xmax(n)-Xmin(n))*rand(n_row,1);endclear n;for n = 1:n_row% define variable power flows% define branch's voltagefor nn = 1:size(A,2)Data(n,nvar+nn)=(Ku(1,nn)+Ku(2,nn)*Data(n,1)+Ku(3,nn)*Data(n,2)+Ku(4,nn)*Data(n,1)*Data(n,2))/(1+Ki(3,Lo(ncase,1))*Data(n,1)+Ki(4,Lo(ncase,1))*Data(n,2)+Ki(5,Lo(ncase,1))*Data(n,1)*Data(n,2));119endclear nn;% define branch's currentsfor nb=1:size(A,2)if nb == Lo(ncase,1)Data(n,nvar+size(A,2)+nb)=(Ki(1,nb)+Ki(2,nb)*Data(n,2))/(1+Ki(3,Lo(ncase,1))*Data(n,1)+Ki(4,Lo(ncase,1))*Data(n,2)+Ki(5,Lo(ncase,1))*Data(n,1)*Data(n,2));elseif nb == Lo(ncase,2)Data(n,nvar+size(A,2)+nb)=(Ki(1,nb)+Ki(2,nb)*Data(n,1))/(1+Ki(3,Lo(ncase,1))*Data(n,1)+Ki(4,Lo(ncase,1))*Data(n,2)+Ki(5,Lo(ncase,1))*Data(n,1)*Data(n,2));elseData(n,nvar+size(A,2)+nb)=(Ki(1,nb)+Ki(2,nb)*Data(n,1)+Ki(3,nb)*Data(n,2)+Ki(4,nb)*Data(n,1)*Data(n,2))/(1+Ki(3,Lo(ncase,1))*Data(n,1)+Ki(4,Lo(ncase,1))*Data(n,2)+Ki(5,Lo(ncase,1))*Data(n,1)*Data(n,2));endend; clear nbData(n,nvar+2*size(A,2)+1:nvar+3*size(A,2))=Data(n,nvar+1:nvar+size(A,2)).*Data(n,nvar+size(A,2)+1:nvar+2*size(A,2));%PvarData(n,nvar+3*size(A,2)+1)=sum(Data(n,nvar+2*size(A,2)+1:nvar+3*size(A,2)));%sum Pvarend%lay nvalue gia tri tot nhatData(:,nvar+3*size(A,2)+2) = 1./Data(:,nvar+3*size(A,2)+1);nvalue =20;Data=sortrows(Data,nvar+3*size(A,2)+1);Data=Data(1:nvalue,:);divpart=Data(1:nvalue,nvar+3*size(A,2)+2)/sum(Data(1:nvalue,nvar+3*size(A,2)+2));test(1)=0;for k=1:length(divpart)test(k+1)=test(k)+divpart(k);endnloop = 200;alpha(1)=1;beta=1000^(1/300);120for k=2:nloopalpha(k)=alpha(k-1)/beta;endgamma = alpha/2;for n = 1:nlooph=nvalue;while h<n_rowt=rand;k=1;while ((t<test(k)) || (t>=test(k+1)))%k=k+1;endt=ceil(nvar*rand);tempo = Data(k,:);tempo(t) = tempo(t)*(1+alpha(n)*2*(rand-0.5));if(tempo(t) < Xmin(t)|| tempo(t) > Xmax(t))% change x1,x2,x3,...
in[Bmin Bmax] !!!continueelsedelta=abs((tempo(1)-Data(k,1))/tempo(1))+abs((tempo(2)Data(k,2))/tempo(2));if delta <= gamma(n)h = h+1;Data(h,:) = tempo;% lay x1, x2, x3,...% define variable power flows% define branch's voltagefor nn = 1:size(A,2)Data(h,nvar+nn)=(Ku(1,nn)+Ku(2,nn)*Data(h,1)+Ku(3,nn)*Data(h,2)+Ku(4,nn)*Data(h,1)*Data(h,2))/(1+Ki(3,Lo(ncase,1))*Data(h,1)+Ki(4,Lo(ncase,1))*Data(h,2)+Ki(5,Lo(ncase,1))*Data(h,1)*Data(h,2));endclear nn;% define branch's currentsfor nb=1:size(A,2)if nb == Lo(ncase,1)Data(h,nvar+size(A,2)+nb)=(Ki(1,nb)+Ki(2,nb)*Data(h,2))/(1+Ki(3,Lo(ncase,1))*Data(h,1)+Ki(4,Lo(ncase,1))*Data(h,2)+Ki(5,Lo(ncase,1))*Data(h,1)*Data(h,2));elseif nb == Lo(ncase,2)121Data(h,nvar+size(A,2)+nb)=(Ki(1,nb)+Ki(2,nb)*Data(h,1))/(1+Ki(3,Lo(ncase,1))*Data(h,1)+Ki(4,Lo(ncase,1))*Data(h,2)+Ki(5,Lo(ncase,1))*Data(h,1)*Data(h,2));elseData(h,nvar+size(A,2)+nb)=(Ki(1,nb)+Ki(2,nb)*Data(h,1)+Ki(3,nb)*Data(h,2)+Ki(4,nb)*Data(h,1)*Data(h,2))/(1+Ki(3,Lo(ncase,1))*Data(h,1)+Ki(4,Lo(ncase,1))*Data(h,2)+Ki(5,Lo(ncase,1))*Data(h,1)*Data(h,2));endend; clear nbData(h,nvar+2*size(A,2)+1:nvar+3*size(A,2))=Data(h,nvar+1:nvar+size(A,2)).*Data(h,nvar+size(A,2)+1:nvar+2*size(A,2));%PvarData(h,nvar+3*size(A,2)+1)=sum(Data(h,nvar+2*size(A,2)+1:nvar+3*size(A,2)));%sum PvarendendendData(:,nvar+3*size(A,2)+2) = 1./Data(:,nvar+3*size(A,2)+1);Data=sortrows(Data,nvar+3*size(A,2)+1);Data=Data(1:nvalue,:);divpart=Data(1:nvalue,nvar+3*size(A,2)+2)/sum(Data(1:nvalue,nvar+3*size(A,2)+2));test(1)=0;for k=1:length(divpart)test(k+1)=test(k)+divpart(k);endclear k;end%Result(ncase,:) = Data(1,:);opti(p,:,ncase) = Data(1,:);clear Data;endendfname = sprintf('ResultPar%d%de%dill.mat', Bus, Times, 1000*eps);save(fname);clear fname;Прил.
13clear;clc;format long e;122tab = fopen('circuit.txt');TAB = fscanf(tab, '%d %d %d %d %d %d %d %d %d %d %d %d', [13 inf]);fclose(tab);A = TAB';NYG = 1; nBranch = (size(A,2)-NYG)/3; nNode = (size(A,1)-NYG)/3;lw = 1;Fz = 18; kF =0.9;% Set fontsizeBUS = [1 2 3 4 5];KZ = [0.8 0.8 0.8 0.8 1; 0.7 0.7 0.7 0.7 1];Lim(1,1:5) = 210;%\Lim(2,1:5) = 225;%/Limit of V(t).Lim(3,1:5) = 0;%\Lim(4,1:5) = [1100 500 1000 400 200]; %/Limit of I(t).Lim(5,1:5) = 0;%\Lim(6,1:5) = [5e4 3e4 4e4 4e4 4e4]; %/Limit of P(t).Lim(7,1:5) = 0;%\Lim(8,1:5) = [800 380 650 260 160]; %/Limit of I1(t).Lim(9,1:5) = 0;%\Lim(10,1:5) = [175 inf inf 55 60]; %/Limit of I2,I0(t).for nf = 5%1:length(BUS)Bus = BUS(nf)fname = sprintf('RealData%d.txt',Bus);tab = fopen(fname);clear fname;colon = char(58);[FG,count1] = fscanf(tab,['%d' colon '%d' '%f' '%f' '%f' '%f' '%f' '%f' '%f''%f' '%f' '%f' '%f' '%f' '%f' '%f' '%f' '%f' '%f' '%f'],[20 inf]);fclose(tab);for p = 1:size(FG,2)%ZA(1:nBranch)=[1e-9;0.001063+j*0.005905;KZ(1,Bus)*FG(18,p)*(FG(15,p)-j*sin(acos(FG(15,p)))); 1.31429-j*8.89675];%ZB(1:nBranch)=[1e-9;0.001063+j*0.005905;KZ(2,Bus)*FG(19,p)*(FG(16,p)-j*sin(acos(FG(16,p)))); 1.31429-j*8.89675];%ZC(1:nBranch) = [1e-9; 0.001063+j*0.005905; FG(20,p)*(FG(17,p)j*sin(acos(FG(17,p)))); 1.31429-j*8.89675];ZA(1:nBranch)=[1e-9;0.001063+j*0.005905;KZ(1,Bus)*FG(18,p)*FG(15,p)-KZ(1,Bus)*FG(18,p)*j*sin(acos(FG(15,p)));1.31429-j*8.89675];ZB(1:nBranch)=[1e-9;0.001063+j*0.005905;KZ(1,Bus)*FG(18,p)*FG(15,p)-KZ(2,Bus)*FG(19,p)*j*sin(acos(FG(16,p)));1.31429-j*8.89675];123ZC(1:nBranch)=[1e-9;0.001063+j*0.005905;KZ(1,Bus)*FG(18,p)*FG(15,p)-FG(20,p)*j*sin(acos(FG(17,p)));1.31429j*8.89675];if NYG == 0ZN = [];elseZN(1:NYG) = 0.0017+j*0.0089;% if the circuit contain Zn!!!endZ=[ZA ZB ZC ZN];Y=1./Z;clear ZA ZB ZC ZN;Y = diag(Y);E = zeros(size(A,2),1);E(1,1) = 220;E(5,1) = 220*exp(-2*pi*i/3);E(9,1) = 220*exp(2*pi*i/3);Z1 = Z;Y1 = diag(1./Z1);% Solve circuitU0 = -(A*Y1*A')\(A*Y1*E);Us = A'*U0;U = Us+E;I = Y1*U;Data(p,1:size(A,2)) = U';Data(p,1+size(A,2):2*size(A,2)) = I.';Data(p,1+2*size(A,2):3*size(A,2)) = U.'.*I.';end%figure%subplot(3,1,1)%plot(abs(Data(:,3)));%hold on;grid on;%plot(abs(Data(:,7)));%plot(abs(Data(:,11)));%subplot(3,1,2)plot_color = [0 0 0];left_color = [0 0 0];right_color = [0 0 0];fig = figure;set(fig,'defaultAxesColorOrder',[plot_color]);subplot(2,1,1)plot(abs(Data(:,16)),'linewidth', lw);124hold on;grid on;ylim([Lim(3,Bus) Lim(4,Bus)]);plot(abs(Data(:,20)),'--','linewidth', lw);plot(abs(Data(:,24)),':','linewidth', lw);legend('Phase A','Phase B','Phase C');title('Currents [A]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );xlabel('24 hours');set(gca,'FontSize', Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);subplot(2,1,2)plot(abs(Data(:,29)+Data(:,33)+Data(:,37)),'linewidth', lw);grid on;ylim([Lim(5,Bus) Lim(6,Bus)]);title('P_v [V.A]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );xlabel('24 hours');set(gca,'FontSize', Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);fname = sprintf('T%d without CS',Bus);print(fname, '-dpng', '-r300'); clear fname;endПрил.
14clear;clc;close allS = [1 1 1; exp(-2*pi*i/3) exp(2*pi*i/3) 1; exp(2*pi*i/3) exp(-2*pi*i/3) 1];Fz = 18; kF =0.9;lw = 1;BUS = [1 2 3 4 5];TIMES = 1;%Good resultLX = [5 1 1 3 4];%Optimal LocalEPS = [0 0.005 0.01];Lim(1,1:5) = 210;%\Lim(2,1:5) = 225;%/Limit of V(t).Lim(3,1:5) = 0;%\Lim(4,1:5) = [1100 800 1000 400 300]; %/Limit of I(t).Lim(5,1:5) = 0;%\Lim(6,1:5) = [2.5e5 1e5 2e5 1e5 7e4]; %/Limit of P(t).Lim(7,1:5) = 0;%\Lim(8,1:5) = [800 380 650 260 200]; %/Limit of I1(t).Lim(9,1:5) = 0;%\Lim(10,1:5) = [175 inf inf inf inf]; %/Limit of I2,I0(t).Lim(11,1:5) = 0;%\125Lim(12,1:5) = [5e4 3e4 4e4 4e4 4e4]; %/Limit of I2,I0(t).for nBus = 5%1:length(BUS)Bus = BUS(nBus); % 1 2 4 5 6;Times = TIMES;lx = LX(Bus);for ie = 2%1:length(EPS)eps = EPS(ie);fname = sprintf('ResultPar%d%de%dill.mat', Bus, Times, 1000*eps);load(fname);clear fname;for n = 1:size(opti,1)sequence(n,1:3) = (S\[opti(n,5);opti(n,9);opti(n,13)]);sequence(n,4:6) = (S\[opti(n,18);opti(n,22);opti(n,26)]);endmsequence = abs(sequence);mopti = abs(opti(:,:,lx));G=[mopti msequence];period = 5;for m=1:15for n = 1:(size(G,1)-period)G(n,1:size(G,2)) = mean(G(n:n+period,1:size(G,2)));end;for m = 1:periodG(n+m,1:size(G,2)) = G(n,1:size(G,2));endendplot_color = [0 0 0];left_color = [0 0 0];right_color = [0 0 0];fig = figure;set(fig,'defaultAxesColorOrder',[plot_color]);subplot(3,1,1)plot(abs(G(:,5)),'LineWidth', lw);hold on;grid on;ylim([Lim(1,Bus) Lim(2,Bus)]);plot(abs(G(:,9)),'--','LineWidth', lw);plot(abs(G(:,13)),':','LineWidth', lw);legend('Phase A','Phase B','Phase C');title('Voltages [V]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );xlabel('24 hours');set(gca,'FontSize', Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);126subplot(3,1,2)plot(abs(G(:,18)),'LineWidth', lw);hold on;grid on;ylim([Lim(3,Bus) Lim(4,Bus)]);plot(abs(G(:,22)),'--','LineWidth', lw);plot(abs(G(:,26)),':','LineWidth', lw);title('Currents [A]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );xlabel('24 hours');set(gca,'FontSize', Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);subplot(3,1,3)plot(G(:,31),'LineWidth', lw);hold on;grid on;ylim([Lim(5,Bus) Lim(6,Bus)]);plot(G(:,35),'--','LineWidth', lw);plot(G(:,39),':','LineWidth', lw);%legend('Phase A','Phase B','Phase C');title('Apparent powers [V.A]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );xlabel('24 hours');set(gca,'FontSize', Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);fname = sprintf('Pvar T%d ',Bus);print(fname, '-dpng', '-r300'); clear fname;fig = figure;set(fig,'defaultAxesColorOrder',[plot_color]);subplot(2,1,1)plot(G(:,18),'LineWidth', lw);hold on;grid on;ylim([Lim(3,Bus) Lim(4,Bus)]);plot(G(:,22),'--','LineWidth', lw);plot(G(:,26),':','LineWidth', lw);legend('Phase A','Phase B','Phase C');title('Currents [A]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );xlabel('24 hours');set(gca,'FontSize', Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);127subplot(2,1,2)hold on;grid on;plot(G(:,42),'LineWidth', lw);ylim([Lim(11,Bus) Lim(12,Bus)]);title('P_v [V.A]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );xlabel('24 hours');set(gca,'FontSize', Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);fname = sprintf('Pvar T%d I and Pv',Bus);print(fname, '-dpng', '-r300'); clear fname;fig = figure;set(fig,'defaultAxesColorOrder',[left_color ; right_color ]);yyaxis leftplot(G(:,47),'LineWidth', lw);hold on;grid on;ylim([Lim(7,Bus) Lim(8,Bus)]);yyaxis rightylim([Lim(9,Bus) Lim(10,Bus)]);plot(G(:,48),'--','LineWidth', lw);plot(G(:,49),':','LineWidth', lw);legend('I1','I2','I0');title('Sequence components of Currents [A]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );xlabel('24 hours');set(gca,'FontSize', Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);fname = sprintf('Pvar T%d I120 ',Bus);print(fname, '-dpng', '-r300'); clear fname;fig = figure;set(fig,'defaultAxesColorOrder',[plot_color]);plot(G(:,1),'--','LineWidth', lw);ylim([0 inf]);hold on;grid on;plot(G(:,2),'LineWidth', lw);LineWidth = lw;legend('X1','X2');title('Reactance values [Ohm]');set(gca,'TitleFontWeight','normal');set(gca,'XTickLabel',[''] );128xlabel('24 hours');set(gca,'FontSize',Fz );set(gca,'TitleFontSizeMultiplier',kF);set(gca,'LabelFontSizeMultiplier',kF);fname = sprintf('Pvar T%d value of X',Bus);print(fname, '-dpng', '-r300'); clear fname;fig = figure;set(fig,'defaultAxesColorOrder',[plot_color]);plot(opti(:,1,lx),'LineWidth', lw);hold on;grid on;plot(opti(:,2,lx),'LineWidth', lw);endend129СПИСОК ЛИТЕРАТУРЫ1.
Grainger J.J., Stevenson W.D. Power System Analysis. New York: McGrawHill, 1994.2. Братцев К.Е. Исследование функциональных возможностей построенияинформационно измерительных систем параметров двухполюсныхэлектрических систем на синусоидальном сигнале фиксированнойчастоты / Братцев К.Е. // Современные информационные технологии.2008.3. Карпов В.Н. Универсальный модуль информационно-измерительнойсистемы для измерения энергетических параметров в потребительскихэнергетических системах апк / Карпов В.Н., Халатов А.Н., Юлдашев З.Ш.,Юлдашев Р.З.
// Известия Санкт-Петербургского государстве.4. Киреев А.О. Информационно-измерительная система для мониторинга ианализа энергопотребления беспроводных сенсорных систем / КиреевА.О. // диссертация на соискание ученой степени кандидата техническихнаук. Пензенский государственный университет. Пенза, 2.5. Герке М. Система электроснабжения для измерительного устройства испособ эксплуатации такой системы / Герке М., Бёель С. // патент наизобретение RUS 2305868 11.04.2003.6. Веселов О.В., Информационно-измерительная система для диагностикиэлектромеханических систем / Веселов О.В., Кононыхина Н.А.,Перепёлкин В.М.