x1=xcorr(x);                                         %求x的自相關函數x1
x1=fft(x1);                                         %求x的一階傅立葉變換
subplot(2,1,1);     
plot(abs(x1))
title('數字基帶信號功率譜密度')              %數字基帶信號畫出功率譜密度圖



yp=xcorr(y);                                         %求x的自相關函數x1
yp=fft(yp);                                         %求x的一階傅立葉變換
     subplot(2,1,2);  
plot(abs(yp))
title('數字基帶信號功率譜密度')              %數字基帶信號畫出功率譜密度圖