MATLAB 产生功率谱密度为0.002的白噪声

可以使用MATLAB中的函数`wgn`产生功率谱密度为0.002的白噪声。具体代码如下：

Fs = 1000; % 采样率
T = 1/Fs; % 采样时间间隔
L = 10000; % 信号长度
t = (0:L-1)*T; % 时间向量
P = 0.002; % 功率谱密度
N = 2*Fs; % FFT点数
x = wgn(L, 1, 10*log10(P), 'linear'); % 产生白噪声
[Pxx, f] = pwelch(x, [], [], N, Fs); % 计算功率谱密度
plot(f, Pxx); % 绘制功率谱密度曲线

解释一下代码：

- `Fs`是采样率，这里设为1000Hz；
- `T`是采样时间间隔，由采样率计算得到；
- `L`是信号长度，这里设为10000；
- `t`是时间向量，用于绘制信号波形；
- `P`是所需的功率谱密度，这里设为0.002；
- `N`是FFT的点数，取2倍采样率；
- `x`是使用`wgn`函数产生的白噪声信号；
- `pwelch`函数用于计算功率谱密度；
- 最后用`plot`函数绘制功率谱密度曲线。

运行以上代码，即可得到功率谱密度为0.002的白噪声信号。