matlab bt法 功率谱估计代码

以下是MATLAB中使用bt函数进行功率谱估计的示例代码：

% 生成样本信号
Fs = 1000; % 采样率
t = 0:1/Fs:1-1/Fs; % 时间向量
x = cos(2*pi*50*t) + randn(size(t)); % 信号

% 计算功率谱密度
[Pxx,f] = bt(x,4,[],Fs);

% 绘制功率谱密度曲线
plot(f,Pxx);
xlabel('频率 (Hz)');
ylabel('功率谱密度');
title('BT法功率谱密度估计');

在上述代码中，`bt`函数的第一个输入参数`x`是输入的信号，后面三个输入参数均为可选参数，其中第四个参数`Fs`为采样率。函数的输出包括功率谱密度`Pxx`和对应的频率向量`f`。

执行上述代码后，将会绘制出信号的功率谱密度估计曲线。