1、（ 2）、修正的周期图法 1、 N=256， M=32 0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5 4 4 . 5 5x 1 04- 1 5- 1 0-50510改进的谱估计 N = 2 5 6 ， M = 3 22、 N=256， M=128 0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5 4 4 . 5 5x 1 04- 2 0- 1 5- 1 0-5051015改进的谱估计 N = 2 5 6 ， M = 1 2 83、 N=512， M=128 0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5 4 4 . 5 5x 1 04- 2 0-
2、 1 5- 1 0-5051015改进的谱估计 N = 5 1 2 ， M = 1 2 8程序： 1.function Per1 N=256； fs=100000； f = 0.02； n = 1 : N； s(n) = sin(2 * pi * f * n)； x = awgn(s, 10)； window1 = hamming(32)； overlap=10； Pxx1,f1=psd(x,N,fs,window1,overlap)； log_Pxx1=10*log10(Pxx1)； n1=0:round(N/2-1)； f1=n1*fs/N； plot(f1,log_Pxx1(n1+1)；
3、 title(改进的谱估计 N=256， M=32)； 2. function Per2 N=256； fs=100000； f = 0.02； n = 1 : N； s(n) = sin(2 * pi * f * n)； x = awgn(s, 10)； window1 = hamming(128)； overlap=10； Pxx1,f1=psd(x,N,fs,window1,overlap)； log_Pxx1=10*log10(Pxx1)； n1=0:round(N/2-1)； f1=n1*fs/N； plot(f1,log_Pxx1(n1+1)； title(改进的谱估计 N=256， M=128)； 3. function Per3 N=512； fs=100000； f = 0.02； n = 1 : N； s(n) = sin(2 * pi * f * n)； x = awgn(s, 10)； window1 = hamming(128)； overlap=10； Pxx1,f1=psd(x,N,fs,window1,overlap)； log_Pxx1=10*log10(Pxx1)； n1=0:round(N/2-1)； f1=n1*fs/N； plot(f1,log_Pxx1(n1+1)； title(改进的谱估计 N=512， M=128)；
