Calculating PSD with a Welch Method Periodogram of a Noisy Sine Wave.m

%Bora Kış 20.04.2019 3.1. Gürültülü Bir Sinüs Dalgasının Welch Metodu Periodogramı İle PSD’sinin Hesaplanması
clc
close all
clear all

t = 0:1/1000:1; %Zaman vektörü
N = length(t); %Örnekleme sayısı
fsig = 50; %Birinci sinyalin frekansı
signal = sin(2*pi*fsig*t); %Birinci sin sinyali
fsig2 = 100; %İkinci sinyalin frekansı
signal2 = sin(2*pi*fsig2*t); %İkinci sin sinyali
noise = 0.25*randn(1,N); %Randn komutu ile oluşturulan rastgele sinyal
x = signal + signal2 + noise; %Sinyallerin toplamı

pxx = pwelch(x,hamming(length(x)),[],64,256); %Welch metodu
subplot(4,1,1);
plot(pxx)
xlabel('Frequency (Hz)');
ylabel('PSD');

pxx2 = pwelch(x,hamming(length(x)),[],128,256); %Welch metodu
subplot(4,1,2);
plot(pxx2);
xlabel('Frequency (Hz)');
ylabel('PSD');

pxx3 = pwelch(x,hamming(length(x)),[],256,256); %Welch metodu
subplot(4,1,3);
plot(pxx3);
xlabel('Frequency (Hz)');
ylabel('PSD');

pxx4 = pwelch(x,hamming(length(x)),[],512,256); %Welch metodu
subplot(4,1,4);
plot(pxx4);
xlabel('Frequency (Hz)');
ylabel('PSD');