经典功率谱估计和现代谱估计方法

信号为两个正弦信号加高斯白噪声，各正弦信号的信噪比均为10dB，长度为N，信号频率分别为f1，f2；初始相位均为0，取f1/fs=0.2，f2/fs取不同数值：0.3、0.25,，fs为采样频率。利用经典功率谱估计和现代功率估计方法进行功率谱估计

%%现代估计
clc;
clear all;
fs=1;%采样频率
N=50;   %信号长度
f1=0.2*fs;
f2=0.3*fs;
M2=15;      %模型阶数
M1=10;
M0=6;
N_0=10;%分段长度
N_1=45;
N_2=50;
n=0:N-1;     
n=n./fs;   %时间序列
x=sin（2*pi*f1*n）+sin（2*pi*f2*n）;
x=awgn（x10）;%加入10dB高斯噪声的信号
[px_1f_1]=pburg（xM2[][]）;%阶数15 Burg法估计
[px_2f_2]=pyulear（xM2[][]）;%阶数15 Yule-walker法估计
[px_11f_11]=pburg（xM1[][]）;%阶数10 Burg法估计
[px_21f_21]=pyulear（xM1[][]）;%阶数10 Yule-walker法估计
[px_10f_10]=pburg（xM0[][]）;%阶数6 Burg法估计
[px_20f_20]=pyulear（xM0[][]）;%阶数6 Yule-walker法估计
figure;
subplot（221）;
plot（nx）;
xlabel（‘t（s）‘）;ylabel（‘x‘）;title（‘原始信号‘）;
subplot（222）;
plot（f_2px_2）;
xlabel（‘频率‘）;ylabel（‘幅度（dB）‘）;title（‘Yule-walker法估计功率谱（长度100，阶数15）‘）;
subplot（223）;
plot（f_1px_1）;
xlabel（‘频率‘）;ylabel（‘幅度（dB）‘）;title（‘Burg法估计功率谱（长度100，阶数15） ‘）;

figure;
subplot（221）;
plot（nx）;
xlabel（‘t（s）‘）;ylabel（‘x‘）;title（‘原始信号‘）;
subplot（222）;
plot（f_21px_21）;
xlabel（‘频率‘）;ylabel（‘幅度（dB）‘）;title（‘Yule-walker法估计功率谱（长度50，阶数10）‘）;
subplot（223）;
plot（f_11px_11）;
xlabel（‘频率‘）;ylabel（‘幅度（dB）‘）;title（‘Burg法估计功率谱（长度50，阶数10） ‘）;
figure;
subplot（221）;
plot（nx）;
xlabel（‘t（s）‘）;ylabel（‘x‘）;title（‘原始信号‘）;
subplot（222）;
plot（f_20px_20）;
xlabel（‘频率‘）;ylabel（‘幅度（dB）‘）;title（‘Yule-walker法估计功率谱（长度100，阶数