signal = 1+0.1*cos(2*pi*f1*t) + 0.7*cos(2*pi*f2*t)，matlab功率谱dBm

在MATLAB中，信号 `signal = 1 + 0.1*cos(2*pi*f1*t) + 0.7*cos(2*pi*f2*t)` 表示一个包含两个正弦波分量的简单时间序列，其中第一个频率是 `f1`，第二个频率是 `f2`。为了计算这个信号的功率谱（PSD, Power Spectral Density），通常会使用 `pwelch` 或 `fft` 函数，并将结果转换为dBm单位。

以下是计算步骤：

1. **创建信号**：

   t = linspace(0, 1, 1000); % 创建时间向量
   f1 = 1; % 第一频率 (Hz)
   f2 = 2; % 第二频率 (Hz)
   signal = 1 + 0.1*cos(2*pi*f1*t) + 0.7*cos(2*pi*f2*t);

2. **计算功率谱**：
- 对于非窗函数，使用 `pwelch` 函数（推荐）：

     [Pxx, f] = pwelch(signal, [], [], [], 'onesided');

- 或者使用离散傅立叶变换（FFT）：

     N = length(signal);
     Pxx = abs(fft(signal)) ./ N; % 去中心化并归一化
     f = (0:N-1) * (1/N) * fs; % 创建频率轴，fs是采样率

3. **转换到dBm**：

   Pxx_dBm = 10*log10(Pxx) + 30; % 以dBm为单位，30是对数基的修正常数

4. **显示结果**：

   plot(f, Pxx_dBm);
   xlabel('Frequency (Hz)');
   ylabel('Power Spectral Density (dBm)');