clc;
clear?all;
close?all;

M?=?1024;?%行
N?=?1;???%列
variance?=?0.1;??????%噪聲方差
P?=?10*log10（variance）;?%以dBW為單位的輸出噪聲強度
y?=?wgn（MNP）;?Py?=?var（y）;
FFT_Y?=?fft（yM）;?%轉換到頻域
y_cor=xcorr（y‘unbiased‘）;?%計算序列的自相關函數
y_corfft=fft（y_cor）;?%自相關函數的傅里葉變換是功率譜密度

n1?=?1:M;?n2?=?1:M*2-1;
figure（1）
subplot（221）
plot（n