ex6_44.m

来自「ex5_1 I-型线性相位FIR滤波器 ex5_2 II-型线性相位FIR」· M 代码 · 共 12 行

nfft = 1024;
fs = 1000;
t = 0:1/fs:1;
xn = sin(2*pi*100*t)+sin(2*pi*300*t)+sin(2*pi*450*t)+randn(size(t));
[P, f] = pburg(xn, 18, nfft, fs);
Pxx = 10 * log10(P);
figure
plot(f, Pxx);
grid on;
xlabel('频率（Hz）');
ylabel('幅度（dB）');