用MATLAB绘制高斯白噪声的功率谱

要绘制高斯白噪声的功率谱，可以使用MATLAB中的fft函数和periodogram函数。下面是一个示例代码：

% 生成高斯白噪声
Fs = 1000; % 采样率
T = 1/Fs; % 采样间隔
L = 1000; % 信号长度
t = (0:L-1)*T; % 时间向量
x = randn(size(t)); % 高斯白噪声

% 计算功率谱
N = length(x); % 信号长度
xdft = fft(x); % 快速傅里叶变换
Px = 1/(Fs*N) * abs(xdft).^2; % 计算功率谱密度
f = Fs*(0:(N/2))/N; % 频率向量

% 绘制功率谱
plot(f, 10*log10(Px(1:N/2+1)))
title('Power Spectral Density')
xlabel('Frequency (Hz)')
ylabel('Power (dB/Hz)')

这段代码首先生成一个长度为1000的高斯白噪声信号x。然后使用FFT计算信号的傅里叶变换，再根据傅里叶变换结果计算功率谱密度Px。最后，使用plot函数绘制频率范围内的功率谱。

注意，上述代码中使用了10*log10函数将功率谱转换为以dB/Hz为单位的对数尺度