clear?all
x=0:1:1000;
Fs=1000;?
nfft=1024;
y=3*square（2*pi*0.01*x50）;???????%矩形波
subplot（321）;plot（xy）;grid;axis（[01000-55]）;
title（‘模擬信號時域波形‘）xlabel（‘時間?t/s‘）ylabel（‘電壓?V/v‘）;
yBW1=abs（fft（ynfft））;?????%模擬信號的?頻域波形
subplot（322）plot（（0:nfft/2-1）/nfft*FsyBW1（1:nfft/2））;
axis（[0250-50400]）;
title（‘模擬信號頻譜‘）xlabel（‘頻率?f/Hz‘）ylabel（‘振幅‘）;

yn=awgn（y15‘measured‘）;?????????%增加噪聲信號
subplot（323）;plot（yn）;grid;axis（[01000-55]）;??
title（‘含噪信號時域波形‘）xlabel（‘時間?t/s‘）ylabel（‘電壓?V/v‘）;
yBW2=abs（fft（ynnfft））;???%含噪信號的?頻域波形???
subplot（324）plot（（0:nfft/2-1）/nfft*FsyBW2（1:nfft/2））;?
axis（[0250-50400]）;
title（‘含噪信號頻譜‘）xlabel（‘頻率?f/Hz‘）ylab