function[ampl,faza]=funkcja1(A,B,C,D)
%% funkcja liczy wyswietla amplitude i faze sygnau 
%% transformaty  fouriera z A*sin(B*x+C)+D
clear all;
ampl=0;
faza=0;
x=[0:99]
for?x
ampl=abs(fft(A*sin(B*x+C)+D));
faza=unwrap(angle(fft(A*sin(B*x+C)+D));
hold on;
subplot(2,1,1);
plot(ampl)
subplot(2,1,2);
plot(faza)
hold off;
%% koniec tej funkcji
