Podstawy cyfrowej automatyki elektroenergetycznej
Wydział Elektryczny Semestr VI
Anna Knap 170739	Ćwiczenie : Przetwarzanie sygnałów przez przekładniki prądowe i napięciowe, efekty dyskretyzacji sygnałów, rola filtrów analogowych.	Termin zajęć: Poniedziałek 13-15
				Ocena:
07.03.2011r.

1. Cel ćwiczenia:

Celem ćwiczenia było wykreślanie transformaty Fouriera różnych sygnałów.

2. Przebieg ćwiczenia:

Zadanie 1:

Zdefiniowanie 4 sygnałów:

clear all;clc;close all;

A=1.0;A0=1.0;

ta=0.1;

w=2*pi*50;

Tp=1/1200;

t=0:Tp:0.2

b=1/3;c=1/6

fp=1200;

Le=length(t)-1;

f=0:fp/Le:fp;

syg1=A*sin(w*t);

syg2=A*sin(w*t)+A0*exp(-t/ta);

syg3=A*sin(w*t)+b*sin(2*pi*150*t)+c*sin(2*pi*300*t);

syg4=A*sin(w*t)-sin(2*pi*1150*t);

Wykreślanie dla kolejnych sygnałów transformaty Fouriera:

Sygnał 1:

y=fft(syg1); %obliczanie DFT ciągu x

m=abs(y)/(Le/2); % widmo amplitudowe

figure(1)

plot(f,m)

Sygnał 2:

y1=fft(syg2); %obliczanie DFT ciągu x

m1=abs(y1/(Le/2)); % widmo amplitudowe

figure(2)

plot(f,m1)

Sygnał 3:

y2=fft(syg3); %obliczanie DFT ciągu x

m2=abs(y2)/(Le/2); % widmo amplitudowe

figure(3)

plot(f,m2)

Sygnał 4:

y3=fft(syg4); %obliczanie DFT ciągu x

m3=abs(y3)/(Le/2); % widmo amplitudowe

figure(4)

plot(f,m3)

Zadanie 2:

Zdefiniowanie sygnału, o skokowej zmianie amplitudy:

A=1.0;A0=1.0;

fp=1200;

Tp=1/1200;

w=2*pi*50;

A=1.0;

A0=1.0;

t=0:Tp:0.2;

t1=0:Tp:0.1;

N=4;

Le=length(t1)-1;f=0:fp/Le:fp;

t2=0.1+Tp:Tp:0.2;

s1=A*sin(w*t1);

s2=N*sin(w*t2);

t=[t1 t2]; s=[s1 s2];

figure(1)

plot(t,s)

Transformata Fouriera tego sygnału przed zmianą amplitudy:

y4=fft(s1); %obliczanie DFT ciągu x

m4=abs(y4)/(Le/2); % widmo amplitudowe

figure(5)

plot(f,m4)

Transformata Fouriera tego sygnału po zmianie amplitudy:

Le1=length(t2)-1;

f=0:fp/Le1:fp;

y5=fft(s2); %obliczanie DFT ciągu x

m5=abs(y5)/(Le1/2); % widmo amplitudowe

figure(6)

plot(f,m5)

Wycięcie fragmentu sygnału do wykreślenia tranformaty wokół momentu przeskoku amplitudy:

sp=s(96:144);

tp=t(96:144);

figure(7)

plot(tp,sp)

Transformata Fouriera wokół momentu przeskoku amplitudy:

Le2=length(tp)-1;

f=0:fp/Le2:fp;

y6=fft(sp); %obliczanie DFT ciągu x

m6=abs(y6)/(Le1/2); % widmo amplitudowe

p6=unwrap(angle(y6))

figure(8)

plot(f,m6)

Zadanie 3:

Wczytanie sygnałów prądowych z pliku ATP i wykreślenie jednego z nich:

x=readpl452;

y=fft(x(:,5));

plot(x(:,5))

Wykreślenie transformaty Fouriera tego sygnału:

t=x(:,1);

y=fft(x(:,5));

Le1=length(t)-1;

f=0:1600/Le1:1600;

m=abs(y);

figure(2)

plot(f,m)

Wykreślenie innego sygnału z tego samego pliku:

figure(10)

plot(x(:,8))

y1=fft(x(:,8));

t=x(:,1);

Le2=length(t)-1;

f2=0:1600/Le2:1600;

m1=abs(y1);

figure(9)

plot(f2,m1)

Oraz jego tranformaty Fouriera:

3. Wnioski:

??????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! :D

---