% İKİ SİNÜZOİDİN TOPLAMINDAN OLUŞAN SİNYALDE FFT ÇİZİMİ

% Sinyalin formüle edilmesi

Fs = 1000;                    % Örnekleme frekansı
T = 1/Fs;                     % Örnekleme periyodu
L = 1000;                     % Sinyalin uzunluğu
t = (0:L-1)*T;                % Time vector
% Sinyalimizi 50 Hz ve 120 Hz' lik iki sinüzoidalin toplamı ise
x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t); 
figure(1)
plot(Fs*t(1:50),x(1:50))
title('Orjinal sinyal')
xlabel('Zaman (milisaniye)')
ylabel('Genlik (Volt)')
y = x + 2*randn(size(x));     % Random gürültü fonksiyonu
figure(2)
plot(Fs*t(1:50),y(1:50))
title('Gürültü ile bozuluma uğramış sinyaller')
xlabel('Zaman (milisaniye)')
ylabel('Genlik (Volt)')

% Sinyalin FFT' sinin alınması
NFFT = 2^nextpow2(L); % L uzunluğundan fazla 2^nin en küçük katı
Y = fft(y,NFFT)/L;  %% NFFT uzunluğunda frekans domenine bilgi taşınıyor
f = Fs/2*linspace(0,1,NFFT/2+1);  %frekans ekseni oluştuur

% Plot single-sided amplitude spectrum.
figure(3);
plot(f,2*abs(Y(1:NFFT/2+1))) 
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

% Sinyalin gürültü eklenmemiş halinin FFT' si
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
X = fft(x,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);

% Plot single-sided amplitude spectrum.
figure(4);
plot(f,2*abs(X(1:NFFT/2+1))) 
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

ind=abs(Y)>(0.3);
denoised_Y=ind.*Y;
A=ifft(denoised_Y,NFFT)/L;
figure(5)
plot(Fs*t(1:50),(A(1:50)))
title('Denoised Signal')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')