Podstawy cyfrowej automatyki elektroenergetycznej
Wydział Elektryczny Semestr VI
Anna Knap 170739	Ćwiczenie : Adaptacyjne algorytmy pomiarowe	Termin zajęć: Poniedziałek 13-15
				Ocena:
23.05.2011r.

1. Cel ćwiczenia:

Celem ćwiczenia było zbadanie zaproponowanie schemat blokowego oraz napisanie algorytmu adaptacyjnego pomiaru amplitudy poprzez adaptację długości okien filtrów do częstotliwości pomiarów.

2. Przebieg ćwiczenia:

Algorytm:

M-plik programu bez adaptacji :

%% adaptacyjne algorytmy pomiarowe fp=1200; f=50; N=fp/f; w=2*pi*50; Au=4; Ai=1; omg=2*pi/N; t=0:1/fp:0.2; m=1; k=1; u=Au*sin(w*t); i=Ai*sin(w*t); %% pomiar czestotliwosci (przez zliczanie impulsów) %dane fp=1200; f1=50; n1=fp/f1; w1=2*pi/n1; d=0; fsd=f1*(1+d); m05=10; a=0; b=0; c1=0; c2=0; nt=1:4*n1; %pomiar for k=2:length(u); if(u(k-1)>0) & (u(k)>0) ,b=0; a=a+1; end

if (u(k-1)<0) & (u(k)<0), a=0; b=b+1; end if (u(k-1)<0) & (u(k)>0); c1=u(k)/(u(k)-u(k-1)); m05=b+1-c1+c2; end if (u(k-1)>0)& (u(k)<0); c2=u(k)/(u(k)-u(k-1)); m05=a+1-c2+c1; end fm(k)=0.5*fp/m05;fs(k)=fsd; end %% ortogonalizacja sinus-cosinus ic1=0;is1=0 for n=25:length(t) k1=1:N%(n); omg=2*pi/N%(n); Lsin=sin((k1+0.5)*omg)*(2/N); Lcos=cos((k1+0.5)*omg)*(2/N); for k=1:N%(n) is=is1+Lsin(k)*i(n-k); ic=ic1+Lcos(k)*i(n-k); is1=is; ic1=ic; end I(n)=sqrt((ic.^2)+(is.^2)); is1=0; ic1=0; end

Przebiegi:

Przebiegi prądów i napięć	Pomiar częstotliwości
Pomiar amplitudy

M-plik programu z adaptacją :

%% adaptacyjne algorytmy pomiarowe fp=1200; f=50; N=fp/f; w=2*pi*50; Au=4; Ai=1; omg=2*pi/N; t=0:1/fp:0.2; m=1; k=1; u=Au*sin(w*t); i=Ai*sin(w*t); %% pomiar czestotliwosci (przez zliczanie impulsów) fp=1200; f1=50; n1=fp/f1; w1=2*pi/n1; d=0; fsd=f1*(1+d); m05=10; a=0; b=0; c1=0; c2=0; nt=1:4*n1; %pomiar for k=2:length(u); if(u(k-1)>0) & (u(k)>0) ,b=0; a=a+1; end

if (u(k-1)<0) & (u(k)<0), a=0; b=b+1; end if (u(k-1)<0) & (u(k)>0); c1=u(k)/(u(k)-u(k-1)); m05=b+1-c1+c2; end if (u(k-1)>0)& (u(k)<0); c2=u(k)/(u(k)-u(k-1)); m05=a+1-c2+c1; end fm(k)=0.5*fp/m05;fs(k)=fsd; end N=round(fp./Fm); %% ortogonalizacja sinus-cosinus ic1=0;is1=0 for n=25:length(t) k1=1:N(n) omg=2*pi/N(n) Lsin=sin((k1+0.5)*omg)*(2/N(n)); Lcos=cos((k1+0.5)*omg)*(2/N(n)); for k=1:N(n) is=is1+Lsin(k)*i(n-k); ic=ic1+Lcos(k)*i(n-k); is1=is; ic1=ic; end I(n)=sqrt((ic.^2)+(is.^2)); is1=0; ic1=0; end

Przebiegi:

Pomiar częstotliwości	Pomiar amplitudy

Wpływ zakłóceń:

a) Stałą nieokresową:

ta=0.1;

i=3*sin(w*t+pi/3)+exp(-t/ta);

Pomiar częstotliwości bez adaptacji:	Pomiar amplitudy bez adaptacji:
Pomiar częstotliwości z adaptacją:	Pomiar amplitudy z adaptacją:

c) Zakłócenia harmoniczne i nieharmoniczne:

- składowa harmoniczna:

u=15*sin(w*t)+sin(150*w*t)

Pomiar częstotliwości bez adaptacji:	Pomiar amplitudy bez adaptacji:
Pomiar częstotliwości z adaptacją:	Pomiar amplitudy z adaptacją:

-składowa stala +1

Pomiar częstotliwości bez adaptacji:	Pomiar amplitudy bez adaptacji:
Pomiar częstotliwości z adaptacją:	Pomiar amplitudy z adaptacją:

d) Odchylenia częstotliwości +- 5Hz

-f=55 Hz

Pomiar częstotliwości bez adaptacji:	Pomiar amplitudy bez adaptacji:
Pomiar częstotliwości z adaptacją:	Pomiar amplitudy z adaptacją:

-f=45 Hz

Pomiar częstotliwości bez adaptacji:	Pomiar amplitudy bez adaptacji:
Przy obniżonej częstotliwości algorytm adaptacyjny nie zadziałał

e) Zmiana częstotliwości próbowania

-fp=1000

Pomiar częstotliwości bez adaptacji:	Pomiar amplitudy bez adaptacji:
Pomiar częstotliwości z adaptacją:	Pomiar amplitudy z adaptacją:

3. Wnioski:

• Dla przebiegow niezak

---