randn('state',0);

Fs = 1024;   
t = 0:1/Fs:1-1/Fs;

A1signal = 3;
F1 = 64;

A2signal = 6;
F2 = 96;

Anoise = 0.0;

x = (A1signal*cos(2*pi*t*F1)) + (A2signal*cos(2*pi*t*F2)) + (Anoise * randn(size(t)));   

Fxx = 2* abs(fft(x));

%[Pxx, f] = pwelch(x,rectwin(1024),0,length(x),Fs);
pwelch(x,rectwin(1024),0,length(x),Fs);

windowName = 'Hamming';
segmentLength = 1024;
overlapPercent = 50;
%calculate the PSD
tensWSHsTest = spectrum.welch(windowName, segmentLength, overlapPercent);
tensionWSPSDDT = psd(tensWSHsTest, detrend(tensionWS), 'Fs', Fs);
tensionWSPSDNDT = psd(tensWSHsTest, tensionWS, 'Fs', Fs);