function[]=funkcja2()
%% funkcja generuje wektory sygnalow: 
%% 2*sin(x)+4*sin(3*x)+sin(2*x)+szum
%% sygnal prostokatny
%% funkcja trojkatna
%% funkcja Gaussa
%% i liczy fft oraz wyswietla amplitude i faze sygnalow

x=[0:0.1:99]
sinus=2*sin(x)+4*sin(3*x)+sin(2*x)+0.5*rand(1);
prostok=square(x);
trojk=sawtooth(x);
GaUss=exp(x.*x/(2*10*10));

four1=fft(sinus);
four2=fft(prostok);
four3=fft(trojk);
four4=fft(GaUss);
ampl1=abs(four1).^2;
faza1=unwrap(angle(four1));
ampl1=abs(four1).^2;
faza2=unwrap(angle(four2));
ampl2=abs(four1).^2;
faza1=unwrap(angle(four1));
ampl1=abs(four1).^2;
faza1=unwrap(angle(four1));
hold on;
subplot(4,2,1);
plot(sinus)
subplot(4,2,2);
plot(ampl)
subplot(4,2,3);
plot(faza)
hold off;
%% koniec tej funkcji
