188

188 Rozdział 14

end

xx=xh(l:N); % dlugosc wektora x musi byc równa N harm=[];

% transformata Fouriera wydzielonej części sygnału X=xh(l:N)

Xh = f f t(xx,N) ;

fN=N/2; % częstotliwość Nyąuista %rozklad amplitudy harmonicznych harm=sqrt(Xh.*conj(Xh) )/N;

% częstotliwość w poszczególnych punktach transformaty fh = fp/N*(0:fN);

% przeniesienie amplitud czest. ujemnych na czest. dodatnie harm(2:fN) = 2*harm(2:fN);

% puste miejsca po ujemnych częstotliwościach w spektrum harm(fN+2:N) = [];

% obliczenie napiec fazowych ula=ea+uN; ulb=eb+uN; ulc=ec+uN;

%% przebiegi czasowe - wykresy subplot(323),plot(T,ulb); title(' Napięcie ulb'); grid on; subplot(324), plot(fh,harm);

title(' Amplitudy harmonicznych w uN'); grid on;

%axis([0 100 0 10]); % skalowanie osi subplot(325),plot(T,X(:,2)); title(' Prąd ilb '); grid on;

%subplot(223),plot(T,uN),title(' Napięcie uN');

subplot(326),plot(X(:,2),ulb);

title(' Punkty pracy (ilb,ulb)'); grid on;

% Wyniki obliczeń

fprintf(fd, '\n Czas symulacji tmax=%f s',tmax) fprintf(fd,’\n\n PARAMETRY ZASTĘPCZE PRZEKLADNIKA w pu') fprintf(fd,'\n    Ln=Lmpu    C0=Cspu    Lr \n')

fprintf(fd,'    %9.6f    %9.6f    %9.6f \n',Ln,CO,Lr)

Xmi=w*Ln; BC=w*C0; XC=1/BC; XCpXmi=XC/Xmi;

fprintf (fd, '\n\n XC=%f, Xmi = %f, XC/Xmi= %f \n\XC, Xmi, XCpXmi) fprintf(fd,'\n Warunki początkowe \n’)

fprintf(fd,¹\n ia0=%f, ib0=%f, ic0=%f, uN0=%f¹,iaO,ibO,icO,uNO)

fprintf(fd,'\n\n Napięcie punktu neutralnego');

fprintf(fd, '\n    t    uN(t) ') ;

nT=length(T) ;

for k=l:nT

fprintf(fd, '\n %8.4f %12.8f' ,T(k) ,uN(k)    );

end

fprintf(fd,'\n\n KONIEC SYMULACJI PRZEPIEC FERROREZONANSOWYCH programem apfsym()') end

fcloseCall'); return function fig = apfsym

%interfejs graficzny do analizy przepiec ferrorezonansowych global C vslider vedit rys format compact

hO = figure('Units','points',    ...