Periodogram of a noisy sine wave and calculation of PSD with a triangular window.m

%Bora KIŞ 21.04.2019 3.3. Gürültülü Bir Sinüs Dalgasının Periodogram ve Üçgen Pencere İle PSD’sinin Hesaplanması
clc
close all
clear all
t = 0:1/1000:2; %Zaman vektörü
N = length(t); %Örnekleme sayısı
fsig = 50; %Birinci sinyalin frekansı
signal = sin(2*pi*fsig*t); %Birinci sin sinyali
fsig2 = 140; %İkinci sinyalin frekansı
signal2 = sin(2*pi*fsig2*t); %İkinci sin sinyali
noise = 2*randn(1,N); %Randn komutu ile oluşturulan rastgele sinyal
x = signal + signal2 + noise; %Sinyallerin toplamı
x = triang(64);

pxx1 = periodogram(x,triang(64),[],64);
subplot(4,1,1);
plot(pxx1);
xlabel('Frequency (Hz)');
ylabel('PSD');

pxx2 = periodogram(x,triang(64),[],128);
subplot(4,1,2);
plot(pxx2);
xlabel('Frequency (Hz)');
ylabel('PSD');

pxx3 = periodogram(x,triang(64),[],256);
subplot(4,1,3);
plot(pxx3);
xlabel('Frequency (Hz)');
ylabel('PSD');

pxx4 = periodogram(x,triang(64),[],512);
subplot(4,1,4);
plot(pxx4);
xlabel('Frequency (Hz)');
ylabel('PSD')